[發(fā)明專利]一種基于不精確牛頓解的預(yù)條件共軛梯度區(qū)域網(wǎng)平差方法

有效

• 申請?zhí)枺?/td>: 201610119376.3
• 申請日: 2016-03-02
• 公開（公告）號: CN105760687B
• 公開（公告）日: 2017-05-17
• 發(fā)明（設(shè)計(jì)）人: 鄭茂騰;熊小東;朱俊鋒
• 申請（專利權(quán)）人: 中國地質(zhì)大學(xué)（武漢）
• 主分類號: G06F19/00
• 分類號: G06F19/00
• 代理公司: 北京輕創(chuàng)知識產(chǎn)權(quán)代理有限公司11212
• 代理人: 陳衛(wèi)

技術(shù)領(lǐng)域

本發(fā)明涉及測繪科學(xué)與技術(shù)領(lǐng)域，具體涉及一種基于不精確牛頓解的預(yù)條件共軛梯度區(qū)域網(wǎng)平差方法，主要應(yīng)用于超大規(guī)模測區(qū)攝影測量等領(lǐng)域。

背景技術(shù)

區(qū)域網(wǎng)平差技術(shù)經(jīng)過幾十年的發(fā)展，其方法和流程已經(jīng)相對成熟，并且在測繪領(lǐng)域得到了廣泛的應(yīng)用。然而隨著科技的加速進(jìn)步，新的傳感器不斷涌現(xiàn)，如航天領(lǐng)域的高分辨率衛(wèi)星，立體測繪衛(wèi)星，航空領(lǐng)域的規(guī)則航空攝影系統(tǒng)，傾斜航空攝影系統(tǒng)，無人機(jī)、飛艇攝影系統(tǒng)等。同時(shí)，全世界范圍內(nèi)三維建模應(yīng)用需求不斷增加，也使得大量地面車載攝影系統(tǒng)，近景攝影系統(tǒng)，普通數(shù)碼相機(jī)，甚至是智能手機(jī)，網(wǎng)絡(luò)圖片庫等采集的影像被用于三維建模。影像數(shù)據(jù)源越來越豐富的同時(shí)，其分辨率也不斷提高，以前數(shù)十米的衛(wèi)星影像如今可以達(dá)到最高0.35米(WorldView-3)，航空影像的分辨率更是進(jìn)入了厘米級時(shí)代。分辨率的增加必然會(huì)帶來數(shù)據(jù)量的增大，攝影時(shí)的航線設(shè)計(jì)也不再滿足傳統(tǒng)的條帶式規(guī)則分布，給相應(yīng)的數(shù)據(jù)處理方法帶來了一定挑戰(zhàn)，受法方程大小的限制，傳統(tǒng)的區(qū)域網(wǎng)平差技術(shù)流程已經(jīng)不能滿足大規(guī)模數(shù)據(jù)處理需求。特別是當(dāng)測區(qū)數(shù)據(jù)大小超過1萬張影像時(shí)，傳統(tǒng)的區(qū)域網(wǎng)平差方法對內(nèi)存的需求急劇增大，即使有少部分圖形工作站內(nèi)存容量足夠大，但大量的內(nèi)存占用使得計(jì)算效率大幅降低，同時(shí)，數(shù)據(jù)處理的硬件成本也隨之增大，上述問題均阻礙了各類新型傳感器數(shù)據(jù)的廣泛應(yīng)用。

發(fā)明內(nèi)容

為了解決上述問題，本發(fā)明在區(qū)域網(wǎng)平差中，引入預(yù)條件共軛梯度法求解大規(guī)模線性方程組(法方程)，避免存儲(chǔ)大規(guī)模法方程系數(shù)矩陣，使得對大規(guī)模測區(qū)(1萬至10萬張影像，下文中統(tǒng)稱此類測區(qū)為大規(guī)模測區(qū))數(shù)據(jù)的區(qū)域網(wǎng)平差成為可能，其流程也更加易于并行化設(shè)計(jì)。預(yù)條件共軛梯度法是一個(gè)迭代求解的過程，每一次迭代都需要遍歷所有的像點(diǎn)觀測值數(shù)據(jù)，每一次區(qū)域網(wǎng)平差迭代內(nèi)部又包括N次預(yù)條件共軛梯度法迭代，N是預(yù)條件共軛梯度法迭代收斂次數(shù)，因而計(jì)算量較傳統(tǒng)的攝影測量區(qū)域網(wǎng)平差流程要大很多，為了減少預(yù)條件共軛梯度法的迭代次數(shù)，本發(fā)明采用一種不精確牛頓解取代法方程的精確解，使得預(yù)條件共軛梯度法可以快速的收斂，從而提高整個(gè)區(qū)域網(wǎng)平差的求解速度。

