[GRASS-user] Re: grass-user Digest, Vol 26, Issue 25

[Richard Chirgwin]

>
> "Unsupported output format (image/geotiff)."
>
> downloading the capabilities file by hand it shows:
>
> <!-- MapServer version 5.0.0 OUTPUT=GIF OUTPUT=PNG OUTPUT=JPEG OUTPUT=WBMP OUTPUT=SVG SUPPORTS=PROJ SUPPORTS=FREETYPE SUPPORTS=WMS_SERVER SUPPORTS=W
> MS_CLIENT SUPPORTS=WFS_SERVER SUPPORTS=WFS_CLIENT SUPPORTS=WCS_SERVER SUPPORTS=GEOS INPUT=EPPL7 INPUT=POSTGIS INPUT=OGR INPUT=GDAL INPUT=SHAPEFILE -
> ->
> ....
>     <GetMap>
>       <Format>image/jpg</Format>
>       <Format>image/gif</Format>
>       <Format>image/png</Format>
>       <Format>image/png; mode=24bit</Format>
>       <Format>image/jpeg</Format>
>       <Format>image/wbmp</Format>
>       <Format>image/tiff</Format>
>       <Format>image/svg+xml</Format>
>
>
> no geotiff support.  set format= to one of the above. 
>
> I tried in a WGS84 lat/lon loc:
>
> g.region n=50:30N s=47:15N w=9E e=13:45E res=0:00:30
> r.in.wms mapserver="$SERVER" layer=TK50 out=TK50_Bavaria format=png -g
>
> it seemed to work, but I only got a blank map. Maybe you will have better luck. Using "-g" I got a "(c) Bayerische ..." in the image, without "-g" I didn't. shrug.
>
> maybe you can get something from v.in.wfs?
> or try QGIS's or GDAL 1.5+'s WMS import tools?
>
>
> Hamish
>   
I've never been able to decide whether "blank raster syndrome" in 
r.in.wms is at the client or the server. In some cases, I guess it's the 
server because things rectify themselves if you come back tomorrow. In 
other cases, adjusting the region parameters - zoom in to a smaller 
area, and use a finer resolution - fixes the problem.

For example, if I use OneEarth on a too-large region at low resolution, 
I'll get a blank map. If I zoom to a smaller region and reset the 
resolution, I will get a successful download.

This query:
r.in.wms output=Trans 
mapserver=http://onearth.jpl.nasa.gov/wms.cgi?request=GetCapabilities 
layers=global_mosaic srs=EPSG:4326 format=geotiff wmsquery=version=1.1.1 
maxcols=1024 maxrows=1024 'wgetoptions=-c -t 5 --user-agent=MSIE5.5' 
method=nearest v=1 --overwrite

Produces a blank raster on this region:
>
> projection: 3 (Latitude-Longitude)
> zone: 0
> datum: wgs84
> ellipsoid: a=6378137 es=0.006694379990141316
> north: 33:36:24.61482S
> south: 34:06:06.11586S
> west: 150:54:42.559128E
> east: 151:24:23.901624E
> nsres: 0:00:10.008433
> ewres: 0:00:10.007542
> rows: 178
> cols: 178
> cells: 31684

Refine the region to this:
> projection: 3 (Latitude-Longitude)
> zone: 0
> datum: wgs84
> ellipsoid: a=6378137 es=0.006694379990141316
> north: 33:50:37.371984S
> south: 33:51:46.767456S
> west: 151:11:49.039152E
> east: 151:12:53.75088E
> nsres: 0:00:00.030002
> ewres: 0:00:00.030001
> rows: 2313
> cols: 2157
> cells: 4989141
And everything works.

Richard