本文主要是介绍supermap idesktop产品postgis数据源支持地方坐标系，希望对大家解决编程问题提供一定的参考价值，需要的开发者们随着小编来一起学习吧！
作者：Neshoir
背景
在supermap idesktop早期版本中postgis数据源不支持地方坐标系，但项目上的矢量数据坐标系均为南京地方坐标系，目前需要用supermap平台，数据需迁移到postgis数据源，针对数据集重设地方坐标系不起作用。原因是不是epsg标准坐标系，postgis不支持。由于postgis支持扩展，所以咱们可以自实现。
数据
- 样例数据zf_jnwg_1.shp文件，带有南京地方坐标系信息。
步骤
- 生成wkt和proj格式的坐标系信息，借助于gdalsrsinfo工具，如截图。
- 在pg库中，将得到的南京地方坐标系的wkt和proj信息插入到spatial_ref_sys表中，定义epsg为9999。
INSERT INTO "spatial_ref_sys" ("srid","auth_name","auth_srid","srtext","proj4text") VALUES (9999,'EPSG',9999,'PROJCS["NJ08_118_50_CM",GEOGCS["GCS_GRS_1980",DATUM["Not_specified_based_on_GRS_1980_ellipsoid",SPHEROID["GRS_1980",6378137,298.257222101]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]], PROJECTION["Transverse_Mercator"],PARAMETER["false_easting",500000],PARAMETER["false_northing",0], PARAMETER["central_meridian",118.833333333],PARAMETER["latitude_of_origin",0], PARAMETER["standard_parallel_1",0],PARAMETER["standard_parallel_2",0],PARAMETER["scale_factor",1], UNIT["Meter",1]]','+proj=tmerc +lat_0=0 +lon_0=118.833333333 +k=1 +x_0=500000 +y_0=0 +ellps=GRS80 +units=m +no_defs ');
- 在idesktop中打开pg数据源，导入该数据集，查看坐标系为4490，查看数据集范围值是投影坐标系下的坐标值。
- 在pg库更新该数据集的几何对象的坐标系为9999。
# 更新几何对象的坐标系
SELECT UpdateGeometrySRID('zf_jnwg','smgeometry',9999);
# 刷新pg，触发事务
# 查询其对象的坐标系，可以从结果中看到已经是9999了。
SELECT st_srid("zf_jnwg".smgeometry) FROM "zf_jnwg";
# 查询数据集在系统表中的坐标系，可以从结果中看到已经是9999了。
SELECT geometry_columns.srid FROM geometry_columns WHERE geometry_columns.f_table_name='zf_jnwg';
- 在idesktop对针对zf_jnwg数据集重设南京地方坐标系，设置后关闭pg数据源和刷新pg库，这一步关键，需要触发smregister表中该数据集坐标系信息的修改。
-
idesktop里查看结果
- pg库查看该数据集的坐标系情况
SELECT smregister.smsrid,smregister.smprojectinfo FROM smregister WHERE smregister.smdatasetname='zf_jnwg';
结论
-
能正常做浏览，编辑，叠加分析，导出shp等操作，结果都是正确，arcmap和qgis也能正常浏览也能识别到其坐标系。但是不保存在生产环境上的复杂操作会有任何风险。
这篇关于supermap idesktop产品postgis数据源支持地方坐标系的文章就介绍到这儿，希望我们推荐的文章对编程师们有所帮助！
