福建省測(cè)繪院，福建 福州 350003

福建省測(cè)繪地理信息局于2014年購(gòu)置了徠卡ALS70-HP激光掃描系統(tǒng)，該系統(tǒng)同時(shí)搭載了徠卡RCD30八千萬(wàn)像素中幅面量測(cè)相機(jī)，航攝時(shí)可同機(jī)獲取高密度高精度的點(diǎn)云數(shù)據(jù)和高分辨率的航攝像片，用于制作全省高精度DEM和DOM。在實(shí)際生產(chǎn)DOM的過(guò)程中，出現(xiàn)了以下三個(gè)較為突出的問(wèn)題：傳統(tǒng)的TerraPhoto軟件制作DOM效率低下且成圖質(zhì)量差、空三加密成果在航攝分區(qū)及加密分區(qū)處存在接邊誤差、原始數(shù)據(jù)資料無(wú)法導(dǎo)入國(guó)產(chǎn)攝影測(cè)量平臺(tái)。文章著重分析了上述問(wèn)題出現(xiàn)的原因，并提出優(yōu)化改進(jìn)方案。

1 傳統(tǒng)作業(yè)方式改進(jìn)

1.1 傳統(tǒng)作業(yè)方式存在的問(wèn)題

傳統(tǒng)的Lidar數(shù)據(jù)生產(chǎn)正射影像是在TerraPhoto軟件中，通過(guò)大量的航片連接點(diǎn)（Tie點(diǎn)）來(lái)計(jì)算并修正航攝相機(jī)安裝角度與IMU系統(tǒng)間的偏差，并統(tǒng)一改正到航片的EO中，再用點(diǎn)云濾波后生成的DEM糾正出正射單片，該過(guò)程無(wú)須外業(yè)控制點(diǎn)。由于TerraPhoto軟件并不是專門(mén)的攝影測(cè)量軟件，它自動(dòng)匹配出的Tie點(diǎn)在數(shù)量和質(zhì)量上與專業(yè)攝影測(cè)量軟件（SSK、INPHO、LPS等）相差較大，特別是在大片植被覆蓋的地方，其自動(dòng)匹配弱的問(wèn)題更加突出，匹配出的點(diǎn)很少，且精度差，幾乎不能用，需要手工量測(cè)大量的Tie點(diǎn)。福建地區(qū)大部分為植被覆蓋區(qū)域，用TerraPhoto軟件生產(chǎn)DOM效率低下，而直接用初始EO及點(diǎn)云濾波后生成的DEM對(duì)原始航片進(jìn)行糾正，發(fā)現(xiàn)局部相鄰正射單片接邊存在少許偏差，說(shuō)明初始EO不是絕對(duì)準(zhǔn)確。

1.2 專業(yè)攝影測(cè)量軟件與點(diǎn)云濾波相結(jié)合的優(yōu)化改進(jìn)方案

替代方案是將航片及初始EO導(dǎo)入SSK中，SSK軟件強(qiáng)大的特征匹配能力使其能在植被茂密區(qū)域匹配到大量分布均勻、可靠的連接點(diǎn)，如圖1所示。

圖1 植被茂密區(qū)域匹配出大量可靠連接點(diǎn)

用初始EO直接進(jìn)行絕對(duì)定向，該過(guò)程同樣無(wú)須外業(yè)控制點(diǎn)，且?guī)缀醪挥檬止ち繙y(cè)連接點(diǎn)，全部由計(jì)算機(jī)自動(dòng)計(jì)算完成。在同等硬件條件下，單個(gè)測(cè)區(qū)的劃分也比TerraPhoto大很多，可以達(dá)到4000片，而TerraPhoto單個(gè)測(cè)區(qū)達(dá)到500片時(shí)，匹配及平差計(jì)算就非常困難。

絕對(duì)定向后輸出新的EO，將其與初始EO進(jìn)行對(duì)比，發(fā)現(xiàn)其存在細(xì)微的變化。用空三后輸出的EO和點(diǎn)云濾波后生成的DEM糾正原始影像，經(jīng)過(guò)檢查發(fā)現(xiàn)相鄰正射單片接邊不存在偏差，且與舊圖件資料DLG或DOM套合情況良好，說(shuō)明空三后輸出的EO是可靠的。前后EO間的細(xì)微變化，正是通過(guò)大量匹配點(diǎn)整區(qū)平差計(jì)算出的初始EO偏差值。

