[Proj] request help with Transverse_Mercator in proj

Rick Graves gravesricharde at yahoo.com
Sun Aug 12 18:41:57 PDT 2012


Hey,

Sorry to trouble you, but I still have a couple of states for which I need to get proj.4 working.  This time, it is Wisconsin.

I am starting with the Western most point in the state.  If you plug the following into maps.google.com, you will see it:

45.64120292600006, -92.88943290699996 

I downloaded the shapefiles from the state web site.  In the shapefile, the minimum X point is this:

294822.7664999999, 578184.1052999999

Here is the metadata from the prj file:

PROJCS["NAD_1983_HARN_Wisconsin_TM",
GEOGCS["GCS_North_American_1983_HARN",
DATUM["D_North_American_1983_HARN",
SPHEROID["GRS_1980",6378137.0,298.257222101]],
PRIMEM["Greenwich",0.0],
UNIT["Degree",0.0174532925199433]],
PROJECTION["Transverse_Mercator"],
PARAMETER["False_Easting",520000.0],
PARAMETER["False_Northing",-4480000.0],
PARAMETER["Central_Meridian",-90.0],
PARAMETER["Scale_Factor",0.9996],
PARAMETER["Latitude_Of_Origin",0.0],
UNIT["Meter",1.0]]

I initialized proj as follows:

proj +proj=tmerc +lat_0=0.0 +lon_0=-90.0 +k=0.9996 +x_0=-4480000.0 +y_0=520000.0

When give proj this: 

-92.88943290699996 45.64120292600006

I get this:

-4705175.18     5578245.14

But I am expecting this:

294822.7664999999, 578184.1052999999

I assume I am not initializing proj correctly.  

If someone could spot my mistake, I would appreciate it.

Thanks,

Rick Graves
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/proj/attachments/20120812/61d17717/attachment.html>


More information about the Proj mailing list