shp文件坐标系转换-PostGIS
问题:shp文件坐标转换
工具:PostGIS、PostgreSQL
本例子使用的是北京54转wgs84
参考文献:
1. srid索引: http://spatialreference.org/ref/epsg/21480/
2. PostGIS官方文档:http://www.postgis.org/docs/reference.html
- 将shp文件导入到Postgresql
1.1 使用命令行方式
shp2pgsql -c 路径\shp数据文件名 新建的数据表名 数据库名|psql -d 数据库名
示例如下:
shp2pgsql -c d:\test.shp test hzmap | psql -d hzmap
1.2使用QGIS
工具在数据库 → SPIT工具 → 将shp文件导入PostgreSQL
如上图所示,主要注意2个地方:
第一个空间参考标识符,即Srid
第二个数据库关系名称,即导入数据库后,新建的表名
2.进行坐标转换
可以在数据库执行如下语句,先查看一下目标的Srid,主要函数为st_srid():
select st_srid(the_geom),* from hehehe;转换是使用如下的语句,主要函数为st_transform(),这块不能直接使用update,因为原来的the_geom字段的srid为21480,不兼容4326的,要所以直接插入到新表当中:
select t.id,st_transform(the_geom,4326) as the_geom into he from hehehe as t;执行完毕后,可使用如下2条语句验证是否成功转换完毕,主要函数为st_asgeojson():
select st_srid(the_geom),* from hehehe;
select st_asgeojson(the_geom) from hehehe;
转载自:https://blog.csdn.net/isowang/article/details/50370756