在正射影像的糾正過(guò)程中，發(fā)現(xiàn)用點(diǎn)云濾波后生成的DEM來(lái)糾正航片，在植被覆蓋茂密區(qū)域糾正出來(lái)的效果并不好。分析其原因，主要是因?yàn)辄c(diǎn)云濾波后是用地面點(diǎn)構(gòu)建的DEM，而DOM在植被區(qū)域顯示出的只是樹(shù)冠部分，并不是地面，當(dāng)樹(shù)很高時(shí)，樹(shù)冠部分的高程與地面的高程相差較大；航攝像片邊緣處的高大地物由于光線入射角的原因，自身存在較大的投影差，再用地面處高程值去糾正樹(shù)冠處的像素會(huì)產(chǎn)生“堆積擠壓”現(xiàn)象，在DOM圖面上就表現(xiàn)為樹(shù)木拉升變形，該現(xiàn)象在遠(yuǎn)離投影中心的航片四周表現(xiàn)得尤為明顯。

福建省大部分為山區(qū)，植被覆蓋率高，用點(diǎn)云濾波方式獲取的DEM糾正的正射單片將會(huì)產(chǎn)生大量的樹(shù)木拉伸現(xiàn)象，最終采取的改進(jìn)方案是將加密成果導(dǎo)入INPHO，用Match-T模塊將匹配參數(shù)設(shè)置成山區(qū)，自動(dòng)匹配出DSM，自動(dòng)濾波，生成DEM。由于Inpho采用基于影像像素的特征匹配算法，其匹配濾波獲取的DEM在山區(qū)不會(huì)緊貼地表，且過(guò)渡均勻、自然，糾正出的正射單片樹(shù)木基本不會(huì)拉伸。

因?yàn)樵贛atch-T模塊中設(shè)置的匹配參數(shù)是山區(qū)，其匹配濾波獲取的DEM在平坦區(qū)域效果不理想，糾正時(shí)房子會(huì)產(chǎn)生明顯的變形，而點(diǎn)云濾波方式獲取的DEM在平坦區(qū)域效果較好，糾正時(shí)房子基本不會(huì)變形，剛好與INPHO匹配方式獲取的DEM互補(bǔ)，所以還需要用點(diǎn)云濾波方式獲取的DEM糾正出一份正射單片。由于測(cè)區(qū)大部分為山地，在DOM修補(bǔ)時(shí)，用INPHO匹配方式獲取的DOM成果作為底圖，碰到房子變形的地方，用點(diǎn)云濾波方式糾正出的正射單片對(duì)房子進(jìn)行修補(bǔ)。另外，在一些比較困難的區(qū)域，如山坡上的密集房屋、立交橋、陡峭山崖等兩種DEM糾正效果都不好的地方，需要將加密成果導(dǎo)入攝影測(cè)量平臺(tái)中，恢復(fù)立體像對(duì)，手工編輯困難區(qū)域的DEM，糾正出正射單片，再進(jìn)行修補(bǔ)。

后續(xù)的生產(chǎn)項(xiàng)目證明，采用上述改進(jìn)方案生產(chǎn)DOM，相比傳統(tǒng)的TerraPhoto生產(chǎn)方式，人工干預(yù)更少，極大地提高了生產(chǎn)效率，且成圖質(zhì)量更有保障，其在山地較多的測(cè)區(qū)效果尤為明顯。

2 航線敷設(shè)及加密分區(qū)改進(jìn)

2.1 航攝分區(qū)及加密分區(qū)間的接邊誤差

