shp文件坐标系转换-PostGIS

问题：shp文件坐标转换
工具：PostGIS、PostgreSQL
本例子使用的是北京54转wgs84
参考文献：
1. srid索引： http://spatialreference.org/ref/epsg/21480/
2. PostGIS官方文档：http://www.postgis.org/docs/reference.html

1. 将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;