[gdal-dev] bad coordinates after transformation

[David Muse]

Hello all,

I wrote some code that reads data from a shapefile, converts it to WGS84
and puts it in a database.  It works fine with some shapefiles, but with
others the lat/lon that I get are really big numbers, like 500000 or
more.

It's possible that there's a bug in my code, but is there a utility that
comes with gdal that I could aim at the shapefile to see if it comes up
with a different result?

Maybe there's something odd with the .prj files too.

I get good results with this .prj file:

PROJCS["NAD_1927_StatePlane_Georgia_West_FIPS_1002",
	GEOGCS["GCS_North_American_1927",
		DATUM["D_North_American_1927",
			SPHEROID["Clarke_1866",6378206.4,294.9786982]
		],
		PRIMEM["Greenwich",0.0],
		UNIT["Degree",0.0174532925199433]
	],
	PROJECTION["Transverse_Mercator"],
	PARAMETER["False_Easting",500000.0],
	PARAMETER["False_Northing",0.0],
	PARAMETER["Central_Meridian",-84.16666666666667],
	PARAMETER["Scale_Factor",0.9999],
	PARAMETER["Latitude_Of_Origin",30.0],
	UNIT["Foot_US",0.3048006096012192]
]

But with this one I get the odd coordinates:

PROJCS["NAD_1983_UTM_Zone_16N",
	GEOGCS["GCS_North_American_1983",
		DATUM["D_North_American_1983",
			SPHEROID["GRS_1980",6378137.0,298.257222101]
		],
		PRIMEM["Greenwich",0.0],
		UNIT["Degree",0.0174532925199433]
	],
	PROJECTION["Transverse_Mercator"],
	PARAMETER["False_Easting",500000.0],
	PARAMETER["False_Northing",0.0],
	PARAMETER["Central_Meridian",-87.0],
	PARAMETER["Scale_Factor",0.9996],
	PARAMETER["Latitude_Of_Origin",0.0],
	UNIT["Meter",1.0]
]

Any ideas?

Thanks in advance,

David Muse
david.muse at firstworks.com


__________________________________________________
D O T E A S Y - "Join the web hosting revolution!"
             http://www.doteasy.com