在生產(chǎn)過(guò)程中發(fā)現(xiàn)，在測(cè)區(qū)的內(nèi)部，正射單片是準(zhǔn)確的；在測(cè)區(qū)的邊緣，每條航線的頭尾幾片及第一條航線和最后一條航線的正射單片不夠準(zhǔn)確，且越靠近測(cè)區(qū)邊緣，精度越差。這是因?yàn)樵谧隹杖用軙r(shí)，分區(qū)是按照航線來(lái)劃分的，分區(qū)間沒(méi)有重疊航線；航飛時(shí)測(cè)區(qū)間重疊只有2～3片，而航線頭尾幾片受飛機(jī)掉頭時(shí)姿態(tài)急劇變化的影響，自身存在較大的角度偏差；在測(cè)區(qū)的邊緣沒(méi)有任何外業(yè)控制點(diǎn)來(lái)約束像片的定向，純粹靠航片初始EO和匹配出的加密點(diǎn)來(lái)進(jìn)行絕對(duì)定向，而自動(dòng)匹配加密點(diǎn)時(shí)在測(cè)區(qū)邊緣缺少足夠的航線間連接點(diǎn)，也缺少傳統(tǒng)空三加密中布設(shè)的同名控制點(diǎn)來(lái)約束測(cè)區(qū)的接邊偏差，就會(huì)造成測(cè)區(qū)邊緣的航片定向不夠準(zhǔn)確，在制作DOM時(shí)就表現(xiàn)為測(cè)區(qū)之間的正射單片接邊出現(xiàn)偏差。

2.2 航線敷設(shè)和加密分區(qū)改進(jìn)

采用的解決方案是在徠卡MissionPro飛行計(jì)劃軟件中，將相鄰航攝分區(qū)間的重疊像片設(shè)置為6～8張，并在空三加密時(shí)將相鄰加密分區(qū)間重疊兩條航線，最后輸出加密成果時(shí)，去掉加密分區(qū)中每條航線頭尾的3～4張航片，并去掉加密分區(qū)的第一條和最后一條航線，如圖2所示。這樣就能保證輸出的每張像片都有足夠的定向點(diǎn)，最終的接邊也不會(huì)產(chǎn)生較大誤差。

圖2 加密分區(qū)策略及輸出定向成果策略

按上述要求對(duì)飛行計(jì)劃進(jìn)行修改，并進(jìn)行加密分區(qū)和定向成果輸出，基本可以消除無(wú)外控條件下分區(qū)間空三成果的接邊差。

3 原始數(shù)據(jù)優(yōu)化改進(jìn)

3.1 原始數(shù)據(jù)無(wú)法導(dǎo)入國(guó)產(chǎn)攝影測(cè)量平臺(tái)的原因

在DOM的生產(chǎn)中，在一些復(fù)雜的地形處，如山坡上的密集房屋、立交橋、陡峭山崖等容易變形的地方，需要在立體模型上對(duì)DEM進(jìn)行特殊編輯，重新糾正出正射單片進(jìn)行修補(bǔ)。但國(guó)內(nèi)大部分的測(cè)圖和DEM編輯工作都是在國(guó)產(chǎn)攝影測(cè)量平臺(tái)上進(jìn)行，如適普、航天遠(yuǎn)景、JX4、吉威DPS等，而Lidar系統(tǒng)在航飛時(shí)的安裝角度是90°，即獲取的原始影像是“橫向”的，無(wú)法在國(guó)產(chǎn)攝影測(cè)量平臺(tái)上生成核線影像來(lái)恢復(fù)立體模型進(jìn)行立體測(cè)圖及DEM編輯。因此，只能利用原始影像構(gòu)建立體像對(duì)進(jìn)行立體測(cè)圖，而原始影像立體像對(duì)視差較大，且無(wú)法將空三后的加密點(diǎn)導(dǎo)入立體模型進(jìn)行輔助匹配，使山區(qū)等紋理單一區(qū)域很難匹配出質(zhì)量高的DSM，也無(wú)法在立體模型上對(duì)易變形區(qū)域DEM進(jìn)行特殊編輯。

目前，大部分的航攝相機(jī)都是采用90°安裝，自帶GPS但沒(méi)有帶慣導(dǎo)，獲取的初始外方位元素只含線元素，沒(méi)有角元素，拍攝出的原始航片為橫向的，在后期處理時(shí)都是將航攝像片向正北方向旋轉(zhuǎn)，即定義一個(gè)方向飛行的像片順時(shí)針轉(zhuǎn)90°，相向飛行的逆時(shí)針轉(zhuǎn)90°。因?yàn)樵谛D(zhuǎn)過(guò)程中并沒(méi)有改變投影中心，所以線元素并沒(méi)有改變，不影響后期的空三使用，只需要重新定義旋轉(zhuǎn)后的相機(jī)參數(shù)即可（將原始相機(jī)參數(shù)中的相幅長(zhǎng)寬交換，按照相機(jī)檢校說(shuō)明文件修改PPA，其他不變），它是通過(guò)初始GPS數(shù)據(jù)輔助外業(yè)控制點(diǎn)來(lái)解算像片的精確外方位元素。

