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
[email protected]
__________________________________________________
D O T E A S Y - "Join the web hosting revolution!"
http://www.doteasy.com
_______________________________________________
gdal-dev mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/gdal-dev
