Ошибки PHP "точка в многоугольнике"

В исходном коде было несколько проблем, закрытие многоугольников устранило одну из них, но код также дал неверные результаты для точек на граничных линиях многоугольника. Оператор if..else в конце функции isWithinBoundary должен выполняться только в том случае, если точка НЕ ​​НАХОДИТСЯ на границе. Поскольку точка на границе фактически не пересекает границу, то количество пересечений всегда будет нечетным для граничной точки, что означает, что этот последний оператор IF всегда будет возвращать FALSE для граничной точки.

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

<?php
    $myPolygon = array('4,3', '4,6', '7,6', '7,3','4,3');

    $test_points = array('0,0','1,1','2,2','3,3','3.99999,3.99999','4,4','5,5','6,6','6.99999,5.99999','7,7');
    echo "The test polygon has the co-ordinates ";
    foreach ($myPolygon as $polypoint){
        echo $polypoint.", ";
    }
    echo "<br/>"; 
    foreach ($test_points as $apoint)
    {
        echo "Point ".$apoint." is ";
        if (!isWithinBoundary($apoint,$myPolygon))
        {
            echo " NOT ";
        }
        echo "inside the test polygon<br />";
    }

    function pointStringToCoordinates($pointString) 
    {
            $coordinates = explode(",", $pointString);
            return array("x" => trim($coordinates[0]), "y" => trim($coordinates[1]));
    }

    function isWithinBoundary($point,$polygon)
    {
        $result =FALSE;
        $point = pointStringToCoordinates($point);
        $vertices = array();
        foreach ($polygon as $vertex) 
        {
            $vertices[] = pointStringToCoordinates($vertex); 
        }
        // Check if the point is inside the polygon or on the boundary
        $intersections = 0; 
        $vertices_count = count($vertices);
        for ($i=1; $i < $vertices_count; $i++) 
        {
            $vertex1 = $vertices[$i-1]; 
            $vertex2 = $vertices[$i];
            if ($vertex1['y'] == $vertex2['y'] and $vertex1['y'] == $point['y'] and $point['x'] > min($vertex1['x'], $vertex2['x']) and $point['x'] < max($vertex1['x'], $vertex2['x'])) 
            { 
                // This point is on an horizontal polygon boundary
                $result = TRUE;
                // set $i = $vertices_count so that loop exits as we have a boundary point
                $i = $vertices_count;
            }
            if ($point['y'] > min($vertex1['y'], $vertex2['y']) and $point['y'] <= max($vertex1['y'], $vertex2['y']) and $point['x'] <= max($vertex1['x'], $vertex2['x']) and $vertex1['y'] != $vertex2['y']) 
            { 
                $xinters = ($point['y'] - $vertex1['y']) * ($vertex2['x'] - $vertex1['x']) / ($vertex2['y'] - $vertex1['y']) + $vertex1['x']; 
                if ($xinters == $point['x']) 
                { // This point is on the polygon boundary (other than horizontal)
                    $result = TRUE;
                    // set $i = $vertices_count so that loop exits as we have a boundary point
                    $i = $vertices_count;
                }
                if ($vertex1['x'] == $vertex2['x'] || $point['x'] <= $xinters) 
                {
                    $intersections++; 
                }
            } 
        }
        // If the number of edges we passed through is even, then it's in the polygon. 
        // Have to check here also to make sure that we haven't already determined that a point is on a boundary line
        if ($intersections % 2 != 0 && $result == FALSE) 
        {
            $result = TRUE;
        } 
        return $result;
    }
?>

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

Что ж, я снова обнаруживаю, что глупо отвечаю на свой собственный глупый вопрос.

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

Итак, это -

$central_park = array('40.768109,-73.981885', '40.800636,-73.958067', '40.796900,-73.949184', '40.764307,-73.972959');

Должно быть это -

$central_park = array('40.768109,-73.981885', '40.800636,-73.958067', '40.796900,-73.949184', '40.764307,-73.972959', '40.764307,-73.972959');

Вот как я сегодня был глуп. Спасибо.

Проблема с вашим кодом заключается в том, что переменная $ result перезаписывается этим кодом.

if ($intersections % 2 != 0) {
    $result = TRUE;
} else {
    $result = FALSE;
}

даже если здесь $ result == TRUE:

if ($xinters == $point['x']) {
    $result = TRUE;
}

В исходном коде было «return», которое было правильным, а не ошибочным.