projecting onto Cartographic coordinates using Gauss Conformal Mercator in python

Question

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:

1. 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?

2. The mapping needs to maintain the areas. Does tha GCM do that?

3. Is Gauss-Krüger the same as Gauss Conformal?

Answer 1

Depending on the context, the Gauss-Krüger and Gauss Conformal are the same or very similar projections (see Gauss–Krüger coordinate system on wiki.

Relevant to your question, there are two important mutually exclusive classes of projections, each having distinctive geometric qualities:

• Conformal (or orthomorphic) projections preserve angular relationships between great circles.
• Equal area (or equivalent) projections preserve areal relationships.

Conformal projections definitely do not maintain correct areal scales.

Perhaps for South Africa, choose an Albers Equal Area Conic projection (see, generally wiki/Albers_projection, and specifically, albers-equal-area-conic-south-africa on spatialreference.org

If you do still wish to use the (conformal) Transverse Mercator (EPSG projection 2047) -- or whichever projection you decide you need -- using Python, perhaps you should ask that separately.