Линейка Яндекса показывает что стороны именно такие и площадь должна быть большей… но координаты я с гугла брал, проверял ставил точку… все совпадает. Мне кажется GPSMapEdit тоже наврал немного… Поставил его демку, но что-то не подружились мы)
liosha, вот реализация Вашей идеи…
import math
from shapely.geometry import Polygon
Lekv = 40075700
Lmer = 19980000
lat = [55.728634, 55.729489, 55.729440, 55.728565, 55.728634]
lon = [37.826026, 37.826097, 37.827862, 37.827780, 37.826026]
nds = [];
for i in range(len(lat)):
mlat = lat[i]*Lekv/360.0
mlon = lon[i]*math.cos(math.radians(lat[i]))*Lmer/360.0
nds.append((mlat, mlon))
polygon = Polygon(nds)
print polygon.area
Результат 5310.72967046
Я уже сомневаюсь в правильности исходных координат, в реальной площади… есть ли какая-нибудь продвинутая считалка у кого-нибудь, может САПР какой?
Получить бы точные исходные данные и площадь для них заведомо правильно посчитанную…