專利名稱：一種逆合成孔徑雷達(dá)的目標(biāo)轉(zhuǎn)速估計方法
技術(shù)領(lǐng)域：
本發(fā)明涉及雷達(dá)技術(shù)領(lǐng)域，更具體地，本發(fā)明涉及一種逆合成孔徑雷達(dá)的目標(biāo)轉(zhuǎn)
速估計方法。
背景技術(shù)：
逆合成孔徑雷達(dá)(ISAR)是一種微波成像技術(shù)，可以提供非合作運(yùn)動目標(biāo)(例如飛機(jī)、艦船)的二維圖像，增強(qiáng)后續(xù)的目標(biāo)識別能力。ISAR通過發(fā)射寬帶信號提供沿雷達(dá)視線方向(距離向)的高分辨能力，通過對雷達(dá)和目標(biāo)之間相對視角的變化進(jìn)行相干積累來提供橫向分辨能力。當(dāng)所成圖像積累的目標(biāo)視角變化較小時，采用高效的FFT處理可獲得目標(biāo)在距離-多普勒(RD)域的圖像。 由于ISAR通常針對非合作目標(biāo)成像，因此，目標(biāo)等效轉(zhuǎn)動的估計(即雷達(dá)和目標(biāo)之間相對視角的變化)非常重要，這是因為存在下述兩個方面的需求。首先，ISAR圖像定標(biāo)的需要；雖然ISAR所成圖像的距離向分辨率取決于雷達(dá)發(fā)射信號的帶寬(其通常是一個已知量)；但是圖像的橫向分辨率取決于相干積累時間內(nèi)相對視角的變化量，對于非合作目標(biāo)成像而言，該變化量通常是未知的。只有對該視角變化量進(jìn)行準(zhǔn)確的估計，才能夠?qū)SAR圖像進(jìn)行橫向定標(biāo)，得到距離-跨距離域的圖像，以用于后續(xù)目標(biāo)識另lj。其次，ISAR成像運(yùn)動補(bǔ)償?shù)男枰划?dāng)成像所積累視角變化較大時，可能出現(xiàn)跨分辨單元徙動(MTRC)問題，引起圖像模糊；同時，當(dāng)存在非勻速的目標(biāo)轉(zhuǎn)動時，也會產(chǎn)生圖像散焦等問題，這些都需要根據(jù)目標(biāo)的轉(zhuǎn)動信息進(jìn)行轉(zhuǎn)動補(bǔ)償，提高成像質(zhì)量。 現(xiàn)有的ISAR轉(zhuǎn)動估計的典型方法主要有軌跡擬合法、圖像質(zhì)量法和基于高階相位項分析的方法和基于光流場分析的方法。軌跡擬合法通過窄帶雷達(dá)提供的目標(biāo)跟蹤數(shù)據(jù)，計算運(yùn)動目標(biāo)相對雷達(dá)的視角變化，通常估計精度較低。圖像質(zhì)量法采用圖像評價函數(shù)來優(yōu)化未知的轉(zhuǎn)角參數(shù)，當(dāng)該圖像評價函數(shù)獲得極值時(此時圖像的聚焦效果最好)，獲得相應(yīng)的轉(zhuǎn)角估計，由于此類方法需要在解空間內(nèi)進(jìn)行反復(fù)搜索，每一次搜索都需要復(fù)雜的成像處理，例如極坐標(biāo)格式成像、巻積逆投影成像，計算量較大。 典型的基于高階相位項的方法諸如特顯點處理(PPP)的方法，該方法將目標(biāo)上的回波信號近似為線性調(diào)頻信號(LFM)，通過估計LFM的起始頻率和調(diào)頻率，并結(jié)合三個散射點所在的距離向位置信息，完成平動補(bǔ)償、轉(zhuǎn)動補(bǔ)償和轉(zhuǎn)角估計，該類方法的主要問題是在實際情況中難以找到穩(wěn)定可靠的高質(zhì)量孤立散射點。對于缺乏特顯點的目標(biāo)，可以引入時頻分析(TFA)的方法提取高階相位項，但TFA類方法對目標(biāo)回波信號的分析都是針對距離單元內(nèi)的回波進(jìn)行的，而散射點回波駐留在同一距離單元內(nèi)的時間有限(特別是對于高分辨雷達(dá))；因而，信號的高階相位系數(shù)通常很小，參數(shù)提取困難，參數(shù)提取及其后的計算量較大，精度較低。

