Определить, находится ли точка внутри полигона - PullRequest
Новая
       7

Определить, находится ли точка внутри полигона

0 голосов
/ 14 октября 2019

Я хотел бы иметь возможность определить, находится ли GPS-положение NMEA в пределах границы. В настоящее время я использую библиотеку geojson .net (https://github.com/GeoJSON-Net/GeoJSON.Net) для преобразования файла геоджона в многоугольник. Объект многоугольника Airspace создает список широт и долгот в airspace.coordinates [0] .coordinates. 


//This void will look for airspace files in /airspace directory.
        private static void _getAirspace()
        {
            string currentDirectory = Directory.GetCurrentDirectory();
            string path = System.IO.Path.Combine(@currentDirectory, "Airspace");
            string[] files = Directory.GetFiles(path, "*.geojson", SearchOption.AllDirectories);

            foreach (string filePath in files)
            {
                string fileContents = System.IO.File.ReadAllText(filePath);
                Polygon airspace = JsonConvert.DeserializeObject<Polygon>(fileContents);


                string json = "{\"coordinates\":[ 5.673479,52.429473 ],\"type\":\"Point\"}";
                Point point = JsonConvert.DeserializeObject<Point>(json);

                if (IsPointInPolygon(point,airspace)){

                }
                else
                {

                }

            }

        }

        public static bool IsPointInPolygon(Point p, Polygon polygon)
        {
            double minX = polygon.Coordinates[0].Coordinates[0].Longitude;
            double maxX = polygon.Coordinates[0].Coordinates[0].Longitude;
            double minY = polygon.Coordinates[0].Coordinates[0].Latitude;
            double maxY = polygon.Coordinates[0].Coordinates[0].Latitude;

            for (int i = 1; i < polygon.Coordinates.Count; i++)
            {               
                minX = Math.Min(polygon.Coordinates[0].Coordinates[i].Longitude, minX);
                maxX = Math.Max(polygon.Coordinates[0].Coordinates[i].Longitude, maxX);
                minY = Math.Min(polygon.Coordinates[0].Coordinates[i].Latitude, minY);
                maxY = Math.Max(polygon.Coordinates[0].Coordinates[i].Latitude, maxY);
            }



            bool inside = false;
            for (int i = 0, j = polygon.Coordinates.Count - 1; i < polygon.Coordinates.Count; j = i++)
            {
                if ((polygon.Coordinates[0].Coordinates[i].Latitude > p.Coordinates.Latitude) != (polygon.Coordinates[0].Coordinates[j].Latitude > p.Coordinates.Latitude) &&
                     p.Coordinates.Longitude < (polygon.Coordinates[0].Coordinates[j].Longitude - polygon.Coordinates[0].Coordinates[i].Longitude) * (p.Coordinates.Latitude - polygon.Coordinates[0].Coordinates[i].Latitude) / (polygon.Coordinates[0].Coordinates[j].Latitude - polygon.Coordinates[0].Coordinates[i].Latitude) + polygon.Coordinates[0].Coordinates[j].Longitude)
                {
                    inside = !inside;
                }
            }

            return inside;
        }

Метод IsPointInPolygon всегда возвращает true, в этом случае я ожидаю, что метод вернет true, когда координаты точки находятся внутри многоугольника.

...