本文介绍了在使用ST_Buffer生成缓冲区时遇到的问题，即地理坐标系下的半径与实际半径不符，导致缓冲区与影像图不匹配。通过在ArcGIS中进行坐标系转换仍无法根本解决。解决方案是利用PostgreSQL的ST_Transform函数，将坐标系从WGS84（4326）转换到墨卡托投影（3857），确保缓冲区计算准确。通过示例展示如何进行转换并应用到ST_Buffer和ST_DWithin操作中。
首先，说一下问题，之前在利用ST_Buffer（上篇博客有说到）做缓冲区生成时，是在地理坐标系的基础上进行生成的，因此会出现生成的半径跟真实半径不匹配的情况甚至相差一千倍以上，并且在生成shp文件后，加载到arcgis时会出现与影像图不匹配，需要根据以下步骤对其进行坐标系转换。
生成的缓冲区shp文件在arcgis中加载时，与影像图坐标不匹配时：
解决方法：对矢量数据加载投影参照坐标系：数据管理工具—>投影和变换à要素à投影—>输入坐标系：地理坐标系；输出坐标系：投影坐标系；
这里可以arcgis种进行坐标系转换，但还是存在一个问题，就是缓冲区半径跟输入的值相差很大。因此这里不能从根本上解决问题。
解决方案，posgresql同样提供了坐标系转换函数ST_Transform（），此函数可以根据坐标系的WKID进行坐标系之间的转换。而ST_Buffer中的缓冲区半径要在墨卡托投影坐标系下，因此这里就可以现将地理坐标系转为墨卡托投影坐标系，然后就可以根据输入的半径进行计算了。ST_DWithin（）也是同样的道理。举例如下，其中4326代表WGS84坐标系（GCS_WGS_1984 是<
最低0.47元/天 解锁文章
扫一扫
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
