1.一種雙站線陣三維成像合成孔徑雷達快速成像方法，它包括以下步驟：?

步驟一、初始化雙站線陣三維成像合成孔徑雷達成像系統參數：?

初始化的成像系統參數包括：發射平臺與接收平臺速度矢量?平臺初始位置矢量?雷達發射電磁波的波數，記做K_c，接收機線陣天線各陣元相對平臺中心位置，記做?其中i為天線各陣元序號，為自然數，i＝0，1...，M，M為接收機線陣天線各陣元總數，雷達發射基帶信號的信號帶寬，記做B，雷達發射信號脈沖寬度，記做T_p，雷達接收波門持續寬度，記做T_d，雷達接收系統的采樣頻率，記做f_s，雷達系統的脈沖重復頻率，記做PRF，及接收系統接收波門相對于發射信號發射波門的延遲，記做T_D；上述參數均為雙站線陣三維成像合成孔徑雷達系統的標準參數，其中，雷達發射電磁波的波數K_c，雷達發射基帶信號的信號帶寬B，雷達發射信號脈沖寬度T_p，雷達接收波門持續寬度T_d，雷達接收系統的采樣頻率f_s，雷達系統的脈沖重復頻率PRF及接收系統接收波門相對于發射信號發射波門的延遲在雙站線陣三維成像合成孔徑雷達設計過程中已經確定；其中，發、收平臺速度矢量?及初始位置矢量?在雙站線陣三維成像合成孔徑雷達觀測方案設計中已經確定；?

步驟二、雙站線陣三維成像合成孔徑雷達原始數據進行距離壓縮：?

采用合成孔徑雷達標準距離壓縮方法對合成孔徑雷達距離向回波數據?進行壓縮，得到距離壓縮后的雙站線陣三維成像合成孔徑雷達數據，記做?

步驟三、獲得雙站線陣三維成像合成孔徑雷達指標集內連續小區域的高度和相應的散射系數最大值，即初始化的成像區域；?

采用雙站線陣三維成像合成孔徑雷達三維反向投影成像方法，對雙站線陣三維成像合成孔徑雷達圖像空間中區域為7×7×H的區域Ω₀＝{(x_i，y_j，h)|i＝0...6，j＝0...6，h∈[-H/2，H/2]}內所有點進行成像處理，H為成像場景的參考高度；得到雙站線陣三維成像合成孔徑雷達圖像空間中區域為Ω₀內的所有點的后向散射系數，也即雙站線陣三維成像合成孔徑雷達指標集為Ξ₀＝{z|z＝h(x_i，y_j)，i＝0...6，j＝0...6}內的所有點沿高度向的后向散射系數的分布函數，記做：σ(x，y，h)；?

對指標集Ξ₀所對應的后向散射系數分布函數沿高度向采用遍歷法進行搜索得到分布函數沿高度向的散射系數最大值及與之對應的最大值處的高度；其中，雙站合成孔徑雷達沿高度向的后向散射系數分布函數σ(x，y，h)最大值用σ_max(x，y)表示，而相應的高度用h(x，y)表示；?

步驟四、沿x向對高度進行擬合預測；?

采用標準的多項式擬合方法對已知的y＝0行的高度數據h(x_i，0)，i＝0…6進行多項式擬合，此處我們選擇擬合的階數為3，即進行三階曲面擬合及預測，得到擬合函數h＝f(x)，把x＝7帶入擬合函數h＝f(x)，得到預測高度h_p(7，0)；?

步驟五、獲取檢測門限；?

根據合成孔徑雷達的后向散射系數的分布規律，選取瑞利——貝塞爾分布；設定虛警概率α或者檢測概率δ，根據虛警概率分布函數或檢測概率分布函數及門限檢測技術計算出檢測門限θ；這里我們已知檢測概率的分布函數即瑞利——貝塞爾函數，通過計算該函數的均值和方差來獲取檢測門限；?

