Найти точку в polygon PHP

у меня есть типичный вопрос с геометрическим типом данных mysql, polygon.

у меня есть полигональные данные, в виде массива широт и долгот, например:

[["x":37.628134,  "y":-77.458334],
["x":37.629867,   "y":-77.449021],
["x":37.62324,    "y":-77.445416],
["x":37.622424,   "y":-77.457819]]

и у меня есть точка (вершина) с координатами широты и долготы, например:

$location = new vertex($_GET["longitude"], $_GET["latitude"]);

теперь я хочу найти, находится ли эта вершина (точка) внутри многоугольника. Как я могу сделать это в php ?

18
mysql php point-in-polygon polygon

4 ответов

это функция, которую я преобразовал с другого языка в PHP:

$vertices_x = array(37.628134, 37.629867, 37.62324, 37.622424);    // x-coordinates of the vertices of the polygon
$vertices_y = array(-77.458334,-77.449021,-77.445416,-77.457819); // y-coordinates of the vertices of the polygon
$points_polygon = count($vertices_x) - 1;  // number vertices - zero-based array
$longitude_x = $_GET["longitude"];  // x-coordinate of the point to test
$latitude_y = $_GET["latitude"];    // y-coordinate of the point to test

if (is_in_polygon($points_polygon, $vertices_x, $vertices_y, $longitude_x, $latitude_y)){
  echo "Is in polygon!";
}
else echo "Is not in polygon";


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;
}

дополнительные: Для дополнительных функций я советую вам использовать polygon.класс php здесь. Создайте класс, используя свои вершины, и вызовите функцию isInside С вашей testpoint в качестве входных данных, чтобы иметь другую функцию, решающую вашу проблему.

37

популярный ответ выше содержит опечатки. В других местах этот код был очищен. Исправленный код выглядит следующим образом:

<?php
/**
  From: http://www.daniweb.com/web-development/php/threads/366489
  Also see http://en.wikipedia.org/wiki/Point_in_polygon
*/
$vertices_x = array(37.628134, 37.629867, 37.62324, 37.622424); // x-coordinates of the vertices of the polygon
$vertices_y = array(-77.458334,-77.449021,-77.445416,-77.457819); // y-coordinates of the vertices of the polygon
$points_polygon = count($vertices_x); // number vertices
$longitude_x = $_GET["longitude"]; // x-coordinate of the point to test
$latitude_y = $_GET["latitude"]; // y-coordinate of the point to test
//// For testing.  This point lies inside the test polygon.
// $longitude_x = 37.62850;
// $latitude_y = -77.4499;

if (is_in_polygon($points_polygon, $vertices_x, $vertices_y, $longitude_x, $latitude_y)){
  echo "Is in polygon!";
}
else echo "Is not in polygon";


function is_in_polygon($points_polygon, $vertices_x, $vertices_y, $longitude_x, $latitude_y)
{
  $i = $j = $c = 0;
  for ($i = 0, $j = $points_polygon-1 ; $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;
}
?>
11

вот возможный алгоритм.

  1. определите новую систему координат с интересующей вас точкой в центре.
  2. в новой системе координат преобразуйте все вершины многоугольника в полярные координаты.
  3. пересеките многоугольник, отслеживая чистое изменение угла, ∆ θ. Всегда используйте наименьшее возможное значение для каждого изменения угла.
  4. Если, как только вы пересекли многоугольник, ваш общий ∆ θ равен 0, то вы находитесь вне полигон. С другой стороны, если это ±2π, то вы внутри.
  5. Если, случайно ∆ θ > 2π или ∆ θ

написание кода остается в качестве упражнения. :)

2

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

вот принятый ответ, измененный для работы с самозакрывающимися полигонами:

$vertices_x = array(37.628134, 37.629867, 37.62324, 37.622424);    // x-coordinates of the vertices of the polygon
$vertices_y = array(-77.458334,-77.449021,-77.445416,-77.457819); // y-coordinates of the vertices of the polygon
$points_polygon = count($vertices_x);  // number vertices = number of points in a self-closing polygon
$longitude_x = $_GET["longitude"];  // x-coordinate of the point to test
$latitude_y = $_GET["latitude"];    // y-coordinate of the point to test

if (is_in_polygon($points_polygon, $vertices_x, $vertices_y, $longitude_x, $latitude_y)){
  echo "Is in polygon!";
}
else echo "Is not in polygon";


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

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

2