I too am interested in trying to use gdal_translate with a WMS getmap request, but using the approach suggested does not seem to preserve the geospatial info:
rsignell@gam:~$ export GDAL_SKIP=WMS rsignell@gam:~$ gdal_translate -of GTiff -co "TFW=YES" "http://geoport.whoi.edu/ thredds/wms/bathy/gom15?LAYERS=topo&ELEVATION=0&TIME=2013-08-21T06%3A26%3A50Z&TR ANSPARENT=true&STYLES=boxfill%2Frainbow&CRS=EPSG%3A4326&COLORSCALERANGE=-50%2C50 &NUMCOLORBANDS=64&LOGSCALE=false&SERVICE=WMS&VERSION=1.1.1&REQUEST=GetMap&EXCEPT IONS=application%2Fvnd.ogc.se_inimage&FORMAT=image%2Fpng&SRS=EPSG%3A4326&BBOX=-7 3.75,37.5,-69.5,41.75&WIDTH=256&HEIGHT=256" foo2.tif Warning 6: A dataset opened by GDALOpenShared should have the same filename (htt p://geoport.whoi.edu/thredds/wms/bathy/gom15?LAYERS=topo&ELEVATION=0&TIME=2013-0 8-21T06%3A26%3A50Z&TRANSPARENT=true&STYLES=boxfill%2Frainbow&CRS=EPSG%3A4326&COL ORSCALERANGE=-250%2C50&NUMCOLORBANDS=64&LOGSCALE=false&SERVICE=WMS&VERSION=1.1.1& REQUEST=GetMap&EXCEPTIONS=application%2Fvnd.ogc.se_inimage&FORMAT=image%2Fpng&SR S=EPSG%3A4326&BBOX=-73.75,37.5,-69.5,41.75&WIDTH=256&HEIGHT=256) and description (/vsimem/http_1/file.dat) Input file size is 256, 256 0...10...20...30...40...50...60...70...80...90...100 - done. rsignell@gam:~$ gdalinfo foo2.tif Driver: GTiff/GeoTIFF Files: foo2.tif Size is 256, 256 Coordinate System is `' Image Structure Metadata: INTERLEAVE=BAND Corner Coordinates: Upper Left ( 0.0, 0.0) Lower Left ( 0.0, 256.0) Upper Right ( 256.0, 0.0) Lower Right ( 256.0, 256.0) Center ( 128.0, 128.0) Band 1 Block=256x32 Type=Byte, ColorInterp=Palette NoData Value=0 The only approach that I got working was passing the WMS request params to gdal_translate again as shown here: http://gis.stackexchange.com/questions/69188/converting-wms-getmap-request-to-geotiff Thanks, Rich On Sun, Apr 14, 2013 at 12:10 PM, Even Rouault <[email protected]> wrote: > Le dimanche 14 avril 2013 13:17:30, Andrea Peri a écrit : >> Even, thx for response >> >> Perhaps I miss the real mean of the parameter -srcwin. >> >> I guess that setting -srcwin 0 0 377 283 for a wms image of >> ...WIDTH=377&HEIGHT=283&... > > Ah, I didn't notice those parameters. Well, they will be ignored by the GDAL > WMS driver. You can check by looking at gdalinfo : > > $ gdalinfo > "WMS:http://www502.regione.toscana.it/cartografia/wmsraster/com.rt.wms.RTmap/wms?map=wmsofc&language=ita&&SERVICE=WMS&VERSION=1.3.0&REQUEST=GetMap&BBOX=1671200.0,4680000.0,1672800.0,4681200.0&CRS=EPSG:3003&WIDTH=377&HEIGHT=283&LAYERS=rt_ofc.10k10&STYLES=&FORMAT=image/jpeg&DPI=96&MAP_RESOLUTION=96&FORMAT_OPTIONS=dpi:96&TRANSPARENT=TRUE"; > Driver: WMS/OGC Web Map Service > Files: none associated > Size is 1073741824, 805306368 > > Actually, as your string is a full WMS GetMap request, you can directly do > this : > > export GDAL_SKIP=WMS > gdal_translate -of GTiff -co "TFW=YES" "http://...."; output_01J31.tif > > Note the removal of "WMS:" prefix in front of the URL, and the GDAL_SKIP=WMS, > since the WMS driver would also try to recognize the string as a description > string of the WMS service, which you don't want. You just want the JPEG driver > to recognize the image resulting of the GetMap query. > _______________________________________________ > gdal-dev mailing list > [email protected] > http://lists.osgeo.org/mailman/listinfo/gdal-dev -- Dr. Richard P. Signell (508) 457-2229 USGS, 384 Woods Hole Rd. Woods Hole, MA 02543-1598 _______________________________________________ gdal-dev mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/gdal-dev
