How to use ST_Intersects in WHERE clause(如何在WHERE子句中使用ST_Intersects)
本文介绍了如何在WHERE子句中使用ST_Intersects的处理方法,对大家解决问题具有一定的参考价值,需要的朋友们下面随着小编来一起学习吧!
问题描述
我想知道如何执行以下操作
WHERE NOT ST_Intersects(dataSource.LIDARDataAsSinglePolygon,geom)
查询:
def exec(self,geometryAsGeoJSONObjectAsString,geom):
query="""
WITH dataSource As(
SELECT DISTINCT
ST_AsGeoJSON(ST_Union(ST_GeomFromGeoJSON(feature->>'geometry'))) As LIDARDataAsSinglePolygonUsingST_UnionAsGeoJSONInEPSG4326,
ST_AsText(ST_Union(ST_GeomFromGeoJSON(feature->>'geometry'))) As LIDARDataAsSinglePolygon
)
FROM
(SELECT json_array_elements('{geometryAsGeoJSONObjectAsString}'::json->'features') AS feature) a WHERE NOT ST_Intersects(dataSource.LIDARDataAsSinglePolygon,geom);
""".format(geometryAsGeoJSONObjectAsString=geometryAsGeoJSONObjectAsString)
print(query)
data = self.connection.query(query,[])
# print(data)
return data
小提琴:
https://dbfiddle.uk/?rdbms=postgres_12&fiddle=8dc0ccc1017a220098022131e8d41063
推荐答案
您最有可能尝试将ST_Intersects
与表达式的别名一起使用,而不是与列或几何图形一起使用。您的逻辑是颠倒的:您需要首先解析和转储要素集合中的几何,然后应用过滤器,而不是反过来:
WITH datasource (feature) AS (
SELECT json_array_elements('{json_string}'::json->'features')
)
SELECT
ST_AsGeoJSON(ST_Union(ST_GeomFromGeoJSON(feature->>'geometry'))),
ST_AsText(ST_Union(ST_GeomFromGeoJSON(feature->>'geometry')))
FROM datasource
WHERE NOT ST_Intersects(ST_GeomFromGeoJSON(feature->>'geometry'),'{a geometry}');
这篇关于如何在WHERE子句中使用ST_Intersects的文章就介绍到这了,希望我们推荐的答案对大家有所帮助,也希望大家多多支持编程学习网!
沃梦达教程
本文标题为:如何在WHERE子句中使用ST_Intersects
猜你喜欢
- 我如何透明地重定向一个Python导入? 2022-01-01
- CTR 中的 AES 如何用于 Python 和 PyCrypto? 2022-01-01
- 如何使用PYSPARK从Spark获得批次行 2022-01-01
- 检查具有纬度和经度的地理点是否在 shapefile 中 2022-01-01
- ";find_element_by_name(';name';)";和&QOOT;FIND_ELEMENT(BY NAME,';NAME';)";之间有什么区别? 2022-01-01
- 使用公司代理使Python3.x Slack(松弛客户端) 2022-01-01
- 计算测试数量的Python单元测试 2022-01-01
- YouTube API v3 返回截断的观看记录 2022-01-01
- 使用 Cython 将 Python 链接到共享库 2022-01-01
- 我如何卸载 PyTorch? 2022-01-01