9

我有一个 shapefile 定义覆盖城市的形状。我还有一组由纬度和经度定义的坐标。

我需要能够确定这些点是否在 shapefile 中的任何形状内。

目前我正在尝试使用 easygis.net 来解决这个问题，但它似乎不起作用。

我相信 shapefile 中的坐标在 UTM 中，但有一个偏移北，我不知道如何更正它，将其转换为纬度/经度，或转换我的纬度/经度对以匹配。

我一直在查看其他库，例如 dotspatial 和 sharpmap，但我没有看到对这个问题的直观答案。

同时，我确信这是一个已经解决的问题。有没有图书馆可以轻松做到这一点？或者我如何将地图的偏移 UTM 转换为纬度/经度，或者我的纬度/经度点到这个偏移 UTM？

4

3 回答 3

8

下面是一些示例代码，介绍如何使用 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;
    }
于 2013-07-17T05:01:43.707 回答
4

或者我如何将地图的偏移 UTM 转换为纬度/经度，或者我的纬度/经度点到这个偏移 UTM？

这需要重新投影您的点（或 shapefile 的数据）。这可以通过Proj4Net完成（除非您的 GIS 已经支持它）。

一旦完成，它就是多边形测试中的一个点。这有很多选择，虽然我经常使用缠绕数法。

于 2013-05-31T17:36:39.530 回答
0

如果您只需要将投影文件转换为地理坐标，那么 Qgis 是我认为最容易使用的工具，它是一款基于 GRASS 的轻量级免费开放的 GIS 软件。使用直接的投影工具。 http://www.qgis.org/

于 2013-06-07T05:19:21.207 回答