Re: [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&amp;x=${x}&amp;y=${y}&amp;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 <cavic...@ferrara.linux.it
<mailto:cavic...@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@lists.osgeo.org <mailto:gdal-dev@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 ]
_______________________________________________
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev