I have shapefile (from the local government gis office of a town in Massachusetts) that contains, in the corresponding
.shp.xml file, the following WKT projection:
PROJCS["NAD_1983_StatePlane_Massachusetts_Mainland_FIPS_2001", 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",200000.0], PARAMETER["False_Northing",750000.0], PARAMETER["Central_Meridian",-71.5], PARAMETER["Standard_Parallel_1",41.71666666666667], PARAMETER["Standard_Parallel_2",42.68333333333333], PARAMETER["Latitude_Of_Origin",41.0], UNIT["Meter",1.0]]
Loading the shapefile into QGIS results in a null CRS ("This layer appears to have no projection specification..."). I've tried converting the above projection into proj4 syntax, but I'm running into some problems caused by the fact that I have no idea what I am doing.
The above projection is awfully close to the standard projection used by MassGIS (
epsg:26986), although "Standard_Parallel_1" and "Standard_Parallel_2" are swapped. Trying to convert this to
proj4 syntax using the Python
osgeo.osr module results in:
ERROR 6: No translation for Lambert_Conformal_Conic to PROJ.4 format is known.
I can replace
Lambert_Conformal_Conic_2SP (as suggested elsewhere), which results in a successful export in
proj4 syntax, but (a) it's not clear to me that's the right thing to do, and in any case (b) the resulting projection still doesn't place the features where they belong.
How do I get this data in the right place?
Edit: I have tried just using
epsg:29686 as well, which didn't work either. The shapefile in question is here.