0

我想将距离的笛卡尔坐标“转换”为米。

我添加了

我尝试使用以下代码：

        return new PagedList<Dto.DrugConsortium.CollectionSiteWithCoordinatesDto>(query.Select(p =>
            new Dto.DrugConsortium.CollectionSiteWithCoordinatesDto { CollectionSite = p, Distance = p.Location.ProjectTo(2855).Distance(myLocation) }).ToList(), pageIndex, pageSize);

我的分机在哪里：

public static class GeometryExtensions
{
    static readonly IGeometryServices _geometryServices = NtsGeometryServices.Instance;
    static readonly ICoordinateSystemServices _coordinateSystemServices
        = new CoordinateSystemServices(
            new CoordinateSystemFactory(),
            new CoordinateTransformationFactory(),
            new Dictionary<int, string>
            {
            // Coordinate systems:

            // (3857 and 4326 included automatically)

            // This coordinate system covers the area of our data.
            // Different data requires a different coordinate system.
            [2855] =
                @"
                PROJCS[""NAD83(HARN) / Washington North"",
                    GEOGCS[""NAD83(HARN)"",
                        DATUM[""NAD83_High_Accuracy_Regional_Network"",
                            SPHEROID[""GRS 1980"",6378137,298.257222101,
                                AUTHORITY[""EPSG"",""7019""]],
                            AUTHORITY[""EPSG"",""6152""]],
                        PRIMEM[""Greenwich"",0,
                            AUTHORITY[""EPSG"",""8901""]],
                        UNIT[""degree"",0.01745329251994328,
                            AUTHORITY[""EPSG"",""9122""]],
                        AUTHORITY[""EPSG"",""4152""]],
                    PROJECTION[""Lambert_Conformal_Conic_2SP""],
                    PARAMETER[""standard_parallel_1"",48.73333333333333],
                    PARAMETER[""standard_parallel_2"",47.5],
                    PARAMETER[""latitude_of_origin"",47],
                    PARAMETER[""central_meridian"",-120.8333333333333],
                    PARAMETER[""false_easting"",500000],
                    PARAMETER[""false_northing"",0],
                    UNIT[""metre"",1,
                        AUTHORITY[""EPSG"",""9001""]],
                    AUTHORITY[""EPSG"",""2855""]]
            "
            });

    public static IGeometry ProjectTo(this IGeometry geometry, int srid)
    {
        var geometryFactory = _geometryServices.CreateGeometryFactory(srid);
        var transformation = _coordinateSystemServices.CreateTransformation(geometry.SRID, srid);

        return GeometryTransform.TransformGeometry(
            geometryFactory,
            geometry,
            transformation.MathTransform);
    }
}

ProjNet4GeoAPI添加了 nuget 包 (v 1.4.1)

但我无法解决GeometryTransform类，它没有找到。

如何解决？

4

1 回答 1

0

First GeometryTransform.TransformGeometry() 是位于 NetTopologySuite.CoordinateSystems.Transformations 命名空间下的方法。你有

“使用 NetTopologySuite.CoordinateSystems.Transformations；”

在您的扩展方法的顶部？

其次，我要补充一点，2855 是西北华盛顿州的坐标系。我会检查您使用的坐标是否在该区域内。

此外，每个坐标系都有自己的测量单位。例如，SRID:4326 对于几何（笛卡尔坐标）以度为单位，对于地理类型以米为单位。

SRID:2855 以英尺为单位

您可以查看此站点以获取可用坐标系的列表。 https://epsg.io/

我是地理类型和转换的新手，但在我的搜索中这是我发现的。

于 2020-09-21T14:46:21.170 回答