曲莉莉，朱豐琪

(山東省國(guó)土測(cè)繪院，山東 濟(jì)南 250013)

隨著《測(cè)繪地理信息事業(yè)“十三五”規(guī)劃》[1]的實(shí)施，高分辨率數(shù)字正射影像圖(digital orthophoto map，DOM)在基礎(chǔ)測(cè)繪、國(guó)土調(diào)查、國(guó)情監(jiān)測(cè)、城市規(guī)劃、地籍測(cè)量等項(xiàng)目中發(fā)揮了重要作用，對(duì)高分辨率的研究顯得十分必要。研究表明[2-5],在航空正射影像生產(chǎn)過(guò)程中，影響正射影像質(zhì)量的重要因素是數(shù)字高程模型(digital elevation model，DEM)編輯和DOM，并提出了多種處理方式方法。 渠甲源[6]曾就高分辨率航空正射影像技術(shù)設(shè)計(jì)難點(diǎn)探討了利用MicroStation軟件分層編輯DEM數(shù)據(jù)，但編輯效率及效果不佳，自動(dòng)化程度較低，需使用大量人工編輯量。目前也無(wú)較為成熟的針對(duì)大面積、多時(shí)相、不同天氣引起的航片色彩差異進(jìn)行整體均衡勻色的處理方法。隨著高科技測(cè)繪技術(shù)的不斷發(fā)展，傳統(tǒng)處理模式已逐漸被先進(jìn)的自動(dòng)化的集群式影像處理系統(tǒng)所替代，與傳統(tǒng)處理模式相比，自動(dòng)化集群式影像處理系統(tǒng)已有數(shù)倍的效率提升，但DEM編輯及DOM勻色這兩個(gè)環(huán)節(jié)的技術(shù)難點(diǎn)是目前數(shù)字正射影像快速生產(chǎn)急需解決的問(wèn)題，將極大地影響自動(dòng)化處理系統(tǒng)的影像生產(chǎn)效率及產(chǎn)品質(zhì)量。

1 主要技術(shù)指標(biāo)及儀器設(shè)備

1.1 主要技術(shù)指標(biāo)

平面坐標(biāo)系采用2000國(guó)家大地坐標(biāo)系(CGCS2000)，投影方式為高斯-克呂格投影，分帶方式為3°分帶，中央經(jīng)線117°。

DOM 成果應(yīng)無(wú)明顯拼接痕跡并保證建筑物等實(shí)體的影像完整，影像色彩應(yīng)接近真實(shí)自然，紋理應(yīng)清晰、無(wú)明顯失真。

1.2 儀器設(shè)備

本文使用的主要儀器設(shè)備及軟件系統(tǒng)主要有PC-6飛機(jī)、DMCIII航攝儀、 Inpho、GEOWAY CIPS、Pixel Factory等。詳見(jiàn)表1。

表1 主要設(shè)備與軟件

2 技術(shù)路線及工藝流程

圖1 航空?qǐng)D像采集和處理流程圖Fig.1 Flowchart of the acquisition and processing of aerial images

3 DEM 優(yōu)化研究

試驗(yàn)證明DEM的質(zhì)量是決定數(shù)字正射影像質(zhì)量的重要因素之一，DEM優(yōu)化質(zhì)量的好壞直接影響DOM精度。 因此，在生產(chǎn)中對(duì)DEM編輯工序必須投入較大精力，在保證DOM的精度質(zhì)量的前提下提高DEM編輯效率。

3.1 DEM難點(diǎn)分析

濱州市地勢(shì)南高北低，大致上由西南向東北傾斜，漸次過(guò)渡到大海。小清河以南的鄒平南部屬丘陵區(qū)，地勢(shì)高峻，其余均為山前傾斜平原，地勢(shì)平緩，海拔高程一般在8～800 m。小清河以北為黃河沖積平原，海拔高程一般在1～20 m，總體上地勢(shì)低平，呈現(xiàn)出濱州地區(qū)整體地貌較平坦但局部地勢(shì)高峻的地貌特征。因此，濱州項(xiàng)目存在以下難點(diǎn)問(wèn)題，影響DEM的精度。

(1)濱州城區(qū)內(nèi)高層建筑物多，航攝影像在航向重疊約60%左右，旁向重疊約30%左右。此類(lèi)重疊度在DSM密集匹配時(shí)高層建筑物會(huì)存在匹配漏洞，導(dǎo)致此類(lèi)建筑物在點(diǎn)云形態(tài)與建筑物形態(tài)匹配度不吻合，在對(duì)建筑物濾波時(shí)會(huì)存在漏濾現(xiàn)象。

