老家县城有个号称风水特别好的楼盘出了设计图，正好手里有县城的DOM和DEM，还有几个控制点坐标。就想着把这规划图叠加上去看看哪栋楼才是风水最好的楼王。这张设计图是手机随意拍下来的，一看上面还标上了红线界址点坐标，于是就有了这个想法，学为所用。
首先分析一下数据。已知项目所在地位于117E和118.5E之间并且很靠近118.5E，设计图Y坐标为47开头6位数字，猜测为118.5°投影不带带号。
然后是确定坐标系，思路是这样子的：现有DOM和DEM数据是西安80坐标系，将界址点展点至图层叠加，看能否和楼盘位置大致对上。位置对上了并且很准的话，可以直接确定设计图也是西安80坐标系；CGCS2000坐标系和80坐标系一般是东西偏差100米左右，南北10米左右，依据这个也能确定是否为CGCS2000坐标系；都不是的话基本可以确定为北京54坐标系了，并且54和80坐标系一般是东西偏差50米左右，国内一般也就这三种。
将红线界址点展点到DOM，发现有部分界址点区域到了右下角红色房子处，而红色房子是可以确定不在楼盘范围。偏差大概是东西50米左右的。所以先确定为设计图为北京54坐标系下的118.5E投影（吐槽一下这都啥年代了还用54坐标系）。
然后就涉及到坐标系转换的操作了，正好手上有该地区的几个控制点，其中还有一个级别还挺高的A级点。数据涉密就不展示了。
这里我们使用平面四参数方式来进行坐标转换，两个以上的控制点可以用最小二乘法计算转换参数
线性代数功底好的同学可以手算，也有很多现成的程序可以计算参数，甚至还可以自个写个程序是吧。
小软件输数据麻烦了点，我们需要54转换至80的参数，所以这里源坐标输入控制点成果中的54坐标系数值，目标控制点输入80坐标系数值。我有五对控制点，总共10个坐标对，20个数字，眼睛都看花了。
需要注意的是计算参数要统一投影方式，源坐标和目标坐标投影带要一样，横坐标去掉带号（控制点成果里坐标一般带着带号）。
计算成果四参数为两个平移量，一个缩放量和一个旋转量。dx,dy,t,k。
现在有了54转换到80的转换参数就可以对红线界址点进行坐标转换了。需要注意的是转换参数针对的是117°E的投影数据，需要将界址点进行换带。用arcgis工具箱里的投影工具，将界址点由54坐标系118.5E投影至54坐标系117E。另外，投影工具需要将输入数据定义一个坐标系，刚刚我们只是展点，还未定义，使用arcgis工具箱里的定义投影工具。118.5°E投影不是arcgis预设的投影方式，还需要新建一个投影坐标系，先随便双击一个54坐标系下的投影坐标系，比如说Beijing_1954_3_Degree_GK_CM_117E。将其中的Central_Meridian改为118.5，并且将名称改为Beijing_1954_1.5_Degree_GK_CM_118.5E。
换带成功后要进行坐标转换。arcgis不支持对图形直接四参数转换，只能将点坐标导出成文本，对数值进行计算，通过四参数公式计算成80坐标系下的坐标值。（后面会写一篇基于arcpy的四参数转换工具的文章，就不用这么麻烦可以直接对图形进行转换了）
数据放到EXCEL中用上面的四参数公式计算一下再展点成图层就是80坐标系下的界址点，叠加到DOM中发现这次范围还是很准的。
然后就是将那张手机拍下来的楼盘设计图配准上去了。
使用配准工具条，添加控制点，将十字丝尽量准确的捕捉到设计图上的界址位置，连接至展点好的界址点图层。因为畸变较大，界址点最好全用上，使用三次多项式方法配准。
最后，用界址点勾勒一个小区范围内的面，裁剪一下就是文章开头的成果图了。
简单的数字化一下，叠加DEM看看哪个是风水最好的楼王。
就是这俩了，不过好像是别墅啊，买不起买不起。
（最后反复看了一下，应该没有涉密数据被我传上去）
1万+
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
