分析圓柱油罐目標合成孔徑雷達圖像的新方法

張月婷*①②④ 丁赤飚①③ 王宏琦①② 胡東輝①②

?

分析圓柱油罐目標合成孔徑雷達圖像的新方法

張月婷丁赤飚王宏琦胡東輝

(中國科學院電子學研究所 北京 100190)(中國科學院空間信息處理與應用系統技術重點實驗室 北京 100190)(微波成像技術國家重點實驗室 北京 100190)(中國科學院研究生院 北京 100190)

該文提出了一種針對油罐目標的合成孔徑雷達(SAR)圖像的分析方法。該方法根據圓柱油罐具有圓型邊緣和圓柱型外形的幾何特點，利用物理光學法(Physical Optics, PO)和增量繞射長度(Incremental Length Diffraction Coefficients, ILDC)理論建立了油罐目標的散射模型，并基于該模型推導了圓柱油罐的散射中心分布特征，同時引入投影映射算法(Mapping and Projection Algorithm, MPA)建立了成像模型，有效利用了單一視角的散射數據進行SAR圖像模擬，仿真結果與實際SAR圖像散射特征吻合，驗證了方法的有效性。同時，利用45組實測數據得到了SAR圖像散射中心的特征，并根據建立的散射模型和成像模型，分析和總結了油罐目標SAR圖像的主要散射特征，為SAR圖像中基于油罐目標的識別領域提供了理論依據。

油罐；SAR圖像；SAR圖像仿真；散射中心；散射特征

1 引言

油罐目標的SAR圖像和它的光學圖像有較大的區別，光學圖像中油罐目標有規則的圓周外形，而在雷達圖像中油罐的圓周輪廓斷斷續續，并且，油罐目標的散射中心的分布與油罐的幾何形狀存在著較大的差別。

然而，目前雷達圖像油罐目標的自動識別多數是采用同光學圖像相同的識別方法，利用油罐目標的橢圓幾何特征進行識別，該方法虛警率較高。油罐目標的有效識別的基礎是針對油罐目標SAR圖像中主要特征的散射機理的有效分析，而SAR圖像仿真是對圖像理解的有效途徑。目前，盡管SAR圖像仿真方法眾多，但針對仿真圖像開展目標特征分析的相關文獻很少。針對油罐目標而言，需要一種根據油罐目標的結構特點建立的SAR圖像特征分析方法以理解油罐目標的SAR圖像特征，為后續的油罐目標識別和信息提取提供依據。本文針對這個問題提出了一種分析油罐目標SAR圖像的方法，通過該方法對油罐目標的散射特征以及罐頂圓周的散射中心分布進行了研究。

本文根據圓柱油罐的圓周邊緣和圓柱外形的幾何特點，利用物理光學法(Physical Optics, PO)和增量繞射長度理論(Incremental Length Diffraction Coefficients, ILDC)建立了油罐目標的近似電磁散射模型，在該散射模型建立的基礎上，引入了映射投影(Mapping and Projection Algorithm, MPA)的SAR成像仿真方法完成了針對油罐目標的SAR圖像仿真，同時，利用得到的仿真SAR圖像與實際數據比對，結合散射模型的構建過程分析和歸納了油罐目標SAR圖像的主要特征，解釋了油罐目標SAR圖像中的顯著細節的產生過程，指明了主要特征對應的散射機理，為油罐目標的識別和信息提取提供了理論的依據。

本文第2節給出了針對油罐目標的SAR圖像仿真方法，包括成像模型、散射模型的分析及具體實現的主要步驟；在此基礎上，第3節計算和仿真了圓柱型油罐的SAR圖像，結合實際數據驗證了仿真圖像的有效性，同時針對散射模型分析和解釋了油罐目標的主要特征；最后，第4節給出了本文研究的結論。

2 圓柱油罐目標SAR圖像仿真方法

本節依次給出了油罐目標的成像模型，散射模型和SAR圖像仿真的主要步驟。

2.1 基于MPA的SAR成像模型

2006年，Feng Xu等在文獻[3]中提出了一種基于映射投影(MPA)的復雜地面場景極化SAR成像模擬算法，該算法采用了根據分辨率劃分網格的思路將目標分解為若干散射單元，有效實現了植被場景和建筑物場景的SAR圖像模擬，該算法利用單角度的散射數據模擬SAR圖像，避免了傳統方法的2維散射數據的計算，節省了計算量。

