[GRASS-user] Fwd: WCS import into GRASS GIS - select region extent only

[Martin Zbinden]

Hi,

I'm all new to this list, but I tried GRASS GIS already for a few weeks.
For my Bachelor-Thesis I want to optimize the erosion risk map of
Switzerland for use in agricultural extension. Although I concentrate
on 3 cantons only (BE, FR, SO),  I have to manage huge amounts of
data. The fact that the raster data I got from the federal office are
overlapping each other (following the defined fieldblocks) doesn't
make things easier, but I've at least found r.patch which can combine
them one by one.

I thought, with Mapserver  (V 6.4.0) I'd found the perfect solution to
all my problems. Now I'm stuck after 5+ hours of trying to import the
data from my local WCS-Service. I have seen the hints on GRASS Wiki
[1] and r.in.gdal starts to download geotiff-tiles from Mapserver. The
problem is, that I don't want all the 12 GB of data, but only the
grass region extent.

The following http-query does what I expect, it makes me a
GeoTIFF-File from my region defined by BBOX=... .
http://localhost/mywcs?SERVICE=WCS&VERSION=1.0.0&REQUEST=GetCoverage&coverage=blw.eros_z&FORMAT=image/tiff&CRS=EPSG:2056&BBOX=2596000,1181000,2597500,1182500&RESX=2&RESY=2


I tried to restrict the download to the  region by defining a BBOX in
this WCS_GDAL xml file. However, this didn't work and isn't supposed
to work [2]. Is there a way to say r.in.gdal to only get the data from
the region extent?

I know there is this new, very clever r.in.wms2 script, which does
exactly what I want -- for WMS. Is there a way to get WMS-Mapserver to
serve RAW-data (erosion in t per annum, elevation data)?  Or maybe
better create a python script to download map by a custom
http-request?

Thanks in advance for any advice! I hope I can get it to work somehow.
Liebe Grüsse
Martin Zbinden

PS: the WCS server given in GRASS wiki seems to be down this time.

[1] http://grasswiki.osgeo.org/wiki/Global_datasets#OGC_WCS_-_Albedo_example
[2] http://lists.osgeo.org/pipermail/gdal-dev/2008-March/016552.html


--
Martin Zbinden
Riedacker 523
3154 Rüschegg Heubach
+41 78 628 28 82