[Gdal-dev] gdal_translate Basic Usage / UTM to latlong

[Michael84]

You use the wrong command line tool...
You should use gdal_warp: "Warp an image into a new coordinate system"

Define the origin cs by -s_srs ... and the resulting cs by -t_srs (
http://www.gdal.org/gdalwarp.html )



John R. wrote:
> 
> Tim Osborn wrote:
>> OK...this is what's going on.....
>> 
>> I'm pretty sure the image is projected in the Wisconsin Transverse
>> Mercator
>> coordinate system WTM83.  You can read about it at
>> http://www.dnr.state.wi.us/maps/gis/wtm8391.html.  They list the
>> parameters
>> there and an email address you can request a prj file.  WTM83 is not the
>> same as regular old UTM NAD83.  The EPSG number code for WTM83 is 3070.
>> 
>> Here are the parameters for WTM83:
> 
> [snip]
> 
> Thanks!
> 
> I requested the prj file (says it is for ArcGIS).  I had done a quick 
> search for the EPSG number but did not find it.  I used the parameters 
> to construct this:
> 
>   -a_srs '+proj=tmerc +lat_0=0 +lon_0=-90 +x_0=520000 +y_0=-4480000 
> +k=0.9996 +ellps=GRS80 +units=m +datum=NAD83 +no_defs'
> 
> The results look pretty close.   Thanks again!
> 
> I just now tried 3070 and I dot get long and lat listed in gdalinfo
> output:
> 
> Driver: GTiff/GeoTIFF
> Size is 10750, 8488
> Coordinate System is:
> PROJCS["unnamed",
>      GEOGCS["NAD83",
>          DATUM["unknown",
>              SPHEROID["unretrievable - using
> WGS84",6378137,298.257223563]],
>          PRIMEM["Greenwich",0],
>          UNIT["degree",0.0174532925199433]],
>      UNIT["metre",1,
>          AUTHORITY["EPSG","9001"]],
>      AUTHORITY["EPSG","3070"]]
> Origin = (659653.722794383182190,352632.358040164515842)
> Pixel Size = (0.329187501713504,-0.329168110366799)
> Metadata:
>    AREA_OR_POINT=Area
>    TIFFTAG_XRESOLUTION=1200
>    TIFFTAG_YRESOLUTION=1200
>    TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)
> Corner Coordinates:
> Upper Left  (  659653.723,  352632.358)
> Lower Left  (  659653.723,  349838.379)
> Upper Right (  663192.488,  352632.358)
> Lower Right (  663192.488,  349838.379)
> Center      (  661423.106,  351235.369)
> Band 1 Block=10750x1 Type=Byte, ColorInterp=Red
> Band 2 Block=10750x1 Type=Byte, ColorInterp=Green
> Band 3 Block=10750x1 Type=Byte, ColorInterp=Blue
> _______________________________________________
> Gdal-dev mailing list
> Gdal-dev at lists.maptools.org
> http://lists.maptools.org/mailman/listinfo/gdal-dev
> 
> 

-- 
View this message in context: http://www.nabble.com/gdal_translate-Basic-Usage---UTM-to-latlong-tf4065198.html#a11600033
Sent from the GDAL - Dev mailing list archive at Nabble.com.