MPA算法的計算過程如下：在成像模擬中把整個場景沿方位向的每個像素間隔劃分為行，對每一行獨立計算散射系數，最后獲得整個場景的散射系數圖，然后與雷達系統沖激響應2維卷積獲得SAR圖像。在具體的實現中，該算法將劃分成網格的各個區域內的散射場強投影到斜距平面，再經過累加得到每個像素內的散射強度值。文獻[3]中對自然目標和規則建筑物目標進行仿真，取得了很好的效果。

根據MPA算法，散射場的計算由下式得到

油罐目標在雷達圖像中多為規則排列的油罐組，對于一個油庫基地而言，在雷達圖像中，它不同于大場景的自然目標，在油罐目標的成像區域，雷達聚焦深度變化微小，本文針對這個特點，將相同距離向弧線的疊加近似為直線疊加，并通過建立兩個坐標系實現散射系數的斜距累加，簡化了成像模型。

考察場景目標尺寸遠遠小于景中心斜距的條件下，雷達波束到達目標不同處的局部入射角近似相同，在這種近似條件下，相應的散射場系數的斜距分布函數滿足下面的等式：

圖1 場景坐標系和成像坐標系

(3)

得到散射矩陣之后，SAR圖像可以通過SAR景中心點目標響應和散射矩陣的2維卷積獲得：

景中心點目標響應的計算方法詳見文獻[5]。

2.2 圓柱型油罐散射模型

油罐按建造材料分為金屬油罐和非金屬油罐，按照形狀分為立式油罐和臥式油罐，本文主要研究臥式油罐，臥式油罐是由圓筒型的油罐壁和罐頂組成的油罐，如圖2所示。

為了建立油罐目標的散射模型，本文針對主要散射結構進行了分析。如圖2所示，在油罐目標與電磁波相互作用時，對散射場起主要貢獻的結構為：油罐側壁以及罐頂部圓周外緣，圓周內緣，罐底圓周。

首先分析油罐壁的散射，通常，罐頂內圓周結構低于外圓周，即：低于，因此油罐外壁和部分油罐內壁產生散射。被照射的油罐內壁貢獻較小，近似用拋物面形(平面與柱面的交線)內壁的散射代替油罐內壁的散射。

圖2 幾何關系圖

(6)

(8)

將式(8)代入式(5)，有：

其中，

(11)

根據駐定相位原理(Principle of Stationary Phase, POSP)，求得：駐相點。

求得散射場為：

(13)

觀察式(13)和式(14)，可以得到以下的散射模型：

物理光學法近似下非金屬油罐和金屬油罐的散射特點是相似的，油罐外壁會產生兩個位于圖2中點和點的散射中心，坐標分別為：，內壁局部的散射會產生一個位于的散射中心，且散射強度小。

以上的計算過程直接建立了油罐壁的散射模型，接下來考慮油罐主要圓周結構的散射模型。

油罐目標主要包含3段主要的圓周邊緣，如前文所述。本文對這些圓周邊緣進行線元建模，應用Incremental Length Diffraction Coefficients，即：ILDC理論，計算它們的散射。在后向散射的情況下，應用ILDC理論推導出每個線元的散射場為：

2.3 SAR圖像仿真的流程

根據前兩節的分析，在SAR圖像模擬前首先設置雷達工作參數和載入油罐模型；主要的參數包括：雷達工作頻率，方位角，入射角，成像分辨率以及油罐的幾何參數和介電常數。油罐模型可以用這樣的方法建立：油罐頂部內外圓周棱邊和底部被照射棱邊分成個線元，根據彈射線理論射線密度大小，線元長度為數量級，每個線元對應一組向量，該向量記錄線元起點終點坐標，線元所屬平面法向量等信息。首先根據雷達參數計算景中心點目標響應，接著計算場景坐標系到成像坐標系的坐標變換矩陣。之后，建立散射矩陣，根據前文建立的散射模型和散射計算方法計算線元的散射和罐壁的散射，根據計算得到的散射場矢量更新散射矩陣，直到所有線元計算完成。在獲得散射矩陣之后，根據地距斜距關系進行幾何校正，最后輸出油罐目標的仿真SAR圖像。

以上建立的仿真油罐目標SAR圖像的算法的具體步驟如下：

(1) 設置雷達參數，建立油罐模型；

(2) 建立場景坐標系和成像坐標系，計算景中心點目標響應和坐標變換矩陣；

(3) 讀取下一個圓周緣結構線元參數；

(4) 重復(3)直到所有線元計算完畢；

(7) 幾何校正，獲取模擬油罐圖像。

