Open araichev opened 4 years ago
I'd love that feature too.
Aha, it's easier than i initially thought:
def utm_epsg_from_latlot(lat, lon):
zone = utm.from_latlon(lat, lon)[2]
return f"326{zone:02d}"
you should make a difference between northern and southern hemisphere:
32633 for zone 33 northern hemisphere 32733 for zone 33 southern hemisphere
So the method could be something like this:
def utm_epsg_from_latlot(lat, lon):
zone = utm.from_latlon(lat, lon)[2]
return f"326{zone:02d}" if lat >= 0 else f"327{zone:02d}"
I'd use that feature.