Feature request: return the EPSG code corresponding to the UTM zone

[araichev]

I'd use that feature.

[rabbl]

I'd love that feature too.

[araichev]

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}"

[rabbl]

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}"

[rico-ci]

Yeah, this is discussed on several platforms (e.g. here). It would be great to have this as part of the package, so +1 on this issue.