濱州市周邊存在大量的工業(yè)園區(qū)，部分工業(yè)園區(qū)廠房較大，廠房?jī)?nèi)露天設(shè)施外型復(fù)雜。對(duì)于人工構(gòu)筑物較大及造型較為復(fù)雜等情況，點(diǎn)云濾波參數(shù)不宜設(shè)置過(guò)大，防止濾波時(shí)地貌特征損失過(guò)多。

(2)濱州市整體地勢(shì)較為平坦，基本以平原為主，但最南邊鄒平市存在丘陵地貌。對(duì)于此類(lèi)多種地貌特征采用同一過(guò)濾參數(shù)效果不佳。

濱州市整體以平原為主，航攝時(shí)間為2020年4月—6月，測(cè)區(qū)內(nèi)存在大面積的小麥種植區(qū)。測(cè)區(qū)北部沿海存在大量鹽田、海水養(yǎng)殖場(chǎng)等，大片紋理類(lèi)似地物在DSM匹配時(shí)會(huì)存在大量的錯(cuò)誤匹配點(diǎn)，這些錯(cuò)誤點(diǎn)在濾波時(shí)影響濾波效果。

(3)集群系統(tǒng)在密集匹配時(shí)可能會(huì)產(chǎn)生噪聲高點(diǎn)或低點(diǎn)，而自動(dòng)濾波只能有效過(guò)濾約80%的非地面點(diǎn)，其余非地面點(diǎn)需人工過(guò)濾。

3.2 DEM生產(chǎn)流程

對(duì)原始航空影像數(shù)據(jù)進(jìn)行空三加密，然后將空三加密成果導(dǎo)入到GEOWAY CIPS、Pixel Factory集群式影像處理系統(tǒng)中，在集群中采用逐像素密集匹配同名點(diǎn)的方式進(jìn)行數(shù)字地表模型(digital surface model，DSM)提取,然后經(jīng)DSM濾波、DSM矢量化、DEM編輯等工序，制作數(shù)字高程模型DEM。

(1)DSM匹配及濾波：利用空三成果，生成立體模型，對(duì)立體模型進(jìn)行核線影像密集匹配生成DSM，對(duì)DSM進(jìn)行植被及建筑物濾波，去除植被及房屋等處的點(diǎn)云。

(2)DSM矢量化：經(jīng)過(guò)濾波的DSM進(jìn)行構(gòu)TIN重采樣，生成DEM過(guò)程數(shù)據(jù)，再生成單片正射影像，用于后續(xù)編輯。

(3)DEM編輯：對(duì)單片正射影像進(jìn)行變形檢查，對(duì)相應(yīng)位置的DEM數(shù)據(jù)進(jìn)行編輯處理。除去濾波過(guò)程未過(guò)濾掉的高程異常值。

3.3 試驗(yàn)結(jié)果分析

在濾波前根據(jù)地貌地物特征，將不同地貌及地物區(qū)域分開(kāi)。如將平地與丘陵地區(qū)分開(kāi)，各自調(diào)整過(guò)濾參數(shù)進(jìn)行過(guò)濾。濾波參數(shù)將最大建筑物長(zhǎng)度、最大地表高差、迭代距離、是否進(jìn)行山脊填補(bǔ)設(shè)為關(guān)鍵參數(shù)值。其中平地過(guò)濾參數(shù)中的建筑物長(zhǎng)度可適當(dāng)放寬，默認(rèn)參數(shù)為60 m，可根據(jù)測(cè)區(qū)實(shí)際建筑長(zhǎng)度及寬度調(diào)整至120 m，丘陵地區(qū)建筑物長(zhǎng)度可適當(dāng)調(diào)小至40 m，即該測(cè)區(qū)如果平地有建筑物長(zhǎng)度超過(guò)120 m，丘陵地區(qū)有超過(guò)40 m的建筑，將不會(huì)進(jìn)行過(guò)濾；最大地表高差默認(rèn)參數(shù)為1.2 m，平地調(diào)整至0.6 m，丘陵調(diào)整至1.4 m，地表高差是判斷地形起伏程度的參數(shù)，平地變化小，丘陵地變化大；迭代距離默認(rèn)參數(shù)為1.8 m，平地可調(diào)整至0.8 m，丘陵調(diào)整至1.2 m，該參數(shù)主要是配合地表高差參數(shù)，形成濾波窗口尺寸；濾波時(shí)有時(shí)會(huì)將比較平坦的山脊過(guò)濾掉。平地?zé)o山脊，不進(jìn)行山脊填補(bǔ)。將濱州市鄒平區(qū)單獨(dú)分開(kāi)，采用丘陵地區(qū)的濾波參數(shù)進(jìn)行過(guò)濾，其余區(qū)域采用平地濾波參數(shù)過(guò)濾。同理將城區(qū)有大量高層建筑的區(qū)域與周邊鄉(xiāng)鎮(zhèn)及農(nóng)村區(qū)域分開(kāi)濾波。

