使用地理数据类型（DbGeography）时查询矩形区域答案

【问题标题】：Query for a rectangular area when using geography data type (DbGeography)使用地理数据类型（DbGeography）时查询矩形区域
【发布时间】：2014-12-06 10:19:06
【问题描述】：

我正在查询 SQL Server 中的表以在 Google 地图上显示结果，并且我正在服务器应用程序上使用 Entity Framework 和 DbGeography 数据类型。

简单地说，我从 Google Maps 实例中获取地图边界，并使用以下方法从中构造一个 DbGeography 实例：

var wkt = string.Format("POLYGON(({0} {1}, {0} {2}, {3} {2}, {3} {1}, {0} {1}))", EastLng, NorthLat, SouthLat, WestLng);
var bounds = DbGeography.FromText(wkt);

然后使用这个查询数据库（作为更大查询的一部分）：

dbset.Where(i => bounds.Intersects(l.GeographicLocation));

这很好用，而且我实际上已经在生产中工作了一段时间。今天我注意到有时有些结果在地图上显示不正确。

经过一番调查，我意识到这是因为地球的曲率；谷歌地图使用地球的矩形投影，但是当我创建一个矩形多边形时，它意味着这样的形状（从 SQL Server Management Studio 复制，从正在运行的查询中复制）：

查看非常大区域的地图（缩小）时效果会增加。但是当我放大时，地球曲率的影响会减小，查询似乎工作正常。

我需要做的是基本上查询数据库中位于 Lat/Lng 边界内的所有点。有没有办法做到这一点？

编辑：

这是我用来创建 DbGeography 实例的完整代码：

public class LatLngBounds
{
    // ...

    public DbGeography ToDbGeography()
    {
        return DbGeographyUtil.CreatePolygon(string.Format("POLYGON(({0} {1}, {0} {2}, {3} {2}, {3} {1}, {0} {1}))", EastLng, NorthLat, SouthLat, WestLng));
    }
}


public static class DbGeographyUtil
{
    public static DbGeography CreatePolygon(string wktString)
    {
        if (string.IsNullOrWhiteSpace(wktString))
            return null;

        var sqlGeography = SqlGeography.STGeomFromText(new SqlChars(wktString), DbGeography.DefaultCoordinateSystemId).MakeValid();

        var invertedSqlGeography = sqlGeography.ReorientObject();
        if (sqlGeography.STArea() > invertedSqlGeography.STArea())
        {
            sqlGeography = invertedSqlGeography;
        }

        return DbSpatialServices.Default.GeographyFromProviderValue(sqlGeography);
    }

}

编辑 2：

上例中来自谷歌地图的坐标是：

东北：Lat=53.6 Lng=92.5
西南：纬度=35.0 Lng=21.2

这是我使用 SQL Server Profiler 捕获的正在运行的查询的一部分：

declare @p4 sys.geography
set @p4=convert(sys.geography,0xE6100000010405000000463E933FAFD64A40070000006E3F354061E412079A894140070000006E3F354061E412079A894140020000808B245740463E933FAFD64A40020000808B245740463E933FAFD64A40070000006E3F354001000000020000000001000000FFFFFFFF0000000003)

SELECT @p4;

（我添加了SELECT 以查看它的外观，并生成了上图）

【问题讨论】：

如果你使用 geography 数据类型，或者 srid 为 4326 的几何，并在地理坐标中定义一个框，你不会有任何问题。我不确定我是否遵循谷歌地图部分——它是非常简单的投影，y 值本质上只是纬度的线性函数，我不明白为什么你会得到你所看到的曲率。给定的纬度线将映射到球形墨卡托（谷歌地图）中的相同 y 坐标。还有其他的投影不是这种情况，例如在英国使用的横轴墨卡托投影，但在这种情况下不是。
@JohnBarça 请看我的编辑。我在数据库中的所有类型都是地理。我想这些线是弯曲的，因为连接东南点和西南点的最短线不平行于赤道，而我从“地理”类型定义中记得（与几何相比）是它考虑了这种情况。
您能否举例说明来自 Google 地图的坐标，这些坐标用于生成上面的图像？
@JohnBarça 请参阅“编辑 2”:)
很抱歉，你是笨拙的，但是你是如何从谷歌地图中得到坐标的。我的理解是谷歌地图使用的是球形墨卡托，单位是米。

标签： sql-server entity-framework google-maps spatial

【解决方案1】：

我设法通过使用 DbGeometry 而不是 DbGeography 来解决这个问题

var boundsAsGeometry = DbGeometry.FromText(bounds.AsText(), bounds.CoordinateSystemId);

dbset.Where(x => boundsAsGeometry.Intersects(DbGeometry.FromText(x.GeographicLocation.AsText(), x.GeographicLocation.CoordinateSystemId)));

请注意，如果您有 GeographicLocation 字段，这将不会使用现有的数据库索引。

【讨论】：