[Proj] UTM Zone projection parameters to pass to pj_transform()

[Eric Miller]

The results look like they're still in radians.  Multiply by 180 / M_PI (or use the RAD_TO_DEG macro).  You can use the pj_is_latlong() function to test if a projection is latlong and then conditionally convert from/to radians.

Eric G. Miller
GIS Analyst
Wildlife and Habitat Data Analysis Branch
Department of Fish and Game

>>> david.jacques at ftw.usda.gov 11/28/2005 12:13:28 PM >>>
The output from listgeo is this...

Corner Coordinates:
Upper Left    ( 360427.280,3278421.510)
Lower Left    ( 360427.280,3270801.510)
Upper Right   ( 367182.280,3278421.510)
Lower Right   ( 367182.280,3270801.510)
Center        ( 363804.780,3274611.510)

The corner coordinates in Decimal Degrees should be

Upper Left    ( 29.628626,-94.441707)
Lower Left    ( 29.559135,-94.441707)
Upper Right   ( 29.628626,-94.371019)
Lower Right   ( 29.559135,-94.371019

What I get is...

-1.648319,0.517104
-1.648302,0.515904
-1.647101,0.517117
-1.647085,0.515917

This prompted me to look at the PJ structure in debug mode (watch)
In my Visual Studio 7 IDE. I saw the following string in the
description...

+	descr	0x0043802c "Universal Transverse Mercator (UTM)	Cyl, Sph
zone= south"	const char *

Funny how I didn't specify "south".

This led me to believe I need to specify more parameters.


David

_______________________________________________
Proj mailing list
Proj at lists.maptools.org 
http://lists.maptools.org/mailman/listinfo/proj