3 數值計算結果和散射特征分析

本節給出了混凝土油罐實驗仿真結果，同時與相應參數的機載雷達圖像進行了對比，并針對45組實際SAR圖像進行了散射特征的分析。

實驗選擇混凝土油罐目標的幾何尺寸為：高17 m，半徑：40.5 m；如圖2所示。雷達工作頻率L波段，電場為水平極化。主要參數見表1。

表1 雷達參數

參數數值 飛行高度5000 m 分辨率5 m×5 m 工作頻率1.2 GHz 飛行速度150 m 過采樣率1.2 波束入射角36°

下面分析和討論主要的數值計算結果。

(1) 景中心點目標響應。

表2 點目標指標

指標距離向方位向 分辨率(m)55 擴展比0.906 0.907 峰值旁瓣比(dB)-13.496-13.375 積分旁瓣比(dB)-10.821 -9.329

(3) 圓周緣結構的散射。

圖4為罐頂面圓周緣散射強度幅值的變化曲線，虛線表示外圓周散射強度隨位置的變化，實線表示內外圓周散射疊加后的散射強度。從圖4中可以看到以下幾個特征：

(a) 外圓周緣的散射主要出現了5個峰值。

(b) 圓周緣散射場出現了閃爍現象。散射強度的波瓣寬度小于0.1度，與實際SAR圖像中邊緣斷斷續續的特征吻合。

(c) 距離電磁波照射的最近點和最遠點的周圍散射強度較大，兩個點位置對應圖2中,兩點。

(4) 仿真SAR圖像。

圖5為仿真得到的油罐目標的SAR圖像，圖6是基于仿真圖像散射中心的分析示意圖。圖6示意了仿真圖像中的散射中心位置。觀察圖5和圖6，油罐目標的散射特征主要有：

(a) 在罐頂圓周圓心位置所處的方位向，3個“月牙”型散射中心沿距離向一字分布，結合數值計算結果，“月牙”中心分別對應圖2中,,3點。3個強散射中心沿距離向排列，由近距到遠距分別對應圖2的,,3點，而圖2中對應的,,,點處并沒有形成較強的散射中心。

(b) 罐頂圓周不連續。

實驗將3個散射中心所處的方位向剖面進行8倍sinc插值，并與實際信號的平均值(實驗選取了45組油罐SAR圖像)進行了對比，如圖7所示。圖7中兩條曲線的主要峰值基本吻合，驗證了本文方法的有效性。

(5) 實際SAR圖像的散射中心的提取。

實驗選取了45組實測數據得到的SAR圖像，利用文獻[9]的方法提取了罐頂圓周緣的散射中心，并計算了每個散射中心與圓心形成向量的幅角，利用輻角信息進行了分析，以10°為一個區間統計了各區間內的散射中心的位置分布信息，表3給出了包含散射中心數目最多的5個角度區間的計算結果。實驗發現，散射中心集中輻角為40°～50°，-50°～-40°，100°～110°，-110°～-100° 4個角度區間內，相應區間內計算得到的4個均值量與圖4的數值計算顯示的峰值位置基本吻合。

圖3 散射強度隨角度變化曲線

圖4 罐頂散射場變化曲線

圖5 仿真油罐SAR圖像

圖6 散射中心分析示意圖

圖7 罐頂中心所處方位向剖面的像素對比

表3實際SAR圖像罐頂散射中心分布

角度區間(°)存在散射中心的油罐目標數散射中心角度均值(°)散射中心角度方差 [-40, -50)23-46.922.45 [100, 110)20106.3510.66 [40, 50)1946.6314.91 [-110, -100)18-105.1718.38 [70, 80)971.114.86 [-120, -110)6-116.525.1

4 結束語

本文提出了一種分析圓柱油罐目標合成孔徑雷達圖像的新方法，該方法根據油罐目標的結構特點建立了散射模型和成像模型，并據此給出了圓柱型油罐目標的SAR圖像的仿真方法。實驗表明，仿真SAR圖像與實際雷達圖像的特性吻合，驗證了模型的有效性。同時，根據建立的散射模型和成像模型，分析和總結了油罐目標SAR圖像的主要散射特征，指明了各個主要特征對應的電磁散射機制。

本文方法為SAR圖像中油罐目標特性的提取工作提供了一種有效的研究途徑，同時，文中得到的油罐目標的主要散射特征為油罐目標的信息提取提供了理論依據。此外，該方法的研究思路：根據目標的幾何特點建立散射模型以分析其SAR圖像主要特征的散射機理，也可以用于其它人工目標的SAR圖像理解工作中。

