Ошибки PHP в Point-in-Polygon

Было несколько проблем с исходным кодом, закрывая многоугольники, зафиксировавшие одну из них, но код также дал неверные результаты для точек на граничных линиях многоугольника. Оператор if..else в конце функции isWithinBoundary должен выполняться только в том случае, если точка IS NOT на границе. Поскольку точка на границе фактически не пересечет границу, то граф пересечений всегда будет нечетным для граничной точки, что означает, что этот окончательный оператор 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; } 

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