李花,秦亮軍
(廣州市城市規劃勘測設計研究院,廣東廣州 510060)
采用衛星遙感方法制作的數字正射影像是獲取基礎空間信息最快速、高效的手段,利用高精度衛星遙感成圖,處理速度快、工藝簡便、圖像清晰,而且精度達到成圖要求。自2001年10月美國DigitalGlobe公司的QuickBird遙感衛星發射升空以來,遙感圖像空間分辨率極高和良好的地物表現能力,廣泛用于土地規劃、城市規劃設計、自然資源管理、環境監測(如空氣污染、油污擴散、洪水災區趨勢分析等)、農作物判別、生態調查、土地利用變化分析、景觀制圖與高精度地圖繪制等各個領域[1]。影像圖集往往也是以衛星影像數據為主要數據源,本文以《廣州市開發區蘿崗區數字影像與規劃圖集》為例,對影像處理的過程及方法進行詳細介紹。
影像圖集是一種集精準、美觀、實用于一體的藝術品,對影像的要求極高。影像處理過程主要是將影像基礎產品處理成為變形小、信息豐富、地理精度高的過程。為了在視覺上改善圖像目視效果,提高分辨率,使目標細部特征更清晰,光譜信息更豐富,灰度范圍更廣,信息量更大,必須首先對基礎產品數據進行數據融合。高分辨率遙感衛星多采用線陣列CCD探測器,按照推掃式掃描成像,同時獲取高分辨率全色和多光譜影像。在遙感成像過程中,由于衛星姿態不穩、傳感器機械原因等因素的影響,在山區特別是地形高差起伏較大的地區,由于中心投影的影響(地面起伏)引起像點位移,使得遙感圖像存在一定的幾何變形,需要對其進行恢復處理。因此,必須對影像進行高精度的正射糾正及幾何糾正。整個處理過程如圖1所示。同時考慮到此次研究前期收集的大部分電子資料都是基于廣州獨立坐標系的,因此最終的成果數據必須全部納入廣州獨立坐標系統下。

圖1 影像處理流程圖
QuickBird高分辨率遙感衛星是目前世界上分辨率最高的商用遙感衛星,傳感器采用多光譜成像(1個全色通道和4個多光譜通道)。
本文研究區域為整個蘿崗區,面積約為393 km2,數據量大(約8 G)。整個區域地勢南北狹長,地勢南低北高,北部的帽峰山是蘿崗區第一高峰,主峰海拔534.9 m。
QuickBird數據是直接購買所得,接收時間為2006年10月~11月,天氣晴朗,數據質量較好。原始數據塊圖情況如圖2所示。

圖2 原始數據塊圖情況
QuickBird多光譜波段分辨率為2.44 m,其地物的分辨能力未能滿足人們日常判讀的需求,為了對觀測目標有一個更加全面、清晰、準確的理解與認知,為了改善解讀的精度、可靠性,將同一區域的多源遙感圖像按統一的坐標系統,通過空間配準和內容復合,生成一幅比單一信息源更準確、更完全、更可靠的新圖像。它提高了影像的空間分解力和清晰度,提高了影像的平面測圖精度、分類精度及可靠性,增強了影像的解譯和動態監測能力,有效提高了遙感影像數據的利用率等[2]。
本文采用IHS變換融合法。先把多光譜影像利用IHS變換(如式1),從RGB(三原色)空間系統變換到IHS(表色系統)空間,同時將單波段的高分辨率圖像經過灰度拉伸,使其灰度的均值與方差和IHS空間中的亮度分量圖像相似,以增大兩影像之間的某種相關性,并抑制光斑效應;然后將拉伸過的全色影像作為新的亮度分量代入到IHS,經過RGB變換(如式2)還原到原始空間。

式中:I表示顏色的亮度;H表示色度;S表示飽和度[3]。
衛星影像在成像的過程中,受到透視投影、攝影軸傾斜、大氣折射、地球曲率及地形起伏等諸多因素影響,致使影像中各像點產生不同程度的幾何變形而失真。正射校正圖像不僅可以消除成像過程中的各種因素導致的影像畸變,而且包含了比線劃圖更豐富和更直觀生動的信息。因而,正射影像制作已經成為各種遙感應用過程中必不可少的基礎工作[4]。

