Find Point in polygon PHP

Question

i have a typical question with the Geometric datatype of mysql, polygon.

I have the polygon data, in the form of an array of latitudes and longitudes, ex:

Answer 1

Above solution is not working as i expect, instead of using the above solution you can prefer below solutions

1. With PHP

   function pointInPolygon($point, $polygon, $pointOnVertex = true) {
       $this->pointOnVertex = $pointOnVertex;
   
       // Transform string coordinates into arrays with x and y values
       $point = $this->pointStringToCoordinates($point);
       $vertices = array(); 
       foreach ($polygon as $vertex) {
           $vertices[] = $this->pointStringToCoordinates($vertex); 
       }
   
       // Check if the lat lng sits exactly on a vertex
       if ($this->pointOnVertex == true and $this->pointOnVertex($point, $vertices) == true) {
           return "vertex";
       }
   
       // Check if the lat lng 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'])) { // Check if point is on an horizontal polygon boundary
               return "boundary";
           }
           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']) { // Check if lat lng is on the polygon boundary (other than horizontal)
                   return "boundary";
               }
               if ($vertex1['x'] == $vertex2['x'] || $point['x'] <= $xinters) {
                   $intersections++; 
               }
           } 
       } 
       // If the number of edges we passed through is odd, then it's in the polygon. 
       if ($intersections % 2 != 0) {
           return "inside";
       } else {
           return "outside";
       }
   }
   
   function pointOnVertex($point, $vertices) {
     foreach($vertices as $vertex) {
         if ($point == $vertex) {
             return true;
         }
     }
   
   }
   
   function pointStringToCoordinates($pointString) {
       $coordinates = explode(" ", $pointString);
       return array("x" => $coordinates[0], "y" => $coordinates[1]);
   }
   // Function to check lat lng
   function check(){
       $points = array("22.367582 70.711816", "21.43567582 72.5811816","22.367582117085913 70.71181669186944","22.275334996986643 70.88614147123701","22.36934302329968 70.77627818998701"); // Array of latlng which you want to find
       $polygon = array(
           "22.367582117085913 70.71181669186944",
           "22.225161442616514 70.65582486840117",
           "22.20736264867434 70.83229276390898",
           "22.18701840565626 70.9867880031668",
           "22.22452581029355 71.0918447658621",
           "22.382709129816103 70.98884793969023",
           "22.40112042636022 70.94078275414336",
           "22.411912121843205 70.7849142238699",
           "22.367582117085913 70.71181669186944"
       );
       // The last lat lng must be the same as the first one's, to "close the loop"
       foreach($points as $key => $point) {
           echo "(Lat Lng) " . ($key+1) . " ($point): " . $this->pointInPolygon($point, $polygon) . "
   ";
       }
   }
   

2. With MySql

CREATE TABLE `TestPoly` (
   `id` int(11) NOT NULL,
   `name` varchar(255) NOT NULL,
   `pol` polygon NOT NULL
 )

SET @g = 'POLYGON((22.367582117085913 70.71181669186944, 22.225161442616514 70.65582486840117, 22.20736264867434 70.83229276390898, 22.18701840565626 70.9867880031668, 22.22452581029355 71.0918447658621, 22.382709129816103 70.98884793969023, 22.40112042636022 70.94078275414336, 22.411912121843205 70.7849142238699, 22.367582117085913 70.71181669186944))';
INSERT INTO TestPoly (pol) VALUES (ST_GeomFromText(@g))

set @p = GeomFromText('POINT(22.4053386588057 70.86240663480157)');
select * FROM TestPoly where ST_Contains(pol, @p);