postgresql - PostgreSQL 数据库中的 ST_Distance 函数计算错误的距离？-6ren

postgresql - PostgreSQL 数据库中的 ST_Distance 函数计算错误的距离？

转载作者：行者123 更新时间：2023-11-29 12:25:31

我们正在开发一个应用程序，该应用程序涉及根据存储在 PostgreSQL 表中的多边形计算给定点的最短距离。

我们正在使用 PostgreSQL 数据库中的 ST_Distance 函数。

当我们将计算出的距离与 Google 地球进行比较时，它们之间存在巨大差异。

当然，Google 和 PostgreSQL 都不可能是错的(或者他们是错的？)，所以我们显然在这里遗漏了一些东西。

知道哪里出了问题吗？

我在下面给出了我们使用 PostgreSQL 测试过的示例查询以及来自 Google 地球的屏幕截图。

SELECT ST_Distance(Place1, Place2) As Place1ToPlace2
    , ST_Distance(Place1, Spot1) As Place1ToSpot1
    , ST_Distance(Place1, Spot2) As Place1ToSpot2
    , ST_Distance(Place2, Spot1) As Place2ToSpot1
    , ST_Distance(Place2, Spot2) As Place2ToSpot2
FROM (SELECT
    ST_PolygonFromText('SRID=4326;POLYGON((-74.0050636293915 40.75265123968514,-74.00500355126653 40.75268991743845,-74.00498169169283 40.75267084386348,-74.00503571044075 40.75263867886528,-74.0050636293915 40.75265123968514))') as Spot1
    ,ST_PolygonFromText('SRID=4326;POLYGON((-74.00503571044075 40.75263867886528,-74.00498225451273 40.75267084385684,-74.00495878551709 40.75265859837483,-74.00501023946696 40.75262521978885,-74.00503571044075 40.75263867886528))') as Spot2
    ,ST_GeogFromText('SRID=4326;POINT(-74.00489 40.752894)') As Place1
    ,ST_GeogFromText('SRID=4326;POINT(-74.004774 40.752846)') As Place2
    ) As foo  ;

它产生以下值:

place1toplace2 |place1tospot1 |place1tospot2 |place2tospot1 |place2tospot2 |
---------------|--------------|--------------|--------------|--------------|
11.152362504   |24.608417285  |25.977083731  |26.004190091  |26.011579435  |

以下是 Google Earth 的屏幕截图:

提前致谢!

以下是从 Google 地球导出的 KML:

<?xml version="1.0" encoding="UTF-8"?>
<kml xmlns="http://www.opengis.net/kml/2.2" xmlns:gx="http://www.google.com/kml/ext/2.2" xmlns:kml="http://www.opengis.net/kml/2.2" xmlns:atom="http://www.w3.org/2005/Atom">
<Document>
    <name>Spot 1.kml</name>
    <Style id="inline">
        <LineStyle>
            <color>ff0000ff</color>
            <width>2</width>
        </LineStyle>
        <PolyStyle>
            <fill>0</fill>
        </PolyStyle>
    </Style>
    <Style id="inline0">
        <LineStyle>
            <color>ff0000ff</color>
            <width>2</width>
        </LineStyle>
        <PolyStyle>
            <fill>0</fill>
        </PolyStyle>
    </Style>
    <StyleMap id="inline1">
        <Pair>
            <key>normal</key>
            <styleUrl>#inline</styleUrl>
        </Pair>
        <Pair>
            <key>highlight</key>
            <styleUrl>#inline0</styleUrl>
        </Pair>
    </StyleMap>
    <Placemark>
        <name>Spot 1</name>
        <styleUrl>#inline1</styleUrl>
        <Polygon>
            <tessellate>1</tessellate>
            <outerBoundaryIs>
                <LinearRing>
                    <coordinates>
                        -74.0050636293915,40.75265123968514,0 -74.00500355126653,40.75268991743845,0 -74.00498169169283,40.75267084386348,0 -74.00503571044075,40.75263867886528,0 -74.0050636293915,40.75265123968514,0 
                    </coordinates>
                </LinearRing>
            </outerBoundaryIs>
        </Polygon>
    </Placemark>
</Document>
</kml>


