![5b50a5abf82da989c83844233ff754f7.png](https://i-blog.csdnimg.cn/blog_migrate/8e2ea1000cdb934a24136d3243683732.jpeg)
![89cf4ff166a99f0bae140e48f3e85237.png](https://i-blog.csdnimg.cn/blog_migrate/7887418759cadf52087fcbf6c951abf0.jpeg)
老家县城有个号称风水特别好的楼盘出了设计图,正好手里有县城的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米左右,国内一般也就这三种。
![bd03a6568ef6108cd235017846a7ade3.png](https://i-blog.csdnimg.cn/blog_migrate/211d4065d1200f87e8aca6e8eeb35b71.jpeg)
![dd5737984b5ab7f701ca7191d1d6b6a4.png](https://i-blog.csdnimg.cn/blog_migrate/bf2dc46119cf4ea3ccf0a48e5d209f99.jpeg)
![7b4b40c9e2956c845fa791893382fe60.png](https://i-blog.csdnimg.cn/blog_migrate/8efa72539fc583af15c251931ad23b2e.jpeg)
将红线界址点展点到DOM,发现有部分界址点区域到了右下角红色房子处,而红色房子是可以确定不在楼盘范围。偏差大概是东西50米左右的。所以先确定为设计图为北京54坐标系下的118.5E投影(吐槽一下这都啥年代了还用54坐标系)。
![88b6797e81be74ebfae3918603c56f40.png](https://i-blog.csdnimg.cn/blog_migrate/4a03268358b373a25d0bc650158e812a.jpeg)
然后就涉及到坐标系转换的操作了,正好手上有该地区的几个控制点,其中还有一个级别还挺高的A级点。数据涉密就不展示了。
这里我们使用平面四参数方式来进行坐标转换,两个以上的控制点可以用最小二乘法计算转换参数
![77c9c163c7504be69dce1d461be55bb9.png](https://i-blog.csdnimg.cn/blog_migrate/e3f099880f501b2cedce77b2701c10a3.png)
线性代数功底好的同学可以手算,也有很多现成的程序可以计算参数,甚至还可以自个写个程序是吧。
![9b38d79a138ffa36c176ffb593461237.png](https://i-blog.csdnimg.cn/blog_migrate/c65e8cf8052ad722198b9d8102a73bb6.jpeg)
![29ab8b7f46eb10377968e6517c1aa0a1.png](https://i-blog.csdnimg.cn/blog_migrate/c528ddcc5dc5497aa201c2cd995c82c3.jpeg)
小软件输数据麻烦了点,我们需要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。
![4d1ffc8afd6a88f37954768485db9699.png](https://i-blog.csdnimg.cn/blog_migrate/c9d945a194733b81b303727364d43082.jpeg)
换带成功后要进行坐标转换。arcgis不支持对图形直接四参数转换,只能将点坐标导出成文本,对数值进行计算,通过四参数公式计算成80坐标系下的坐标值。(后面会写一篇基于arcpy的四参数转换工具的文章,就不用这么麻烦可以直接对图形进行转换了)
数据放到EXCEL中用上面的四参数公式计算一下再展点成图层就是80坐标系下的界址点,叠加到DOM中发现这次范围还是很准的。
![04048932154ae005f54241ee8f90b372.png](https://i-blog.csdnimg.cn/blog_migrate/edaf49d698f32ddc00f8de1da73776b3.jpeg)
然后就是将那张手机拍下来的楼盘设计图配准上去了。
使用配准工具条,添加控制点,将十字丝尽量准确的捕捉到设计图上的界址位置,连接至展点好的界址点图层。因为畸变较大,界址点最好全用上,使用三次多项式方法配准。
![2ec6e4026791b8c4fa6437a77f507ec0.png](https://i-blog.csdnimg.cn/blog_migrate/3cffce1b0f5ba4959d25c413c5ba61c9.jpeg)
最后,用界址点勾勒一个小区范围内的面,裁剪一下就是文章开头的成果图了。
简单的数字化一下,叠加DEM看看哪个是风水最好的楼王。
![2f413ca30fea98c358fb553d6733dba6.png](https://i-blog.csdnimg.cn/blog_migrate/ddac69c506916f81b65ba6e6420c8715.jpeg)
就是这俩了,不过好像是别墅啊,买不起买不起。
(最后反复看了一下,应该没有涉密数据被我传上去)