I want to run various algorithms in Shapely on some geographic data that I have. It is acceptable to have some measure of error, as long as it is "reasonable". e.g. find closest pair.
I would like to project the points from latitude-longitude to cartesian space such that
euclidan_distance(P(p0), P(p1)) ~= geodesic_distance(p0, p1)
I was thinking of using PyProj but it may be overkill, and it doesn't seem the easiest thing to do.
y = latitude * 110574 x = longitude * 11320 * cos(radians(latitude))
Results of naive approach:
I put in a set of coordinates nearby, and got relative errors of 0.5%-15.3%. 15% seems excessive.