新聞中心

EEPW首頁 > 模擬技術(shù) > 設計應用 > 旋轉(zhuǎn)目標干涉逆合成孔徑三維成像技術(shù)

旋轉(zhuǎn)目標干涉逆合成孔徑三維成像技術(shù)

作者: 時間:2011-11-24 來源:網(wǎng)絡 收藏

孔徑技術(shù)

本文根據(jù)微波成像原理,導出了采用孔徑(INISAR)處理獲取第三維高程信息的INISAR技術(shù).這種技術(shù)通過鑒別由俯仰角微小差異的兩副天線分別獲取的兩幅相位相干的二維圖像對應散射中心的相位變化來確定散射中心在豎直維的位置;散射中心的橫坐標、縱坐標及幅度由一幅二維圖像給出.模擬計算與飛機模型實驗測量均得到滿意結(jié)果.
  關(guān)鍵詞:雷達目標;微波成像;測量

INISAR 3-D Imaging Technique for Rotating Targets

XIAO Zhi-he,DAI Chao-ming,CHAO Zeng-ming,XU Xiao-jian
(Beijing Institute of Environmental Features,Beijing 100854,China)

  Abstract:Based on the principle of microwave imaging of the rotating targets,an INISAR (ISAR Interferometry) 3-D imaging technique of rotating targets is described in this paper.The vertical position of scattering centers are given by comparing their phase difference of two phase coherent 2-D images which are obtained by using two antennas with a little change in pitch angle.The down-range and cross-range position and amplitude of scattering centers of the 3-D image are given by one 2-D image.The simulated and measured results are all satisfactory.
  Key words:radar targets;microwave imaging;measurement

一、引  言
  雷達目標二維成像的理論和測量技術(shù)已經(jīng)取得了巨大的發(fā)展和廣泛應用,并成為電磁散射實驗室、目標特性及其控制研究、微波遙感等領(lǐng)域不可缺少的標志性技術(shù).作為二維成像的推廣,能獲取復雜目標散射中心三維分布的技術(shù),由于要求獲取巨大數(shù)量的測量數(shù)據(jù),從目前的試驗手段而言可能需要幾十個小時,實際難以實現(xiàn).因此,完全意義上的三維成像一直未能在實驗室中實現(xiàn).如何采用新的可行的技術(shù)途徑獲得散射中心高程坐標是人們一直試圖解決的問題.
  圖樣可以用來捕捉目標表面相對參考平面的位移和變形,這在光學領(lǐng)域已有悠久歷史.隨著合成孔徑雷達(SAR)在微波遙感領(lǐng)域的發(fā)展和應用,到九十年代初,人們提出并研究了干涉合成孔徑雷達(INSAR)[1~4].INSAR法即通過測量在高度方向上位置有微小差別的兩副天線的相位差或是多軌道(可以非同時)情況下測得的相位差,提取有關(guān)地形的第三維高程信息.
  旋轉(zhuǎn)目標二維成像基于圓周孔徑的ISAR成像,通過寬帶脈沖壓縮得到徑向高分辨;通過圓周合成孔徑得到橫向高分辨.本文根據(jù)旋轉(zhuǎn)目標微波成像原理,導出了采用干涉孔徑(INISAR)處理獲取第三維高程信息的INISAR三維成像技術(shù).這種技術(shù)通過鑒別由俯仰角微小差異的兩副天線所獲取的相位相干的兩幅二維圖像對應散射中心的相位變化來確定散射中心在豎直維的位置,從而得到目標散射中心三維坐標.這種技術(shù)適用于在同一徑向、橫向分辨單元內(nèi)只有一個強散射源的情況.

二、旋轉(zhuǎn)目標二維成像原理
  圖1為轉(zhuǎn)臺逆合成孔徑成像的幾何關(guān)系,其中X—軸為橫向,Y—軸為徑向.設目標散射中心空域分布為g(x,y),則相參雷達的接收數(shù)據(jù)(譜域)為:

旋轉(zhuǎn)目標干涉逆合成孔徑三維成像技術(shù)

圖1 轉(zhuǎn)臺成像幾何圖

G(K,ψ)=∫∫g(x,y)ej2πK(ycosψ+xsinψ)dxdy (1)

式中:K=2/λ,為空間頻率,ψ為雷達視線與Y軸(徑向坐標)夾角.
  (ycosψ+xsinψ)為g(x,y)在雷達視線上的投影;G(k,ψ)和g(x,y)互為傅立葉變換關(guān)系.
  在有限的頻帶和有限的觀察角范圍內(nèi),G(K,ψ)的自變量取值形式如圖2所示,為極坐標格柵數(shù)據(jù).

