matlab坐标转换求参数,基于Matlab的七参数空间坐标转换研究与实现

2014年 12月 第 12卷第 6期 地 理 空 间 信 息 ( ()SI ^⋯11iAI I ()l{MA I 1()N Dec．，2014 Vo1．12,No．6 doi：10．3969／j．issrt．1672—4623．2014．06．022 基于 Matlab的七参数空问坐标转换研究与实现 刘 平 ，段 志强 ，谢 超 (1．湖北省测绘成果档案馆 ，湖北 武汉 430074) 摘 要 ：分析了七参数空间坐标转换模型和Madab语言在矩阵运算方面的优势。在此基础上开发实现了七参数空间坐标转换程 序。简要介绍 了利用 Madab语言开发此程序的一点体会 。 关键词 ：Matlab；七参数 ；空间坐标转换 中图分类号 ：1'226 文献标志码 ：}{ 文意编号 ：I672—4623(20I4)06一()()62—03 1 七参数空间坐标转换基本原理 坐标 系统 间的差异主要 来 自于坐标 系统的定 义差 ， 即原 点位 置、坐标轴 向 的定 向和 尺度 的定义 差。坐标 转换时应优先考虑 坐标 系统定义差异的转换 ，即首先 应完成相似 变换。在相似变换的基础上，再考虑对剩 余误 差进行 拟合 ，使 精度较 低 的坐标框 架点 附合 到精 度较 高的 坐标 系统 的框 架点 坐标 ，使统 一后 的坐标 系 框架点坐标具有较好的一致性。 七参数空间坐标转换 ，实际上是 2个不同基准面 Ro Y2 =(1+m)R0 剖 COSgyCOS COS sinez — siney — cOssxsinez+ sinexsineyCOSsz COS~xcOSSz+ sinexsinersinez sinexCOS~y sinexsinez+ COSExsineyCOSEz —sinexcOs￡z+ COSsxsineysinez cOssxCOSEy 一 般 ‰ ￡ 、 为微 小转角 ，可取 ： COSex COSCy= coscz 1 sincx = ，siney=。y， sincz=6z sinexsiney= sinexsinez= sineysinez： 0 于是有 ： I =(1+m)l一￡ ￡ 1 1+[藿] 式 (1)即为布 尔沙 一沃 尔夫七参 数空 间坐标转换 模 型。此模 型的转换步骤是首先作三轴旋转 ，再统 一 尺度 ，最后再平移 。如果能求解 出这 7个参数值 ， 则所有待转换 点的坐标转换也就 可以求得 了。公式 中含有 7个转换参数 ，为 了求得 7个转换参数 ，至 少 需 要 3个 公 共 点 (也 称 为 重 合 点 )， 当 多于 3个 公共点时，可按最小二乘法求得 7个参数的最或是 值 。 设 a1= 1， =a1ex，a3=a1sy，a4=a1sZ，可取 ： 收稿日期 ：2i)14-01—26 1 0 0 Xl 0 一 0 1 0 g zl 0 一 0 0 1 Zl 一 0 ％ 则可写 出如 下形式的误 一 1 ． 7 I ‘吧J已知 值 差 方程 ： 0 0 x1 0 一 f Y1 1 0 0 一j 0 1 z1f 一 Xlf 0 换值 △ △ 口1 口2 口4 AXo ZXYo △ 以1 以2 3 4 + 知 值 第 12卷第 6期 冀n 平等 ：基于 Matlab的七参数空J葡坐标转换研究与实现 ·63· 改写成矩阵形式为 ： V=一B ×6X + (5) 式中， =(△ ，△y0，△Zo，al，口2， ， 4) ，为 待托 最转换参数 向量 ； ：( ，， ， ) 为改正 数 向量 ；工= ( ，y2l，Z2