将未知 EPSG 转换为 EPSG 4326答案

【问题标题】：Convert Unknown ESPG to ESPG 4326将未知 EPSG 转换为 EPSG 4326
【发布时间】：2021-11-10 22:18:30
【问题描述】：

在 postgres 中运行空间查询时，出现以下错误：

psycopg2.errors.InternalError_: lwgeom_intersection: GEOS Error: TopologyException: Input geom 1 is invalid: Self-intersection at or near point 676597.26377722982 3673373.5709498846 at 676597.26377722982 3673373.5709498846

我不知道返回的是什么投影。输入是全部（据我所知）ESPG 4326，但这些显然不是在该投影中起作用的点。如何破译要转换的内容以便找到错误？

【问题讨论】：

您可以添加引发错误的查询吗？您是否尝试过ST_SRID 来检查您的几何图形的 SRS？ postgis.net/docs/ST_SRID.html
你能检查你的输入几何是否有效吗？

标签： postgresql postgis psycopg2 spatial-query

【解决方案1】：

错误被 psycopg2 误报告为几何错误，而实际上它是 SQL 错误。 SQL 已更改为使用返回为simple_envelope（以前为envelope）的简化多边形，但查询在CASE 语句的某些分支中查找envelope。实际错误是缺少列，但 psycopg2 返回了带有神秘坐标的错误几何错误。

【讨论】：