旋轉(zhuǎn)目標干涉逆合成孔徑三維成像技術(shù)

圖2 極坐標格柵數(shù)據(jù)

  根據(jù)式(1)得到目標二維散射中心分布的估計值旋轉(zhuǎn)目標干涉逆合成孔徑三維成像技術(shù)(x,y)的濾波——逆投影算法:

旋轉(zhuǎn)目標干涉逆合成孔徑三維成像技術(shù) (2)

三、INISAR三維成像原理
  INSAR三維成像的幾何關(guān)系見圖3所示,其中,X—為橫向;Y—為徑向;Z—為豎向.設體目標散射中心的直角坐標為g(x,y,z),對應的球坐標為g′(ρ,φ,θ),坐標之間的關(guān)系為:

旋轉(zhuǎn)目標干涉逆合成孔徑三維成像技術(shù)

圖3 三維成像的幾何關(guān)系圖

旋轉(zhuǎn)目標干涉逆合成孔徑三維成像技術(shù) (3)

并記g的方向矢量為:

g=(sinθcosφ,sinθsinφ,cosθ).

  設INISAR的兩副接收天線分別為M1和M2,均指向目標坐標原點;其入射線均位于包含Z軸在內(nèi)的入射內(nèi)平面P內(nèi),P與XOY平面交線的方向角為φ1.設M1的方向矢量為

M1=(sinθ1cosφ1,sinθ1sinφ1,cosθ1) (4)

M2的方向矢量為:

M2=(sinθ2cosφ1,sinθ2sinφ1,cosθ2) (5)

則目標散射中心g在M1和M2上的投影分別為:

R1=ρg.M1=xsinθ1cosφ1+ysinθ1sinφ1+zcosθ1
R2=ρg.M2=xsinθ2cosφ1+ysinθ2sinφ1+zcosθ2

  設目標ISAR成像的旋轉(zhuǎn)面為XOY平面,且M1在XOY平面上,M2與XOY平面有夾角Δθ,則得:θ1=90°,θ2=90°-Δθ.記ψ為M1(M2相同)與Y軸的夾角,則ψ=90°-φ1,ψ∈(ψ1,ψ2)為二維成像合成孔徑角范圍.則得:

R1=ycosψ+xsinψ (6)
R2=(ycosψ+xsinψ)cosΔθ+zsinΔθ≈ycosψ+xsinψ+zsinΔθ (7)

  則散射中心g在M1所得圖像I1和M2所得圖像I2上的程差為:

ΔR=R2-R1≈zsinΔθ (8)

  對應的傳輸相位為:

旋轉(zhuǎn)目標干涉逆合成孔徑三維成像技術(shù) (9)

  式中所示的傳輸相位將反映在圖像I1和圖像I2對應像素點相應散射中心的相位差上.因此通過鑒別圖像I1和圖像I2對應像素點相應散射中心的相位差,則由式(9)可以得到其豎向坐標為:

z=(λ.Θ)/(4π.sinΔθ) (10)

  為避免相位模糊,Θ的取值范圍應為-π<Θ<π,設|z|<z0,則有

Δθ(rad)<λ/(4z0) 或
Δθ(deg)<14λ/z0=4.2/(f0z0) (11)

  式中f0單位:GHz;z0單位:m.
  在INISAR三維成像中,目標散射中心分布的橫向(x)坐標、徑向(y)坐標及散射強度由天線M1所得二維成像結(jié)果給出,散射中心的豎向(z)坐標由式(10)給出.
  由式(10)可知,散射中心豎向坐標的誤差主要由角度值Δθ和相位差值Θ決定.其中Δθ是一固定值,可以通過一定的技術(shù)手段得到較準確的計量.由Δθ引起的誤差為:

旋轉(zhuǎn)目標干涉逆合成孔徑三維成像技術(shù) (12)

  設f0=10GHz,z0=1m,則由式(12)得,Δθ<0.42°.取Δθ=0.4°.設Δθ的誤差為0.01°,則Δzθ=-2.5%,z的絕對誤差小于2.5cm.可見由Δθ引起的誤差為恒定相對誤差,隨z坐標的增大,其絕對誤差將增大.因此,在滿足式(11)的條件下Δθ的值應盡量取大些.當頻率高時,由于Δθ的值取得小,Δzθ會增大.
  散射中心豎向坐標的另一主要誤差來源是相位量Θ的測量誤差,其相對誤差為:

