[發明專利]一種基于廣義S變換的微弱目標檢測方法

無效

• 申請號: 200810147693.1
• 申請日: 2008-11-27
• 公開（公告）號: CN101493934A
• 公開（公告）日: 2009-07-29
• 發明（設計）人: 彭真明;雍楊;張杰;甄莉
• 申請（專利權）人: 電子科技大學
• 主分類號: G06T7/00
• 分類號: G06T7/00
• 代理公司: 暫無信息
• 代理人: 暫無信息

1.一種基于廣義S變換的微弱目標檢測方法，包括如下步驟：

步驟1：對一幅大小為M×N的原始圖像I(x′，y′)，進行二維廣義S變換，M為圖像寬度，N為圖像高度，(x′，y′)表示原始圖像I(x′，y′)的空間坐標，x′＝1、2、…、M，y′＝1、2、…、N，具體包括以下幾個步驟：

步驟1-1：對原始圖像I(x′，y′)進行快速傅立葉變換：I(x′，y′)→F(α，β)，F(α，β)為圖像的快速傅立葉變換結果，α，β為頻率域變量；

步驟1-2：對所有空間頻率點(k_x，k_y)，k_x＝1、2、…、M，k_x表示x方向上的空間頻率點，k_y＝1、2、…、N，k_y表示y方向上的空間頻率點，進行步驟1-3-步驟1-6的操作：

步驟1-3：在空間頻率點(k_x，k_y)，對高斯局部化窗函數進行傅立葉變換，得到它的傅立葉頻譜W(α，β)，高斯局部化窗函數中的μ和η為取值范圍為正數的調節參數，其定義為σ_x＝μ/k_x，σ_y＝η/k_y，σ_x、σ_y是尺度因子；

步驟1-4：移動傅立葉頻譜：F(α，β)→F(α+k_x，β+k_y)；

步驟1-5：計算F(α+k_x，β+k_y)與W(α，β)的點積，表示為

步驟1-6：計算的傅立葉反變換得到原始圖像I(x′，y′)的廣義S變換結果S(x，y，k_x，k_y)；

步驟2：根據步驟1所得的原始圖像I(x′，y′)的廣義S變換結果S(x，y，k_x，k_y)對圖像中的微弱目標進行檢測，固定表示y方向上的空間頻率點k_y的值，使廣義S變換結果S(x，y，k_x，k_y)從四維降為三維數據，計算沿y方向每一個取值下所對應的垂直切面圖像(x，k_x)的頻率成分能量之和P_x(i)的平均值作為閾值T₁，其中i＝1、2、…、M；比較P_x(i)與T₁的大小，其中P_x(i)大于T₁的行就是微弱目標所在行的位置；

步驟3：根據步驟1所得的原始圖像I(x′，y′)的廣義S變換結果S(x，y，k_x，k_y)對圖像中的微弱目標進行檢測，固定表示x方向上的空間頻率點k_x的值，使廣義S變換結果S(x，y，k_x，k_y)從四維降為另一組三維數據，計算沿x方向每一個取值下所對應的垂直切面圖像(y，k_y)的頻率成分能量之和P_y(j)的平均值作為閾值T₂，其中j＝1、2、…、N；比較P_y(j)與T₂的大小，其中P_y(j)大于T₂的列就是微弱目標所在列的位置；

步驟4：因S變換用一個窗函數與原始圖像I(x′，y′)相乘，相當于截取了原始圖像I(x′，y′)在(x，y)位置的局部信息，變換后的(x，y)位置能反映原始圖像I(x′，y′)的空間位置信息；所以步驟2、3所確定的微弱目標所在行和列的位置就是該目標在原始圖像I(x′，y′)中的行列位置。