Postgres使用距离连接查询返回不正确的结果

Question

我创建了一个表，其中包含位置及其空间坐标（纬度、经度（度数）），以postgres表示。我用于相同用途的命令有：

我在postgres中插入了空间坐标（经度、纬度）

我创建了一个表，其中包含位置及其空间坐标（纬度、经度（度数）），以postgres表示。我用于相同用途的命令有：

CREATE TABLE spatialTest(
  name character varying(100),
  the_geo geography
);
\copy spatialTest(name,the_geo) FROM 'testSpatial.csv' DELIMITERS E'\t' CSV;

testSpatial。csv包含以下值：

A   SRID=4326;POLYGON((0.178773 -127.681841,0.178711 -127.681962,0.179125 -127.682083,0.179176 -127.682006,0.179153 -127.681986,0.179143 -127.681962,0.179147 -127.681935,0.179166 -127.681913,0.179195 -127.681897,0.179244 -127.681886,0.179284 -127.681887,0.179336 -127.681904,0.179464 -127.681757,0.179489 -127.681736,0.179429 -127.681680,0.179370 -127.681516,0.179221 -127.681331,0.179184 -127.681185,0.179051 -127.681264,0.178822 -127.681499,0.178761 -127.681698,0.178796 -127.681703,0.178839 -127.681721,0.178857 -127.681736,0.178861 -127.681740,0.178871 -127.681756,0.178873 -127.681782,0.178859 -127.681809,0.178843 -127.681825,0.178812 -127.681839,0.178773 -127.681841))
B   SRID=4326;POINT(0.628912 -127.700922)

现在，我想找到彼此距离在50公里以内的所有空间位置。为此，我使用了以下命令：

select s1.name, s2.name  from spatialTest s1, 
     spatialTest s2  where ST_DWithin(s1.the_geo, s2.the_geo, 50000);

然而，令我惊讶的是，我发现尽管A和B彼此之间的距离大于50公里（精确地说是50.0995公里），但使用Chris Venese的大地测量公式发现的​ （以经纬度计算点和线段之间的距离）），但结果由postgres返回。有人能帮我弄清楚我到底错在哪里吗。

我使用的是PostgreSQL 9.6devel和Postgis版本，我使用的是：POSTGIS="2.2.1 r14555"

Answer 1

正如其他人所说，这是你的负担。。

NOTICE:  Coordinate values were coerced into range [-180 -90, 180 90] for GEOGRAPHY
CONTEXT:  COPY spatialtest, line 2, column the_geo: "SRID=4326;POINT(0.628912 -127.700922)"

很可能，您正在尝试加载（lat， long），而PostgreSQL需要（long， lat）。但是，我们可以绕过输入列的顺序。我们要做的就是，

>

创建一个临时表asgeometry`。

CREATE TEMP TABLE foo(name text, the_geo geometry);

切断SRID位（EKWT），在导入时失去SRID。

sed -i -e's/SRID=4326;//' ./foo.csv

加载输入。

\copy foo(name,the_geo) FROM 'foo.csv' DELIMITERS E'\t' CSV;

复制到主表中。

INSERT INTO spatialTest(name, the_geo)
SELECT
  name,
  ST_Union(
    ST_SetSRID(
      ST_MakePoint(st_y(geom), st_x(geom))
    ,4326)
  )::geography
FROM foo
CROSS JOIN ST_DumpPoints(foo.the_geo) AS t
GROUP BY name
ORDER BY name;

作为旁注：不要使用varchar作为名称。只需将其保留为TEXT。除非你有很好的理由设定角色限制。

然后跑你的距离...

SELECT ST_Distance(
  (SELECT the_geo FROM spatialtest WHERE name='A'), (SELECT the_geo FROM spatialtest WHERE name='B')
);
  st_distance   
----------------
 49740.52214618
(1 row)