旋轉(zhuǎn)目標干涉逆合成孔徑三維成像技術(shù) (13)

可見,由相位量Θ的測量誤差引起的散射中心豎向坐標的相對誤差與其豎向坐標的量值成反比,但有固定的最大絕對誤差.高精度相干測量雷達的相位測量誤差小于2°,即ΔΘ=2°,設f0=10GHz,Δθ=0.4°,z0=1.0m,則Δzθ=1.2%;絕對誤差小于1.2cm.

四、模擬計算
  為便于計算,以金屬球組合目標為例進行模擬計算.金屬球后向散射用級數(shù)解求得.球心不在坐標中心時,按徑向距離加入傳輸相位.
  假設5個球的組合目標,球的半徑及其分布位置設定見下表.先分別計算每一金屬球的級數(shù)解,并根據(jù)其球心位置調(diào)整其相位;將5個球的散射貢獻矢量相加,逐一生成極坐標格柵的雷達模擬數(shù)據(jù).設成像頻帶范圍8~12GHz;合成孔徑角范圍24°;兩天線俯仰角差Δθ=0.4°.分別計算θ=90°和θ=90°-0.4°時的兩幅二維成像.根據(jù)式(10)計算得到每一像素的豎向坐標,畫成立體圖.如圖4所示.強散射源的三維位置(鏡面)及幅度見表1.從表中可見:橫向、徑向位置誤差小于1cm,豎向位置誤差小于2cm,均在一個分辨力范圍內(nèi).

旋轉(zhuǎn)目標干涉逆合成孔徑三維成像技術(shù)
旋轉(zhuǎn)目標干涉逆合成孔徑三維成像技術(shù)

圖4 五個球的組合目標三維成像模擬計算結(jié)果

表1

球半徑(m)及RCS(dBm2)球心位置預設(m)球鏡面位置預設(m)球鏡面位置測量值(m)
序號半徑預估RCS測量RCS橫向徑向豎向橫向徑向豎向橫向徑向豎向
10.05-21.05-21.30.800.800.50.800.750.50.800.760.51
20.04-22.22-23.00.400.400.00.400.360.00.400.360.00
30.03-25.41-26.00.000.00-0.20.00-0.03-0.20.00-0.02-0.20
40.02-30.93-29.3-0.40-0.40-0.4-0.40-0.42-0.4-0.40-0.42-0.40
50.01-33.43-34.4-0.80-0.80-0.6-0.80-0.81-0.6-0.80-0.80-0.62
五、微波暗室實驗測試
  根據(jù)前述三維成像測量技術(shù),用飛機模型進行了實驗測量.測試中心頻率9.25GHz,帶寬1.8GHz,水平極化,俯仰角差Δθ=-0.4°.測試結(jié)果見圖5所示.如圖中所示,P2為飛機的面天線;P1為壓氣機,距天線1.09m;P3、P5為進氣口,相距0.15m,距天線0.435m;P4為座椅,距天線0.41m,比飛機軸心高約0.1m;P6為座艙前沿,距天線約0.29m;P7、P8為導彈架,相距0.61m,距天線0.86m,比飛機軸心低0.17m.

旋轉(zhuǎn)目標干涉逆合成孔徑三維成像技術(shù)
旋轉(zhuǎn)目標干涉逆合成孔徑三維成像技術(shù)

圖5 某飛機縮比模型(頭部方位)三維成像測量結(jié)果

  在測量中,飛機軸心在豎向軸的位置約為-4cm.分析圖5所示結(jié)果,飛機模型散射中心的相對位置關(guān)系吻合良好,誤差在一個分辨單元之內(nèi).表明了旋轉(zhuǎn)目標干涉逆合成孔徑三維成像技術(shù)的可行性.

六、結(jié)  論
  本文提出了一種簡單易行的三維成像測量技術(shù),該技術(shù)能得到目標強散射源散射中心的三維坐標.模擬計算與飛機縮比模型暗室實驗測量表明該技術(shù)切實可行.三維成像測試技術(shù)的建立為目標散射中心診斷、多散射源建模與應用等提供了新的技術(shù)手段.

合成孔徑雷達相關(guān)文章:合成孔徑雷達原理


評論


相關(guān)推薦

技術(shù)專區(qū)

關(guān)閉