使用逆函数将几何转换为坐标会导致混合SRID错误

Question

我正在使用带postig扩展名的postgres。 我将坐标放入以下postgis函数中：

select ST_SetSRID(ST_point(18.533935546875, 48.23199134320962),4674) 

结果，我得到的几何类型为：01010000204212000000000000B08832402D0570E4B11D4840

当我想做逆向动作时

select st_x(st_transform(way,4326)), st_y(st_transform(way,4326)) from planet_osm_point where way ='01010000204212000000000000B08832402D0570E4B11D4840'

我得到错误：

ERROR:  Operation on mixed SRID geometries
********** Error **********
ERROR: Operation on mixed SRID geometries
SQL state: XX000

可能是什么问题？

Answer 1

错误消息中对该问题进行了某种描述：您试图在planet_osm_point.way找到相同的几何，但是它们的SRID与4674不同。实际上，这是这样的：

SELECT ST_SetSRID(ST_Point(0, 0), 4674) = ST_SetSRID(ST_Point(0, 0), 4326);

因此，您需要选择一个通用的SRID，然后将其中一个转换为另一个

SELECT ST_SetSRID(ST_Point(0, 0), 4674) =
                      ST_Transform(ST_SetSRID(ST_Point(0, 0), 4326), 4674);

除非我不建议使用=(geometry, geometry)运算符，否则它需要相同的几何，这在转换后对于浮点错误不是很常见。

更好的办法是找到所有的way使用误差距离的小幅度内的几何形状ST_DWithin ：

select st_x(st_transform(way,4326)), st_y(st_transform(way,4326))
from planet_osm_point
where ST_DWithin(
         way,
         ST_Transform(
           ST_SetSRID(ST_Point(18.533935546875, 48.23199134320962), 4674),
           <way_SRID>),
         1e-6);

在这里，它将点转换为way的SRID（您从未指定过，因此请更新<way_SRID> ），然后在1e-6度距离（约10厘米）内找到任何东西。 如果您想收集更多附近的way点，请增加它。