濾波參數(shù)過(guò)濾效果成功率約在80%左右，濾波后仍有較多的高程異常值需要處理。針對(duì)這種情況，為減小人工干預(yù)工作量，將平地區(qū)域與丘陵區(qū)域分開(kāi)。平地經(jīng)濾波后，采用比較大的平滑窗口，對(duì)整體地形進(jìn)行平滑處理，可獲得較好的平滑結(jié)果，同時(shí)能避免一些地物輕微變形,見(jiàn)圖2。丘陵地區(qū)采用較小的平滑窗口，以免平滑損失過(guò)多的地形信息，從而影響后續(xù)正射影像糾正精度。針對(duì)不同地貌特征采用合理的平滑優(yōu)化，可極大程度減少異常高程值引起的地物變形，節(jié)省大量人工編輯工作量。

圖2 建筑物變形修正Fig.2 Building deformation correction

全自動(dòng)化處理時(shí)不可避免會(huì)存留部分粗差，采用人機(jī)交互的方式對(duì)DEM過(guò)程數(shù)據(jù)進(jìn)行編輯，解決建筑物扭曲變形、地形地貌失真等問(wèn)題。設(shè)置平地參數(shù)處理后的地區(qū)有可能造成山頭缺失，需要人工添加匹配點(diǎn)；大面積水域可手工沿水域的邊緣選擇編輯區(qū)域，輸入水面高程，用平面擬合算法進(jìn)行平滑處理；將異常錯(cuò)誤的點(diǎn)云手動(dòng)人工刪除，手動(dòng)添加缺失的特征點(diǎn)，使高程網(wǎng)準(zhǔn)確貼合地面。

4 DOM勻色優(yōu)化研究

4.1 DOM難點(diǎn)問(wèn)題分析

濱州測(cè)區(qū)航攝分區(qū)6個(gè)，原始航攝數(shù)據(jù)量大小為3.4 TB。航攝時(shí)間跨度為2020年4月中旬至6月底。測(cè)區(qū)南北狹長(zhǎng)，整體為平原地貌，且北臨渤海灣。濱州地區(qū)4月至5月份以晴朗天氣為主，6月份天氣逐漸向多云多霧變化。以上情況會(huì)導(dǎo)致下列問(wèn)題。

(1)因航攝時(shí)間跨度較長(zhǎng)，地表植被覆蓋生長(zhǎng)情況差異明顯，造成各架次影像色調(diào)嚴(yán)重不一致，在影像鑲嵌時(shí)不同架次接邊處相同地物、植被存在不同的色彩及紋理特征，影響影像接邊質(zhì)量及效果。

(2)測(cè)區(qū)北部臨渤海海域，存在大面積鹽池及海水養(yǎng)殖場(chǎng)，水面整體反光嚴(yán)重，局部存在因強(qiáng)反光產(chǎn)生的曝光現(xiàn)象，影響水體圖面效果。

(3)部分架次影像存在輕度霧氣，整體視覺(jué)通透性差，影像噪點(diǎn)多，地物紋理表現(xiàn)力差。

4.2 DOM制作流程

基于制作完成的DEM成果和空三加密成果，進(jìn)行單片正射糾正、勻光勻色、鑲嵌、圖面修復(fù)等工序，制作DOM。

(1)單片糾正：利用DEM成果和空三加密成果對(duì)濱州原始航攝影像進(jìn)行逐片微分正射校正。

(2)勻光勻色：調(diào)整濱州測(cè)區(qū)較為合適的勻色模板，對(duì)因跨時(shí)相、跨架次、天氣等引起的顏色差異基于糾正后的單片正射影像采用勻色模板進(jìn)行勻光勻色處理，使所有單片影像顏色趨于一致。

(3)影像鑲嵌：對(duì)勻色后的正射影像自動(dòng)化鑲嵌處理，對(duì)鑲嵌線切割各類(lèi)線狀地物、建筑物造成明顯錯(cuò)位等問(wèn)題進(jìn)行人工鑲嵌線編輯，使影像各要素完整、拼接痕跡不明顯、紋理清晰、城區(qū)高層建筑物投影方向一致，避免建(構(gòu))筑物相互壓蓋、交錯(cuò)、反向等現(xiàn)象。

(4)圖面修復(fù)：對(duì)標(biāo)準(zhǔn)分幅影像成果進(jìn)行人工檢查，對(duì)因DEM編輯不當(dāng)造成的影像變形問(wèn)題以及其他各類(lèi)圖面問(wèn)題進(jìn)行修復(fù)。

4.3 試驗(yàn)結(jié)果分析