本發(fā)明所要解決的技術(shù)問題是提供一種基于不精確牛頓解的預(yù)條件共軛梯度區(qū)域網(wǎng)平差方法，能夠解決現(xiàn)有技術(shù)的不足。

本發(fā)明解決上述技術(shù)問題的技術(shù)方案如下：

一種基于不精確牛頓解的預(yù)條件共軛梯度區(qū)域網(wǎng)平差方法，包括：

S1、導(dǎo)入?yún)^(qū)域網(wǎng)平差計(jì)算需要的原始數(shù)據(jù)，且將導(dǎo)入的原始數(shù)據(jù)進(jìn)行時(shí)空基準(zhǔn)統(tǒng)一，得到初始數(shù)據(jù)，所述原始數(shù)據(jù)至少包括初始內(nèi)外方位元素?cái)?shù)據(jù)以及點(diǎn)位數(shù)據(jù)；

S2、進(jìn)入?yún)^(qū)域網(wǎng)平差迭代流程：分別讀取每一個(gè)像點(diǎn)對應(yīng)的初始數(shù)據(jù)，計(jì)算每一個(gè)像點(diǎn)對應(yīng)的法方程常數(shù)項(xiàng)向量c的一個(gè)分量，并將各個(gè)法方程常數(shù)項(xiàng)向量c的分量累加，得到完整的法方程常數(shù)項(xiàng)向量c作為步驟S4中殘差向量s的初值；

S3、分別讀取每一個(gè)像點(diǎn)對應(yīng)的初始數(shù)據(jù)，計(jì)算每一個(gè)像點(diǎn)對應(yīng)的預(yù)條件矩陣M的分量，并將該分量求逆并乘以方程常數(shù)項(xiàng)向量c得到矩陣-向量積M^-1c的一個(gè)分量，將各個(gè)矩陣-向量積M^-1c的分量累加，得到完整的矩陣-向量積M^-1c作為步驟S4中方向向量d的初值；

S4、進(jìn)入預(yù)條件共軛梯度法迭代流程：對未知數(shù)改正數(shù)向量u、殘差向量s以及方向向量d進(jìn)行初始化；

S5、分別讀取每一個(gè)像點(diǎn)對應(yīng)的初始數(shù)據(jù)，計(jì)算每一個(gè)像點(diǎn)對應(yīng)的法方程系數(shù)矩陣B的與方向向量d的矩陣-向量積Bd分量，將各矩陣-向量積Bd分量累加，得到完整的矩陣向量積Bd；

S6、根據(jù)預(yù)條件共軛梯度算法，以及本次迭代中的未知數(shù)改正數(shù)向量u，殘差向量s，方向向量d，預(yù)條件矩陣M以及法方程系數(shù)矩陣B，計(jì)算新的未知數(shù)改正數(shù)向量u、新的殘差向量s以及新的方向向量d；

S7、根據(jù)不精確牛頓解法，計(jì)算本次迭代的強(qiáng)制序列系數(shù)η；

S8、判斷此次預(yù)條件共軛梯度迭代是否符合預(yù)定收斂條件，若符合，則結(jié)束預(yù)條件共軛梯度迭代，輸出新的未知數(shù)改正數(shù)向量u，執(zhí)行步驟S9；否則，返回步驟S5；

S9、判斷此次區(qū)域網(wǎng)平差迭代是否符合預(yù)定收斂條件，若符合，則結(jié)束區(qū)域網(wǎng)平差迭代，執(zhí)行步驟S10，否則，返回步驟S2；

S10、根據(jù)步驟S8中輸出的新的未知數(shù)改正數(shù)向量u，更新所有的未知數(shù)數(shù)值，并輸出所有的未知數(shù)數(shù)值。