The shapefiles and boundaries I work with are on WGS 84 referenced to SA:
GEOGCS["WGS 84", DATUM["World Geodetic System 1984", SPHEROID["WGS 84", 6378137.0, 298.257223563, AUTHORITY["EPSG","7030"]], AUTHORITY["EPSG","6326"]], PRIMEM["Greenwich", 0.0, AUTHORITY["EPSG","8901"]], UNIT["degree", 0.017453292519943295], AXIS["Geodetic longitude", EAST], AXIS["Geodetic latitude", NORTH], AUTHORITY["EPSG","4326"]]
I need to project them onto a plane to make area calculations and intersect shapes. I'm stuck with a few questions:
- local official projection is defined by Gauss Conformal Mercator. How do I project files with this map definition in Python? here are a number of definitions: http://spatialreference.org/ref/epsg/2047/, as well as http://epsg.io/?q=south%20africa.
Which package/wrapper/function do I choose to do the heavy-lifiting?
I read that planar projection in SA are done in a variety of ways- Which of the Lo/Wg do I choose for the projection to cover the entire span of the country?
The mapping needs to maintain the areas. Does tha GCM do that?
Is Gauss-Krüger the same as Gauss Conformal?