Определение зоны UTM (для преобразования) из долготы / широты

Я пишу программу, которая ожидает количество точек широты / долготы, и конвертирую их внутренне в UTM, чтобы производить некоторые вычисления в метрах.

Диапазон самих точек широты / долготы довольно невелик - около 200 м x 200 м.Можно полагаться, что они почти всегда находятся в пределах одной зоны UTM (если только вам не повезет и вы не пересечете границу зоны).

Однако зона, в которой находятся широта и долгота, не ограничена. Однажды программа может быть запущена для людей в Австралии (и ох, сколько зон находится даже в одном штате, и сколько боли это уже причинило мне ...), а в другой день - для людей в Мексике.

У меня вопрос - есть ли способ определить, в какой зоне находится конкретная долгота / широта, чтобы ее можно было передать в библиотеку преобразования (в настоящее время я использую proj4, а также пакет R rgdal ).

Мой язык - R, но ответа не должно быть - может быть, это просто простое вычисление, или, может быть, я могу встроить системный вызов в исполняемый файл proj .

ура.