發(fā)明內(nèi)容
為克服現(xiàn)有的ISAR圖像橫向定標(biāo)方法計算量大、精度低以及強(qiáng)散射中心提取困
4難的缺陷，本發(fā)明提出一種基于相鄰圖像旋轉(zhuǎn)相關(guān)的目標(biāo)轉(zhuǎn)速估計方法。
根據(jù)本發(fā)明的一個方面，提供一種逆合成孔徑雷達(dá)的目標(biāo)轉(zhuǎn)速估計方法，包括 步驟10)、根據(jù)接收的回波數(shù)據(jù)，獲取兩幅逆合成孔徑雷達(dá)的距離多普勒圖像； 步驟20)、通過平移相關(guān)來對齊所述兩幅距離多普勒圖像，并獲取所述兩幅距離多
普勒圖像中的目標(biāo)區(qū)域圖像； 步驟30)、對兩幅所述目標(biāo)區(qū)域距離多普勒圖像進(jìn)行旋轉(zhuǎn)相關(guān)，根據(jù)旋轉(zhuǎn)相關(guān)系數(shù)，估計目標(biāo)的轉(zhuǎn)速。其中，步驟IO)中，對所述回波數(shù)據(jù)進(jìn)行滑窗成像，獲取兩幅距離-多普勒圖像序
列；其中，使用滑窗寬度確定每幅距離多普勒圖像積累的脈沖數(shù)，滑窗步長確定兩幅圖像之間的視角差。 其中，步驟20)還包括通過平移相關(guān)方法消除目標(biāo)在距離多普勒平面上的位置偏移，對齊所述兩幅距離多普勒圖像。 其中，步驟20)中還包括裁剪目標(biāo)所在的區(qū)域的距離多普勒圖像的邊緣區(qū)域。其中，步驟30)中，所述兩幅目標(biāo)區(qū)域距離多普勒圖像之間的對應(yīng)關(guān)系為
碌2〕
isin(a7。
cos (《
=R
'浙其中，(Xc(tm》，Yc(U))和(Xc(U，Yc(U)分別表示相同散射中心在兩幅目標(biāo)區(qū)
域距離多普勒圖像中的位置，edif = "。(tm2-tj表示所述兩幅圖像之間的視角差，nr和na分別為距離多普勒圖像的距離向和方位向尺度因子，R( e dif)限定兩幅圖像之間的旋轉(zhuǎn)關(guān)系，其中，R(&小SR(d)S-'=
77r 、"7。
sin(&,) cos(~) 其中，步驟30)中，所述兩幅目標(biāo)區(qū)域距離多普勒圖像的旋轉(zhuǎn)操作可通過如下三個一維的頻域巻積插值實現(xiàn)，
1
A tan1
浙
1
"。 〃
1
tan
《
其中，步驟30)中，所述兩幅目標(biāo)區(qū)域距離多普勒圖像之間的相關(guān)系數(shù)定義為
i:ix,u(x,力)/2(x,力 其中，fjx， y)和f2(x， y)為兩幅圖像的幅度，Trat ()表示所述旋轉(zhuǎn)操作。 其中，步驟30)中，所述旋轉(zhuǎn)相關(guān)系數(shù)的極大值點對應(yīng)估計的目標(biāo)轉(zhuǎn)速，獲得目標(biāo)
轉(zhuǎn)速和旋轉(zhuǎn)方向的估計。
5
其中，上述方法進(jìn)一步包括 步驟40)、根據(jù)估計的目標(biāo)轉(zhuǎn)速對距離_多普勒圖像進(jìn)行橫向定標(biāo)，所述橫向標(biāo)度
因子為化^^T，其中，A為載波波長，^為脈沖重復(fù)頻率，M為圖像積累的脈沖數(shù)目，"。為目標(biāo)轉(zhuǎn)速。 通過應(yīng)用本發(fā)明，對ISAR系統(tǒng)接收的回波數(shù)據(jù)進(jìn)行分段成像得到目標(biāo)轉(zhuǎn)速的估計，不需要額外的系統(tǒng)硬件；通過三級的一維巻積插值過程實現(xiàn)高質(zhì)量的圖像旋轉(zhuǎn)，可以極大地減少計算量，提高轉(zhuǎn)速估計的速度；同時，該方法對目標(biāo)散射特性沒有具體的要求，應(yīng)用范圍較大。