首先利用GEOWAY CIPS 對(duì)測(cè)區(qū)所有單片先進(jìn)行第一次勻色處理，使每張影像的色彩色調(diào)趨于一致。然后將勻色后的單片導(dǎo)入到Pixel Factory中進(jìn)行鑲嵌和二次勻色，提高了影像成果的質(zhì)量和顯示效果，方法如下。

(1)合理選擇勻色模板。時(shí)相或影像質(zhì)量一致或接近的為一組，采用同一勻色模板；存在輕度霧氣的數(shù)據(jù)選出幾張典型影像，微調(diào)勻色模板測(cè)試勻色效果直至影像地物紋理清晰、層次分明、色彩均衡一致、直方圖基本呈正態(tài)分布為止，然后用調(diào)試好的勻色模板對(duì)所有輕度霧氣影像進(jìn)行勻色；大面積水域反光或有曝光現(xiàn)象的影像降低勻色模板亮度后進(jìn)行勻色。

(2)在Pixel Factory中編輯鑲嵌線，使其沿塊狀或線狀地物的邊沿，避免切割同一地物，同時(shí)設(shè)置合理的羽化參數(shù)，過(guò)渡鑲嵌效果，此步驟可進(jìn)一步優(yōu)化基于第一步勻色后的影像勻色效果；對(duì)于大面積鹽池和水域反光，盡量將鑲嵌線從鹽池或養(yǎng)殖場(chǎng)的堤岸處過(guò)渡，避免鑲嵌線切割水面，造成陰陽(yáng)水面。

(3)鑲嵌線編輯完畢后將整個(gè)濱州測(cè)區(qū)導(dǎo)出一個(gè)低分辨率、數(shù)據(jù)量小的鑲嵌影像模板，在Photoshop圖像處理軟件對(duì)其調(diào)整色彩色調(diào)，然后將其重新導(dǎo)入至Pixel Factory中，將其顏色以地理映射的方式應(yīng)用于全局，以小帶大的方式獲得整個(gè)測(cè)區(qū)最佳勻色效果,見(jiàn)圖3。

圖3 測(cè)區(qū)勻色前和勻色后Fig.3 Before and after image uniformity of testing zone

鑲嵌線沿街道、公路、河流線狀地物進(jìn)行編輯，盡量避免鑲嵌線從線狀地物中間穿過(guò)分割建筑物、山脊等，使影像鑲嵌后色彩自然合理、地物完整。

4 結(jié)論

本文采用較先進(jìn)的技術(shù)方案，利用集群式影像處理系統(tǒng)，針對(duì)正射影像生產(chǎn)流程中對(duì)精度及效果影響較大的DEM編輯及影像勻光勻色關(guān)鍵性技術(shù)進(jìn)行研究，得出了如下結(jié)論：

(1)DEM編輯處理中采用按地貌類(lèi)別分區(qū)濾波、合理使用平滑優(yōu)化參數(shù)、人機(jī)交互編輯等方法，在保證DEM質(zhì)量的同時(shí)最大化自動(dòng)化處理，減少了大量人工干預(yù)環(huán)節(jié)，節(jié)省人力成本和時(shí)間。處理后的DEM成果平面精度及高程精度滿足后續(xù)DOM生產(chǎn)要求，即無(wú)明顯高程異常值，地形表述與實(shí)際地形一致，為后續(xù)生產(chǎn)提供了質(zhì)量較優(yōu)的DEM數(shù)據(jù)。

(2) DOM處理中的影像勻光勻色環(huán)節(jié)按照影像質(zhì)量合理使用不同的勻色模板進(jìn)行初次勻色，根據(jù)測(cè)區(qū)地物合理編輯鑲嵌線及羽化效果，采用地理映射的方式對(duì)整測(cè)區(qū)進(jìn)行二次勻色。勻色后的影像地物細(xì)節(jié)清晰，反差適中，層次分明，色彩均衡，鑲嵌線兩側(cè)和相鄰影像色調(diào)應(yīng)基本一致。

通過(guò)以上試驗(yàn)方法制作完成DOM，平面精度及色彩紋理質(zhì)量均達(dá)到了預(yù)期的0.2 m分辨率的標(biāo)準(zhǔn)要求。試驗(yàn)方法合理有效，基于目前市面最先進(jìn)的集群式影像處理系統(tǒng)，進(jìn)一步提升生產(chǎn)了DEM編輯自動(dòng)化程度、減少了人工干預(yù)；勻光勻色消除因原始影像引起的質(zhì)量問(wèn)題，形成一套快速規(guī)?；瘮?shù)字正射影像處理方案，能有效指導(dǎo)后續(xù)此類(lèi)項(xiàng)目的生產(chǎn)工作。