遙感圖像道路信息提取的仿真分析研究

宋偉 賀偉 鹿旭升 崔紅婕

摘? 要：隨著遙感技術的飛速發展，高分辨率遙感圖像變得極易獲得，而提取遙感圖像中道路信息更是具有極大的現實意義。從遙感圖像中提取的道路信息是結合直方圖閾值分割技術和數學形態學研究出一種半自動的遙感圖像道路提取方法。通過觀察遙感圖像的灰度直方圖，選擇出合適的閾值并分割出初始道路信息。初始道路會包含大量的顆粒噪聲且道路會出現間斷等情況。利用中值濾波器能夠濾除原始道路圖像中的顆粒噪聲，利用數學形態學的方法可以消除物體間的粘連和填充道路孔洞。通過仿真分析，驗證了研究的遙感圖像道路提取方法和有效性。

關鍵詞：遙感圖像;閾值分割;數學形態學;道路提取

中圖分類號：TP751.1 文獻標識碼：A 文章編號：2095-2945（2019）25-0018-04

Abstract： With the rapid development of remote sensing technology， high-resolution remote sensing images become easy to obtain， and the extraction of road information from remote sensing images is of great practical significance. The road information extracted from remote sensing image is combined with histogram threshold segmentation technology and mathematical morphology to study a semi-automatic road extraction method from remote sensing image. By observing the gray histogram of remote sensing image， the appropriate threshold is selected and the initial road information is segmented. The initial road will contain a lot of particle noise and the road will be interrupted and so on. The median filter can be used to filter the particle noise in the original road image， and the mathematical morphology method can be used to eliminate the adhesion between objects and fill the road holes. Through the simulation analysis， the road extraction method and effectiveness of the remote sensing image are verified.

Keywords： remote sensing image; threshold segmentation; mathematical morphology; road extraction

1 概述

遙感技術具有采集數據速度快、實時性高、探測范圍大、不受時間和地域條件的限制等優勢。在近二十年的飛速發展中，遙感技術已廣泛應用于軍事偵察、海洋、地理、測繪、農業、水利、環境保護等方面，在國防建設和科學研究等方面發揮著重要的作用。隨著我國基礎設施的飛速發展，尤其是在加快信息化的進程中，我國的遙感技術得到了飛速的發展。在現代遙感技術應用中，由于我國高分辨率衛星遙感圖像發展速度快，應用遙感技術對地觀測的能力不斷提高，高分辨率遙感圖像的研究已成為研究熱點。早期的高分辨率遙感衛星主要應用于軍事，隨著遙感信息的圖像成本逐漸降低，遙感技術逐漸應用于商業。

在遙感圖像提供的信息中，道路信息無疑是非常重要的。隨著我國的經濟飛速發展，人民對于美好生活的標準逐漸提高，準確獲取實時交通信息和道路流量狀況是人民和交通部門最基本的需求。早期的道路信息獲取主要需要人工現場勘查，這需要投入大量的人力、物力，而且采集周期長。因此，如何利用遙感道路圖像來快速的提取道路顯得尤為重要。

2 遙感圖像道路特征

不同的物體往往有不同的特征，了解道路的基本特征有助于我們設計實驗從而設計遙感道路提取算法。遙感道路圖像的特征可分為幾何特征、輻射特征、拓撲特征、功能特征和上下文特征。

（1）幾何特征：道路一般呈長條形，而且道路的長度大于寬度。在一定范圍內，道路會呈現直線形狀，在彎道中，曲率有一定的限制。

（2）輻射特征：道路的灰度值與其他區域的灰度值有明顯的區別，所以圖像中的道路有明顯的邊緣。

（3）拓撲特性：道路往往呈現為網狀結構，一般是貫通的、連續的。

（4）功能特征：道路分為快速路、主干路、次干路和支路。主干路是城市道路的骨架，主要功能是交通運輸;次干路兼有服務功能;支路起疏散交通的作用。

（5）上下文特征：上下文特征分為局部和全局，局部上下文特征指的是道路兩旁的信息，如樹木。全局上下文特征指的是道路所處的區域，如鄉村、田野、城區等。

3 圖像處理基礎理論

3.1 灰度直方圖

灰度圖像的直方圖（簡稱灰度直方圖）是一種表示數字圖像中各級灰度值與其出現頻率關系的函數，可以看作是像素灰度值出現的概率函數圖。灰度直方圖的橫軸表示像素的灰度級別，范圍為0～255，共256個灰度級別;縱坐標表示各個灰度級別出現的概率。若灰度圖像的灰度級范圍為[0，L-1]，則灰度圖像的直方圖定義為

式中，rk是第k級灰度值，nk是灰度值為rk的像素個數，h（rk）是灰度圖像直方圖的離散函數。

在觀察彩色遙感圖像的灰度直方圖時，首先把彩色圖像轉化為灰度圖像，如圖1所示，然后再顯示灰度圖像的直方圖，如圖2所示。

3.2 閾值分割技術

