浅谈MIKE前处理中投影坐标处理问题

MIKE 中投影坐标一直是个问题，尤其对 2d 里的科氏力影响很大， 由于我们现获取基础资料都是 CAD 格式，在 GIS 里转 shp 后我们会发现很多是地方坐标，对于这种情况，小编也是无能无力，只有想办法 让 CAD 提供方去转换。现讲下国内普遍坐标的投影定义方法吧。

状况一：X 坐标为 6 位数，Y 坐标是 7 位数

一般从业主手里拿来的 CAD 图都是有坐标的（一般如下图坐标）， 但是我们不知道这个 CAD 图地理坐标是什么？这样我们就必须问下 CAD 提供方，求知地理坐标（beijing1954），转成 shp 后导入 GIS，这 种 X 坐标为 6 位数，Y 坐标是 7 位数（图 1），在国内基都属于高斯投 影，此时我们也可以通过谷歌、天地图去搜下该区域所在地的经纬度，通过经纬度（经度 119°）来确定中央经线（图 2）,然后查 GIS 坐标 系统中距离（经度 119°）较近的几投影坐标（图 3），不断进行定义投影，查看与在线地图匹配情况，如果匹配，则可考虑定义的投影坐标校准，本次定义 Beijing_1954_3_Degree_GK_CM_120E，可以相互匹 配（图 4）。此时认为本次投影坐标较为准确，可直接用于 2d 科氏力计算。

坐标位数图

经纬度查询图

GIS 中邻近投影坐标系

数据匹配图

状况二：X 坐标为 8 位数，Y 坐标是 7 位数

如果此时业主提供的 CAD 图转 shp 后，发现 X 坐标为 8 位 数，Y 坐标是 7 位数（图 5），此时中央经线是 X 坐标前两位数乘 以 3，投影就应该选 Beijing_1954_3_Degree_GK_Zone_40。

坐标位数图

坐标定义图

关于坐标小编亦数外行，但是 MIKE 软件模型中坐标一直是个很难的问题，基础数据的坐标会对计算结果及精度都会有一定的影响，尤其大区域湖面、跨区域海面等，坐标直接影响科氏力的作用；同时坐标不统一完全可以导致 FLOOD 耦合模型因11和21文件位置不匹配而无法计算。希望群友有更多好的方法可以共同交流学习分享。