查找 Lat Long position 是否在 esri shapefile 中

Question

我正在尝试在一个简单的控制台应用程序中使用 Net Topology Suite (v 1.13.2) 来查找包含道路信息的形状文件中的点。

我已经加载了形状文件并将数据存储在一个类型化的列表中，并且我检查了数据是否符合预期..

当我遍历列表时，我想查找给定点是否在形状文件中的任何形状内。

听起来很简单！

我的代码看起来像；

    private static Feature FindPoint(double lat, double lon)
    {
        Coordinate c = new Coordinate(lat, lon);

        IGeometry g = factory.CreateGeometry(Geometry.DefaultFactory.CreatePoint(c));

        foreach(Feature f in Features)
        {
            if (f.Geometry.Overlaps(g))
                return f;
            if (f.Geometry.EnvelopeInternal.Contains(c))
                return f;
            if (f.Geometry.Boundary.Contains(g))
                return f;
            if (f.Geometry.Contains(g))
                return f;
        }
        return null;
    }

这些陈述或我尝试过的任何其他陈述都没有返回任何指示该点位于任何形状内的指示。

我正在尝试从文件中的一个形状中挑选的经纬度。 所以它应该在那里。

关于我哪里出错的任何想法？

@Habib，我已经检查过了，几何类型是 LineString，但是信封被定义为多边形。

我查看了您发送的链接，下面是新代码；

DistanceOp dop = new DistanceOp(f.Geometry,g);
var np = dop.NearestPoints();
var d = dop.Distance();

尝试按照那里的指示进行操作，但我得到的关于我选择的点与线之间的距离的答案是 71.236662957979718。

71.236662957979718 什么？ 厘米，米，度？？？

无论如何，很奇怪，因为我在其中一条线上选择了一个点，所以我希望答案为 0。

Answer 1

在创建坐标时将经度放在纬度之前让它对我有用。

更正的代码：

private static Feature FindPoint(double lat, double lon)
{
    Coordinate c = new Coordinate(lon, lat);
    Geometry g = Geometry.DefaultFactory.CreatePoint(c);

    foreach(Feature f in Features)
    {
        if (f.Geometry.Contains(g))
            return f;
    }
    return null;
}

我正在使用 NetTopologySuite.IO.Esri 在这里找到： https://github.com/NetTopologySuite/NetTopologySuite.IO 。