<?xml version="1.0" encoding="UTF-8"?>
<kml xmlns="http://www.opengis.net/kml/2.2" xmlns:gx="http://www.google.com/kml/ext/2.2" xmlns:kml="http://www.opengis.net/kml/2.2" xmlns:atom="http://www.w3.org/2005/Atom">
<Document>
    <name>Spot 2.kml</name>
    <Style id="inline">
        <LineStyle>
            <color>ff0000ff</color>
            <width>2</width>
        </LineStyle>
        <PolyStyle>
            <fill>0</fill>
        </PolyStyle>
    </Style>
    <StyleMap id="inline0">
        <Pair>
            <key>normal</key>
            <styleUrl>#inline1</styleUrl>
        </Pair>
        <Pair>
            <key>highlight</key>
            <styleUrl>#inline</styleUrl>
        </Pair>
    </StyleMap>
    <Style id="inline1">
        <LineStyle>
            <color>ff0000ff</color>
            <width>2</width>
        </LineStyle>
        <PolyStyle>
            <fill>0</fill>
        </PolyStyle>
    </Style>
    <Placemark>
        <name>Spot 2</name>
        <styleUrl>#inline0</styleUrl>
        <Polygon>
            <tessellate>1</tessellate>
            <outerBoundaryIs>
                <LinearRing>
                    <coordinates>
                        -74.00503571044075,40.75263867886528,0 -74.00498225451273,40.75267084385684,0 -74.00495878551709,40.75265859837483,0 -74.00501023946696,40.75262521978885,0 -74.00503571044075,40.75263867886528,0 
                    </coordinates>
                </LinearRing>
            </outerBoundaryIs>
        </Polygon>
    </Placemark>
</Document>
</kml>

!编辑 2!

我又尝试了一件事情。我从 Spot1 中提取了视觉上看起来最接近 Place1 的点并将其标记为 Place3。

当我们在 map 上看时，Place3 看起来比 Place2 更接近 Place1，但是当查询时，到 Place3 的距离给出了更高的值。

我已经检查了以下查询:

SELECT ST_Distance(Place1, Place2) As Place1ToPlace2
,ST_Distance(Place1, Place3) As Place1ToPlace3
FROM (SELECT
    ST_GeogFromText('SRID=4326;POINT(-74.00489 40.752894)') As Place1
    ,ST_GeogFromText('SRID=4326;POINT(-74.004774 40.752846)') As Place2
    ,ST_GeogFromText('SRID=4326;POINT(-74.00500355126653 40.75268991743845)') As Place3
    ) As foo;

结果如下:

place1toplace2 |place1toplace3 |
---------------|---------------|
11.152362504   |24.608417285   |

在 map 上时: Place1ToPlace3

以下是 Place1、Place2 和 Place3 的 KML(我已经从 KML 中删除了与样式相关的标签)

<?xml version="1.0" encoding="UTF-8"?>
<kml xmlns="http://www.opengis.net/kml/2.2" xmlns:gx="http://www.google.com/kml/ext/2.2" xmlns:kml="http://www.opengis.net/kml/2.2" xmlns:atom="http://www.w3.org/2005/Atom">
<Document>
    <name>Distance Comparison.kml</name>
    <Folder>
        <name>Distance Comparison</name>
        <open>1</open>
        <Style>
            <ListStyle>
                <listItemType>check</listItemType>
                <bgColor>00ffffff</bgColor>
                <maxSnippetLines>2</maxSnippetLines>
            </ListStyle>
        </Style>
        <Placemark>
            <name>Place 1 -  40°45&apos;9.78&quot;N  74° 0&apos;18.07&quot;W (40.752894, -74.00489)</name>
            <open>1</open>
            <LookAt>
                <longitude>-74.00500758183839</longitude>
                <latitude>40.75269419172616</latitude>
                <altitude>0</altitude>
                <heading>-0.0008536233435993688</heading>
                <tilt>29.8433509629012</tilt>
                <range>47.16429940085073</range>
                <gx:altitudeMode>relativeToSeaFloor</gx:altitudeMode>
            </LookAt>
            <styleUrl>#msn_1</styleUrl>
            <Point>
                <gx:drawOrder>1</gx:drawOrder>
                <coordinates>-74.00501944444444,40.75271666666666,0</coordinates>
            </Point>
        </Placemark>
        <Placemark>
            <name>Place 2 - 40°45&apos;9.54&quot;N 74° 0&apos;17.70&quot;W (40.752846, -74.004774)</name>
            <open>1</open>
            <LookAt>
                <longitude>-74.00500758183839</longitude>
                <latitude>40.75269419172616</latitude>
                <altitude>0</altitude>
                <heading>-0.0008536233435993688</heading>
                <tilt>29.8433509629012</tilt>
                <range>47.16429940085073</range>
                <gx:altitudeMode>relativeToSeaFloor</gx:altitudeMode>
            </LookAt>
            <styleUrl>#msn_2</styleUrl>
            <Point>
                <gx:drawOrder>1</gx:drawOrder>
                <coordinates>-74.00491666666667,40.75265,0</coordinates>
            </Point>
        </Placemark>
        <Placemark>
            <name>Place 3 -  40°45&apos;9.68&quot;N  74° 0&apos;18.01&quot;W (40.75268991743845 -74.00500355126653)</name>
            <open>1</open>
            <styleUrl>#msn_3</styleUrl>
            <Point>
                <coordinates>-74.00500277777778,40.75268888888889,0</coordinates>
            </Point>
        </Placemark>
    </Folder>
