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;