Как вычислить площадь многоугольника на земной поверхности с помощью Python?

Название в основном говорит само за себя. Мне нужно вычислить площадь внутри многоугольника на поверхности Земли с помощью Python. "ГИС" аромат, то нужно выбрать единицу измерения для вашего района и найти подходящую проекцию, которая сохраняет территорию (не все делают). Поскольку ты говорят о вычислении произвольный многоугольник, я бы использовал что-то вроде азимутального ламберта Равноплощадочная проекция. Установить начало / центр проекции, которая будет центр вашего многоугольника, проект многоугольник в новую координату системе, затем рассчитайте площадь, используя стандартные планарные методы.

Итак, как мне сделать это в Python?