最小二乘法求解平面坐标转换四参数

本文链接：https://blog.csdn.net/dibowei2069/article/details/106039613

本文深入解析了四参数的转换公式，通过构造条件方程并利用最小二乘原理，详细介绍了如何计算未知数X，进而求解m和α，为地理坐标转换提供了数学依据。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

四参数的转换公式为：

$\begin{bmatrix} x_{2}\\ y_{2} \end{bmatrix}= \begin{bmatrix} \Delta x\\ \Delta y \end{bmatrix}+ m\begin{bmatrix} cos\alpha &-sin\alpha \\ sin\alpha & cos\alpha \end{bmatrix}\begin{bmatrix} x_{1}\\ y_{1} \end{bmatrix}$

构造条件方程：

$\begin{bmatrix} x_{2}\\ y_{2} \end{bmatrix} - \begin{bmatrix} x_{1}\\ y_{1} \end{bmatrix}= \begin{bmatrix} \Delta x\\ \Delta y \end{bmatrix}+ \begin{bmatrix} mcos\alpha &-msin\alpha \\ msin\alpha & mcos\alpha \end{bmatrix}\begin{bmatrix} x_{1}\\ y_{1} \end{bmatrix}-\begin{bmatrix} 1 &0 \\ 0 & 1 \end{bmatrix}\begin{bmatrix} x_{1}\\ y_{1} \end{bmatrix}\Rightarrow$

$\begin{bmatrix} x_{2}\\ y_{2} \end{bmatrix} - \begin{bmatrix} x_{1}\\ y_{1} \end{bmatrix}= \begin{bmatrix} \Delta x\\ \Delta y \end{bmatrix}+ \begin{bmatrix} mcos\alpha -1 &-msin\alpha \\ msin\alpha & mcos\alpha -1\end{bmatrix}\begin{bmatrix} x_{1}\\ y_{1} \end{bmatrix}$

设 $a$ = mcosa-1, b=msina,有:

$\begin{bmatrix} x_{2}\\ y_{2} \end{bmatrix} - \begin{bmatrix} x_{1}\\y_{1} \end{bmatrix} =\begin{bmatrix} \Delta x\\ \Delta y\ \end{bmatrix} + \begin{bmatrix} a & -b\\ b & a \end{bmatrix}\begin{bmatrix} x_{1}\\ y_{1} \end{bmatrix} \Rightarrow \begin{bmatrix} x_{2}\\ y_{2} \end{bmatrix} - \begin{bmatrix} x_{1}\\y_{1} \end{bmatrix} =\begin{bmatrix} \Delta x\\ \Delta y \end{bmatrix} + \begin{bmatrix} x_{1} & -y_{1}\\ y_{1} & x_{1} \end{bmatrix}\begin{bmatrix} a\\b \end{bmatrix}$

等价表达：

$\begin{bmatrix} x_{2}\\ y_{2} \end{bmatrix} - \begin{bmatrix} x_{1}\\y_{1} \end{bmatrix} = \begin{bmatrix} 1&0&x_{1}&-y_{1}\\0&1&y_{1}&x_{1} \end{bmatrix}\begin{bmatrix} \Delta x \\ \Delta y \\a\\b\end{bmatrix}$

其中， $Li=\begin{bmatrix} x_{2i}\\ y_{2i} \end{bmatrix} - \begin{bmatrix} x_{1i}\\y_{1i} \end{bmatrix}$ ， $Bi=\begin{bmatrix} 1 &0& x_{1i} &-y_{1i} \\ 0 &1 &y_{1i} & x_{1i} \end{bmatrix}$ , $X=\begin{bmatrix} \Delta x\\ \Delta y \\ a \\ b\end{bmatrix}$

即：

$L=\begin{bmatrix} x_{21}-x_{11}\\ y_{21} -y_{11}\\...\\ x_{2i}-x_{1i}\\ y_{2i} -y_{1i}\\ ...\\ x_{2n}-x_{1n}\\ y_{2n} -y_{1n} \end{bmatrix}$ , $B=\begin{bmatrix} 1 &0& x_{11} &-y_{11} \\ 0 &1 &y_{11} & x_{11} \\ &&... \\ 1 &0& x_{1i} &-y_{1i} \\ 0 &1 &y_{1i} & x_{1i} \\ &&...\\ 1 &0& x_{1n} &-y_{1n} \\ 0 &1 &y_{1n} & x_{1n} \\ \end{bmatrix}$ , $X=\begin{bmatrix} \Delta x\\ \Delta y \\ a \\ b\end{bmatrix}$ ， $P=\begin{bmatrix} 1& 0 &0 &0 \\ 0 & 1 & 0 &0 \\ 0 & 0 & 1&0 \\ 0 &0 &0 & 1 \end{bmatrix}$

根据最小二乘原理计算X：

$X=(B^TPB)^-^1B^TPL$

解得X后，可知 $a$ 和 $b$ ，可解得m和α：

$\left\{\begin{matrix}a = m\cos \alpha-1\\ b = m\sin \alpha \end{matrix}\right \Rightarrow \left\{\begin{matrix} \sin \alpha/\cos \alpha = b/(a +1) \\ (a+1)^2 + b^2 = m^2(\cos ^2 \alpha +\sin^2 \alpha) \end{matrix}\right. \Rightarrow \left\{\begin{matrix} \alpha = \arctan[b/(a+1)] \\ m = \sqrt{(a+1)^2+b^2} \end{matrix}\right.$