python实现四参数七参数坐标转换

四参和七参转换，转换坐标没有什么问题，我附录了一些坐标信息检验，不过七参数的数值有些问题，不知道是哪里有些错误，希望大家指正，下面是完整代码：

在这里插入代码片
```import numpy as np
from numpy import linalgdef four():def four_number(x1,x2,x3,x4,xx1,xx2,xx3):"""由三个共同点求两个不同平面直角坐标系的坐标转换"""def caculate_four(x4,old,new):"""计算新坐标"""        mid=new-oldA_=np.array([[1,0,x4[0],-x4[1]],[0,1,x4[1],x4[0]]])old_=np.array([[x4[0]],[x4[1]]])xxx=linalg.lstsq(A,mid,rcond=-1)v=np.dot(A_,xxx[0])new_=old_+vreturn new_A=np.array([[1.0,0,x1[0],-x1[1]],[0,1.0,x1[1], x1[0]],[1.0,0,x2[0],-x2[1]],[0,1.0,x2[1], x2[0]],[1.0,0,x3[0],-x3[1]],[0,1.0,x3[1], x3[0]]])old=np.array([[x1[0]],[x1[1]],[x2[0]],[x2[1]],[x3[0]],[x3[1]]])new=np.array([[xx1[0]],[xx1[1]],[xx2[0]],[xx2[1]],[xx3[0]],[xx3[1]]])return caculate_four(x4,old,new)x1 = [4031675.025,552545.805]xx1 = [4031634.029,552496.906]x2 = [4032992.528,547359.488]xx2 = [4032951.494,547310.569]x3 = [4026900.438,534348.966]xx3 = [4026859.289,534300.075]x4 = [4027240.551,541727.683]print("平面直角新坐标为\n",four_number(x1,x2,x3,x4,xx1,xx2,xx3))##    四参数转换数据
##    旧：a 4031675.025,552545.805
##        b 4032992.528,547359.488
##        c 4026900.438,534348.966
##
##    新：a 4031634.029,552496.906
##        b 4032992.528,547359.488
##        c 4026900.438,534348.966
##
##    待求旧：4027240.551,541727.683
##    待求新：4027199.461,541678.802
four()def seven():def seven_number(x1,x2,x3,x4,xx1,xx2,xx3):"""由三个共同点求两个不同空间直角坐标系的坐标转换"""def caculate_seven(x4,old,new):"""计算新坐标"""mid=new-oldA_=np.array([[1,0,0,x4[0],0,-x4[2],x4[1]],[0,1,0,x4[1],x4[2],0,-x4[0]],[0,0,1,x4[2],-x4[1],x4[0],0]])old_=np.array([[x4[0]],[x4[1]],[x4[2]]])xxx=linalg.lstsq(A,mid,rcond=-1)v=np.dot(A_,xxx[0])new_=old_+vreturn new_A=np.array([[1,0,0,x1[0],0,-x1[2],x1[1]],[0,1,0,x1[1],x1[2],0,-x1[0]],[0,0,1,x1[2],-x1[1],x1[0],0],[1,0,0,x2[0],0,-x2[2],x2[1]],[0,1,0,x2[1],x2[2],0,-x2[0]],[0,0,1,x2[2],-x2[1],x2[0],0],[1,0,0,x3[0],0,-x3[2],x3[1]],[0,1,0,x3[1],x3[2],0,-x3[0]],[0,0,1,x3[2],-x3[1],x3[0],0]])old=np.array([[x1[0]],[x1[1]],[x1[2]],[x2[0]],[x2[1]],[x2[2]],[x3[0]],[x3[1]],[x3[2]]])new=np.array([[xx1[0]],[xx1[1]],[xx1[2]],[xx2[0]],[xx2[1]],[xx2[2]],[xx3[0]],[xx3[1]],[xx3[2]]])return caculate_seven(x4,old,new)#共同点坐标x1=[3516952.028,519078.188,36.393]x2=[3516556.925,519081.738,31.720]x3=[3516569.509,518700.605,31.097]xx1=[3516904.452,519026.990,34.565]xx2=[3516509.353,519030.538,29.893]xx3=[3516521.939,518649.409,29.272]#待转换坐标x4=[3516587.828,518636.550,30.246]print("\n空间直角新坐标为\n",seven_number(x1,x2,x3,x4,xx1,xx2,xx3))##    七参数转换数据
##    旧：a 3516952.028,519078.188,36.393
##        b 3516556.925,519081.738,31.720
##        c 3516569.509,518700.605,31.097
##
##    新：a 3516904.452,519026.990,34.565
##        b 3516509.353,519030.538,29.893
##        c 3516521.939,518649.409,29.272
##
##    待求旧：d 3516587.828,518636.550,30.246
##    待求新：d 3516540.258,518585.355,28.421seven()

python实现四参数七参数坐标转换相关推荐

中海达ihand30手簿使用说明_如何使用RTK手簿求坐标转换参数（四参数/七参数）...
参数计算是用于两个坐标系统之间的转换关系,GPS 接收机输出的数据是CGCS2000经纬度坐标,但是有些时候根据甲方的要求,需要转化到施工测量坐标,这个时候就需要软件进行坐标转换参数的计算和设置,而转 ...
南方h5手簿求转换参数_如何使用RTK手簿求坐标转换参数（四参数/七参数）
参数计算是用于两个坐标系统之间的转换关系,GPS 接收机输出的数据是CGCS2000经纬度坐标,但是有些时候根据甲方的要求,需要转化到施工测量坐标,这个时候就需要软件进行坐标转换参数的计算和设置,而转 ...
坐标转换四参数/七参数/正形变换 ∈ C# 编程笔记
更新日期:2020/4/7 文章目录 2.平面坐标转换四参数/正形变换 2.1 直接法 2.2 平差法 2.3 正形变换法附简化版正形变换法代码 2.4 七参数坐标转换 [注1]其中的代码也许 ...
ArcGIS栅格数据四参数七参数投影变换
对于不同坐标系的投影变化,在arcgis的Toolbox中可以直接调用工具,且支持矢量和栅格的投影变换.然而,在平常的项目开发过程中,通常需要把地方坐标系和cgcs2000坐标系进行坐标转换.由于转换 ...
四参数坐标转换c++_写给测绘新手，四参数与七参数坐标转换含义及区别
坐标转换时RTK技术里不可缺少的重要部分.不同坐标系,其椭球体的长半径,短半径和扁率是不同的.比如我们常用的四种坐标系北京54.西安80.WGS84.CGCS2000所对应的椭球体,它们的椭球体参数就 ...
什么是七参数？什么是四参数？有什么不同
上一节我们说到了各种坐标系问题,今天我们先来看看什么是测绘中的七参数和四参数? 1.四参数两个不同的二维平面直角坐标系之间转换通常使用四参数模型,四参数适合小范围测区(5 km以内)的空间坐标转换. ...
四参数坐标转换c++_GPSRTK坐标转换及四参数、七参数适用条件
工程测量仪器已由经纬仪.全站仪过渡到GNSS(全球卫星导航系统),特别是公路行业,GPS-RTK作为GNSS的一种应用目前已十分普及.现阶段GPS-RTK以WGS-84 坐标系统为主流,所发布的星历参 ...
四参数与七参数坐标转换含义及区别
目录一.前言二.四参数转换三.七参数转换四.区别一.前言不同坐标系,其椭球体的长半径,短半径和扁率是不同的.比如我们常用的四种坐标系北京54.西安80.WGS84.CGCS2000所对应的 ...
坐标转换流程与公式七参数四参数
7参数转换整体流程(不同椭球之间): 平面直角坐标系-->大地坐标系-->空间直角坐标系转换(To) 空间直角坐标系-->大地坐标系-->平面直角坐标系举个栗子,比如从BJ ...

python实现四参数七参数坐标转换

python实现四参数七参数坐标转换相关推荐

最新文章

热门文章