postgis 通过 st_dwithin 检索一定距离范围内的 元素

psql=‘SELECT x.geom,x.小区名称 FROM orcal网元geom工参lte铜陵 as x,(SELECT * FROM orcal网元geom工参lte铜陵 where 小区名称 like ‘%%TL联通大厦%%’ or 小区名称 like ‘%%TL枞阳旗山路与银塘路交叉口%%’) as b WHERE ST_DWithin(ST_Transform(x.geom, 2415),ST_Transform(b.geom, 2415),1000)’

c=gpd.read_postgis(psql,con=engine,geom_col=‘geom’)
print©

1、需要把4326格式转换为2415
ST_Transform(x.geom, 2415)
2、检索一定范围内的网元元素
ST_DWithin(ST_Transform(x.geom, 2415),ST_Transform(b.geom, 2415),1000)

你可能感兴趣的:(gis,python,pandas,sql,python)