閾值分割是圖像分割方法中較為簡單的一種，因其實現較為簡單且性能穩定而成為應用最廣泛的分割方法。設圖像的灰度直方圖的灰度級范圍為[0，L-1]，灰度級為k時，像素為nk，則該圖像的總像素數N可表示為

灰度級出現的概率為

目標與背景有較大差異時，圖像的灰度直方圖往往會出現兩個雙峰的情況，此時圖像可大概分為兩個部分。在此情況下，從背景中提取目標首先要確定閾值。在灰度直方圖中，左側的山峰對應背景，右側的山峰對應目標，所以要選取兩個山峰之間的谷底對應的灰度值T作為閾值。在進行直方圖分割時，判決條件如下：

（1）灰度值Ti大于或等于閾值T時，判決為目標;

（2）灰度值小于閾值T時，判決為背景。

即對于圖像f（x，y），基于閾值的圖像分割方法可以定義為

也可以定義為

其中，g（x，y）為分割出的二值圖像。

3.3 非線性平滑濾波圖像增強方法——中值濾波法

中值濾波是一種去除噪聲的方法，是非線性濾波器，能夠有效地保留圖像的邊緣信息。中值濾波的原理是利用中值濾波窗口，取窗口所覆蓋區域的所有灰度值中的中間值，用中間值代替所覆蓋區域中心像素點的灰度值。

中值濾波器可以定義為

式中，f（x，y）表示圖像中位于（x，y）點的灰度值，g（x，y）表示濾波窗口所覆蓋區域為A的中值濾波結果。

中值濾波法大致有三種實現方法：

（1）直接選擇特定窗口濾波;

（2）先使用小窗口濾波，再使用大窗口濾波;

（3）先使用一維濾波，再使用二維濾波。

利用中值濾波法的過程可簡單敘述如下，先選定窗口的形狀并根據窗口形狀確定窗口的中心像素在圖像上的重合方式，然后利用窗口進行掃描并且將掃描到灰度值按大小進行排列，然后找到中間灰度值，最后把中間灰度值賦值給窗口中心像素。

3.4 灰度形態學的基本運算

3.4.1 灰度腐蝕

假設灰度圖像用A表示，結構元素用B表示，則用結構元素B對灰度圖像A進行灰度腐蝕的過程為：

（1）將結構元素B的原點重疊在灰度圖像A的中心元素上。

（2）依次用灰度圖像A的中心元素減去B的各個元素并將計算結果放在對應位置上。

（3）將結構元素B的原點移動到與原始圖像A的中心元素相鄰的8個元素上進行（2）的相同操作，得到8個結果。

（4）取步驟（2）和步驟（3）所得到的9個結果的最小值，此最小值即為腐蝕結果。

按以上步驟依次計算，即可得到圖像A的腐蝕結果。

灰度腐蝕的運算是逐像素點進行的，求某像素點的腐蝕結果即求在該像素點和其附近的8個像素點與結構元素中的對應位置進行相減運算，運算對得到9個結果，其中的最小值即為該像素點的灰度腐蝕結果。

3.4.2 灰度膨脹

假設灰度圖像用A表示，結構元素用B表示，則用結構元素B對灰度圖像A進行腐蝕的過程如下：

（1）將結構元素B的原點重疊在灰度圖像A的中心元素上。

（2）依次用灰度圖像A的中心元素加上結構元素B的各個元素并進行相加計算，并將計算結果放到對應位置上。

（3）將結構元素的原點移動到與原始圖像的中心元素相鄰的8個元素上進行（2）的相同操作，得到8個結果。

（4）取步驟（2）和步驟（3）所得到的9個結果的最大值，此最大值即為膨脹結果。

（5）按以上步驟依次計算，即可得到圖像A的腐蝕結果。

灰度膨脹的運算是逐像素點進行的，求某像素點的膨脹結果即求在該像素點和其附近的8個像素點與結構元素中的對應位置進行相加運算，運算對得到9個結果，其中的最大值即為該像素點的灰度腐蝕結果。

3.4.3 灰度開運算

假設結構元素用b表示，灰度圖像用f表示，則用結構元素b對灰度圖像f進行灰度開運算可表示為

灰度開運算的過程為：在灰度圖像下方滑動結構元素，在每個像素點計算結構元素上的最高點，最高點的集合即為開運算的結果。

灰度開運算能夠去除灰度圖像中的較小的明亮細節，保持整體灰度不變的作用。

3.4.4 灰度閉運算

假設結構元素用b表示，灰度圖像用f表示，則用結構元素b對灰度圖像f進行灰度閉運算可表示為

灰度閉運算的過程為：在灰度圖像下方滑動結構元素，在每個像素點計算結構元素上的最低點，最低點的集合即為閉運算的結果。

灰度閉運算能夠去除圖像中比較暗的細節部分，能夠保持明亮區域。

4 基于閾值分割和形態學的道路提取方法

4.1 道路提取流程和步驟

本文研究的是直方圖和數學形態學結合的方法，流程圖如3所示。