3.2 導(dǎo)入國(guó)產(chǎn)攝影測(cè)量平臺(tái)對(duì)原始影像和初始EO進(jìn)行處理

徠卡ALS70激光掃描系統(tǒng)自帶高精度GPS和慣性導(dǎo)航設(shè)備，其獲取的初始EO已包含線元素和角元素，且相對(duì)比較精準(zhǔn)，不需要外業(yè)控制點(diǎn)來(lái)做后期的空三，只需要利用設(shè)備獲取的初始EO自動(dòng)匹配出加密點(diǎn)進(jìn)行絕對(duì)定向，以消除航攝相機(jī)安裝角度與IMU系統(tǒng)間的微小系統(tǒng)差。旋轉(zhuǎn)過(guò)程中雖然沒(méi)有改變線元素，但明顯改變了角元素，因此需要知道在航片旋轉(zhuǎn)過(guò)程中角元素發(fā)生了怎樣的變化，并總結(jié)出規(guī)則，以進(jìn)行相應(yīng)的修改。

選擇相向飛行的兩條航線作為試驗(yàn)區(qū)域，并按要求布設(shè)外業(yè)控制點(diǎn)。先不旋轉(zhuǎn)，在SSK中，導(dǎo)入初始EO進(jìn)行影像匹配，再直接用初始EO進(jìn)行絕對(duì)定向，并輸出最終EO1；再將像片向正北方向旋轉(zhuǎn)，并記錄旋轉(zhuǎn)方向，在SSK中重新定義旋轉(zhuǎn)后的相機(jī)參數(shù)，量測(cè)外業(yè)控制點(diǎn)，直接用控制點(diǎn)進(jìn)行絕對(duì)定向，輸出最終EO2。將EO1和EO2進(jìn)行比對(duì)，可以發(fā)現(xiàn)航片旋轉(zhuǎn)后，除了相機(jī)參數(shù)需要重新定義，外方位元素中只有KAPPA角發(fā)生了加減90°的變化，其他項(xiàng)并未發(fā)生明顯變化。

采用同樣的方法對(duì)不同方向飛行的航線做實(shí)驗(yàn)，并分析總結(jié)出航片旋轉(zhuǎn)時(shí)EO的變化規(guī)律，最終結(jié)果如表1所示。

表1 不同飛行方向下的像片旋轉(zhuǎn)及初始EO修改規(guī)則

批量將航片及初始EO按上述規(guī)則進(jìn)行處理，并重新定義相機(jī)參數(shù)，做空三處理，就可以將空三成果導(dǎo)入國(guó)產(chǎn)攝影測(cè)量軟件中，生成核線影像來(lái)恢復(fù)立體像對(duì)，進(jìn)行立體測(cè)圖及DEM編輯等工作。

4 結(jié)束語(yǔ)

文章系統(tǒng)地分析了基于Lidar數(shù)據(jù)快速生產(chǎn)DOM中出現(xiàn)的典型問(wèn)題，并提出了改進(jìn)方案。通過(guò)改進(jìn)的工藝流程，短時(shí)間內(nèi)完成了福建省12.4萬(wàn)km2無(wú)外業(yè)控制點(diǎn)空三加密及1∶5000 DOM生產(chǎn)任務(wù)，并順利通過(guò)了質(zhì)檢部門(mén)的驗(yàn)收，且該DOM成果已經(jīng)被當(dāng)作福建省多批次衛(wèi)星影像生產(chǎn)項(xiàng)目（包括第三次全國(guó)土地調(diào)查和福建省年度地理國(guó)情監(jiān)測(cè)）的參考資料來(lái)使用，說(shuō)明該改進(jìn)工藝是可靠的，尤其是在山地較多的測(cè)區(qū)值得推廣。