8

我现在正在做一个室内地图导航应用程序，我想做的是在建筑物中建立一个地图点数据库。

我使用的所有坐标都取自谷歌地图（这意味着 EPSG 是 3857）。我现在需要做的是以米为单位的距离以及以米为单位使用 D_Within

当我尝试提取2点之间的距离时：

SELECT ST_DISTANCE(
ST_GeomFromText('POINT(' || StartLon || ' ' || StartLat || ')',3857),
ST_GeomFromText('POINT(' || EndLon || ' ' || EndLat || ')',3857))
FROM i3_building.floordata;

对于前 2 行：

Start: 103.776047 1.292149; End: 103.77607 1.292212 (3 meters away)
Start: 103.776070 1.292212; End: 103.77554 1.292406 (50 meters away)

给出的结果是：

2.59422435413724e-005
4.11096095831604e-005

即使它们在 rad 中，第二个结果也只是第一个结果的两倍。所以这让我很困惑。然后我尝试将其输出为米：

SELECT ST_DISTANCE(
    ST_GeographyFromText('POINT(' || StartLon || ' ' || StartLat || ')'),
    ST_GeographyFromText('POINT(' || EndLon || ' ' || EndLat || ')'))
FROM i3_building.floordata;

为相同的行给出的结果是：

2.872546829
4.572207435

这也不是我所期望的。我对 PostGis 和 SRID 不是很熟悉，所以这个问题可能看起来很简单，但请帮帮我，我不知道@@

4

2 回答 2

19

您的坐标参考系 (CRS) 是 4326，纬度/经度。这是与 Google 地图混淆的常见原因：3857 是 Google 地图用于其图块的 CRS，并且是基于球形地球的投影米。添加到 Google 地图的矢量源（KML 数据、GPS 转储等）往往采用纬度/经度 4326，以度为单位进行测量并即时转换。

如果您想要两个纬度/经度点之间的米距离，请使用ST_Distance_Sphere。例如，对于您的第一组点，

SELECT ST_Distance_Sphere(ST_MakePoint(103.776047, 1.292149),ST_MakePoint(103.77607, 1.292212));

这给出了 7.457 米。根据相同的查询，您的第二组点彼此相距 62.74 米。

请注意，还有ST_Distance_Spheroid，它采用第三个参数，即测量球体，即地球形状的近似值。这可能会更准确，但在小距离上可能并不重要。

ST_Distance 以投影坐标给出距离，这可能就是为什么插入纬度/经度值会得到奇怪的结果。

编辑：如评论中所述，从 Postgis 2.2 开始，此函数重命名为ST_DistanceSphere

于 2014-07-08T06:04:24.933 回答
0

这是伪代码

SELECT ST_Distance_Sphere(ST_MakePoint(lng1, lat1),ST_MakePoint(lng2, lat2));

此外，对于球体距离新功能“ST_DistanceSphere”使用以下

SELECT ST_DistanceSphere(ST_MakePoint(lng1,lat1),ST_MakePoint(lng2,lat2))
于 2021-02-02T19:56:50.947 回答