Площадь по геокоординатам

[Vago]

Сообщение от Merkator

Нет, они уже даны в порядке обхода. А что?

А, ну раз так, то всё нормально.

[Vago]

Сообщение от Merkator

Vago, мне кажется, что Z - не просто выпуклость Земли. Насколько я понял систему ECEF, то просто из центра Земли определенным образом отложена декартова система координат.

Да. Так оно и есть. Центр координат - где-то в центре Земли, OX - к точке пересечения Гринвичского меридиана и экватора, OY - к точке 0 СШ; 90 ВД; OZ - в сторону Северного полюса.

Сообщение от Merkator

И если я откажусь от Z, то скажем не в Финском заливе, а где-нибудь в Атлантике программа перестанет работать. (координата Z будет сильно меняться, а какая-нибудь Y наоборот не будет)

Тоже, в принципе, верно. Больше всего - на полюсах.

Так что, самую грамотную идею тут подал уважаемый alexBlack

Добавляем

def MoveBase( data ):
	new_data = []
	deX = data[0][0]
	deY = data[0][1]
	deZ = data[0][2]
	for a in data:
		new_data.append( (a[0]-deX, a[1]-deY, a[2]-deZ) )

	return new_data

Вызовы:

dataECEF = Geo2ECEF( data )
print "Data in ECEF format:", dataECEF
ndata = MoveBase( dataECEF )

(дальше площадь - как у Вас)

Для

data = [
  (0., 0.),
  (1., 0.),
  (1., 1.),
  (1., 0.),
]

("квадрат" на экваторе со стороной примерной 111км) получаем площадь 12306.9 км^2, что похоже на правду

[Merkator]

Спасибо огромное!
Вы все мне очень помогли. Отзывы добавил.