[发明专利]楼房密集区隐蔽点的坐标测量方法

权利要求书

1.一种楼房密集区隐蔽点的坐标测量方法，其特征是：包括

步骤一、布设像控点：基于获取的包含测区的米级遥感影像选取至少六个像控点，且所述像控点均选择在楼房的顶部；依次对选择的像控点进行编号，同时对应记录每个像控点的像方坐标(u_i,v_i)；

步骤二、地面控制点数据采集：依据步骤一选取的像控点，使用多台GPS RTK同时测得像控点对应的地面控制点的物方坐标(X_i,Y_i)，并进行记录；

步骤三、建立像方坐标和物方坐标的数学关系式：基于步骤一获取的像控点的像方坐标(u_i,v_i)和步骤二获取的地面控制点的物方坐标(X_i,Y_i)，利用最小二乘算法构建如下数学关系式：

$u+k_1{\mathrm{ur}}^2=\frac{L_{1}X+L_{2}Y+L_{3}Z+L_4}{L_{9}X+L_{10}Y+L_{11}Z+1}$

$v+k_1{\mathrm{vr}}^2=\frac{L_{5}X+L_{6}Y+L_{7}Z+L_8}{L_{9}X+L_{10}Y+L_{11}Z+1}$

并变形换算成如下的误差方程式：

M＝BL+P

其中：

$B=-\frac{1}{A}\left[\begin{array}{cccccccccccc}X& Y& Z& 1& 0& 0& 0& 0& uX& uY& uZ& {\mathrm{Aur}}^2\\ 0& 0& 0& 0& X& Y& Z& 1& vX& vY& vZ& {\mathrm{Avr}}^2\end{array}\right]$

L＝[L₁ L₂ L₃ L₄ L₅ L₆ L₇ L₈ L₉ L₁₀ L₁₁ k₁]^T

$P=-\frac{1}{A}\left[\begin{array}{c}u\\ v\end{array}\right]$

A＝L₉X+L₁₀Y+L₁₁Z+1

解求出L系数；

步骤四、解算隐蔽点的物方坐标：基于步骤三建立的像方坐标和物方坐标的数学关系式及解求的L系数，换算得到隐蔽点的物方坐标的误差方程式：

V＝NS+Q

其中：

S＝[X Y Z]^T

$Q=-\frac{1}{A}\left[\begin{array}{c}{L}_4+u\\ L_8+v\end{array}\right]$

u′＝u+k₁ur²

v′＝v+k₁vr²

并代入隐蔽点的像方坐标，解算出隐蔽点对应的物方坐标。