[gdal-dev] frmt_wms_googlemaps_tms.xml help ...
Travis Kirstine
traviskirstine at gmail.com
Thu Dec 6 06:25:48 PST 2012
Mario,
>From what I can tell the output of the gdalinfo is fine, I'm guessing that
it is the correct mercator / dms coords of the bounding box?
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)
I would try the gdalwarp using the -s_srs flag
gdalwarp -s_srs EPSG:900913 -t_srs "+proj=utm +zone=32 +datum=wgs84"
tempfile.tif tempfile-warp.tif
or
gdalwarp -s_srs EPSG:3857 -t_srs "+proj=utm +zone=32 +datum=wgs84"
tempfile.tif tempfile-warp.tif
Regards
On 6 December 2012 04:15, Mario Cavicchi <cavicchi at ferrara.linux.it> wrote:
> 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}<http://mt.google.com/vt/lyrs=s&x=$%7Bx%7D&y=$%7By%7D&z=$%7Bz%7D>
> </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 <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<http://members.ferrara.linux.it/cavicchi>>
>> ]
>> ______________________________**___________________
>> gdal-dev mailing list
>> gdal-dev at lists.osgeo.org <mailto:gdal-dev at lists.osgeo.**org<gdal-dev at lists.osgeo.org>
>> >
>> http://lists.osgeo.org/__**mailman/listinfo/gdal-dev<http://lists.osgeo.org/__mailman/listinfo/gdal-dev>
>> <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<http://members.ferrara.linux.it/cavicchi>]
> ______________________________**_________________
> gdal-dev mailing list
> gdal-dev at lists.osgeo.org
> http://lists.osgeo.org/**mailman/listinfo/gdal-dev<http://lists.osgeo.org/mailman/listinfo/gdal-dev>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20121206/ae90e648/attachment.html>
More information about the gdal-dev
mailing list