圖1為根據(jù)本發(fā)明的實施例的飛機(jī)目標(biāo)模型； 圖2為根據(jù)本發(fā)明的實施例的目標(biāo)轉(zhuǎn)速估計流程圖； 圖3為基于滑窗的圖像序列獲取示意圖； 圖4為逆合成孔徑雷達(dá)信號處理流程示意圖； 圖5為基于旋轉(zhuǎn)體模型的逆合成孔徑雷達(dá)成像示意圖； 圖6為頻率的巻積插值實現(xiàn)過程示意圖； 圖7為逆時針旋轉(zhuǎn)的目標(biāo)轉(zhuǎn)速搜索結(jié)果示意圖； 圖8為順時針方向旋轉(zhuǎn)的目標(biāo)搜索結(jié)果示意圖； 圖9為采用1024脈沖的未定標(biāo)距離多普勒圖像示意圖； 圖10為采用1024脈沖的定標(biāo)距離多普勒圖像示意圖； 圖11為采用1024脈沖的極坐標(biāo)格式定標(biāo)圖像示意圖。
具體實施例方式
下面結(jié)合附圖和具體實施例，對本發(fā)明提供的一種逆合成孔徑雷達(dá)的目標(biāo)轉(zhuǎn)速估計方法作進(jìn)一步說明。 在根據(jù)本發(fā)明的一個實施例中，本發(fā)明所述的方法針對我國某實驗ISAR系統(tǒng)提供的數(shù)據(jù)進(jìn)行處理，在本實施例中，圖1給出的飛機(jī)目標(biāo)模型，作為本實施例中定標(biāo)的對象。通常，本領(lǐng)域內(nèi)普通技術(shù)人員可以認(rèn)識到本發(fā)明不限于圖中的飛機(jī)目標(biāo)，還包括諸如艦船等其它非合作目標(biāo)。所述實驗ISAR系統(tǒng)，雷達(dá)載波頻率為5. 52GHz ;發(fā)射線性調(diào)頻信號，帶寬400MHz ;采用Dechirp方式獲取脈沖壓縮信號，并經(jīng)正交I/Q雙通道采樣，采樣頻率10MHz ;等效脈沖重復(fù)頻率200Hz (經(jīng)過了脈沖抽取)；獲取2048個脈沖進(jìn)行處理。
圖2為根據(jù)本發(fā)明實施例的目標(biāo)轉(zhuǎn)速估計方法總體流程圖。如圖2所示，首先，對接收到的原始回波數(shù)據(jù)進(jìn)行平動補(bǔ)償，獲取平動補(bǔ)償后的脈壓數(shù)據(jù)；對所述脈壓數(shù)據(jù)使用滑窗方式成像，獲取兩幅距離-多普勒(RD)圖像；在根據(jù)本發(fā)明實施例所述的方法中，該兩幅RD圖像可以有脈沖重疊，也可以沒有脈沖重疊；然后，通過平移相關(guān)方法對齊兩幅RD圖像，并選取目標(biāo)所在的距離單元中心區(qū)域作為等效旋轉(zhuǎn)中心；定義需要搜索的目標(biāo)轉(zhuǎn)速范圍，并根據(jù)已知參數(shù)得到相應(yīng)的RD圖像橫向尺度因子和兩幅圖像之間的視角差；旋轉(zhuǎn)其中的一幅圖像(順時針或者逆時針方向)，將旋轉(zhuǎn)圖像和另一幅圖像作相關(guān)；根據(jù)相關(guān)系數(shù)的
6最大值確定目標(biāo)的旋轉(zhuǎn)速度和旋轉(zhuǎn)方向，并對傳統(tǒng)RD圖像橫向定標(biāo)或進(jìn)行PFA/CBP成像。 在下文中，將對上述方法進(jìn)行詳細(xì)描述。
平動補(bǔ)償及滑窗成像 首先選擇目標(biāo)平穩(wěn)運(yùn)動段的數(shù)據(jù)，如圖3所示，對接收到的脈沖采用滑窗方式獲 取兩幅RD圖像?；诨安介L和滑窗大小，兩幅圖像之間可以有脈沖重疊，滑窗大小是 指每幅圖像中相干積累的脈沖數(shù)。在本實施例中，設(shè)定滑窗寬度為1024脈沖，滑窗步長為 1024脈沖，平分所有脈沖獲取兩幅RD圖像。 圖4為逆合成孔徑雷達(dá)信號處理流程示意圖，首先通過平動補(bǔ)償去除目標(biāo)相對于 雷達(dá)的平動成分，得到圖5所示的轉(zhuǎn)臺目標(biāo)運(yùn)動模型。在圖4中，對目標(biāo)的平動補(bǔ)償通常分 兩步進(jìn)行，包括回波包絡(luò)對齊和平動初相校正。較穩(wěn)健的包絡(luò)對齊方法主要有基于相鄰脈 沖積累的回波相關(guān)法、全局最小熵法等；較穩(wěn)健的平動初相校正方法主要有多普勒質(zhì)心跟 蹤(DopplerCentroid Tracking, DCT)法及其改進(jìn)算法，如基于圓移位處理的多普勒質(zhì)心跟 蹤(Circular Shifting based DCT， CS-DCT)法等。 在平面波照射下，通過處理目標(biāo)的寬帶回波獲得沿雷達(dá)視線方向(距離向)的高
分辨；通過對回波相干積累處理回波的多普勒頻率信息，形成跨距離向高分辨。 圖5為基于旋轉(zhuǎn)體模型的逆合成孔徑雷達(dá)成像示意圖，如圖5所示，假設(shè)在較短的
時間間隔內(nèi)經(jīng)過平動補(bǔ)償之后目標(biāo)以"。勻速轉(zhuǎn)動，則目標(biāo)上一點(x。，y。)到雷達(dá)天線相位
中心(APC)的瞬時距離可表達(dá)為
)=[汙+《-2r。r。cos (《+ "人)]"2《r。 — [x。cos ("人)-凡sin ("人)](1);
其中，(r。， e 。)為該散射中心位置的極坐標(biāo)，ra為恒定項。 對接收數(shù)據(jù)進(jìn)行距離_多普勒(RD)成像，則在某脈沖時刻tm，該散射中心在離散 RD圖像上的位置可表達(dá)如下Z= _乂 =《'+ (x。 cos (w人)-;;。sin( 人))/ ^ (2);化,)=尋4^^ = [x。 sin) +少。cos(w人)]/7。 (3);
乂 欲 」 其中，X。 = _ra/ nr為恒定項，fs為采樣頻率，fr為脈沖重復(fù)頻率，c為電磁波傳播 速度，近似為300000000米/秒，A為波長，M為RD圖像積累的脈沖數(shù)，nr和na分別為 RD圖像的距離向和跨距離向尺度因子(scalingfactor， SF)。 上述表達(dá)式中，脈沖壓縮采用匹配濾波方式實現(xiàn)，此時距離向和跨距離向尺度因 子分別為，
c W o、
7r = 77"77。 _-
2乂 2勘0 當(dāng)采用De-chirp處理進(jìn)行脈沖壓縮時，表達(dá)式(2)和(3)依然成立，此時距離向 和跨距離向尺度因子分別為，仏=左上=^^=^ (5);
r 7V2y 2層 2Mw。
7其中，Y為線性調(diào)頻信號的調(diào)頻率，N為一次脈沖的采樣點數(shù)；B為LFM信號帶寬； Tp為脈沖持續(xù)時間。 將表達(dá)式(2)和(3)寫成矩陣形式，可以得到

X。 少。
其中，S和R(tm)分別為尺度矩陣和旋轉(zhuǎn)矩陣 IK
(6);
s =R(0 =
(8); 對于勻速旋轉(zhuǎn)目標(biāo)，根據(jù)表達(dá)式(6)中的映射關(guān)系，在U和U時刻分別獲得的兩 幅RD圖像存在如下的關(guān)系
m2,其中，e
之間的旋轉(zhuǎn)關(guān)系。
二SR(^)R-'(C
.義w
.化,).
雄
AW
化l)
(9);
'。(tm2-tml)表示兩幅圖像之間的視角差，R(edif)定義了兩幅圖像R(~) = SR(,m2)R-'(d)S-'=
cos f (9,7。
sin《
浙,
(10); 表達(dá)式(9)表明，將tml時刻所稱的圖像經(jīng)過相應(yīng)的旋轉(zhuǎn)可以讓其展現(xiàn)出和tm2所 成圖像相同的目標(biāo)姿態(tài)。對圖像進(jìn)行旋轉(zhuǎn)需要進(jìn)行插值，為了降低圖像旋轉(zhuǎn)所需的計算量， 將上述旋轉(zhuǎn)矩陣分解成如下的形式
1
0
1

、i乂
1
07。
sin(《',)1
1 -衛(wèi)tan 0 1
《
浙
(11); S卩，通過上述的矩陣分解，可以將該傳統(tǒng)圖像旋轉(zhuǎn)過程所需的二維插值以3個連 續(xù)的一維插值實現(xiàn)，且該一維插值相當(dāng)于信號延遲的過程，可以采用如圖6所示的頻域方 法快速實現(xiàn)，在圖6中，對輸入信號進(jìn)行N點的FFT變化，并與頻域因子相乘，進(jìn)而進(jìn)行IFFT 變換，獲取延時輸出。
圖像預(yù)處理 對于所得的兩幅圖像，需要做一些預(yù)處理操作，才可以進(jìn)行目標(biāo)轉(zhuǎn)速的搜索。根據(jù) 算法的特點，該預(yù)處理過程包含以下兩步 (1)、由于兩幅RD圖像形成過程中平動補(bǔ)償效果的不同，目標(biāo)在RD平面上的位置 會有偏移；一般來說，由于該偏移量很小，可以通過二維平移相關(guān)的過程搜索該偏移量，對 齊兩幅圖像。 (2)、由于經(jīng)過RD成像之后，目標(biāo)僅占據(jù)RD圖像中央的很小一部分區(qū)域，因此通過
8對圖像裁剪，去除一些圖像邊緣區(qū)域，降低所需的計算量。
目標(biāo)轉(zhuǎn)速搜索 由于實際目標(biāo)的旋轉(zhuǎn)中心是未知的，因此可將目標(biāo)回波所占距離單元的中心設(shè)為 目標(biāo)的等效旋轉(zhuǎn)中心。同時根據(jù)計算能力，可以設(shè)定一定的轉(zhuǎn)速搜索范圍，根據(jù)轉(zhuǎn)速就可以 得到相應(yīng)的圖像橫向尺度因子和圖像間視角差，實現(xiàn)對其中一幅圖像的旋轉(zhuǎn)，并將旋轉(zhuǎn)圖 像和另外一幅圖像相關(guān)。圖像相關(guān)因子的計算方式如下 其中，fjx， y)和f2(x， y)為兩幅圖像的幅度，Trot ()表示(11)中的旋轉(zhuǎn)變換。
因此，得到的旋轉(zhuǎn)相關(guān)系數(shù)公式(12)為目標(biāo)轉(zhuǎn)速的函數(shù)，當(dāng)該相關(guān)因子取得最大 值時，可以估計得到實際的目標(biāo)轉(zhuǎn)速。 由于目標(biāo)旋轉(zhuǎn)信息是完全未知的，因此對于每一次搜索，都需要采用順時針和逆 時針兩次旋轉(zhuǎn)方式進(jìn)行，分別如圖7和8所示。當(dāng)然，由于兩幅RD圖像已經(jīng)獲得，可以根據(jù) 人工判斷旋轉(zhuǎn)方向，將搜索的計算量降低一倍。 在本案例中，由于順時針旋轉(zhuǎn)搜索的最大相關(guān)系數(shù)較大，因此判斷順時針方向為 正確的旋轉(zhuǎn)方向。如圖7所示，所獲得的目標(biāo)轉(zhuǎn)速搜索結(jié)果大約為0.0103rad./sec.。
圖像橫向定標(biāo) 根據(jù)所估計的目標(biāo)轉(zhuǎn)速，采用極坐標(biāo)格式成像算法(PFA)抑制可能出現(xiàn)的越分辨 單元徙動(MTRC)現(xiàn)象，獲得良好聚焦的、定標(biāo)的I SAR二維圖像輸出。在本例中，圖9是沒 有目標(biāo)轉(zhuǎn)動信息時，采用距離-多普勒成像算法獲得的未定標(biāo)的圖像，顯然無法獲知目標(biāo)
的確切二維形狀信息。圖io是根據(jù)所估計的目標(biāo)轉(zhuǎn)速，對距離-多普勒像橫向定標(biāo)結(jié)
果；圖11是根據(jù)所估計的目標(biāo)轉(zhuǎn)速，采用PFA成像所提供的具有良好聚焦效果且已定標(biāo)的
目標(biāo)ISAR圖像，對比圖9，經(jīng)過定標(biāo)的圖像能較準(zhǔn)確反映目標(biāo)的形狀信息。 最后應(yīng)說明的是，以上實施例僅用以說明本發(fā)明的技術(shù)方案而非對其限制，并且
在應(yīng)用上可以延伸到其他的修改、變化、應(yīng)用和實施例，同時認(rèn)為所有這樣的修改、變化、應(yīng)
用、實施例都在本發(fā)明的精神和范圍內(nèi)。
(12)
權(quán)利要求
一種逆合成孔徑雷達(dá)的目標(biāo)轉(zhuǎn)速估計方法，包括步驟10)、根據(jù)接收的回波數(shù)據(jù)，獲取兩幅逆合成孔徑雷達(dá)的距離多普勒圖像；步驟20)、通過平移相關(guān)來對齊所述兩幅距離多普勒圖像，并獲取所述兩幅距離多普勒圖像中的目標(biāo)區(qū)域圖像；步驟30)、對兩幅所述目標(biāo)區(qū)域距離多普勒圖像進(jìn)行旋轉(zhuǎn)相關(guān)，根據(jù)旋轉(zhuǎn)相關(guān)系數(shù)，估計目標(biāo)的轉(zhuǎn)速。
2. 權(quán)利要求l的方法，其中，步驟IO)中，對所述回波數(shù)據(jù)進(jìn)行滑窗成像，獲取兩幅距離_多普勒圖像序列；其中，使用滑窗寬度確定每幅距離多普勒圖像積累的脈沖數(shù)，滑窗步長確定兩幅圖像之間的視角差。
3. 權(quán)利要求1的方法，其中，步驟20)還包括通過平移相關(guān)方法消除目標(biāo)在距離多普勒平面上的位置偏移，對齊所述兩幅距離多普勒圖像。
4. 權(quán)利要求1的方法，其中，步驟20)中還包括裁剪目標(biāo)所在的區(qū)域的距離多普勒圖像的邊緣區(qū)域。
5. 權(quán)利要求l的方法，其中，步驟30)中，所述兩幅目標(biāo)區(qū)域距離多普勒圖像之間的對應(yīng)關(guān)系為<formula>formula see original document page 2</formula>其中，(Xc(tm》，Yc(U))和(Xc(U，Yc(U)分別表示相同散射中心在兩幅目標(biāo)區(qū)域距離多普勒圖像中的位置，edif = "。(tm2-tj表示所述兩幅圖像之間的視角差，nr和n,分別為距離多普勒圖像的距離向和方位向尺度因子，R( e dif)限定兩幅圖像之間的旋轉(zhuǎn)關(guān)系，其中，<formula>formula see original document page 2</formula>
6.權(quán)利要求5的方法，其中，步驟30)中，所述兩幅目標(biāo)區(qū)域距離多普勒圖像的旋轉(zhuǎn)操作可通過如下三個一維的頻域巻積插值實現(xiàn)，<formula>formula see original document page 2</formula>
7.權(quán)利要求6的方法，其中，步驟30)中，所述兩幅目標(biāo)區(qū)域距離多普勒圖像之間的相關(guān)系數(shù)定義為<formula>formula see original document page 2</formula>其中，fjx， y)和f2(x， y)為兩幅圖像的幅度，Trot()表示所述旋轉(zhuǎn)操作。
8. 權(quán)利要求7的方法，其中，步驟30)中，所述旋轉(zhuǎn)相關(guān)系數(shù)的極大值點對應(yīng)估計的目標(biāo)轉(zhuǎn)速，獲得目標(biāo)轉(zhuǎn)速和旋轉(zhuǎn)方向的估計。
9. 權(quán)利要求1的方法，進(jìn)一步包括步驟40)、根據(jù)估計的目標(biāo)轉(zhuǎn)速對距離-多普勒圖像進(jìn)行橫向定標(biāo)，所述橫向標(biāo)度因子為77。 =^T，其中，入為載波波長，^為脈沖重復(fù)頻率，M為圖像積累的脈沖數(shù)目，"。為目標(biāo)轉(zhuǎn)速。
全文摘要
本發(fā)明涉及雷達(dá)技術(shù)領(lǐng)域，提供一種逆合成孔徑雷達(dá)的目標(biāo)轉(zhuǎn)速估計方法，包括根據(jù)接收的回波數(shù)據(jù)，獲取兩幅逆合成孔徑雷達(dá)的距離多普勒圖像；通過平移相關(guān)對齊所述兩幅距離多普勒圖像，并獲取所述兩幅距離多普勒圖像中的目標(biāo)區(qū)域圖像；對兩幅所述目標(biāo)區(qū)域距離多普勒圖像進(jìn)行旋轉(zhuǎn)相關(guān)，根據(jù)旋轉(zhuǎn)相關(guān)系數(shù)，估計目標(biāo)的轉(zhuǎn)速。本發(fā)明通過三級一維卷積插值過程實現(xiàn)圖像旋轉(zhuǎn)，可以減少計算量，提高轉(zhuǎn)速估計的速度，同時，該方法對目標(biāo)散射特性沒有具體的要求，應(yīng)用范圍較大。
文檔編號G01S13/58GK101738615SQ20081022660
公開日2010年6月16日 申請日期2008年11月17日 優(yōu)先權(quán)日2008年11月17日
發(fā)明者葉春茂, 左渝, 彭應(yīng)寧, 王秀壇, 許稼 申請人:清華大學(xué)