Рассчитайте географическое расстояние между двумя точками, используя NetTopologySuite с EF Core и postgis - PullRequest
Новая
       69

Рассчитайте географическое расстояние между двумя точками, используя NetTopologySuite с EF Core и postgis

0 голосов
/ 01 июля 2019

Я использую образ докера kartoza для запуска сервера postgres с postgis. У меня есть база данных, которую я использую, используя приложение ASP.NET Core и Enity Framework Core. База данных содержит таблицу с именем Park, которая представлена ​​следующей сущностью: 

[Table("Park")]
public class Park
{
     [Key]
     public int Id { get; set; }

     [Column(TypeName = "geography (point)")]
     public Point Location { get; set; }
}

Я создаю DbContext, как указано в документации: 

protected override void OnModelCreating(ModelBuilder modelBuilder)
{
    modelBuilder.HasPostgresExtension("postgis");

    modelBuilder.Entity<Park>()
        .Property(p => p.Location)
        .HasColumnType("geography (point)");
}

Я заполняю базу данных следующим кодом: 

if (!Parks.Any())
{
    var park = new Park
    {
        Location = new Point(48.8566, 2.3522)
        {
            SRID = 4326
        }
    };

    Parks.Add(park);
    this.SaveChanges();
}

Теперь, чтобы проверить мой код, я вычисляю расстояние между этим парком и некоторой точкой:

enter image description here

Результат соответствует геометрическому расстоянию между двумя точками в декартовой системе координат, а не географическому расстоянию. Можно ли получить эквивалент следующего запроса:

rap

1 Ответ

1 голос
/ 02 июля 2019

Как уже упоминалось здесь , для расчета географического расстояния необходима проекция системы координат.ProjNet4GeoAPI используется для выполнения проекции с использованием метода расширения: 

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);
    }
}

Теперь результат кажется более правильным:

enter image description here

PS: В дополнение к пакету слепков ProjNet4GeoAPI для метода расширения необходим также NetTopologySuite.

...