Расчет долготы центра над 180-м меридианом

Question

Расчет долготы центра над 180-м меридианом

Я вычисляю простое среднее значение пар широта / долгота, чтобы получить центральную точку (проецируемую на плоскую плоскость), и я изо всех сил пытаюсь разработать алгоритм для правильного учета 180-го меридиана. Я разрабатываю это на PHP.

Я вычисляю простое среднее значение, складываю все долготы и делю на количество. Однако в ситуации, когда есть, например, точки -175 и 175, он возвращает среднее значение 0, тогда как мне нужно значение 180.

В одном из подходов я попытался добавить 360 ко всем отрицательным значениям, взять среднее и затем вычесть 360. Это устраняет проблему с 180-м меридианом, но создает ту же проблему с основным меридианом (0 меридиан).

$neg_lng_flag = false;
//calculate centerpoint
$i = 0;
while($i < $vertice_count) {
$lat_sum += $vertice_obj[$i]->lat;
//if lng is negative, add 360
if($vertice_obj[$i]->lng < 0) {
$lng_sum += $vertice_obj[$i]->lng + 360;
$neg_lng_flag = true;
} else {
$lng_sum += $vertice_obj[$i]->lng;
}
$i++;
}
$avg_lat = round($lat_sum / $vertice_count,2);
$avg_lng = round($lng_sum / $vertice_count,2);

if($neg_lng_flag) {
$avg_lng = $avg_lng - 360;
if($avg_lng < -180) {
$avg_lng = $avg_lng + 180;
}
}

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

0

geography php

Решение

Другие решения

Других решений пока нет …

Источник

Accepted Answer

Я смог найти надежное решение для моих целей. В моем случае мне просто нужно было вычислить минимальную и максимальную координаты долготы, и если разница превысила 180, я знал, что нужно сдвинуть отрицательные значения на +360, а затем вычислить среднее значение. Если среднее было больше 180, мне нужно было вычесть 360, а если меньше -180, мне нужно было добавить 360.

//COMPUTE LONGITUDAL MIDPOINT
function getLngMidpoint($vertice_obj) {
//get min and max vals
$i = 0;
while($i < count($vertice_obj)) {
$min = $vertice_obj[$i]->lng;
$max = $vertice_obj[$i]->lng;

if($vertice_obj[i]->lng > $max) {
$max = $vertice_obj[$i]->lng;
}
if($vertice_obj[$i]->lng < $min) {
$min = $vertice_obj[$i]->lng;
}

$i += 1;
}

$shift = 0;
//check if distance between min and max > 180. If so, need to shift
if(($max - $min) > 180) {
//shift all lng by 180
$shift = 360;
}

$i = 0;
$sum_lng = 0;
while($i < count($vertice_obj)) {
if($vertice_obj[$i] < 0) {
$sum_lng += $vertice_obj[$i]->lng + $shift;
} else {
$sum_lng += $vertice_obj[$i]->lng;
}
$i += 1;
}

$avg_lng = $sum_lng / count($vertice_obj);
if($avg_lng > 180) {
$avg_lng = $avg_lng - 360;
}
if($avg_lng < -180) {
$avg_lng = $avg_lng + 360;
}

return $avg_lng;
}

0