WMS: Difference between revisions

From GRASS-Wiki
Jump to navigation Jump to search
mNo edit summary
mNo edit summary
Line 2: Line 2:


{{wikipedia|Web_Map_Service|desc=Web Map Service}} (WMS) data can be downloaded and imported into GRASS using the {{cmd|r.in.wms}} module.
{{wikipedia|Web_Map_Service|desc=Web Map Service}} (WMS) data can be downloaded and imported into GRASS using the {{cmd|r.in.wms}} module.
Mini tutorial for downloading [[SRTM]] elevation and LANDSAT imagery from [http://onearth.jpl.nasa.gov NASA's OnEarth WMS server].
Mini tutorial for downloading [[SRTM]] elevation and [[LANDSAT]] imagery from [http://onearth.jpl.nasa.gov NASA's OnEarth WMS server].


Create a NZTM location (EPSG:2193) with the location wizard.
Create a NZTM location (EPSG:2193) with the location wizard.

Revision as of 15:38, 26 December 2011

Web_Map_Service (WMS) data can be downloaded and imported into GRASS using the r.in.wms module. Mini tutorial for downloading SRTM elevation and LANDSAT imagery from NASA's OnEarth WMS server.

Create a NZTM location (EPSG:2193) with the location wizard.

# Zoom to AKL region
g.region n=5936130 s=5900340 \
         w=1740030 e=1776030 res=30

If downloading for your own area, try to keep the region rows and columns smaller than the maxrow,maxcol options given to r.in.wms. Please don't abuse the server, just download small chunks.

### Download and process SRTM elevation:

r.in.wms output=akl_srtm \
  mapserver="http://onearth.jpl.nasa.gov/wms.cgi" \
  -c layers=worldwind_dem format=geotiff style=short_int \
  maxcols=1600 maxrows=1600

g.rename akl_srtm.1,akl_srtm
r.colors akl_srtm color=srtm

### Download and process LANDSAT:

r.in.wms output=akl_landsat \
  mapserver="http://onearth.jpl.nasa.gov/wms.cgi" \
  -c layers=global_mosaic_base format=geotiff style=visual \
  maxcols=1600 maxrows=1600

i.landsat.rgb -p r=akl_landsat.red g=akl_landsat.green b=akl_landsat.blue
d.rgb r=akl_landsat.red g=akl_landsat.green b=akl_landsat.blue
r.composite r=akl_landsat.red g=akl_landsat.green b=akl_landsat.blue \
   output=akl_landsat --overwrite

### Shaeded relief map
r.shaded.relief akl_strm
d.shadedmap relief=akl_strm.shade drape=akl_landsat bright=70

### 3D view
nviz akl_srtm color=akl_landsat

In the NVIZ Visualization → Raster Surface menu set the fine raster grid size to 1. Set the vertical exaggeration to 1.0 or 2.0