</Document>
</kml>

最佳答案

您在 PostGIS 中的功能似乎与它们在 Google 中的位置不同...这是我通过可视化得到的结果。

现在还有两个问题。一、geometry and geography的区别.二、距离计算。

geographies 上的 st_distance 返回以米为单位的值。这就是为什么你的第一个查询的输出基本上是你所期望的，给定特征的位置(因为你的点是地理，Postgres 为你将多边形转换为地理)。但是，geometries 上的 st_distance 返回以基础投影为单位的值，在本例中为度数。这解释了您的第二个查询的输出。

转到 distance calculation .在这里尽量简洁，但影响您计算的最大因素是基础投影。根据我的理解，地理类型更适合较大的区域——在适合您正在分析的区域的投影中使用几何图形，您会得到更准确的结果。 “suitable projection”取决于您正在查看世界的哪些部分以及您的目的是什么。

此外，不要假设 Google 的数据是正确的 - 我进行了快速搜索，但没有找到任何关于他们如何计算距离的信息。如果他们使用不同的方法，您会得到不同的答案...

编辑这是带有显示在屏幕截图中的坐标的查询。注意坐标值的差异!另外请注意，您屏幕截图中显示的度分秒与括号中的十进制度值不匹配!

SELECT ST_Distance(Place1, Place2) As Place1ToPlace2
    , ST_Distance(Place1, Spot1) As Place1ToSpot1
    , ST_Distance(Place1, Spot2) As Place1ToSpot2
    , ST_Distance(Place2, Spot1) As Place2ToSpot1
    , ST_Distance(Place2, Spot2) As Place2ToSpot2
FROM (SELECT
    ST_PolygonFromText('SRID=4326;POLYGON((-74.0050636293915 40.75265123968514,-74.00500355126653 40.75268991743845,-74.00498169169283 40.75267084386348,-74.00503571044075 40.75263867886528,-74.0050636293915 40.75265123968514))') as Spot1
    ,ST_PolygonFromText('SRID=4326;POLYGON((-74.00503571044075 40.75263867886528,-74.00498225451273 40.75267084385684,-74.00495878551709 40.75265859837483,-74.00501023946696 40.75262521978885,-74.00503571044075 40.75263867886528))') as Spot2
    ,ST_GeogFromText('SRID=4326;POINT( -74.005019 40.752717)') As Place1
    ,ST_GeogFromText('SRID=4326;POINT(-74.004917 40.752650)') As Place2
    ) As foo  ;

返回的值为 11.38223433、3.27827391、5.99175215、5.93327383、3.65564537。它们与 Google 的结果相差不到 cm。

关于postgresql - PostgreSQL 数据库中的 ST_Distance 函数计算错误的距离？，我们在Stack Overflow上找到一个类似的问题： https://stackoverflow.com/questions/40504193/

文章推荐： mysql - 如何使用其程序恢复删除的 mysql 数据库？

文章推荐： python - Django authenticate() 重置密码

文章推荐： python - 使用 idx 自动增量将数据框插入 postgresql sqlalchemy

