亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        視頻合成孔徑雷達技術發(fā)展現(xiàn)狀綜述

        2024-11-26 00:00:00顏上取付耀張文鵬楊威余若峰張法桐
        系統(tǒng)工程與電子技術 2024年8期

        摘 要: 視頻合成孔徑雷達 (video synthetic aperture radar, ViSAR) 由于能夠連續(xù)觀測目標區(qū)域的變化,在許多領域得到了廣泛的應用。首先梳理了國內(nèi)外ViSAR相關的研究動態(tài),分析了ViSAR高幀率成像以及動目標陰影形成的原理;然后系統(tǒng)地闡述了國內(nèi)外ViSAR系統(tǒng)及其處理技術的研究進展,對成像算法、運動補償算法以及動目標檢測與跟蹤技術等方面的研究進展進行了梳理和總結;最后總結了目前ViSAR相關的發(fā)展,并對未來ViSAR技術的潛在研究方向進行了展望。

        關鍵詞: 合成孔徑雷達; 高幀率成像; 視頻合成孔徑雷達; 動目標檢測

        中圖分類號: TN 95

        文獻標志碼: A

        DOI:10.12305/j.issn.1001-506X.2024.08.13

        Review of the development status for ViSAR techniques

        YAN Shangqu, FU Yaowen*, ZHANG Wenpeng, YANG Wei, YU Ruofeng, ZHANG Fatong

        (College of Electronic Science and Technology, National University of Defense Technology, Changsha 410073, China)

        Abstract: Video synthetic aperture radar (ViSAR) has been widely used in many fields due to its ability to continuously observe changes in the target area. Firstly, the research trends related to ViSAR at home and abroad are reviewed, and the principles of high frame rate imaging and shadow formation of moving targets in ViSAR are analyzed. Secondly, the research progress of ViSAR systems and their processing technologies of both domestical and international is systematically elaborated, and the research progress in imaging algorithms, motion compensation algorithms, and moving target detection and tracking technologies is summarized. Finally, the current developments related to ViSAR are summarized, and potential research directions for future ViSAR technology are discussed.

        Keywords: synthetic aperture radar (SAR); high frame rate imaging; video synthetic aperture radar (ViSAR); moving target detection

        0 引 言

        20世紀50年代,美國數(shù)學家Carl Wiley研究阿特拉斯液體洲際彈道導彈時發(fā)現(xiàn)如果對多普勒頻移進行處理,能改善雷達方位向的分辨率。根據(jù)這一原理,可利用雷達的運動來合成等效長的方位向長孔徑,并獲得目標區(qū)域的二維地表圖像。這一思想后來被發(fā)展成合成孔徑雷達(synthetic aperture radar, SAR)技術[1。SAR技術具備全天候、全天時、高分辨率以及可穿透云霧等優(yōu)點[2-3,經(jīng)過多年的發(fā)展,現(xiàn)已作為一種主動的航天、航空遙感手段,被廣泛地應用于地形測繪[4、資源普查5、目標偵察6以及戰(zhàn)場態(tài)勢感知等民用或軍用領域7。但傳統(tǒng)SAR僅能獲取目標成像區(qū)域的靜態(tài)二維(距離和方位)平面信息,無法得到更豐富的場景信息。因此,研究人員通過對SAR的深入研究,提出不同的方案,來獲取更多的場景信息。其中,干涉SAR(interferometric SAR, InSAR)可獲取目標成像區(qū)域的高度信息,InSAR通過SAR復圖干涉得到干涉相位,來獲取目標的高度信息,從而實現(xiàn)對成像區(qū)域內(nèi)高程信息的檢測[8。極化SAR(polarimetric SAR, PolSAR)利用多種極化方式獲取成像區(qū)域內(nèi)多種散射信息,使SAR圖像由近似“黑白”的圖像向 “彩色”圖像發(fā)展[9。視頻SAR(video SAR, ViSAR)將視頻顯示技術與SAR成像技術相結合,通過對目標成像區(qū)域的凝視獲取該區(qū)域的連續(xù)幀圖像,并且最終以“視頻”的形式展示,使SAR圖像由“靜態(tài)”向“動態(tài)”發(fā)展,是目前SAR技術中最為前沿的課題之一[10-15。

        美國桑迪亞國家實驗室(Sandia National Laboratory, SNL)于2003年首次提出了ViSAR的概念。ViSAR可通過實時調(diào)整雷達波束來對地面目標區(qū)域進行長時間凝視,同時對接收的回波信號進行持續(xù)成像,來形成連續(xù)的高幀率圖像,并且將SAR技術與視頻顯示技術進行結合,高效地展示出地面的變化情況,同時利用觀察者對運動目標的敏感性提升對場景動態(tài)信息的分析能力。當ViSAR的成像分辨率與幀率分別達到0.2 m和5 Hz時,觀察者能更為準確地對地面動目標進行實時檢測與跟蹤[16-17。作為一種新型的對地探測技術,ViSAR大致可分為兩類,一類是高頻ViSAR;另一類是非高頻ViSAR。高頻ViSAR技術通過提高工作頻率來提高成像幀率,其工作頻段位于太赫茲頻段,最大的優(yōu)勢在于成像延遲和合成孔徑積累時間較短,易于實現(xiàn)對成像場景的持續(xù)實時監(jiān)測,并且高頻ViSAR的雷達天線尺寸較小,有利于搭載至無人機等小型平臺。但由于高頻ViSAR的工作頻率較高,受大氣衰減影響較大,其工作距離較短,并且其極高的實時性也對處理算法的效率提出了一定的要求。非高頻ViSAR技術采用孔徑重疊模式提高成像幀率,其最大優(yōu)勢在于不需要較高的工作頻率,可適用于目前大多數(shù)的SAR系統(tǒng)。此外,非高頻ViSAR的相關算法以及系統(tǒng)較為成熟,能提供較長的工作距離。但選擇使用孔徑重疊模式,務必需要較長的積累時間,造成一定的成像延時,實時性較差,可能無法滿足對成像場景的實時監(jiān)測[18-19。

        本文主要圍繞這兩類ViSAR的相關技術展開討論,梳理了國內(nèi)外相關的研究動態(tài),對兩類ViSAR系統(tǒng)的高幀率成像以及動目標陰影形成原理進行分析;對成像算法、運動補償算法以及動目標檢測與跟蹤技術等方面的研究進展進行了梳理和總結,并且在此基礎上對該技術未來的發(fā)展趨勢進行了展望。

        1 ViSAR系統(tǒng)原理

        圖1為經(jīng)典條帶SAR成像幾何的示意圖,載機以固定高度H、飛行速度V沿直線軌跡進行勻速飛行,Pr為雷達天線相位中心 (antenna phase center, APC),Pn為測繪帶中任意一個點目標。根據(jù)SAR成像理論可知,條帶SAR距離向分辨率20可表示為

        ρr=c2Br(1)

        式中:Br為距離向帶寬;c為光速。同時,方位向分辨率可表示為

        ρa=VBa(2)

        式中:Ba為多普勒帶寬。

        在ViSAR體制中,因需要對目標成像區(qū)域進行長時間持續(xù)觀測,并要以動態(tài)的形式反映成像區(qū)域內(nèi)的運動狀態(tài),ViSAR系統(tǒng)通常工作于聚束模式下[21。在聚束模式下,為了滿足對地面目標區(qū)域的持續(xù)觀測需求,在載機平臺飛行的過程中,需實時地調(diào)整雷達波束,使其持續(xù)指向目標成像區(qū)域,獲得持續(xù)觀測所需的回波信息以進行成像。當照射時間增加時,合成孔徑積累角也會對應地增大,這樣可形成更大的多普勒帶寬Ba,從而達到比條帶模式距離向分辨率更高的方位向分辨率。因此,在目前的研究中,為了便于對幀率進行分析和對目標成像區(qū)域進行監(jiān)測,ViSAR成像大都采用圓周軌跡聚束模式[22-25

        圖2為圓周軌跡聚束模式示意圖,平臺以固定的高度沿圓周軌跡飛行,同時不斷地調(diào)整雷達波束的指向來保證飛行過程中成像區(qū)域始終在雷達的波束范圍之內(nèi)。ViSAR的距離向分辨率[26可表示為

        ρr=c2Brcos φ(3)

        同時,在聚束模式中,ViSAR的方位向分辨率[27可表示為

        ρa=λ2Δθcos θsq(4)

        式中:Δθ為方位向積累角;λ為波長;θsq為斜視角。一般來說,為了使SAR圖像不出現(xiàn)壓縮展寬等現(xiàn)象,需要保持ρra,即保持距離向和方位向分辨率一致。因此,合成孔徑積累角也需要滿足以下條件:

        Δθ=Brλc=Brfc(5)

        式中:fc為雷達工作頻率。

        1.1 高頻ViSAR基本原理

        根據(jù)式(4)可知,當雷達的工作頻率fc增大時,ViSAR僅需要較小的Δθ就可以達到較高的地面方位向分辨率。對于速度恒定為V的載機平臺來說,對應的合成孔徑時間Ts也會相對應減小。當忽略成像所需的時間時,每幀SAR圖像的合成孔徑時間Ts

        Ts=Raλ2ρav(6)

        同時,ViSAR的獨立幀率Fs為Ts的倒數(shù),即:

        Fs=2ρavcRafc(7)

        式中:Ra為雷達距離成像中心的斜距。獨立幀率的定義為ViSAR系統(tǒng)每秒鐘可獲得的SAR幀圖像的個數(shù),即單幀SAR圖像的合成孔徑時間Ts的倒數(shù)。對一個ViSAR系統(tǒng)而言, Fs由地面方位向分辨率ρa、載機飛行速度V、雷達工作頻率fc以及Ra共同決定。其中,ρa與Ra受成像應用的限制,通常無法改變。如需提高ViSAR系統(tǒng)幀率,只能提高載機飛行速度V或雷達工作頻率fc。通常采用提高雷達工作頻率fc的方式來達到高幀率的需求。SAR成像必須經(jīng)過一個合成孔徑時間Ts才能形成所需的方位向分辨率,那么對應地,成像就會有Ts的延遲時間。根據(jù)獨立幀率Fs和合成孔徑時間Ts的關系可知,獨立幀率越高,合成孔徑時間越短,成像延遲時間越低,系統(tǒng)的實時性也越高。同時,根據(jù)美國國防部高級研究計劃局(Defense Advanced Research Projects Agency, DARPA)和SNL的研究表明,ViSAR的成像分辨率以及幀率需分別到達0.2 m和5 Hz,才能通過肉眼實現(xiàn)對地面目標的實時跟蹤與監(jiān)測[28-29。

        在此,以成像分辨率0.2 m、作用距離2 km、載機飛行速40 m/s為例,根據(jù)式(7)計算fc與Fs之間的關系,計算結果如圖3所示。從圖3的計算結果可知,ViSAR的獨立幀率要達到5 Hz,fc則需要在187.4 GHz以上,這個頻段已經(jīng)屬于太赫茲波段,相對于Ku與Ka波段的雷達來說,該波段更易受到大氣衰減的影響。

        根據(jù)圖4所示的不同頻段下的大氣衰減曲線可知,在越高的工作頻率下,大氣衰減是越大的,但在94 GHz、220 GHz、340 GHz、580 GHz等頻段內(nèi),大氣對電磁波的吸收較弱,因此這幾個頻帶被稱為“大氣窗口”[30。結合圖3的計算結果以及目前的硬件水平來看,為了滿足0.2 m的分辨率以及5 Hz的獨立幀率,高頻ViSAR系統(tǒng)的雷達工作頻率應在220 GHz以及340 GHz附近。

        1.2 非高頻ViSAR基本原理

        對于X波段、K波段以及Ka波段的非高頻雷達系統(tǒng),為了實現(xiàn)5 Hz的幀率,需要通過孔徑重疊模式來提高幀率,以達到ViSAR系統(tǒng)的要求。

        圖5為孔徑重疊模式的原理示意圖,從圖中可知,在相鄰幀圖像之間,有部分的數(shù)據(jù)是重疊的,只需更新小部分的數(shù)據(jù)就可得到下一幀的圖像,從而提高幀率。使用孔徑重疊模式得到的ViSAR系統(tǒng)的幀率稱之為交疊幀率,其表達式為

        Fl=2ρavcRa(1-γ)fc(8)

        式中:γ(0≤γ<1)表示兩幀孔徑之間的重疊率。除第一幀SAR圖像之外,形成一幀圖像所需的時間為Tl=(1-γ)Ts。

        由式(8)可知,在一個ViSAR系統(tǒng)的其他參數(shù)不變的情況下,可以通過提高孔徑重疊率γ提高系統(tǒng)成像的幀率。在理論上,雷達每接收一個脈沖均可形成一幀新的SAR圖像。因此,脈沖重復頻率(pulse repetition frequency, PRF)的數(shù)值是孔徑重疊模式下ViSAR可達到的最大交疊幀率的數(shù)值。另外,采用孔徑重疊模式成像還可以增加相鄰幀圖像之間的相干性[31,讓相鄰圖像之間有更強的連續(xù)性,有助于觀察者的判斷以及后續(xù)的處理。但是,過高的孔徑重疊率γ會導致相鄰幀圖像的變化較小,造成較高的信息冗余。因此,選擇合適的孔徑重疊率γ也是非高頻ViSAR系統(tǒng)設計中的重要一環(huán)[32

        1.3 動目標陰影形成原理

        ViSAR圖像序列的動目標陰影可用于地面動目標檢測與跟蹤,本節(jié)將詳細描述ViSAR圖像陰影成像的原理。在ViSAR體制中,雷達的工作頻率較高,目標較小的運動就會使目標成像發(fā)生較大的偏移和散焦。其次,動目標自身的遮擋會導致目標在地面真實位置的回波能量較弱,在成像結果中會出現(xiàn)真實反映動目標位置以及狀態(tài)信息的陰影。

        圖6為動目標陰影形成原理示意圖,在合成孔徑時間Ts內(nèi),車輛從A點運動到B點,A點為起始時刻車尾的位置,B點為結束時刻車頭的位置。在整個Ts內(nèi),車輛的陰影區(qū)域于A點開始,至B點終止。O1為動目標沿運動方向的陰影強度剖面(shadow intensity profile, SIP)的起點,一般在背景為均勻介質的模型中,陰影的強弱與雷達波束照射的時間成正比。rs表示O1距離SIP起點的距離。r表示SIP上任意一點P到O1的距離。動目標開始被遮擋時的時刻ts可表示為

        ts=r-rs-lvt(9)

        式中:l表示動目標在運動方向上的長度;vt表示動目標的運動速度。同樣,結束遮擋的時刻te可表示為

        te=r-rsvt(10)

        根據(jù)式(9)和式(10)可知,SIP上任意一點P被遮擋的時間t為:t=te-ts=l/vt。由于合成孔徑時間Ts是確定的,因此需要對P點的SIP進行修正,保證一幀圖像內(nèi)的有效遮擋時間在[0,Ts]內(nèi),其修正函數(shù)的表達式為

        C(t)=t, 0≤t≤Ts

        Ts, Ts<t(11)

        設地面成像區(qū)域的強度為1,則可通過式(11)計算SIP:

        S(P)=Ts-C(te-ts)Ts=1-C(te-ts)Ts(12)

        從式(12)中可知,當C(te-ts)與Ts的比值越大時,S(P)越小,即陰影越明顯。根據(jù)Ts=Raλ/2ρav可知,波長越短,雷達載頻就越大,Ts越小,陰影成像的效果越好。此外,目標的運動速度vt越小,陰影成像的效果也越明顯,陰影檢測也越容易。

        2 ViSAR系統(tǒng)研究進展

        2.1 國外研究進展

        SNL于2009年公布了其ViSAR系統(tǒng)的專利。在該專利中,SNL采用孔徑重疊模式來進行ViSAR系統(tǒng)的設計與成像,并且SNL開發(fā)了先進的ViSAR算法,同時擁有多項ViSAR模式的專利,實現(xiàn)了革命性的情報獲取方式[31。同時,SNL在其官網(wǎng)上公布了他們所研發(fā)的ViSAR系統(tǒng)的三段成像結果[33,圖7為公布數(shù)據(jù)中其中一段SAR視頻的截圖。

        圖7為一空軍基地大門的ViSAR成像結果,其成像結果包含了運動車輛通過大門的場景。從圖7(b)中的紅框放大細節(jié)可明顯看出,當運動目標停止移動時,目標的反射能量呈現(xiàn)在陰影頂部。當目標開始移動時,陰影就會再次出現(xiàn)在其真實位置。同時,圖7中位于道路中的陰影全部為移動車輛的真實位置,這也給基于陰影的動目標檢測奠定了理論基礎。此外,圖7也清楚地顯示了基地大門周圍物體的精細細節(jié),展示了ViSAR成像的高分辨以及高幀率能力。

        加拿大國防研究及發(fā)展中心(Defence Research and Development Canada, DRDC)自2007年起開始對ViSAR進行立項,并在2010年對其展開了大量的研究工作,把其研究重心放在如何提高ViSAR成像的圖像質量中。在文獻[14]中,DRDC基于X波段的雷達系統(tǒng),通過孔徑重疊模式來實現(xiàn)ViSAR成像,并且該中心通過對相鄰幀圖像的累積求和來減少散斑以及填充陰影,以達到增加后續(xù)幀圖像對比度的目的。

        DRDC通過圓周軌跡聚束模式共成像了53幀圖像。如圖8所示,第10幀圖像是在圍繞目標成像區(qū)域旋轉約10°后累積形成的,第53幀圖像是在旋轉約55°后累積形成的。從3個成像結果的對比中可明顯得出,經(jīng)過多幀累積之后的SAR圖像噪聲更少、目標的輪廓更加清晰。這也證實了相鄰幀圖像的累積求和能提高成像分辨率。

        2012年,美國國防部高級研究計劃署(Defense Advanced Research Projects Agency,DARPA)的ViSAR項目開始實施,該項目旨在開發(fā)出具有穿透云層及煙塵、同時具有高分辨率和高幀率的ViSAR成像系統(tǒng)[12,29。在短短數(shù)年之間,DARPA聯(lián)合雷神、貝宜系統(tǒng)有限公司等多家美國國防公司設計出了231.5~235 GHz的ViSAR系統(tǒng),并在2018年年底公布了其ViSAR系統(tǒng)成像的結果[34

        從圖9(a)的ViSAR系統(tǒng)結構示意圖可知,DARPA所設計的235 GHz“L”型一發(fā)四收ViSAR系統(tǒng)能在豎直方向上實現(xiàn)干涉,獲得成像區(qū)域的高度信息,同時在水平方向上也可實現(xiàn)動目標的檢測。圖9(b)和圖9(c)為DARPA在德克薩斯州米德洛錫安進行機載實驗的ViSAR成像結果中的兩幀圖像,從成像結果可以清晰地看出各目標以及其陰影的圖像特征。DARPA的ViSAR項目成功地證明了235 GHz雷達系統(tǒng)的可行性,并展示了以低圖像延遲和高SAR幀率生成實時精細分辨率的SAR圖像的能力。

        表1為SNL與DARPA兩個實驗室研發(fā)的ViSAR系統(tǒng)的具體參數(shù)對比,從表中的對比也可以得出高頻ViSAR系統(tǒng)和非高頻ViSAR系統(tǒng)之間的差異。

        美國通用原子航空系統(tǒng)公司(General Atomics Aero-nautical Systems Incorporated, GA-ASI)于2013年開始對圓周軌跡聚束模式下的非高頻ViSAR系統(tǒng)進行研究,同時提出采用反向投影(back projection, BP)算法結合移位寄存器來實現(xiàn)孔徑的循環(huán)位移以及成像處理,最終的交疊幀率達到了4 Hz[35。

        德國弗勞恩霍夫(Fraunhofer-Gesellschaft, FhG)協(xié)會于2013年開展太赫茲關鍵器件研究,并研制了型號為MIRANDA-300的太赫茲雷達成像系統(tǒng)。該系統(tǒng)的工作頻率為300 GHz,其分辨率達到了3.7 mm[36。2016年,F(xiàn)hG研制出了新一代MIRANDA-300 ViSAR成像系統(tǒng),并進行了成像實驗,其雷達系統(tǒng)與成像結果如圖10所示,成像結果也很好地驗證了太赫茲ViSAR系統(tǒng)的高分辨成像能力[37-40

        2017年開始,部分亞洲國家也開始對ViSAR進行仿真驗證與研究。韓國延世大學對多輸入多輸出調(diào)頻連續(xù)波(multiple input and multiple output frequency modulated continuous wave, MIMO-FMCW)ViSAR系統(tǒng)進行了研究,并設計出了多輸入多輸出調(diào)頻連續(xù)波ViSAR(multiple input and multiple output ViSAR, MIMO ViSAR)系統(tǒng)[41,達到了ViSAR體制成像幀率的要求[42。

        本節(jié)對國外ViSAR系統(tǒng)的發(fā)展路徑作了詳細分析,現(xiàn)對其進行總結概括,如表2所示。

        2.2 國內(nèi)研究進展

        ViSAR的概念最早由國外提出,但國內(nèi)各研究機構對其極為重視,系統(tǒng)理論的研究與機載實驗也在相繼開展中,并且已經(jīng)取得了一些實質性進展。

        2014年,文獻[43]設計了一款雷達工作頻率范圍在336.6~343.8 GHz的FMCW-SAR系統(tǒng),其成像結果在方位向上可達到較高的分辨率。圖11為FMCW-SAR系統(tǒng)示意圖,該系統(tǒng)總共由四大模塊組成:信號發(fā)生器、chirp源、收發(fā)器以及信號處理單元。該系統(tǒng)以340 GHz的中心頻率、7.2 GHz的帶寬進行FMCW模式的工作,并且其成像結果驗證了太赫茲雷達系統(tǒng)的成像能力。同時,還對太赫茲雷達的三維成像技術以及太赫茲ViSAR的成像技術進行了相關研究。

        2015年,文獻[44]首次提出了星載ViSAR模式。星載ViSAR系統(tǒng)工作于聚束模式與孔徑重疊模式下,該系統(tǒng)對采集到的回波數(shù)據(jù)按照90%的重疊率進行分割,接著對所有數(shù)據(jù)進行成像處理,最終成功地得到了連續(xù)的視頻幀圖像序列。一年后,文獻[45]在文獻[44]星載ViSAR模式的基礎上設計了基于雙站模式的天基ViSAR系統(tǒng),同時詳細研究了ViSAR成像性能與各參數(shù)之間的關系。該ViSAR系統(tǒng)采用靜止軌道上的SAR衛(wèi)星作為信號發(fā)射源,利用低軌上的SAR順序接收成像區(qū)域的回波信號,這樣能使系統(tǒng)具備較強的反輻射攻擊能力,同時能適用于較多的工作模式。文獻[45]的仿真結果表明,天基ViSAR系統(tǒng)可以較為正確地反映目標的運動情況,能為后續(xù)的基于星載ViSAR的動目標檢測提供重要依據(jù)。

        2016年,文獻[46]通過合理地設計和對緊湊器件的應用,成功地設計了機載FMCW-SAR系統(tǒng)?;谠摳咝阅艿腇MCW-SAR系統(tǒng),通過對原始數(shù)據(jù)進行預處理、擴展有限狀態(tài)自動機(finite state automata, FSA)并行實現(xiàn)以及對圖像數(shù)據(jù)的再處理,實現(xiàn)了ViSAR成像。在FMCW-SAR系統(tǒng)的處理上還進行了準確的耗時測試,并進行了多次機載實驗,獲得了最終的FMCW-SAR視頻成像結果。圖12(a)為FMCW-SAR系統(tǒng)的處理板,在FMCW-SAR視頻的處理過程中,F(xiàn)PGA讀取存儲卡中的原始數(shù)據(jù),并將回波數(shù)據(jù)傳輸?shù)接米鲾?shù)據(jù)緩沖區(qū)的DDR3sdram中。在對幀數(shù)據(jù)進行緩沖后,數(shù)字信號處理芯片通過擴展的FSA處理相應的原始數(shù)據(jù)。處理后的圖像數(shù)據(jù)通過網(wǎng)絡傳輸?shù)接嬎銠C上。最后,在計算機上顯示相應的幀圖像。下一幀的數(shù)據(jù)也以同樣的方式進行處理和傳輸。圖12(b)為FMCW-SAR系統(tǒng)處理后的ViSAR成像結果,其成像場景為800 m×3 000 m的城鎮(zhèn)。從成像結果中可明顯看出,成像區(qū)域道路清晰而連續(xù),同時建筑結構比較集中,耕地面積不同的農(nóng)田也比較容易區(qū)分[47。

        文獻[48-49]分別在2016年詳細分析了條帶式高頻ViSAR與孔徑重疊模式ViSAR在不同平臺速度下的成像幀率、成像分辨率與雷達工作頻率的關系,為后續(xù)開發(fā)ViSAR系統(tǒng)的參數(shù)設計提供了一定參考。

        2018年,文獻[50]設計了一種220 GHz的ViSAR系統(tǒng),該系統(tǒng)采用一發(fā)一收的天線布局,發(fā)射信號為線性調(diào)頻連續(xù)波信號。該系統(tǒng)可達到最大為4.8 GHz的信號帶寬,信號的輸出功率約為20 mW,并可通過成像試驗實現(xiàn)5 Hz的獨立幀率,成像分辨率約為3.2 cm。圖13為該ViSAR系統(tǒng)的示意框圖[50。

        2021年,文獻[51]根據(jù)ViSAR系統(tǒng)的技術指標分析報告計算出了ViSAR系統(tǒng)中數(shù)字處理器件各個模塊的主要參數(shù),并依據(jù)參數(shù)設計出了ViSAR數(shù)字處理組件的PCB板。圖14為所設計的ViSAR數(shù)字處理組件的實物圖。通過測試,該系統(tǒng)可以實現(xiàn)2路SAR信號的產(chǎn)生、4路中頻信號的采集和數(shù)據(jù)預處理,并且可以與多個外置模塊(如慣導、伺服和頻綜等)協(xié)同完成整個合成ViSAR系統(tǒng)的控制、數(shù)據(jù)采集和狀態(tài)監(jiān)控等功能[51

        本節(jié)對國內(nèi)ViSAR系統(tǒng)的發(fā)展作了詳細分析,現(xiàn)對其進行總結概括,如表3所示。

        3 ViSAR處理技術研究進展

        ViSAR的出現(xiàn)從某種意義上來講就是為了更加精確地進行地面機動目標的實時監(jiān)測與跟蹤,并且由于其合成孔徑時間短,有利于動目標的陰影成像,因此基于陰影的動目標檢測是ViSAR的最大應用之一。此外,如何效率更高、更穩(wěn)定地進行高幀率持續(xù)成像,以及對ViSAR系統(tǒng)進行運動補償,也是ViSAR的研究重點。本節(jié)將對ViSAR的成像算法、運動補償以及動目標檢測這3個方面來詳細分析ViSAR處理技術的研究現(xiàn)狀。

        3.1 成像算法研究進展

        目前,聚束模式下的ViSAR成像算法主要有兩大類:時域成像算法以及頻域成像算法。時域成像算法中最為基礎的是反向投影(back projection, BP)算法[35。BP算法雖成像精度高、成像處理過程簡單,但其計算量較大,實時性較差。因此,文獻[52]提出了一種用于ViSAR成像的快速BP(fast BP, FBP)算法。該算法將成像場景進行劃分,分為興趣區(qū)域(region of interest, ROI)和一般區(qū)域(general region, GR),在不同的區(qū)域使用不同的孔徑長度進行成像,重點提升ROI的分辨率,并通過遞歸思想降低ViSAR成像的計算復雜度[52。圖15為迭代FBP算法的說明。

        此外,Bishop等[10利用GPU中單精度比雙精度運算更快的特性,對BP算法中的距離公式進行了二階Taylor展開,并使用Taylor三階系數(shù)估計成像誤差,提高了成像算法的計算效率[10

        根據(jù)文獻[53]的推導,差分范圍ΔR相對于點(x,y,z)=(0,0,0)的Taylor展開式為

        ΔR=ΔR(0)+ΔR(1)+ΔR(2)+ΔR(3)+…(13)

        式中:ΔR(0)=0;ΔR(1)=-xxar-yyar-zzar。當(xa,ya,za)為雷達平臺的實際坐標,r為雷達天線到點目標的距離時,ΔR(2)可表達為

        ΔR(2)=-xyxayar3-yzyazar3-xzzaxar3+x2

        (ya2+za2)r3+y22·(xa2+za2)r3+z22·(ya2+xa2)r3(14)

        最后,用式(13)和式(14)代替平方根,推導出三階項以估計相位誤差。根據(jù)文獻[10]所推導的ΔR(3),ΔR(3)是由常數(shù)項、線性項、二次項以及三次項組成,根據(jù)ΔR(3)所確定的散焦相位誤差為

        ?=4πfc·xy22·xa((xa2-2(ya2+za2))r5(15)

        圖16為文獻[10]在ViSAR系統(tǒng)上使用不同方法計算的實測成像結果:① 是使用平方根計算形成的,在圖像上方的邊緣所成像的道路幾乎為直線;② 為ΔR線性近似(一階項)成像的結果,圖像上方的道路明顯變得扭曲;③ 為使用二階項進行計算所成的圖像,圖像上方的道路也較為筆直,與①難以區(qū)分。這驗證了在圖形處理器上使用ΔR的二階近似能提高數(shù)據(jù)的計算效率,并能提升在大規(guī)模場景下實時處理ViSAR數(shù)據(jù)的能力。

        針對天基ViSAR系統(tǒng)的成像問題,文獻[45]提出一種大視角下重疊孔徑模式的并行運算的濾波BP(filtering BP, FBP)重建算法。該算法首先將回波信號劃分成若干個子孔徑,并且為了實現(xiàn)子孔徑疊加,幀周期需為子孔徑長度的整數(shù)倍。接著,在單個子孔徑內(nèi)進行距離向壓縮和插值,之后進行回波相位補償,最后對每個子孔徑的成像進行并行計算,得到由子孔徑合成的單幀全分辨率圖像[45,54。最后由該算法得到的ViSAR成像結果正確地反映了地面目標的運動情況,能為基于天基ViSAR的運動目標檢測提供依據(jù)[55。

        但上述的成像結果會有主瓣展寬、旁瓣較高的問題。因此,Hawley等[56利用ViSAR相鄰幀數(shù)據(jù)的重疊特性,提出通過對子孔徑圖像進行疊加成像來進行旁瓣抑制的算法。在孔徑重疊模式中,生成ViSAR連續(xù)幀圖像時會存在歷史相位重疊的問題,造成一定的重復處理。Hawley等先使用BP算法從不重疊的子孔徑中形成過采樣的SAR圖像,在得到子孔徑圖像之后再相干地求和產(chǎn)生全分辨率的圖像,這樣可以控制旁瓣的窗口。在圖像求和時對子孔徑進行加權,可以為每個子孔徑產(chǎn)生兩個或更多的加權圖像,這些圖像可以通過線性組合來近似任何期望的孔徑加權。在該方法中,實現(xiàn)了與相位歷史數(shù)據(jù)加權幾乎相同的旁瓣控制,并為每一幀形成一個新的圖像,節(jié)省了計算成本。圖17和圖18分別為連續(xù)子孔徑相干求和、BP成像與加權亞孔徑組合方法成像的比較結果。從兩圖的對比可以看出,圖17中的連續(xù)子孔徑相干求和的結果有著明顯的交叉范圍偽影。圖18(a)為對原始相位歷史進行加權并使用BP算法來形成使用整個加權相位歷史的圖像的結果,圖18(b)為18個子孔徑圖像相干求和的結果,由圖可以清楚地看到周期性的旁瓣偽影,圖18(c)為二次窗口近似進行相干求和的結果,該方法使用方位角上6個子孔徑中的每一個子孔徑的3個加權子孔徑圖像。除了峰值響應的兩側都有一個小的旁瓣偽影之外,很難與圖18(a)進行區(qū)分。

        極坐標格式算法(polar format algorithm, PFA)是ViSAR頻域成像算法中最有代表性的一種,由SNL實驗室的Walker提出[57。頻域算法基于二維匹配濾波的原理,通過對回波信號的快速傅里葉變換(fast Fourier transform, FFT),將其轉化為頻域內(nèi)的復振幅和相位信息,再通過相位補償、插值等處理方法,實現(xiàn)快速成像。PFA則利用平面波假設原理,將信號在波數(shù)域內(nèi)的分布近似為簡單的極坐標分布,再使其頻譜信號重采樣為直角分布,通過二維FFT得到最終的SAR圖像。同時,PFA更適用于聚束模式的SAR成像[58。

        美國SNL的Doerry等[59提出了一種基于一維Chirp-Z變換的極坐標實時成像算法,該算法先通過硬件改變系統(tǒng)的采樣率和采樣時間,避免插值運算,實現(xiàn)實時的運動補償,之后采用Chirp-Z變換對補償后的信號進行成像處理,這樣進一步提高了成像算法的效率,并且適用于高速運動場景。

        在此研究的基礎上,2017年,Zuo等[60通過自主設計的太赫茲調(diào)頻連續(xù)波SAR(frequency modulated continuous wave SAR, FMCW-SAR)系統(tǒng)采集了實驗數(shù)據(jù),并提出了一種基于PFA的ViSAR成像算法,驗證了太赫茲ViSAR的高幀率成像能力[60。文獻[60]在傳統(tǒng)PFA成像算法中的笛卡爾坐標系中預定義了一個二維譜矩陣,然后通過對極坐標系中分布的子孔徑的二維譜進行插值運算,可以將其填充到預定義的矩陣中,即可計算出子孔徑二維譜的最小封閉矩形。最后,對基于最小封閉矩形的預定義矩陣得到的二維譜進行二維FFT,得到笛卡爾坐標系上固定方向的ViSAR幀圖像,這樣就可避免重復的插值計算。圖19為該算法的過程示意圖。

        頻域成像算法中的圖像插值運算會導致成像的效率不高,因此文獻[61]提出了一種僅需二維FFT、無需圖像插值就能實現(xiàn)的ViSAR成像算法。該成像算法的流程為:① 對回波信號進行距離向預處理(去調(diào)頻與去斜處理);② 根據(jù)分辨率匹配的條件,計算分辨率匹配角以及對應的幀孔徑角,并以此計算最長共同子序列(longest common subsequence, LCS)下的幀波數(shù)譜分布,同時生成插值索引表和權值系數(shù)表;③ 按孔徑重疊率和幀孔徑角劃分合適的幀孔徑數(shù)據(jù),并依據(jù)插值表得到LCS下的幀波數(shù)譜;④ 對幀波數(shù)譜進行二維線性Chirp-Z變換(two-dimensional Chirp-Z transform, 2D-CZT),得到單幀的地面輸出坐標系(ground output coordinate system, GOCS)圖像;⑤ 對下一幀數(shù)據(jù)重復上述步驟,得到連續(xù)幀圖像。傳統(tǒng)與改進后ViSAR成像的PFA過程如圖19所示。

        Jiang等[62針對PFA算法中的二維插值影響計算速度這一缺點,提出了一種基于Chirp Scaling的太赫茲ViSAR成像算法。該算法利用了太赫茲ViSAR的合成小角特性和線性調(diào)頻的固有特性,采用兩步Chirp Scaling算法代替PFA中的二維插值,以獲得類似的聚焦效果,成像效率更高。

        英國倫敦大學的Pu等[63為了減少ViSAR的采樣量,將圖像處理表述為一個低Tubal-Rank張量恢復問題。該文獻所提出的算法是對矩陣感知問題中正交矩陣追蹤算法的一種改進,能夠有效重建目標張量,而不需要任何關于預先指定或預先估計的張量秩值的先驗信息。為此,該算法需先迭代估計目標張量的基張量和權值張量,利用線性映射觀測張量與估計張量之間的殘差,并將殘差作為停止條件。最終,證明了該算法在采樣效率方面優(yōu)于現(xiàn)有的ViSAR成像算法。從聯(lián)合采樣和計算效率的角度來看,與其他低秩張量恢復算法相比,該算法是最佳的選擇,并從理論上證明了所提方法的收斂性和正確性。

        3.2 高幀率運動補償算法研究進展

        在實際的成像實驗中,由于載機平臺時常會受到大氣環(huán)境等因素的影響,雷達APC總會偏離理想位置[64。其次,由于ViSAR的工作頻率較高,微小的載機平臺位置誤差也會引起較大的相位誤差,必須考慮ViSAR系統(tǒng)的運動補償[65。本節(jié)對已經(jīng)公開的ViSAR運動補償文獻進行分析及總結。

        Linnehan等[66針對圓周軌跡聚束模式下的ViSAR系統(tǒng)難以聚焦的問題進行了研究,提出通過偽點散射體方法進行相位梯度的估計,解決了載機飛行帶來的場景旋轉和坐標變換的問題,在固定方向取得了良好的聚焦效果,并在圖像被編譯成SAR視頻之前對其進行聚焦。該方法的第一步是對復雜圖像數(shù)據(jù)進行線性相位校正,使偽相位歷史(pseudo-phase history, PPH)域中的能量以域軸為中心。頻譜居中后,將圖像進行旋轉,以便將水平軸定向,并可以應用相位梯度自聚焦(phase gradient autofocusing, PGA)成像算法。同時,將這個圖像方向指定為遠距離像(far range up, FRU)。接著,使用零填充的FFT將圖像轉換到PPH域。然后,使用簡單的最近鄰方法分別旋轉插值頻譜的實部和虛部,再通過逆FFT將旋轉后的數(shù)據(jù)轉換回圖像域,從而不丟失PGA算法所需的相位信息。針對圖像水平方向上的方位角模糊,應用了文獻[27]所描述的傳統(tǒng)PGA算法。經(jīng)過多次迭代后,自動聚焦向量被充分估計,以減輕點目標的模糊響應,此時在PPH域可以使用雙線性插值將聚焦的FRU圖像幅度旋轉回基本方向(cardinal direction up, CDU)。這樣,聚焦的SAR幀的像素將與原始的反向投影過程進行匹配,并可以編譯成一個ViSAR圖像序列。

        針對大斜視情況下的高頻ViSAR振動補償問題,文獻[67]提出了一種適用于斜視模式的振動相位誤差補償成像算法。該算法首先在時域中校正信號的大距離徒動,之后再在頻域中校正因多普勒Keystone變換引起的距離徒動,從而實現(xiàn)對因載機平臺振動引起的相位誤差的運動補償。該算法能在斜視角為30°時對成像場景進行運動補償,并且成像質量較好。

        文獻[68]分析了因載機平臺引起的ViSAR振動的誤差閾值與振動頻率的關系,并提出了針對ViSAR平臺的高頻振動的補償方案。為了減少高頻微小振動的影響,該方案采用了隔振技術,通過在振動源和其他元器件之間的振動傳輸通道上加裝彈性元件以減少振動的影響。接下來,利用平臺上的多種傳感器來測量載機平臺的各種運動數(shù)據(jù),計算出平臺的運動誤差并進行粗補償。最后,基于回波數(shù)據(jù),采用PGA算法進行精補償,完成對相位誤差的估計和補償,同時從復圖像域出發(fā),將散焦的圖像聚焦和加窗,完成對相位梯度的估計與補償。該方案能以較高精度實現(xiàn)ViSAR的運動補償,同時能隔絕大部分的高頻振動。仿真和實際實驗驗證了該補償方案的有效性。

        針對太赫茲ViSAR的運動補償問題,文獻[30]提出了一種適用于高頻ViSAR的運動補償算法,此算法可大幅減少計算復雜度,更加符合ViSAR系統(tǒng)的實時性。該算法在二維重采樣之后進行空變誤差的補償,相較于其他的算法能極大地提高補償算法的計算效率,同時該算法可在二維重采樣時與PFA算法相結合并進行成像,進一步提高了運動補償與成像算法的效率。

        此外,文獻[69]研究了一種ViSAR地面移動目標聚焦的框架,可提供更為清晰的動目標陰影。該框架首先使用BP成像算法獲取SAR視頻,然后使用基于深度學習跟蹤器跟蹤多個移動目標的陰影,并最終通過所提出的移動目標BP(moving target BP, mt-BP)算法對移動目標進行重新聚焦。在該框架下,該算法可以對多個移動目標進行綜合檢測、跟蹤、成像,大大提高了ViSAR系統(tǒng)的動目標檢測能力[69。

        3.3 動目標檢測與跟蹤技術研究進展

        在ViSAR處理技術中,動目標的檢測與跟蹤是最重要的應用場景。在常規(guī)的SAR動目標檢測和跟蹤中,基于SAR的地面動目標檢測(SAR-based ground moving target indication, SAR-GMTI)技術的解決方法是解決SAR動目標檢測與跟蹤的主流方案,但是SAR-GMTI技術存在最小檢測速度限制以及不能檢測靜止目標等問題,此外還對檢測目標的雷達散射截面(radar cross section, RCS)有著一定的要求,因此在實際的戰(zhàn)場情況中有著較大缺陷。研究人員發(fā)現(xiàn),在ViSAR體制中,多普勒頻移會使ViSAR圖像序列中的運動目標在實際位置上留下陰影,且相鄰的SAR幀圖像之間具有很強的相關性。因此,較多的研究人員提出了將ViSAR圖像序列中動目標的陰影作為特征來進行動目標的檢測與跟蹤。

        Jahangir[70提出了一種將動目標陰影作為背景分布中的異常值來進行動目標跟蹤的新穎方法,該方法在較弱的動目標陰影成像中的檢測結果非常優(yōu)異。經(jīng)過實際驗證,該方法可在最遠10 km的距離內(nèi)通過陰影探測到移動的目標。

        美國SNL實驗室的Raynal等[71詳細地分析了ViSAR的靜態(tài)目標以及動目標陰影形成的過程,同時從數(shù)學角度提供了ViSAR的任意運動的非特定目標的陰影尺寸和強度的理論公式。此外,也利用中值濾波算法對 SAR視頻中的相干斑噪聲進行空間濾波,以提高陰影雜波比,并且給出了濾波后計算檢測概率和虛警概率的公式[72

        文獻[73]基于美國SNL實驗室公布的ViSAR數(shù)據(jù)提出了一種低虛警率的ViSAR動目標檢測與跟蹤方法。該方法利用了ViSAR圖像序列的空間與時間的特征,通過圖像序列的數(shù)據(jù)匹配與相干斑濾波進行ViSAR圖像序列的預處理,再利用背景差分以及對稱差分進行動目標陰影的檢測處理與跟蹤,同時利用圖像序列的空間域信息抑制虛警率,實現(xiàn)ViSAR的動目標檢測與跟蹤,最終利用該方法成功檢測動目標的正確率達到了85%。

        文獻[74]采用孔徑重疊模式生成了TerraSAR-X星載SAR收錄的荷蘭某地的ViSAR幀圖像,并以此作為數(shù)據(jù)集,針對星載ViSAR提出了一種ViSAR動目標恒虛警率(constant 1 alarm rate, CFAR)檢測方法。該方法首先對ViSAR每幀的回波數(shù)據(jù)進行成像、幾何校正與圖像匹配等預處理,再利用生成的參考圖像計算ViSAR幀圖像的單幀變化因子,根據(jù)每幀的變化因子確定單幀檢測的虛警概率與檢測門限。緊接著,對每一幀圖像進行檢測,同時根據(jù)潛在運動目標的速度、運動時間、圖像分辨率等確定檢測門限,得到最后的運動目標檢測結果。在一定的條件下,檢測概率可達90%以上。

        為了對低RCS的動目標進行檢測,文獻[75]提出了一種新的低RCS動目標檢測方法,并將其成功應用到了K波段的ViSAR系統(tǒng)中。該方法首先對ViSAR圖像序列進行非相干平均以減少散斑噪聲,接著對圖像序列進行集合校正與圖像配準,然后選擇一組相鄰的圖像計算一個表示靜態(tài)場景的參考圖像,之后將目標區(qū)域的每個SAR圖像與參考圖像分割成比例圖像,抑制圖像的靜態(tài)部分,突出正在發(fā)生變化的部分。最后,采用CFAR檢測方法對動目標陰影進行檢測。為了提高對噪聲和剩余雜波的魯棒性,該方法還采用了聚類算法,將可疑目標分成唯一的聚類,并刪除了那些與動目標陰影不匹配的聚類,提高了動目標檢測的精度。

        文獻[76]利用動目標陰影在ViSAR圖像序列中的強相關性,提出了一種利用背景補償和背景分割來進行動目標檢測的方法。該方法將美國Sandia實驗室的ViSAR數(shù)據(jù)作為數(shù)據(jù)集,先對每幀圖像進行尺度不變特征變換(scale-invariant feature transform, SIFT),并利用隨機抽樣一致性(random sample consensus, RANSAC)算法進行背景補償,之后再采用CattePM(Perona-Malik, PM)擴散模型抑制相干斑噪聲。最后,通過灰度熵提取出動目標的陰影,獲得僅包含陰影的二值圖像,對相鄰幀圖像背景建模并進行差分,結合三幀間差分法分離出多個動目標陰影,實現(xiàn)多目標的運動檢測。該算法最終的精度達到856%。文獻[77]通過建立ViSAR仿真模型,分析了動目標的回波特性,并采用圖像統(tǒng)計技術,成功地檢測出每幀SAR圖像中的運動目標。

        文獻[78]針對目標運動時產(chǎn)生的模糊陰影的問題,提出了一種基于關節(jié)角化相關濾波器(joint kernelized correlation filter, JKCF)的運動目標跟蹤框架。該框架基于JKCF的特征訓練,將ViSAR圖像序列中的陰影與相應的距離-多普勒譜中的能量相結合,對目標進行跟蹤。針對跟蹤漂移和塌陷的問題,采用交互處理方法,基于置信度評估,增強目標定位和特征更新。通過融合初始化和特征更新策略,可以顯著提高跟蹤的成功率和精度。

        上述的傳統(tǒng)目標檢測方法較易受到圖像中復雜背景的干擾,并且有人工提取特征的主觀性。因此,國內(nèi)外較多的研究人員將深度學習的方法應用到ViSAR的動目標檢測中。文獻[79]首次提出了一種基于卷積神經(jīng)網(wǎng)絡(convolution neural network, CNN)的陰影檢測與跟蹤的動目標檢測方法,該方法基于SAR圖像中動目標的陰影特征,對一系列從時間和空間子孔徑SAR圖像中提取出的潛在候選動目標進行CNN跟蹤分類。通過對該算法的仿真實驗和性能分析,驗證了該算法的有效性。

        文獻[80]借助美國Sandia實驗室的ViSAR數(shù)據(jù),提出了一種基于滑窗密度聚類、雙向長短期記憶(bi-long short term memory, Bi-LSTM)網(wǎng)絡與快速區(qū)域CNN(faster region-based CNN, Faster R-CNN)相結合的方法,對動目標陰影進行檢測。該方法利用滑動密度聚類算法和Bi-LSTM網(wǎng)絡減少漏警,并用Faster R-CNN進行動目標陰影的追蹤。表4為基于實測數(shù)據(jù)的檢測性能對比,可以看出文獻[80]所提方法具有較優(yōu)秀的性能,能實現(xiàn)ViSAR動目標檢測與追蹤。文獻[81]基于動目標的能量和陰影的雙域(信號域與圖像域)聯(lián)合檢測思想,分別通過Faster R-CNN和航跡關聯(lián)兩種技術實現(xiàn)了ViSAR動目標檢測與跟蹤,給出了機載實測數(shù)據(jù)處理結果,并進行了詳細分析[81

        針對ViSAR動目標陰影特征不明顯、配準難度大與虛警概率高等問題,文獻[82]提出一種基于改進的Faster R-CNN的ViSAR動目標檢測方法,并通過其課題組自制的Mini-SAR系統(tǒng)獲取的實測ViSAR數(shù)據(jù)驗證了該方法的有效性。該方法在Faster R-CNN的最后一層(卷積池化層)引入特征金字塔結構,以提高對ViSAR圖像特征的利用。接著,使用K-means聚類方法對ViSAR圖像序列進行預處理,得到一組合理的ancho box設置數(shù)值,這能有效地提高該算法針對ViSAR運動目標的檢測率。圖20所示為文獻[82]提出的改進的Faster R-CNN算法的流程圖,最終的檢測結果的正確率達到了87.2%,漏檢率為12.8%。

        此外,文獻[83]將CVPR2019中精度最高的算法“SiamMask”應用于ViSAR圖像的動目標跟蹤。實驗結果表明,SiamMask算法在處理測試數(shù)據(jù)方面的性能優(yōu)于其他傳統(tǒng)算法。同時,由于動目標的陰影不會散焦或移位,SiamMask算法同時可以用于動目標的跟蹤。

        文獻[84]提出了一種新的深度學習網(wǎng)絡 “ShadowDeNet”來對ViSAR進行動目標陰影檢測。該網(wǎng)絡利用了5種主要的框架來保證其優(yōu)越的檢測性能:① 直方圖均衡陰影增強(histogram equalization shadow enhancement, HESE)用于增強陰影顯著性,便于特征提取;② Transformer自注意機制(transformer self-attention mechanism, TSAM)用于聚焦感興趣區(qū)域,抑制雜波干擾;③ 形狀變形自適應學習(shape deformation adaptive learning, SDAL)用于學習移動變形目標陰影以克服運動速度變化;④ 語義引導的自適應學習(semantic guided anchor-adaptive learning, SGAAL)用于優(yōu)化陰影的位置和形狀;⑤ 在線硬挖掘(online hard-example mining, OHEM)用于選擇典型的困難負樣本以提高背景識別能力。同時,進行了廣泛的消融實驗,以確認上述每個模塊的有效性。實驗結果表明,與其他5種方法(Faster R-CNN、特征金字塔網(wǎng)絡(feature pyramid network, FPN)、YOLO(you only look once, YoLO)v3、RetinaNet、CenterNet)相比,Sha-dowDeNet具有最先進的性能,并且ShadowDeNet只是在一個可接受的范圍內(nèi)犧牲了檢測速度。

        4 結束語

        ViSAR可以生成傳統(tǒng)SAR體制沒有的高幀率、高分辨率的連續(xù)SAR圖像序列,并且可將其以更直觀的動態(tài)形式展示出來,獲得成像區(qū)域內(nèi)目標的時間維度信息,提高了對場景的動態(tài)監(jiān)測能力。因此,ViSAR在動目標檢測與跟蹤領域有著強大的應用潛力,是目前各國大力研究的重點。本文首先分析了ViSAR系統(tǒng)的基本原理;然后綜述了國內(nèi)外ViSAR系統(tǒng)與相關處理技術的研究進展。根據(jù)目前國內(nèi)外ViSAR公開的報道以及文獻資料來看,ViSAR在以下幾個方面還需繼續(xù)深入研究。

        (1) 太赫茲ViSAR系統(tǒng)的研究

        因太赫茲波段具備較高的多普勒帶寬,基于太赫茲的ViSAR系統(tǒng)無需子孔徑成像合成,只需很短的合成孔徑時間,就能進行高分辨率、高幀率的全動態(tài)ViSAR圖像序列生成。但由于太赫茲波段的特性也決定了太赫茲ViSAR系統(tǒng)的發(fā)展需要研發(fā)更高靈敏度、更穩(wěn)定的太赫茲頻段器件,提升太赫茲ViSAR系統(tǒng)的可行性;然后需設計更加緊湊合理的太赫茲ViSAR系統(tǒng)結構,使之能更為方便地放置在無人機等輕便平臺。

        (2) ViSAR處理技術的研究

        成像技術、運動補償算法以及動目標檢測與跟蹤技術是目前ViSAR處理技術的核心。為了提升ViSAR成像算法的普適性,還需在ViSAR的大場景成像、曲線模式成像以及大斜視角成像等方面進行深入研究;然后,研發(fā)適合高幀頻的ViSAR運動補償自聚焦算法也是支撐ViSAR技術蓬勃發(fā)展的重點;另外,ViSAR體制可將SAR-GMTI技術與基于動目標陰影的檢測技術進行結合,對地面的非穩(wěn)態(tài)運動目標進行檢測,更加精確地進行地面機動目標的實時監(jiān)測與跟蹤。同時,需加強對ViSAR連續(xù)幀圖像中不同速度的動目標形成的陰影的分析,并結合人工智能算法對動目標進行檢測與跟蹤。

        (3) ViSAR技術拓展研究

        ViSAR系統(tǒng)因需提供連續(xù)的視頻幀圖像,目前只適用于對固定區(qū)域的連續(xù)監(jiān)測上,并應工作于聚束模式下。ViSAR體制雖可提供距離、方位、時間的三維信息,但為了提高觀測精度,可將干涉測量技術或PolSAR技術等與ViSAR技術結合在一起,使ViSAR技術具備更全面的優(yōu)勢。

        參考文獻

        [1]LOVE A. In memory of Carl A[J]. IEEE Antennas and Propagation Society Newsletter, 1985, 27(3): 17-18.

        [2]MOREIRA A, PRATS-IRAOLA P, YOUNIS M, et al. A tutorial on synthetic aperture radar[J]. IEEE Geoscience and Remote Sensing Magazine, 2013, 1(1): 6-43.

        [3]保錚, 邢孟道, 王彤. 雷達成像技術[M]. 北京: 電子工業(yè)出版社, 2005.

        BAO Z, XING M D, WANG T. Radar imaging technology[M]. Beijing: Publishing House of Electronics Industry, 2005.

        [4]李濤, 唐新明, 高小明, 等. SAR衛(wèi)星業(yè)務化地形測繪能力分析與展望[J]. 測繪學報, 2021, 50(7): 891-904.

        LI T, TANG X M, GAO X M, et al. Analysis and outlook of the operational topographic surveying and mapping capability of the SAR satellites[J]Acta Geodaetica et Cartographica Sinica, 2021, 50(7): 891-904.

        [5]方臣, 胡飛, 陳曦, 等. 自然資源遙感應用研究進展[J]. 資源環(huán)境與工程, 2019, 33(4): 563-569.

        FANG C, HU F, CHEN X, et al. Advances in application of remote sensing to natural resources[J]. Resources Environment amp; Engineering, 2019, 33(4): 563-569.

        [6]黃國滿, 程春泉, 趙爭, 等. 機載SAR遙感測圖技術及應用[J]. 測繪科學, 2019, 44(6): 105-113.

        HUANG G M, CHENG C Q, ZHAO Z, et al. Mapping technology and application of airborne SAR[J]. Science of Surveying and Mapping, 2019, 44(6): 105-113.

        [7]范西昆, 何麗娜. 美軍地面戰(zhàn)場態(tài)勢感知系統(tǒng)發(fā)展啟示[J]. 現(xiàn)代雷達, 2018, 40(5): 5-7, 12.

        FAN X K, HE L N. The development of U. S ground battlefield situational awareness system and its implications[J]. Modern Radar, 2018, 40(5): 5-7, 12.

        [8]明婧. 多基線干涉SAR高精度高程反演方法研究[D]. 成都: 電子科技大學, 2020.

        MING J. A research of high-precision elevation reconstruction method on multi-baseline interferometric SAR[D]. Chengdu: University of Electronic Science and Technology of China, 2020.

        [9]楊祥立. 多極化SAR圖像相干斑濾波與變化檢測研究[D]. 武漢: 武漢大學, 2020.

        YANG X L. Speckle reduction and change detection multi-polarization SAR images[D]. Wuhan: Wuhan University, 2020.

        [10]BISHOP E, LINNEHAN R, DOERRY A. Video-SAR using higher order Taylor terms for differential range[C]∥Proc.of the IEEE Radar Conference, 2016.

        [11]PALM S, WAHLEN A, STANKO S, et al. Real-time onboard processing and ground based monitoring of FMCW-SAR videos[C]∥Proc.of the 10th European Conference on Synthetic Aperture Radar, 2014.

        [12]HUANG X J, DING J X, GUO Q H. Unsupervised Image Registration for VIDEO SAR[J]. IEEE Journal of Selected Topics in Applied Earth Observation and Remote Sensing, 2021, 14(1): 1075-1083.

        [13]WALLACE H B. DARPA MMW system programs and how they drive compound semiconductor technology needs[C]∥Proc.of the IEEE Compound Semiconductor Integrated Circuit Symposium, 2015.

        [14]DAMINI A, BALAJI B, PARRY C, et al. A videoSAR mode for the X-band wideband experimental airborne radar[C]∥Proc.of the SPIE-the International Society for Optical Engineering, 2010: 135-145.

        [15]DAMINI A, MANTLE V, DAVIDSON G. A new approach to coherent change detection in VideoSAR imagery using stack averaged coherence[C]∥Proc.of the IEEE Radar Conference, 2013.

        [16]WELLS L, SORENSEN K, DOERRY A, et al. Developments in SAR and IFSAR systems and technologies at sandia national laboratories[C]∥Proc.of the IEEE Aerospace Conference, 2003: 1085-1095.

        [17]WALLACE H B. Development of a video SAR for FMV through clouds[C]∥Proc.of the SPIE Open Architecture/Open Business Model Net-Centric Systems and Defense Transformation, 2015: 64-65.

        [18]肖忠源, 張振華, 馬曉萌. 視頻SAR發(fā)展現(xiàn)狀研究[C]∥第三屆航天電子戰(zhàn)略研究論壇, 2017: 37-40.

        XIAO Z Y, ZHANG Z H, MA X M. Research on the development of video SAR[C]∥Proc.of the 3rd Aerospace Electronics Strategy Research Forum, 2017: 37 -40.

        [19]CERUTTI-MAORI D, KLARE J, BRENNER A R, et al. Wide-area traffic monitoring with the SAR/GMTI system PAMIR[J]. IEEE Trans.on Geoscience and Remote Sensing, 2008, 46(10): 3019-3030.

        [20]ROSEN P A, HENSLEY S, JOUGHIN I R, et al. Synthetic aperture radar interferometry[J]. Proceedings of the IEEE, 2000, 88(3): 333-382.

        [21]YAN H, MAO X H, ZHANG J D, et al. Frame rate analysis of video synthetic aperture radar (ViSAR)[C]∥Proc.of the IEEE International Symposium on Antennas and Propagation, 2016: 446-447.

        [22]邢孟道, 林浩, 陳濺來, 等. 多平臺合成孔徑雷達成像算法綜述[J]. 雷達學報, 2019, 8(6): 732-757.

        XING M D, LIN H, CHEN J L, et al. A review of imaging algorithms in multi-platform-borne synthetic aperture radar[J]. Journal of Radars, 2019, 8(6): 732-757.

        [23]范文娜. 基于電磁模型的大斜視SAR成像研究[D]. 西安: 西安電子科技大學, 2019.

        FAN W N. Research on highly squinted SAR imaging based on electromagnetic model[D]. Xi’an: Xidian University, 2019.

        [24]HU R Z, MIN R, PI Y M. A video-SAR imaging technique for aspect-dependent scattering in wide angle[J]. IEEE Sensors Journal, 2017, 17(12): 3677-3688.

        [25]胡睿智. 視頻合成孔徑雷達成像理論與關鍵技術研究[D]. 成都: 電子科技大學, 2018.

        HU R Z. Research on imaging theory and key technology of video synthetic aperture radar[D]. Chengdu: University of Electronic Science and Technology of China, 2018.

        [26]SOUMEKH M. Synthetic aperture radar signal processing[M]. New York: Wiley, 1999.

        [27]JAKOWATZ C V J, WAHL D E, EICHEL P H, et al. Spotlight-mode synthetic aperture radar: a signal processing approach: a signal processing approach[M]. New York: Springer, 2012.

        [28]HARMONY D W, BICKEL D L, MARTINEZ A. A velocity independent continuous tracking radar concept[R]. Albuquerque: Sandia National Lab, 2011.

        [29]WALLACE H, GORMAN J, MALONEY P. Video synthetic aperture radar (ViSAR)[R]. Arlington: Defense Advanced Research Projects Agency, 2012.

        [30]左峰. 視頻合成孔徑雷達成像算法研究[D]. 成都: 電子科技大學, 2019.

        ZUO F. Research on video synthetic aperture radar imaging algorithm[D]. Chengdu: University of Electronic Science and Technology of China, 2019.

        [31]BIELEK T P, THOMPSON D G, WALKER B C. Synthetic aperture design for increased SAR image Rate[P]. U.S.: Patent 7498968, 2009-03-03.

        [32]GU C F, CHANG W G. An efficient geometric distortion correction method for SAR video formation[C]∥Proc.of the IEEE International Conference on Modern Circuits and Systems Technologies, 2016.

        [33]SANDIA NATIONAL LABORATORIES. Eubank gate and traffic VideoSAR[EB/OL]. [2023-01-10]. https:∥www.sandia.gov/radar/video/index.html.

        [34]KIM S H, FAN R, DOMINSKI F. ViSAR: a 235 GHz radar for airborne applications[C]∥Proc.of the IEEE Radar Confe-rence, 2018: 1549-1554.

        [35]MILLER J, BISHOP E, DOERRY A. An application of back projection for video SAR image formation exploiting a subaperature circular shift register[C]∥Proc.of SPIE-the International Society for Optical Engineering, 2013: 66-79.

        [36]STANKO S, CARIS M, WAHLEN A, et al. Millimeter resolution with radar at lower terahertz[C]∥Proc.of the IEEE International Radar Symposium, 2013: 235-238.

        [37]CARIS M, STANKO S, PALM S, et al. 300 GHz radar for high resolution SAR and ISAR applications[C]∥Proc.of the IEEE International Radar Symposium, 2015: 577-580.

        [38]STANKO S, PALM S, SOMMER R, et al. Millimeter resolution SAR imaging of infrastructure in the lower THz region using MIRANDA-300[C]∥Proc.of the IEEE European Microwave Conference, 2016: 1505-1508.

        [39]PALM S, SOMMER R, CARIS M, et al. Ultra-high resolution SAR in lower terahertz domain for applications in mobile mapping[C]∥Proc.of the IEEE German Microwave Confe-rence, 2016: 205-208.

        [40]PALM S, SOMMER R, TESSMANN A, et al. Ultra-high resolution imaging of facades and vertical infrastructure by carborne SAR and airborne CSAR[J]. ISPRS Annals of Photogrammetry, Remote Sensing and Spatial Information Sciences, 2019: 129-136.

        [41]KIM S, YU J, JEON S Y, et al. Signal processing for a multiple-input, multiple-output (MIMO) video synthetic aperture radar (SAR) with beat frequency division frequency-modulated continuous wave (FMCW)[J]. Remote Sensing, 2017, 9(5): 491.

        [42]KIM S. SAR video generation of MIMO video SAR with beat frequency division FMCW[C]∥Proc.of the IEEE International Conference on Signal Processing and Communication Systems, 2017.

        [43]ZHANG B, PI Y M, LI J. Terahertz imaging radar with inverse aperture synthesis techniques: system structure, signal processing, and experiment results[J]. IEEE Sensors Journal, 2014, 15(1): 290-299.

        [44]ZHAO S T, CHEN J, YANG W, et al. Image formation method for spaceborne video SAR[C]∥Proc.of the IEEE Asia-Pacific Conference on Synthetic Aperture Radar, 2015: 148-151.

        [45]梁健, 張潤寧, 包敏鳳. 天基視頻SAR系統(tǒng)設計及成像算法研究[J]. 中國空間科學技術, 2016, 36(6): 22-28.

        LIANG J, ZHANG R N, BAO M F. Research on spaceborne video SAR system design and image formation algorithm[J]. Chinese Space Science and Technology, 2016, 36(6): 22-28.

        [46]GU C F, CHANG W G, LI X Y, et al. The formation of high-resolution FMCW SAR video[C]∥Proc.of the IEEE Progress in Electromagnetic Research Symposium, 2016: 496-499.

        [47]GU C F, CHANG W G. An efficient geometric distortion correction method for SAR video formation[C]∥Proc.of the IEEE International Conference on Modern Circuits and Systems Technologies, 2016.

        [48]宋曉燊, 禹衛(wèi)東. 條帶式VideoSAR參數(shù)依賴關系的推導及應用[J]. 中國科學院大學學報, 2016, 33(1): 121-127.

        SONG X S, YU W D. Derivation and application of stripmap VideoSAR parameter relations[J]. Journal of University of Chinese Academy of Sciences, 2016, 33(1): 121-127.

        [49]YAN H, MAO X H, ZHANG J D, et al. Frame rate analysis of video synthetic aperture radar (ViSAR)[C]∥Proc.of the IEEE International Symposium on Antennas and Propagation, 2016: 446-447.

        [50]李大圣, 吳福偉, 孫俊, 等. 一種220 GHz波段太赫茲合成孔徑成像雷達[J]. 微波學報, 2018, 34(4): 55-59.

        LI D S, WU F W, SUN J, et al. Terahertz synthetic aperture imaging radar based on 220 GHz band[J]. Journal of Microwave Science, 2018, 34(4): 55-59.

        [51]鄭南南. 視頻SAR信號產(chǎn)生與采集處理組件設計[D]. 西安: 西安電子科技大學, 2021.

        ZHENG N N. Design of video SAR Signal Generation and acquisition module[D]. Xi’an: Xidian University, 2021.

        [52]SONG X S, YU W D. Processing video-SAR data with the fast back projection method[J]. IEEE Trans.on Aerospace and Electronic Systems, 2016, 52(6): 2838-2848.

        [53]RIGLING B D, MOSES R L. Taylor expansion of the differential range for monostatic SAR[J]. IEEE Trans.on Aerospace and Electronic Systems, 2005, 41(1): 60-64.

        [54]MOSES R L, ASH J N. An autoregressive formulation for SAR back projection imaging[J]. IEEE Trans.on Aerospace and Electronic Systems, 2011, 47(4): 2860-2873.

        [55]梁健, 張潤寧. GEO-LEO雙站視頻SAR系統(tǒng)若干問題研究[J]. 現(xiàn)代雷達, 2017, 39(3): 17-20.

        LIANG J, ZHANG R N. A study on key technologies of spaceborne video SAR system based on GEO-LEO bistatic model[J]. Modern Radar, 2017, 39(3): 17-20.

        [56]HAWLEY R W, GARBER W L. Aperture weighting technique for video synthetic aperture radar[C]∥Proc.of SPIE-the International Society for Optical Engineering, 2011: 67-73.

        [57]WALKER J L. Range-Doppler imaging of rotating objects[J]. IEEE Trans.on Aerospace and Electronic Systems, 1980, AES-16(1): 23-52.

        [58]ZHU D Y, ZHU Z D. Range resampling in the polar format algorithm for spotlight SAR image formation using the Chirp-Z-transform[J]. IEEE Trans.on Signal Processing, 2007, 55(3): 1011-1023.

        [59]DOERRY A W. Forming rotated SAR images by real-time motion compensation[R]. Albuquerque: Sandia National Lab, 2012.

        [60]ZUO F, LI J. A persistent imaging method for video-SAR in terahertz band[C]∥Proc.of the IEEE International Applied Computational Electromagnetics Society Symposium, 2017.

        [61]HU R Z, MIN R, PI Y M. Interpolation-free algorithm for persistent multi-frame imaging of video-SAR[J]. IET Radar, Sonar amp; Navigation, 2017, 11(6): 978-986.

        [62]JIANG J W, LI Y W, ZHENG Q B. A THz video SAR imaging algorithm based on chirp scaling[C]∥Proc.of the CIE International Conference on Radar, 2021: 656-660.

        [63]PU W, WU J J, HUANG Y L, et al. ORTP: a video SAR imaging algorithm based on low-tubal-rank tensor recovery[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2022, 15(1): 1293-1308.

        [64]KANTOR J M. Minimum entropy autofocus correction of residual range cell migration[C]∥Proc.of the IEEE Radar Conference, 2017: 11-16.

        [65]JIA G W, CHANG W G, ZHANG Q L, et al. The analysis and realization of motion compensation for circular synthetic aperture radar data[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2016, 9(7): 3060-3071.

        [66]LINNEHAN R, MILLER J, BISHOP E, et al. An autofocus technique for video-SAR[C]∥Proc.of the 20th Conference on Algorithms for Synthetic Aperture Radar Imagery, 2013: 56-65.

        [67]孫偉, 孫進平, 張遠, 等. 大斜視直升機載太赫茲ViSAR振動補償成像算法[J]. 北京航空航天大學學報, 2016, 42(12): 2755-2761.

        SUN W, SUN J P, ZHANG Y, et al. High squint helicopter-borne terahertz ViSAR vibration compensation imaging algorithm[J]. Journal of Beijing University of Aeronautics and Astronautics, 2016, 42(12): 2755-2761.

        [68]趙雨露, 張群英, 李超, 等. 視頻合成孔徑雷達振動誤差分析及補償方案研究[J]. 雷達學報, 2015, 4(2): 230-239.

        ZHAO Y L, ZHANG Q Y, LI C, et al. Vibration error analysis and motion compensation of video synthetic aperture radar[J]. Journal of Radars, 2015, 4(2): 230-239.

        [69]YANG X Q, SHI J, ZHOU Y Y, et al. Ground moving target tracking and refocusing using shadow in video-SAR[J]. Remote Sensing, 2020, 12(18): 3083.

        [70]JAHANGIR M. Moving target detection for synthetic aperture radar via shadow detection[C]∥Proc.of the IET International Conference on Radar Systems, 2007.

        [71]RAYNAL A M, BICKEL D L, DOERRY A W. Stationary and moving target shadow characteristics in synthetic aperture radar[C]∥Proc.of the 18th Conference on Radar Sensor Technology, 2014: 413-427.

        [72]RAYNAL A M, DOERRY A W, MILLER J A, et al. Shadow probability of detection and 1 alarm for median-filtered SAR imagery[R]. Albuquerque: Sandia National Lab, 2014.

        [73]何志華, 陳興, 于春銳, 等. 一種穩(wěn)健的視頻SAR動目標陰影檢測與跟蹤處理方法[J]. 電子與信息學報, 2022, 44(11): 3882-3890.

        HE Z H, CHEN X, YU C R, et al. A robust moving target shadow detection and tracking method for VideoSAR[J]. Journal of Electronics and Information Science, 2022, 44(11): 3882-3890.

        [74]李錦偉, 張升, 李財品. 一種星載視頻SAR動目標恒虛警檢測方法[C]∥第六屆高分辨率對地觀測學術年會, 2019: 58-80.

        LI J W, ZHANG S, LI C P. A CFAR detection method for spaceborne video SAR moving targets[C]∥Proc.of the 6th Annual Conference of High Resolution Earth Observation, 2019: 58 -80.

        [75]WANG H, CHEN Z S, ZHENG S C. Preliminary research of low-RCS moving target detection based on Ka-band video SAR[J]. IEEE Geoscience and Remote Sensing Letters, 2017, 14(6): 811-815.

        [76]聊蕾, 左瀟麗, 云濤, 等. 基于圖像序列的VideoSAR動目標檢測方法[J]. 雷達科學與技術, 2016, 14(6): 563-567, 573.

        LIAO L, ZUO X L, YUN T, et al. An approach to detect moving target in VideoSAR imagery sequence[J]. Radar Science and Technology, 2016, 14(6): 563-567, 573.

        [77]GOU L T, ZHU D Y, LI Y. A novel moving target detection method for VideoSAR[C]∥Proc.of the IEEE International Applied Computational Electromagnetics Society Symposium, 2019.

        [78]ZHONG C, DING J S, ZHANG Y H. Video SAR moving target tracking using joint kernelized correlation filter[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2022, 15: 1481-1493.

        [79]ZHANG Y, YANG S Y, LI H B, et al. Shadow tracking of moving target based on CNN for video SAR system[C]∥Proc.of the IEEE International Geoscience and Remote Sensing Symposium, 2018: 4399-4402.

        [80]丁金閃. 視頻SAR成像與動目標陰影檢測技術[J]. 雷達學報, 2020, 9(2): 321-334.

        DING J S. Focusing algorithms and moving target detection based on video SAR[J]. Journal of Radars, 2020, 9(2): 321-334.

        [81]丁金閃, 仲超, 溫利武, 等. 視頻合成孔徑雷達雙域聯(lián)合運動目標檢測方法[J]. 雷達學報, 2022, 11(3): 313-323.

        DING J S, ZHONG C, WEN L W, et al. Joint detection of moving target in video synthetic aperture radar[J]. Journal of Radars, 2022, 11(3): 313-323.

        [82]閆賀, 黃佳, 李睿安, 等. 基于改進快速區(qū)域卷積神經(jīng)網(wǎng)絡的視頻SAR運動目標檢測算法研究[J]. 電子與信息學報, 2021, 43(3): 615-622.

        YAN H, HUANG J, LI R A, et al. Research on video SAR moving target detection algorithm based on improved faster region-based CNN[J]. Journal of Electronics and Information Technology, 2021, 43(3): 615-622.

        [83]LIANG Z H, LIANG C J, ZHANG Y, et al. Tracking of moving target based on SiamMask for video SAR system[C]∥Proc.of the IEEE International Conference on Signal, Information and Data Processing, 2019.

        [84]BAO J Y, ZHANG X L, ZHANG T W, et al. ShadowDeNet: a moving target shadow detection network for video SAR[J]. Remote Sensing, 2022, 14(2): 320.

        作者簡介

        顏上?。?997—),男,博士研究生,主要研究方向為雷達智能信息處理。

        付耀文(1976—),男,研究員,博士,主要研究方向為信息融合、雷達信號處理與目標識別。

        張文鵬(1989—),男,講師,博士,主要研究方向為雷達信號處理、時頻分析、信號分解。

        楊 威(1985—),男,副教授,博士,主要研究方向為信息融合、多目標追蹤。

        余若峰(1993—),男,博士研究生,主要研究方向為智能感知與處理。

        張法桐(1999—),男,碩士研究生,主要研究方向為微小型無人機載SAR成像及運動補償。

        无码国产精品色午夜| 久久精品中文字幕大胸| 久久精品国产www456c0m | 欧美自拍视频在线| 日本第一区二区三区视频| 中文字幕综合一区二区三区| 夜夜爽妓女8888888视频| 国产成人精品一区二区视频| 麻豆国产成人AV网| 女女同女同一区二区三区| 欧美人与禽z0zo牲伦交| 日本强好片久久久久久aaa| 国产精品久久一区性色a| 青青草小视频在线观看| 欧美人牲交| 最近高清中文在线字幕观看| 资源在线观看视频一区二区| 中文字幕中文字幕在线中二区| 国产va免费精品高清在线观看| 日韩毛片在线看| 久久久精品国产亚洲av网| 玖玖色玖玖草玖玖爱在线精品视频 | 久久爱91精品国产一区| 电驱蚊液可以插一晚上吗| 欧美性猛交xxxx富婆| 亚洲mv国产精品mv日本mv| 亚洲国产综合精品一区最新| 婷婷五月深深久久精品| 亚洲黄色av一区二区三区| 少妇扒开毛茸茸的b自慰| 日韩AV不卡六区七区| 日本一区二区久久精品亚洲中文无| 免费人成视频网站在在线| 高清破外女出血av毛片| 国产成人免费一区二区三区| 亚洲av天堂在线免费观看| 欧美精品videosse精子| 无码视频一区二区三区在线观看| 女优免费中文字幕在线| 97人妻精品一区二区三区男同| 精产国品一二三产区m553麻豆|