PostGis距离计算

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

我使用的所有坐标均取自 Google 地图（这意味着 EPSG 为 3857）。我现在需要做的是找到以米为单位的距离并使用以米为单位的 D_Within

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

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

尽管它们以拉德为单位，但第二个结果仅是第一个结果的两倍。所以这让我很困惑。 然后我尝试将其输出为米：

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

相同行给出的结果是：

2.872546829
4.572207435

这也不是我所期望的。 我对 PostGis 和 SRID 不太熟悉，所以这个问题可能看起来很简单，但请帮助我，我被困了没有@@

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

如果您想要两个纬度/经度点之间的距离（以米为单位），请使用ST_距离_球体。例如，对于您的第一组点，

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

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

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

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

EDIT:正如评论中指出的，从 Postgis 2.2 开始，该函数被重命名ST_距离球体