Richard Chirgwin wrote:
> I am trying to access layers on the Geosciences Australia
> OGC Web Map server:
> http://www.ga.gov.au/wms/getmap?dataset=national&request=getCapabilities
....
> Is it a bug in Grass-GIS or at the other end?

I'll be kind and say GRASS because /technically/ the other end conforms to the 
WMS spec. The problem was that the other end does not accept queries in the 
form of HTTP POST data. Yes, a query server which doesn't accept queries. 
[EMAIL PROTECTED]

Anyway, in the past I tried to auto-detect that and switch to the other method, 
but the error mode the remote server sends seems to vary for each, so it is 
hard to do that well. I've now added a new flag to r.in.wms (-g) to force using 
the GET method. With that I can fetch data from Geosciences Australia's WMS:

r.in.wms -l -g mapserv="http://www.ga.gov.au/wms/getmap?dataset=national&";

# search for Gold:
r.in.wms mapserv="http://www.ga.gov.au/wms/getmap?dataset=national&"; \
   -g layers=au_occdd output=au_gold format=tiff


I also discovered that Geosciences Australia's WMS server does not ship 
GeoTiffs, which is r.in.wms's default. So you'll have to use format=<other>.
An issue with format=tiff is that the alpha (transparency) layer that comes 
from Geosci. Aust. runs from 1-255 instead of 0-255, and so no-data is treated 
as some-data and the resulting data looks a bit funky. That is easily fixed 
with r.mapcalc or r.null after import or in the case of the Gold, a simple call 
to r.colors:
  echo "1023 yellow" | r.colors au_gold col=rules

format=gif and png are really slow as the individual tiles use their own custom 
8bit color palettes and when merged together the palette from the first tile is 
used and the rest of the tiles look nuts. So r.in.gdalwarp goes the long way to 
ensure good results, which is very very slow. gdal_merge.py is probably a much 
better method for the future.

format=jpeg now works (fixed in latest commit; it has no alpha layer), but the 
results are subject to lossy JPEG compression/blurring.


yikes.


If you get the contents of the scripts/r.in.wms/ dir from develbranch_6 and put 
replace 6.3.0's versions and re-make that dir you should be able get it to work 
(proof attached :)
http://trac.osgeo.org/grass/browser/grass/branches/develbranch_6/scripts/r.in.wms



calmly,
Hamish


ps- sorry if this email is wrapping at about char 90, yahoo as upgraded their 
system & removed the "wrap at 72" compose setting, I suppose to make things 
simpler for us.



      
____________________________________________________________________________________
Be a better friend, newshound, and 
know-it-all with Yahoo! Mobile.  Try it now.  
http://mobile.yahoo.com/;_ylt=Ahu06i62sR8HDtDypao8Wcj9tAcJ

<<attachment: au_gold.png>>

_______________________________________________
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Reply via email to