圖3 DEM數據
進行影像正射校正前必須制作整個區域的DEM。首先利用廣州市城市規劃勘測設計研究院1∶500基礎地形圖數據提取已賦高程值的等高線及高程點數據分塊構造不規則三角網(TIN),內插得到格網DEM數據,并轉換到GRID格式,得到正射校正用的DEM數據如圖3所示。由于基礎地形數據量非常大,因此整個DEM的生成過程必須分塊進行,因此后期必須進行圖像拼接。圖像拼接處理是將具有地理參考的若干相鄰圖像合并成一幅圖像或者一組圖像[5]。在進行正射校正前必須對廣州獨立坐標系下的DEM數據進行坐標轉換,將其轉為UTM坐標。
為提高DOM的精度,此次控制點的采集我們選用野外GPS采點,作業中采取GPS聯測的方法,每景數據控制點采集28個點左右,以便分析后進行取舍。利用配套軟件進行GPS控制點坐標的平差分析,分析結果表明控制點精度滿足要求。然后利用Erdas軟件正射校正模塊,采用共線方程糾正法進行正射糾正,生成整個蘿崗區的DOM。
在衛星遙感的各種應用中,特別是在專題圖集應用中還需要疊加現狀及規劃要求的情況下,為了有效地使用衛星影像,幾何糾正是必不可少的步驟。影像幾何糾正的方法有許多種,主要可分成兩大類:數學擬合模型,不管其物理意義,直接以地面控制點對影像進行數學變換,本法一般適用于平坦地區,屬于較不嚴密解法;幾何成像模型,考慮成像時造成影像變形的物理意義,再利用這些物理條件來構建成像幾何模型,屬于嚴密解法[6]。
幾何成像模型是最嚴格的影像糾正模型,但算法比較復雜,對輔助數據(如星歷參數、DEM數據等)要求較高[6]。本次研究利用的輔助數據主要為廣州市城市規劃勘測設計研究院的1∶500基礎地形圖數據,完全符合幾何成像模型的要求。
此次研究采用了幾何成像模型,選取1∶500地形圖中336個控制點結合高精度DEM成果數據進行幾何精糾正,控制點分布如圖4所示。
糾正前后衛星影像的拼接效果如圖5所示。在ArcGIS中同時加載糾正后的衛星影像和1∶500基礎地形數據進行檢驗,其匹配效果如圖6所示。由以上兩圖明顯看出基本滿足無縫拼接的同時,與現狀地形的套合情況也是比較好的。

圖4 蘿崗區幾何糾正控制點分布圖

圖5 航帶間的拼接效果對比

圖6 與1∶500地形圖疊加效果

圖7 精度檢驗點分布圖
將處理完成的遙感圖像與1∶500地形圖疊加,并在影像圖上均勻選取了94個點位,與1∶500地形圖上相應點位坐標進行對比計算其點位誤差。其選點分布情況如圖7所示,點位覆蓋北部的九龍鎮與南部各個街道,較差中誤差為0.965 m。完全滿足圖集疊加的各種現狀要素、規劃要素的要求,其統計的精度也比較高。
[1]BESSETTES V,DESACHY J.Using directional variance for urban area analysis on QuickBird 5 panchromatic images[J].Image and Signal Processing for Remote Sensing,2002(4):288~296
[2]賈永紅,李德仁,孫家柄.多源遙感影像數據融合[J].遙感技術與應用,2000,5(1):41~44
[3]王建梅,李德仁.QuickBird全色與多光譜數據融合方法用于土地覆蓋分類中的比較研究[J].測繪通報,2005,(10):37~43
[4]欒慶祖,劉慧平,肖志強.遙感影像的正射校正方法比較[J].遙感技術與應用,2007,22(6):743 ~748
[5]黨安榮,王曉棟,陳曉峰,等.ERDAS IMAG IN E遙感圖像處理方法[M].北京:清華大學出版社,2003,86~90
[6]馬洪超,趙向東.基于地形隨機場模型的遙感圖像幾何糾正[J].測繪學報,2006,35(3):251~254