时间戳

首页 > 解决方案 > 为什么相同几何的 ST_Equals 返回 false？

sql - 为什么相同几何的 ST_Equals 返回 false？

问题描述

为什么相同对象的ST_Equals（一个来自几何列，另一个来自从同一几何列提取的 geojson）返回 false？

前任。：

select st_equals(geom, st_setsrid(st_geomfromgeojson(st_asgeojson(geom)), 2180)) from city

柱几何的定义是：

ALTER TABLE city ADD COLUMN geom geometry(MultiPolygon,2180);

编辑：

使用 st_asgeojson(geom, 999999999) 查询也不起作用

select st_equals(geom, st_setsrid(st_geomfromgeojson(st_asgeojson(geom, 999999999)), 2180)) from city

编辑2：

可重现的例子：

SELECT st_equals(geom, st_setsrid(st_geomfromgeojson(st_asgeojson(geom, 999999999)), 2180))
FROM
  (SELECT st_geomfromewkb(decode('AQYAACCECAAABAAAAAEDAAAAAQAAAAsAAAAPqGkZxtsfQR+Tq/icNydBNwTIEZsUIEF/uKr0AfIm
QTsvQxUcvB9Bdw1O9C7IJkEVnpoEBD0fQcFK5/MexSZBwMHT9KXEHkFNt8zxIdwmQezy8+91nx5B
JW6L+tkeJ0G68RXtlYkeQdnY3QSabSdB7MjG/EsBH0FuiKD6qpQnQczdMBNOrB9BK783BvmeJ0GO
URITcx4gQbSMcQaieSdBD6hpGcbbH0Efk6v4nDcnQQEDAAAAAQAAAAUAAABERDkK5GcfQdNUQuem
ZCZBkJbI/psQH0FUMIjfsikmQQDtYP3jBR9BVr0n5jpcJkHS02H+iw0fQQL74OyGjyZBREQ5CuRn
H0HTVELnpmQmQQEDAAAAAQAAAAUAAAAkHKMesHYgQTy3RfkU7iZBi9UrHiJzIEG7SVjr0oMmQR/r
kBTaKSBBH3xW6YJ0JkHnHCARmg8gQbVkpO76nCZBJByjHrB2IEE8t0X5FO4mQQEDAAAAAQAAAAUA
AAA8b3IhIIwgQcCpIQq1bidBINh3JxC6IEFwnmYBGSwnQYN91h44eCBB6PT+/2AhJ0HeIjUWYDYg
QQi4hQfNWidBPG9yISCMIEHAqSEKtW4nQQ==', 'base64')) geom) a

标签： sqlpostgresqlgispostgis

解决方案

st_asGeoJson 默认将坐标小数减少到 9 位。如果源几何具有更多数字，则转换后的几何因此不同。

几何图形也使用双精度保存，有效数字为 15。如果您的数据（输入）较少，剩余的数字仍然被设置，但可以被视为噪声......但是这种噪声被导出/导入回来，并且可能与输入略有不同。

最简单的解决方法是将坐标捕捉到网格，其大小是您的输入精度。您应该在源几何体上进行捕捉（将数据添加到数据库时）以及读取 geojson 时

st_equals(ST_SnaptoGrid(geom,5),ST_SnaptoGrid(st_setsrid(st_geomfromgeojson(st_asgeojson(ST_SnaptoGrid(geom,5))), 2180),5))

推荐阅读