具體操作步驟為：（1）將彩色遙感道路圖像轉換成灰度圖像。（2）顯示灰度直方圖，觀察灰度直方圖并選擇合適的閾值。（3）根據所選擇的閾值進行圖像分割，得到大致道路圖像。（4）選擇合適的中值濾波模板濾除初始道路圖像中的小顆粒噪聲。（5）利用形態學處理中值濾波后的圖像，粘連道路并去除非道路信息。（6）圖像疊加。本設計流程中，選擇合適的閾值是進行圖像分割的關鍵，而且選擇合適的中值濾波模板和選擇合適的結構元素也同樣重要。選擇閾值要觀察灰度直方圖或借助數學方法，選擇中值濾波模板和選擇結構元素需要借助實驗進行觀察實驗后的效果。

4.2 遙感圖像道路信息提取的具體方案

4.2.1 彩色遙感圖像灰度化

將彩色圖像轉化為灰度圖像通常有三種處理方法，分別為最大值法、平均值法和加權平均值法。假設彩色圖像用f（x，y）表示，灰度圖像用g（x，y）表示，每個像素的分量為R、G、B。

（1）最大值法用公式可以表示為

（2）平均值法用公式可以表示為

（3）加權平均值法可以表示為

4.2.2 閾值分割提取初始道路

觀察圖4中原始圖像的灰度直方圖，選擇合適的閾值，把初始道路提取出來。灰度直方圖如圖5所示。

由圖5所示的灰度直方圖可以觀察得到：該灰度圖像的直方圖有兩個明顯的山峰，左邊的山峰較高而且山比較寬，右邊的山峰較低而且山較窄。結合灰度圖像可以發現，圖像中有兩個明顯的特征，分別為道路和道路兩旁的植被，而且道路兩旁的植被在該灰度圖像中所占的區域比較大。由此可以推斷出，左邊的山峰對應圖像中的植被，右邊的山峰對應圖像的道路。

閾值應當選取為兩個山峰之間的山谷，經過觀察圖5選擇閾值為T=110，提取灰度值大于110的圖像即為初始道路圖像。初始道路圖像如圖6所示。

4.2.3 中值濾波去除顆粒噪聲

由于受光照、反射、陰影等客觀因素的影響，衛星在獲取圖像的過程中會出現同一類別的區域（如植被）的顏色不統一，這就會造成灰度值的差別。所以在閾值分割后的初始道路中會出現大量的顆粒噪聲。

中值濾波器能夠很好地消除圖像中的顆粒噪聲。選擇合適的中值濾波模板是進行中值濾波的關鍵，中值濾波器的模板選擇通常借助實驗的方法。

4.2.4 數學形態學處理

雖然初始道路圖像經過中值濾波消除了大量的顆粒噪聲，但是在圖像中依然存在大量的道路間斷和道路孔洞。由數學形態學理論分析可知，開運算能夠消除兩個物體間的粘連，閉運算能夠連接間斷道路和填充孔洞。

4.2.5 最終道路提取效果

將開運算處理后的結果與原始彩色遙感圖像疊加，道路的顏色為白色，得到最終的遙感圖像道路信息提取效果如圖7所示。

5 結論

本文所研究的遙感道路提取方法的主要步驟有：圖像灰度化、閾值分割、中值濾波、開運算和閉運算和圖像疊加。將彩色圖像灰度化能夠加快圖像處理的速度，而且中值濾波器能夠濾除圖像中大量的顆粒噪聲。在研究中，閾值的選取是非常重要的，選擇爭取的閾值是后續研究的基礎。在后續的中值濾波器的設計中，設置合適的中值濾波模板也同樣重要。本文通過MATLAB圖像處理平臺，驗證了本文研究的道路提取方法的有效性。

參考文獻：

[1]李華勝.基于遙感影像的道路信息提取方法研究[D].內蒙古工業大學，2015.

[2]李潤生，曹聞.基于視覺顯著性特征的遙感影像道路網提取方法[J].計算機系統應用，2014，23（8）：114-118.

[3]闕昊懿.基于遙感圖像的城市道路及車輛信息采集研究[D].湘潭大學，2014.

[4]李俊山，等.數字圖像處理（第三版）[M].北京：清華大學出版社，2017：1-278.

[5]Wijesingha J.S.J.， Kumara R.W.D.M.， Kajanthan P.， et al. Automatic Road Feature Extraction from High Resolution Satellite? Images Using LVQ Neural Networks[J]. Asian Journal of Geoinformation， 2013，13（1）：30-36.

[6]Alex D.M.. Robust and Efficient Method to Extract Roads from Satellite Images[J]. International Journal of Latest Trends in Engineering and Technology（IJTET）， 2013，2（3）：26-29.

[7]Miao Zelang， Wang Bin， Shi Wenzhong， et al. A Method for Accurate Road Centerline Extraction From a Classified Image[J]. IEEE JOURNAL OF SELECTED TOPICS IN APPLIED EARTH OBSERVATIONS AND REMOTE SENSING，2014，7（12）： 4762-4771.