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

[John R.]

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