Как определить, находится ли точка в форме в шейп-файле?

Question

У меня есть shapefile, определяющий формы для покрытия города. У меня также есть набор координат, определяемый широтой и долготой.

Мне нужно определить, находятся ли эти точки в пределах любой формы в шейп-файле.

В настоящее время я пытаюсь использовать easygis.net, чтобы понять это, но он, похоже, не работает.

Я считаю, что координаты в шейп-файле находятся в UTM, но имеют смещение на север, и я не знаю, как исправить это, перевести его на lat/long или преобразовать мои пары lat/long для соответствия.

Я смотрел на другие библиотеки, такие как dotspatial и sharpmap, но я не вижу интуитивного ответа на проблему.

В то же время, я уверен, что это проблема, которая была решена. Есть ли библиотеки, которые делают это легко? Или как мне преобразовать сдвиг UTM карты в lat/long, или мой lat/long указывает на это смещение UTM?

Answer 1

Или как преобразовать смещение UTM карты в lat/long, или мои lat/long указывает на это смещение UTM?

Для этого требуется перепрофилирование ваших точек (или данных шейп-файла). Это можно сделать через Proj4Net (если ваша ГИС уже не поддерживает его).

Как только это будет сделано, это точка в полигоне. Есть many options for this, хотя я часто использую winding number method.

Answer 2

Если вам нужно только преобразовать проецируемый файл в географическую координату, то Qgis - это самый простой инструмент, который, на мой взгляд, будет бесплатным и открытым. С прямыми проекционными инструментами. http://www.qgis.org/

Answer 3

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

public bool PointInShape() { 
     // Load a shapefile. If the shapefile is already using your custom projection, we don't need to change it. 
     Shapefile wierdShapefile = Shapefile.OpenFile("C:\\MyShapefile.shp"); 

     // Note, if your shapefile with custom projection has a .prj file, then we don't need to mess with defining the projection. 
     // If not, we can define the projection as follows: 

     // First get a ProjectionInfo class for the normal UTM projection 
     ProjectionInfo pInfo = DotSpatial.Projections.KnownCoordinateSystems.Projected.UtmNad1983.NAD1983UTMZone10N; 

     // Next modify the pINfo with your custom False Northing 
     pInfo.FalseNorthing = 400000; 

     wierdShapefile.Projection = pInfo; 

     // Reproject the strange shapefile so that it is in latitude/longitude coordinates 
     wierdShapefile.Reproject(DotSpatial.Projections.KnownCoordinateSystems.Geographic.World.WGS1984); 

     // Define the WGS84 Lat Lon point to test 
     Coordinate test = new Coordinate(-120, 40); 

     foreach (Feature f in wierdShapefile.Features) { 
      Polygon pg = f.BasicGeometry as Polygon; 
      if (pg != null) 
      { 
       if (pg.Contains(new Point(test))) 
       { 
        // If the point is inside one of the polygon features 
        return true; 
       } 
      } 
      else { 
       // If you have a multi-part polygon then this should also handle holes I think 
       MultiPolygon polygons = f.BasicGeometry as MultiPolygon; 
       if (polygons.Contains(new Point(test))) { 
        return true; 
       } 
      } 
     } 

     return false; 
    }