地理坐标系、投影坐标系及对应WKIDW在生成缓冲区时的应用

首先，说一下问题，之前在利用ST_Buffer（上篇博客有说到）做缓冲区生成时，是在地理坐标系的基础上进行生成的，因此会出现生成的半径跟真实半径不匹配的情况甚至相差一千倍以上，并且在生成shp文件后，加载到arcgis时会出现与影像图不匹配，需要根据以下步骤对其进行坐标系转换。

 生成的缓冲区shp文件在arcgis中加载时，与影像图坐标不匹配时：

解决方法：对矢量数据加载投影参照坐标系：数据管理工具—>投影和变换à要素à投影—>输入坐标系：地理坐标系；输出坐标系：投影坐标系；

这里可以arcgis种进行坐标系转换，但还是存在一个问题，就是缓冲区半径跟输入的值相差很大。因此这里不能从根本上解决问题。

解决方案，posgresql同样提供了坐标系转换函数ST_Transform（），此函数可以根据坐标系的WKID进行坐标系之间的转换。而ST_Buffer中的缓冲区半径要在墨卡托投影坐标系下，因此这里就可以现将地理坐标系转为墨卡托投影坐标系，然后就可以根据输入的半径进行计算了。ST_DWithin（）也是同样的道理。举例如下，其中4326代表WGS84坐标系（GCS_WGS_1984 是WGS1984），3857代表墨卡托投影坐标系（102100 WGS_1984_web_mercator_auxiliary_sphere）：

1、ST_Transform(ST_GeomFromText('POLYGON((743238 2967416,743238 2967450,743265 2967450,743265.625 2967416,743238 2967416))',4326),3857)

2、update poly set oneGradeBuffer =ST_Buffer(ST_Transform(ST_GeomFromText('LINESTRING(-122.323002 47.6455388,-122.3229709 47.6434522)'),4326),3857 ),50 , 'endcap=round join=round')where objectid=243

3、select * from poly where ST_DWithin(ST_Transform(geom,3857), (ST_Transform((select geom from poly where objectid=1)),3857), 50)=true