I'm using Shapely to do some GIS analysis in Python. I have a
Polygon and a
Point, both in EPSG 3857 (web Mercator). I want to check whether the point is within a given distance (in metres) of the polygon. The data is guaranteed to be within the UK.
I can see how to do this in Shapely terms by buffering the polygon by the distance and then checking whether the result contains the point. But I don't think I can do this in EPSG 3857 because the scale is different in different directions (in the UK, one degree of latitude is getting on for twice as large as one degree of longitude). I need to first do a coordinate transform into something that has the same scale in every direction, at least locally.
The sort of distances I'm doing this over are less than 1km, and the polygons I start with are considerably smaller than that. Approximations are okay; I'd be happy with wrong answers up to 5% of the distance I'm checking for.
What is the projection property that I'm looking for? Can someone suggest a suitable projection please?