根據步驟三得到的后向散射系數最大值σ_max(x，y)，選取y＝0行的σ_max(x_i，0)，i＝0…6，分別計算這一組值的中值和方差σ_med、s_σ，利用式子?得到均值幅度β，按照3dB帶寬檢測準則，?即為分布函數的概率均值、s_σ為其方差；?

設定檢測概率為0.8，利用門限檢測技術，我們得到檢測門限θ；?

步驟六、搜索雙站線陣三維成像合成孔徑雷達圖像點的真實高度；?

采用雙站線陣三維成像合成孔徑雷達三維反向投影成像方法，對三維成像空間中的點(7，0，h_p)進行成像，得到點(7，0，h_p)的后向散射系數σ(7，0，h_p)；?

當σ(7，0，h_p)＜θ，執行上行操作，上行操作定義為采用雙站線陣三維成像合成孔徑雷達三維反向投影成像方法，在點(x，y，h)已經成像的情況下，對位于成像空間中的點(x，y，h+kΔh)，k≥1進行成像；因此先對三維成像空間中的點(7，0，h_p+Δh)進行成像，得到該點的后向散射系數σ(7，0，h_p+Δh)，其中Δh為雙站線陣三維成像合成孔徑雷達沿?高度向的分辨率；判斷σ(7，0，h_p+Δh)是否大于門限θ，如果σ(7，0，h_p+Δh)＞θ，采用雙站線陣三維成像合成孔徑雷達三維反向投影成像方法，對三維成像空間中的點(7，0，h_p+2Δh)進行成像，得到該點的后向散射系數σ(7，0，h_p+2Δh)，當σ(7，0，h_p+2Δh)＜σ(7，0，h_p+Δh)，跳出本次搜索，且取σ_max(7，0)＝σ(7，0，h_p+Δh)，h(7，0)＝h+Δh，當σ(7，0，h_p+2Δh)＞σ(7，0，h_p+Δh)，對點(7，0，h_p+3Δh)進行成像，得到該點的后向散射系數σ(7，0，h_p+3Δh)，當σ(7，0，h_p+3Δh)＜σ(7，0，h_p+2Δh)，跳出本次搜索，且取σ_max(7，0)＝σ(7，0，h_p+2Δh)，h(7，0)＝h+2Δh，當σ(7，0，h_p+2Δh)＜σ(7，0，h_p+Δh)，繼續本次的上行操作，直至得到滿足條件σ(7，0，h_p+(k-1)Δh)＜σ(7，0，h_p+kΔh)＞σ(7，0，h_p+(k+1)Δh)的k值，此時令σ_max(7，0)＝σ(7，0，h_p+kΔh)，h(7，0)＝h+kΔh；當σ(7，0，h_p+Δh)＜θ，進行下行操作，下行操作定義為采用雙站線陣三維成像合成孔徑雷達三維反向投影成像方法，在點(x，y，h)已經成像的情況下，對位于成像空間中的點(x，y，h+kΔh)，k≤-1進行成像；因此對點(7，0，h_p-Δh)進行成像，得到該點的后向散射系數σ(7，0，h_p-Δh)，如果σ(7，0，h_p-Δh)＞θ，對點(7，0，h_p-2Δh)進行成像，得到該點的后向散射系數σ(7，0，h_p-2Δh)，如果σ(7，0，h_p-2Δh)＜σ(7，0，h_p-Δh)，跳出本次搜索，且令σ_max(7，0)＝σ(7，0，h_p-Δh)，h(7，0)＝h-Δh，當σ(7，0，h_p-2Δh)＞σ(7，0，h_p-Δh)，對點(7，0，h_p-3Δh)進行成像，得到該點的后向散射系數σ(7，0，h_p-3Δh)，當σ(7，0，h_p-3Δh)＜σ(7，0，h_p-2Δh)，跳出本次搜索，且取σ_max(7，0)＝σ(7，0，h_p-2Δh)，h(7，0)＝h-2Δh，當σ(7，0，h_p-3Δh)＞σ(7，0，h_p-2Δh)繼續本次的下行操作，直至得到滿足條件σ(7，0，h_p-(k-1)Δh)＜σ(7，0，h_p-kΔh)＞σ(7，0，h_p-(k+1)Δh)的k值，此時令σ_max(7，0)＝σ(7，0，h_p-kΔh)，h(7，0)＝h-kΔh??；其中k＝2…ΔhH/2；當σ(7，0，h_p-Δh)＜θ，轉入k＝2處的上行操作，當σ(7，0，h_p+2Δh)＜θ，轉入k＝2處的下行操作，這樣上下交替搜索，直至搜索到滿足條件的后向散射系數及對應的高度值；?當沒有搜索到滿足條件的后向散射系數，取后向散射系數最大處的高度值作為成像的最終高度；?

