Алгоритм Point in Polygon дает неверные результаты для отрицательных точек

Я пытаюсь проверить, находится ли широта, долгота в многоугольнике или нет. Вот мой массив:

$vertices_x : 
Array
(
    [0] => -32.581189
    [1] => -38.785885
    [2] => -39.26384
    [3] => -34.919383
    [4] => -32.284464
)

$vertices_y:
Array
(
    [0] => 170.643905
    [1] => 170.424179
    [2] => -178.15004
    [3] => -176.524063
    [4] => -178.325821
)

$longitude_x : 173.5385
$latitude_y : -34.472
$points_polygon = count($vertices_x) - 1;

Я использую функцию ниже для проверки:

 function is_in_polygon($points_polygon, $vertices_x, $vertices_y, $longitude_x, $latitude_y) {
        $i = $j = $c = 0;
        for ($i = 0, $j = $points_polygon; $i < $points_polygon; $j = $i++) {
            if ((($vertices_y[$i] > $latitude_y != ($vertices_y[$j] > $latitude_y)) &&
                    ($longitude_x < ($vertices_x[$j] - $vertices_x[$i]) * ($latitude_y - $vertices_y[$i]) / ($vertices_y[$j] - $vertices_y[$i]) + $vertices_x[$i])))
                $c = !$c;
        }
        return $c;
    }

И эта функция всегда дает мне 0 (не в полигоне), но если вы проверите, то моя точка $longitude_x : 173.5385 , $latitude_y : -34.472 находится в этой области полигона.

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


php polygon point-in-polygon
person Community    schedule 19.04.2020    source источник
comment
Как вы можете видеть здесь, этот алгоритм имеет проблемы с отрицательными значениями. Альтернативы и решения доступны в этой теме   -  person Claudio    schedule 19.04.2020
comment
эта ветка тоже бесполезна   -  person    schedule 19.04.2020

Ответы (1)


arrow_upward
2
arrow_downward

Я не думаю, что ваша точка (173,5385,-34,472) находится в многоугольнике.

Ваше значение x намного больше, чем самая большая вершина x. Вы можете увидеть это, если вы построите это. Вы можете видеть по оранжевой точке, что даже если вы просто перепутали широту и долготу, она все равно не в полигоне.

введите здесь описание изображения

РЕДАКТИРОВАТЬ:

При первом осмотре вашего нового многоугольника и точки кажется, что этот, скорее всего, содержится внутри многоугольника (x_min ‹ x ‹ x_max и y_min ‹ y ‹ y_max).

-36.236432, 176.467563
-37.936530, 172.688266
-39.801068, 177.895786
-35.345287,-177.446011
-34.625208,-177.907437

точка:

(-37.0278,176.6158)

Тем не менее, построение этого графика снова показывает, что точка находится за пределами многоугольника:

введите здесь описание изображения

...и увеличение...

введите здесь описание изображения

Я построил их на питоне с помощью matplotlib и рекомендую вам делать что-то подобное при отладке подобных вещей. Если вы хотите использовать php, html и т. д., вы можете использовать многоугольники svg. или html canvas.

person compuphys    schedule 21.04.2020
comment
А как насчет этого многоугольника: -36.236432,176.467563|-37.93653,172.688266|-39.801068,177.895786|-35.345287,-177.446011|-34.625208,-177.907437 и точка -37.0278,176.6158 - person ; 22.04.2020
comment
Спасибо за ваш ответ, но если вы заглянете в карты Google, то эта точка будет показана в новом многоугольнике. Я взял этот полигон только с карт Google. - person ; 22.04.2020
comment
Ах, мы говорим о широте/долготе на самом деле в трехмерном смысле ... Я думал, что это просто имена ваших переменных ... это совсем другая игра с мячом. Помогает ли этот пост: gis.stackexchange.com/a/46720? - person compuphys; 22.04.2020