[GRASS-user] r.in.wms issues

[Markus Neteler]

On Tue, Jan 27, 2009 at 10:50 PM, Adam Dershowitz
<adershowitz at exponent.com> wrote:
> I am trying to use r.in.wms to get elevation data from wms.jpl.nasa.gov, but
> I have been running into two different problems.  I have been using the gui
> on a mac with 6.3, although I also tried 6.4 RC2 and got similar results.
>  The files do download but then I can't seem to work with them properly.  I
> have tried geotiff and png.  I have tried a bunch of other combinations as
> well.
>
>
> The first problem seems to be with how temporary files are being named so
> that patching is not working properly.  Here is the command that is being
> generated by the gui:
>
> r.in.wms output=elev_us_ned2 mapserver=wms.jpl.nasa.gov/wms.cgi
> layers=us_ned format=png wmsquery=version=1.1.1 maxcols=1024 maxrows=1024
> {wgetoptions=-c -t 5 -nv} {curloptions=-C - --retry 5 -s -S} method=cubic
> v=1
>
> This is followed by this long list of errors:
> Calculating tiles
> Requesting 2 tiles.
> Downloading tiles
> Tile already downloaded
> Tile already downloaded
> All tiles downloaded successfully
> Creating output file that is 972P x 558L.
> Processing input file /Users/dersh/grass/wms_download/elev_us_ned2__0.png.
> Using band 2 of source image as alpha.
> 0...10...20...30...40...50...60...70...80...90...100 - done.
> Invalid map <elev_us_ned2_tile_0_tmp.alpha>
> Parse error
> Raster map <elev_us_ned2_tile_0.1> not found
> Invalid map <elev_us_ned2_tile_0_tmp.alpha>
...

I tried in 6.4 and got (copying your command line):

r.in.wms output=elev_us_ned2 mapserver="wms.jpl.nasa.gov/wms.cgi" \
   layers=us_ned format=png wmsquery=version=1.1.1 maxcols=1024 \
   maxrows=1024 "wgetoptions=-c -t 5 -nv" curloptions="-C - --retry 5 -s -S" \
   method=cubic v=1
Calculating tiles
Requesting 1 tiles.
Downloading tiles
Downloading data
2009-01-27 22:53:36 URL:http://wms.jpl.nasa.gov/wms.cgi [53504] ->
"/home/neteler/grassdata/wms_download/elev_us_ned2__0.png" [1]
All tiles downloaded successfully
Creating output file that is 650P x 488L.
Processing input file /home/neteler/grassdata/wms_download/elev_us_ned2__0.png.
Using band 2 of source image as alpha.
0...10...20...30...40...50...60...70...80...90...100 - done.
Rename raster <elev_us_ned2_tile_0.1> to <elev_us_ned2.1>
Rename raster <elev_us_ned2_tile_0.2> to <elev_us_ned2.2>
ERROR: Raster map <elev_us_ned2> not found in current mapset
ERROR: Raster map <elev_us_ned2> not found in current mapset
ERROR: Raster map <elev_us_ned2> not found in current mapset
ERROR: Raster map <elev_us_ned2> not found in current mapset
ERROR: Raster map <elev_us_ned2> not found in current mapset
ERROR: Raster map <elev_us_ned2> not found in current mapset

Slightly better/different but eventually failing, too. :(

Markus