[gdal-dev] frmt_wms_googlemaps_tms.xml help ...

[Mario Cavicchi]

The contents of the frmt_wms_googlemaps_tms.xml:

<GDAL_WMS>
     <Service name="TMS">

<ServerUrl>http://mt.google.com/vt/lyrs=s&x=${x}&y=${y}&z=${z}</ServerUrl>
     </Service>
     <DataWindow>
         <UpperLeftX>-20037508.34</UpperLeftX>
         <UpperLeftY>20037508.34</UpperLeftY>
         <LowerRightX>20037508.34</LowerRightX>
         <LowerRightY>-20037508.34</LowerRightY>
         <TileLevel>20</TileLevel>
         <TileCountX>1</TileCountX>
         <TileCountY>1</TileCountY>
         <YOrigin>top</YOrigin>
     </DataWindow>
     <Projection>EPSG:900913</Projection>
     <BlockSizeX>256</BlockSizeX>
     <BlockSizeY>256</BlockSizeY>
     <BandsCount>3</BandsCount>
     <MaxConnections>5</MaxConnections>
     <Cache />
</GDAL_WMS>

the gdalinfo:

Driver: GTiff/GeoTIFF
Files: tempfile.tif
Size is 363, 380
Coordinate System is:
PROJCS["Google Maps Global Mercator",
     GEOGCS["WGS 84",
         DATUM["WGS_1984",
             SPHEROID["WGS 84",6378137,298.257223563,
                 AUTHORITY["EPSG","7030"]],
             AUTHORITY["EPSG","6326"]],
         PRIMEM["Greenwich",0],
         UNIT["degree",0.0174532925199433],
         AUTHORITY["EPSG","4326"]],
     PROJECTION["Mercator_1SP"],
     PARAMETER["central_meridian",0],
     PARAMETER["scale_factor",1],
     PARAMETER["false_easting",0],
     PARAMETER["false_northing",0],
     UNIT["metre",1,
         AUTHORITY["EPSG","9001"]]]
Origin = (1299378.124538742238656,5588829.177167181856930)
Pixel Size = (0.298582141697407,-0.298582141697407)
Metadata:
   AREA_OR_POINT=Area
Image Structure Metadata:
   INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  ( 1299378.125, 5588829.177) ( 11d40'21.04"E, 44d59' 3.39"N)
Lower Left  ( 1299378.125, 5588715.716) ( 11d40'21.04"E, 44d59' 0.79"N)
Upper Right ( 1299486.510, 5588829.177) ( 11d40'24.55"E, 44d59' 3.39"N)
Lower Right ( 1299486.510, 5588715.716) ( 11d40'24.55"E, 44d59' 0.79"N)
Center      ( 1299432.317, 5588772.447) ( 11d40'22.80"E, 44d59' 2.09"N)
Band 1 Block=363x7 Type=Byte, ColorInterp=Red
Band 2 Block=363x7 Type=Byte, ColorInterp=Green
Band 3 Block=363x7 Type=Byte, ColorInterp=Blue


Regards,
Mario


On 12/05/2012 09:51 PM, Travis Kirstine wrote:
> what do the contents of the frmt_wms_googlemaps_tms.xml and the output
> of gdalinfo tempfile.tif look like
>
>
> On 5 December 2012 12:48, Mario Cavicchi <cavicchi at ferrara.linux.it
> <mailto:cavicchi at ferrara.linux.it>> wrote:
>
>     Hi folks,
>     Maybe this question is old, but I don't find an answer on Internet.
>     If I perform this command:
>
>     gdal_translate -outsize 50% 50% -of GTiff -projwin
>     1299378.2102987265 5588829.1661406181 1299486.6170906540
>     5588715.6742054483 frmt_wms_googlemaps_tms.xml tempfile.tif
>
>     the image obtained is perfect in term of area of interest, but the
>     georeference is wrong if I open the image with ENVI.
>     I think that it is caused for the "Google Maps Global Mercator", but
>     I don't know how I can convert the projection to (for example) an
>     UTM projection or something else. I tried to reproject using gdalwarp:
>
>     gdalwarp -t_srs "+proj=utm +zone=32 +datum=wgs84" tempfile.tif
>     tempfile-warp.tif
>
>     but, still wrong georeference.
>
>     Someone can help me?
>
>     --
>     Ciao,
>              Mario
>
>     [ Web Info: http://members.ferrara.linux.__it/cavicchi
>     <http://members.ferrara.linux.it/cavicchi> ]
>     _________________________________________________
>     gdal-dev mailing list
>     gdal-dev at lists.osgeo.org <mailto:gdal-dev at lists.osgeo.org>
>     http://lists.osgeo.org/__mailman/listinfo/gdal-dev
>     <http://lists.osgeo.org/mailman/listinfo/gdal-dev>
>
>


-- 
Ciao,
         Mario

[ Web Info: http://members.ferrara.linux.it/cavicchi ]