如何在EF中使用NpSql计算Geography类型的ST_Distance？

你当前的问题核心有两点：

• 直接调用Point.Distance()计算的是平面笛卡尔距离，而SRID 4326是WGS84经纬度坐标系，这种计算方式无法得到真实的地表球面距离。
• 在EF+NpSql场景下，你使用的是Geometry类型，而要计算地理距离需要切换为Geography类型（对应PostGIS的geography类型）。

修正方案 #

1. 实体属性配置（EF Core） #

确保你的实体中地理坐标属性被配置为geography类型，而非默认的geometry：

protected override void OnModelCreating(ModelBuilder modelBuilder)
{
    modelBuilder.Entity<YourEntity>()
        .Property(e => e.Location)
        .HasColumnType("geography (point)") // 指定为地理类型
        .HasSRID(4326);
}

2. EF查询中计算地理距离 #

使用NpSql提供的EF函数StDistance来调用PostGIS的原生地理距离计算：

// 目标点：注意地理类型坐标顺序是【纬度，经度】
var targetPoint = NetTopologySuite.Geometries.GeographyPoint.Create(51.25, 12.25, 4326);

// 查询计算每个实体到目标点的距离（单位：米）
var distanceResults = dbContext.YourEntities
    .Select(e => new {
        EntityId = e.Id,
        DistanceInMeters = EF.Functions.StDistance(e.Location, targetPoint)
    })
    .ToList();

3. 内存中计算球面距离 #

如果需要在内存中直接计算两个地理点的距离，使用GeographyPoint类型：

// 地理类型坐标顺序：纬度在前，经度在后
var point1 = NetTopologySuite.Geometries.GeographyPoint.Create(51.25, 12.25, 4326);
var point2 = NetTopologySuite.Geometries.GeographyPoint.Create(51.25, 12.50, 4326);

// 输出球面距离，单位为米
Console.WriteLine(point1.Distance(point2));

关键注意事项 #

• 坐标顺序：Geography类型的坐标是【纬度，经度】，和Geometry类型的【经度，纬度】相反，不要搞混。
• 单位：SRID 4326下StDistance返回的距离单位是米，如需公里可除以1000转换。
• 类型区分：PostGIS的geometry用于平面几何计算，geography用于球面地理计算，必须对应正确的类型才能得到准确的地理距离。

内容的提问来源于stack exchange，提问作者Igor Gavrilov