當σ(7，0，h_p)＞θ，采用雙站線陣三維成像合成孔徑雷達三維反向投影成像方法，對點(7，0，h_p+Δh)進行成像，得到該點的后向散射系數σ(7，0，h_p+Δh)，對三維成像空間中的點(7，0，h_p-Δh)進行成像，得到該點的后向散射系數σ(7，0，h_p-Δh)，Δh意義同上；當σ(7，0，h_p-Δh)＜σ(7，0，h_p)＞σ(7，0，h_p+Δh)，跳出本次搜索，且取σ_max(7，0)＝σ(7，0，h_p)，h(7，0)＝h；當σ(7，0，h_p-Δh)＜σ(7，0，h_p)＜σ(7，0，h_p+Δh)，進入上面敘述的上行操作，直至獲得滿足條件的后向散射系數及對應的高度值；當σ(7，0，h_p-Δh)＞σ(7，0，h_p)＞σ(7，0，h_p+Δh)，進入上面敘述的下行操作，直至獲得滿足條件的后向散射系數及對應的高度值；當σ(7，0，h_p-Δh)＞σ(7，0，h_p)＜σ(7，0，h_p+Δh)，進入全高度搜索，并取最大散射系數及其對應的高度作為滿足條件的散射系數及高度；當在執行上行操作或下行操作時沒有跳出搜索，則同樣取最大散射系數及其對應的高度作為滿足條件的散射系數及高度；?

在經過上述的操作步驟后，我們得到用曲線擬合及預測方法的最大后向散射系數σ_max(7，0)及對應的高度h(7，0)；?

步驟七、基于滑窗的整行預測；?

取高度值的窗長為7，即沿x向滑動得到一組新的高度值h(x_i，0)，i＝1…7，重復運用步驟五與步驟六相同的預測方法，得到滿足條件的最大后向散射系數σ_max(8，0)及對應的高度h(8，0)，這樣把窗位置依次進行滑動，最后得到y＝0行的所有滿足條件的散射系數和對應的高度；?

步驟八、獲取全場景預測的初始數據；?

按照步驟六和步驟七的方法，獲取y＝1…6行的所有滿足條件的散射系數和對應的高度；?

步驟九、進行全場景滑窗預測；?

同樣取高度值的窗長為7，此時取沿y向滑動得到的高度值h(x_i，y_j)，其中?i＝0…128，j＝0…128；重復運用步驟五與步驟六相同的預測方法，得到所有滿足條件的最大后向散射系數σ_max(x_i，y_j)及對應的高度h(x_i，y_j)；?

步驟十、獲取雙站線陣三維成像合成孔徑雷達最終成像；?

為了減小預測成像的誤差，對經過上述步驟得到的高度值h(x_i，y_j)組成的矩陣H進行窗為7×7的二維中值濾波，二維中值濾波在MATLAB7.0中的函數為medfilt2()，得到雙站線陣三維成像合成孔徑雷達最終成像的高度矩陣結果H及其對應的后向散射系數矩陣σ；?

經過上述操作，就可以得到雙站線陣三維成像合成孔徑雷達快速成像。?