I am trying to convert some coordinates in my shape file to GPS coordinates. Here is the prj file:
PROJCS["NAD_1983_StatePlane_North_Carolina_FIPS_3200_Feet", 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["Lambert_Conformal_Conic"], PARAMETER["False_Easting",2000000.002616666], PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",-79.0], PARAMETER["Standard_Parallel_1",34.33333333333334], PARAMETER["Standard_Parallel_2",36.16666666666666], PARAMETER["Latitude_Of_Origin",33.75],UNIT["Foot_US",0.3048006096012192]] Based on the PROJCS line, it looks like it would fit this projection: http://www.spatialreference.org/ref/epsg/2264/
EDIT: Looking through my *.prj file further, this looks to be a better match (http://spatialreference.org/ref/esri/102719/) but still isn't working when I plug in init='esri:102719
Here is my code to convert into lat/long coordinates:
in_proj = Proj(init='epsg:2264') out_proj = Proj(init='epsg:4326') new_points = [] for point in shapes[0].points: print point new_point = transform(in_proj, out_proj, point[0], point[1]) print new_point I'm not getting what I'm expecting (~36,~-78). Instead I'm getting (~-62,~39). Any idea what I could be doing wrong? Thanks so much.
أكثر...
PROJCS["NAD_1983_StatePlane_North_Carolina_FIPS_3200_Feet", 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["Lambert_Conformal_Conic"], PARAMETER["False_Easting",2000000.002616666], PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",-79.0], PARAMETER["Standard_Parallel_1",34.33333333333334], PARAMETER["Standard_Parallel_2",36.16666666666666], PARAMETER["Latitude_Of_Origin",33.75],UNIT["Foot_US",0.3048006096012192]] Based on the PROJCS line, it looks like it would fit this projection: http://www.spatialreference.org/ref/epsg/2264/
EDIT: Looking through my *.prj file further, this looks to be a better match (http://spatialreference.org/ref/esri/102719/) but still isn't working when I plug in init='esri:102719
Here is my code to convert into lat/long coordinates:
in_proj = Proj(init='epsg:2264') out_proj = Proj(init='epsg:4326') new_points = [] for point in shapes[0].points: print point new_point = transform(in_proj, out_proj, point[0], point[1]) print new_point I'm not getting what I'm expecting (~36,~-78). Instead I'm getting (~-62,~39). Any idea what I could be doing wrong? Thanks so much.
أكثر...