[1] 夏昕. SAR圖像典型目標特征提取與識別方法研究[D]. [碩士論文], 四川大學, 2006, 第1章.

Xia Xin. Study of feature extraction and recognition of typical targets in SAR images[D]. [Master dissertation], Sichuan University, 2006, Ch. 1.

[2] Wang Yuanxun and Ling Hao. Radar signature prediction using moment method codes via a frequency extrapolation technique[J]., 1999, 47(6): 1008-1014.

[3] Xu Feng and Jin Ya-qiu. Imaging simulation of polarimetric SAR for a comprehensive terrain scene using the mapping and projection algorithm[J]., 2006, 44(11): 3219-3234.

[4] Franceschetti Giorgio and Schirinzi Gilda. A SAR processor based on two-dimensional FFT codes[J]., 1990, 26(2): 356-366.

[5] Cumming Ian G and Wong Frank H. Digital Processing of Synthetic Aperture Radar Data: Algorithms and Implement- tation[M]. London: Artech House, 2005, Ch. 5.

[6] Ling H, Chou R, and Lee S W. Shooting and bouncing rays: calculating RCS of an arbitrary cavity[J]., 1989, 37(2): 194-205.

[7] Shore Roben A and Yaghjian Anhur D. Incremental diffraction coefficients for planar aurfaces[J].1988, 36(1): 55-70.

[8] 王新賽, 張天序, 桑農. 實孔徑雷達景象模擬研究[J]. 紅外與毫米波學報, 2002, 21(1): 54-58.

Wang Xin-sai, Zhang Tian-xu, and Sang Nong. Simulation of real aperture radar scene[J]., 2002, 21(1): 54-58.

[9] Rihaczek August W and Hershkowitz Stephen J. Theory and Practice of Radar Target Identification[M]. London: Artech House, 2000, Ch. 2.

A New Algorithm for Analyzing the SAR Image of the Cylindrical Oil Tank

Zhang Yue-tingDing Chi-biaoWang Hong-qiHu Dong-hui

(Institute of Electronics Chinese Academy of Sciences, Beijing 100190, China)(Key Laboratory of Spatial Information Processing and Application System Technology, Chinese Academy of Sciences, Beijing 100190, China)(The National Key Laboratory of Microwave Imaging Technology, Beijing 100190, China)(Graduate University of Chinese Academy of Sciences, Beijing 100190, China)

In this paper, a new approach for the understanding of the SAR image of the cylinder oil tank is presented. In the approach, the Physical Optics (PO) and the theory of Incremental Length Diffraction Coefficients (ILDC) are employed to set up the scattering model of the oil tank according to the circle edges and cylinder geometry futures, and the position of the main scattering centers are deduced. The SAR image of the tank is simulated through the scattering model and the imaging model based on Mapping and Projection Algorithm (MPA). The real SAR data and the simulation test the validity of the approach. Furthermore, the scattering centers are detected with 45 groups of the measured data. According the scattering model and the imaging model in the approach proposed, the main scattering characteristics of the cylindrical oil tank are interpreted and summarized. It is hoped that the work in this paper could provide some useful guidelines for the identification of the oil tank in the SAR images.

Oil tank; SAR image; SAR image simulation; Scattering center; Scattering characteristics

TN959

A

2095-283X(2012)02-0190-06

10.3724/SP.J.1300.2012.20020

2012-03-29收到，2012-06-07改回；2012-06-13網絡優先出版

國家自然科學基金(40871209)資助課題

張月婷 zhangyueting06@mails.gucas.ac.cn

張月婷(1983-)，女，中國科學院電子學研究所在讀博士生，主要從事SAR圖像仿真和SAR圖像理解領域的研究。

丁赤飚(1969-)，男，研究員，博士生導師，現任中國科學院電子學研究所副所長，主要從事合成孔徑雷達、遙感信息處理和應用系統等領域的研究工作，先后主持多項國家863重點項目和國家級遙感衛星地面系統工程建設項目，曾獲國家科技進步一等獎、二等獎各一項。

王宏琦(1964-)，男，研究員，博士生導師，長期從事信息處理及相關技術的科研工作，曾經參與的研究工作主要涉及地下探測雷達、遙感圖像典型目標探測技術和遙感數據地面處理系統等方面。目前研究領域主要為遙感圖像處理及目標識別。

胡東輝(1970-)，男，副研究員，主要從事合成孔徑雷達成像處理領域的研究。