algorithm - 距离(B)+ 距离(A-B)
A是不同元素的序列，B是A的子序列，A-B是A中的所有元素，但不是B中的所有元素距离(A) = 总和|a(i)-a(i+1)|从 i=1 到 n-1找到一个子序列 B 使得 Dist(B)+Dist(
r - 许多矩阵对之间的相似性/距离
我想通过计算每对中所有(多维)点集之间距离的平均值来量化组相似性。我可以很容易地手动为每对组手动完成此操作，如下所示: library(dplyr) library(tibble) library(
OpenXML 距离、大小单位
在 OpenXML 中用于指定大小或 X、Y 坐标的度量单位是什么？ (介绍)。将那些与像素匹配是否有意义，如果是这样，那些如何转换为像素？ graphicFrame.Transform = new
r - 如何替换过渡层中的值？ (距离)
我想知道是否有人可以帮助我替换过渡层中的值。如果我尝试: transitionlayer[transitionlayer >= 0.14] = 0.14 : comparison (5) is
Firebase - 按自定义功能排序(距离)
我在 firebase 中有一个列表，其中包括地理位置(经度和纬度)，并且我想获得距给定坐标最近的 10 个位置。我正在从 MySQL 过渡，在那里我将计算 SELECT 中的距离, 并在 ORDE
Python根据2个GPS坐标计算速度、距离、方向
如何在 Python 中根据 2 个 GPS 坐标计算速度、距离和方向(度)？每个点都有纬度、经度和时间。我在这篇文章中找到了半正矢距离计算: Calculate distance between
java - 距离出租车几何形状
关闭。这个问题需要多问focused 。目前不接受答案。想要改进此问题吗？更新问题，使其仅关注一个问题 editing this post . 已关闭 6 年前。 Improve this ques
python - 标记曲线之间的最大偏差/距离
我只想使用 matplotlib 标记两条曲线之间发生最大偏差的位置。请帮助我。垂直距离适用于 Kolmogorov–Smirnov test import numpy as np %matplot
linux - 查找重复行之间的平均时间/距离
我有一个包含数万行重复项的文件。我想根据行号找到重复项之间的平均时间/距离。例如:(其中第一列是行号) 1 string1 2 string2 3 string2 4 string1 5 strin
使用公式速度=距离/时间计算时间
用公式speed=distance/time计算时间但时间总是0我的输入是 distance=10 和 speed=5 我的输出必须 = 2 #include int main() { in
字符串相似度 -> Levenshtein 距离
我正在使用 Levenshtein 算法来查找两个字符串之间的相似性。这是我正在制作的程序的一个非常重要的部分，因此它需要有效。问题是该算法没有发现以下示例相似: CONAIR AIRCON 算法给出
mysql - 距离+关键词搜索方案
对于一个房地产网站，我需要实现一个允许搜索文本和距离的搜索机制。当 lat 和 lon 记录在单独的列中时，在 MySQL 表上进行距离计算很容易，但房子往往有 LOT true/false 属性。
iphone - UIPanGestureRecognizer 距离
是否可以在触发前更改 UIPanGestureRecognizer 的距离？目前的实现似乎在触发前有 5-10 像素的距离余量，我想降低它如果可能的话。原因是我将 UIPanGestureRecog
3d - 计算两个网格之间的 Hausdorff 距离
我试图找到两个网格之间的偏差。例如在 3d 空间中定义的两组点之间的差异，我计划使用一些 3d 可视化工具来可视化距离，例如QT3d 或一些基于开放式 gl 的库。我有两组网格，基本上是两个 .ST
excel - 找出哪些细胞具有最小的 levenshtein 距离
所以，我有这个函数可以快速返回两个字符串之间的 Levenshtein 距离: Function Levenshtein(ByVal string1 As String, ByVal string2
OCR:加权 Levenshtein 距离
我正在尝试用字典创建一个光学字符识别系统。事实上，我还没有实现字典=) 我听说有一些基于 Levenstein 距离的简单指标，这些指标考虑了不同符号之间的不同距离。例如。 'N' 和 'H' 彼此
gis - 带有经纬度SRID的PostGIS中的真实(大圆)距离？
我在PostGIS数据库(-4326)中使用经纬度/经度SRID。我想以一种有效的方式找到最接近给定点的点。我试图做一个 ORDER BY ST_Distance(point, ST_GeomF
r - 沿线串查找坐标 x 距离
我想从线串的一端开始提取沿线串已知距离处的点的坐标。例如: library(sf) path % group_by(L1) %>% summarise(do_union =
r - 确定基于序列(距离)的聚类的理想聚类数
我已经编写了这些用于聚类基于序列的数据的函数: library(TraMineR) library(cluster) clustering <- function(data){ data <- s
iphone - 是否可以设置线之间的 UILabel 距离？
是否可以设置 UILabel 的行之间的距离，因为我有一个 UILabel 包含 3 行，并且换行模式是自动换行？最佳答案如果您指的是“前导”，它指的是类型行之间的间隙 - 您无法在 UILabe

行者123

个人简介

我是一名优秀的程序员,十分优秀！

作者热门文章

滴滴打车优惠券免费领取

全站热门文章

首页

博学

6Ren·AI

商城

postgresql - PostgreSQL 数据库中的 ST_Distance 函数计算错误的距离？