CN107064895A - 星载合成孔径雷达图像的模糊抑制方法和装置 - Google Patents
星载合成孔径雷达图像的模糊抑制方法和装置 Download PDFInfo
- Publication number
- CN107064895A CN107064895A CN201710134700.3A CN201710134700A CN107064895A CN 107064895 A CN107064895 A CN 107064895A CN 201710134700 A CN201710134700 A CN 201710134700A CN 107064895 A CN107064895 A CN 107064895A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mfrac
- image
- eta
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/9021—SAR image post-processing techniques
- G01S13/9027—Pattern recognition for feature extraction
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Radar Systems Or Details Thereof (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
本公开涉及一种星载合成孔径雷达图像的模糊抑制方法和装置。所述方法包括:对星载合成孔径雷达图像进行方位模糊重建,确定出所述星载合成孔径雷达图像的模糊图像;根据与所述星载合成孔径雷达图像中的数据点对应的所述模糊图像的数据点,判断所述星载合成孔径雷达图像中的数据点是否满足预定抑制条件;当所述星载合成孔径雷达图像中的数据点满足所述预定抑制条件时,利用所述模糊图像,对所述星载合成孔径雷达图像中、满足所述预定抑制条件的数据点进行模糊抑制。这样,能够有效去除星载合成孔径雷达图像中的方位模糊,提高星载合成孔径雷达图像的清晰度。
Description
技术领域
本公开涉及信号处理技术领域,具体地,涉及一种星载合成孔径雷达图像的模糊抑制方法和装置。
背景技术
星载合成孔径雷达是一种可以对地球进行观测的高分辨率遥感探测雷达。由于星载合成孔径雷达卫星可以不受雨雾、暗夜等能见度低的气象条件的限制,因此,可以通过它实现全天时、全天候、高分辨率、宽幅对地观测。正因为这样,星仔合成孔径雷达在农业、海洋、灾害监测、3D绘图等领域具有广阔的应用前景。
星载合成孔径雷达具有聚束模式、条带模式和滑动聚束模式三种工作模式,其中,条带模式是它最主要的工作模式。然而,星载合成孔径雷达本身的工作机理决定了雷达图像中可能存在方位模糊的现象,这将会明显降低星载合成孔径雷达图像的质量,后续的与星载合成孔径雷达图像相关的应用也会受到一定程度的影响。
相关技术中,通常采用理想滤波器法和维纳滤波器法等方法对星载合成孔径雷达图像进行模糊抑制,但这些方法均难以将模糊信号抑制干净,甚至会损失图像分辨率和信噪比。
发明内容
本公开的目的是针对现有技术中的方位模糊信号抑制不彻底的问题,提供一种星载合成孔径雷达图像的模糊抑制方法和装置。
为了实现上述目的,本公开提供一种星载合成孔径雷达图像的模糊抑制方法,所述方法包括:
对星载合成孔径雷达图像进行方位模糊重建,确定出所述星载合成孔径雷达图像的模糊图像;
根据与所述星载合成孔径雷达图像中的数据点对应的所述模糊图像的数据点,判断所述星载合成孔径雷达图像中、满足所述预定抑制条件的数据点是否满足预定抑制条件;
当所述星载合成孔径雷达图像中的数据点满足所述预定抑制条件时,利用所述模糊图像,对所述星载合成孔径雷达图像中的数据点进行模糊抑制。
可选地,所述模糊图像包括负频率方位模糊图像;
所述对星载合成孔径雷达图像进行方位模糊重建,确定出所述星载合成孔径雷达图像的模糊图像,包括:
对所述星载合成孔径雷达图像s0(τ,η)按列做快速傅里叶变换,得到图像S1(τ,fη),其中,j为列索引,Nr为图像矩阵总列数,fs为信号采样率,i为行索引,Na为图像矩阵总行数,fp为脉冲重复频率,τ为距离向时间,η为方位向时间,fη为方位向频率;
对所述图像S1(τ,fη)按列进行方位相位补偿和天线方向图补偿,每一列数据乘以传递函数H1,得到图像S2(τ,fη):
S2(τ,fη)=H1·S1(τ,fη)
其中,H1为传递函数,f0为信号载频,R0为斜距,c为光速,Rref为场景中心斜距,λ为波长,Vr为等效速度,fηc为多普勒中心频率,D(fη-fp)为徙动因子,Wa(fη-fηc-fp)为卫星方位向天线方向图的频域表达式,wa(τ)为方向图的时域表达式;
对所述图像S2(τ,fη)按行做快速傅里叶变换,得到图像S3(fτ,fη);
对所述图像S3(fτ,fη)按行进行一致模糊距离徙动校正和方位向调制,得到图像S4(fτ,fη);
对所述图像S4(fτ,fη)按行做快速逆傅里叶变换,得到图像S5(τ,fη);
对所述图像S5(τ,fη)按行进行补余模糊距离徙动校正,得到图像S6(τ,fη);
对所述图像S6(τ,fη)按行做快速傅里叶变换,得到图像S7(fτ,fη);
对所述图像S7(fτ,fη)按行进行脉冲压缩和二次距离压缩误差校正,得到图像S8(fτ,fη);
对所述图像S8(fτ,fη)按行做快速逆傅里叶变换,得到图像S9(τ,fη);
对所述图像S9(τ,fη)补偿残余误差和进行方位散焦处理,得到图像S10(τ,fη);
对所述图像S10(τ,fη)按列做快速逆傅里叶变换,得到所述负频率方位模糊图像。
可选地,所述模糊图像包括正频率方位模糊图像;
所述对星载合成孔径雷达图像进行方位模糊重建,确定出所述星载合成孔径雷达图像的模糊图像,包括:
对所述星载合成孔径雷达图像s0(τ,η)按列做快速傅里叶变换,得到图像S11(τ,fη),其中,j为列索引,Nr为图像矩阵总列数,fs为信号采样率,i为行索引,Na为图像矩阵总行数,fp为脉冲重复频率,τ为距离向时间η为方位向时间,fη为方位向频率;
对所述图像S11(τ,fη)按列进行方位相位补偿和天线方向图补偿,得到图像S12(τ,fη):
S12(τ,fη)=H1·S11(τ,fη)
其中,H1为传递函数,f0为信号载频,R0为斜距,c为光速,Rref为场景中心斜距,λ为波长,Vr为等效速度,fηc为多普勒中心频率,D(fη+fp)为徙动因子,Wa(fη-fηc+fp)为卫星方位向天线方向图的频域表达式,wa(τ)为方向图的时域表达式;
对所述图像S12(τ,fη)按行做快速傅里叶变换,得到图像S13(fτ,fη);
对所述图像S13(fτ,fη)按行进行一致模糊距离徙动校正和方位向调制,得到图像S14(fτ,fη);
对所述图像S14(fτ,fη)按行做快速逆傅里叶变换,得到图像S15(τ,fη);
对所述图像S15(τ,fη)按行进行补余模糊距离徙动校正,得到图像S16(τ,fη);
对所述图像S16(τ,fη)按行做快速傅里叶变换,得到图像S17(fτ,fη);
对所述图像S17(fτ,fη)按行进行脉冲压缩和二次距离压缩误差校正,得到图像S18(fτ,fη);
对所述图像S18(fτ,fη)按行做快速逆傅里叶变换,得到图像S19(τ,fη);
对所述图像S19(τ,fη)补偿残余误差和进行方位散焦处理,得到图像S20(τ,fη);
对所述图像S20(τ,fη)按列做快速逆傅里叶变换,得到所述正频率方位模糊图像。
可选地,所述根据与所述星载合成孔径雷达图像中的数据点对应的所述模糊图像的数据点,判断所述星载合成孔径雷达图像中的数据点是否满足预定抑制条件,包括:
判断所述模糊图像的数据点的绝对值是否大于所述星载合成孔径雷达图像中对应的数据点的绝对值的预定倍数,其中,所述预定倍数大于0且小于1;
判断所述方位向时间是否满足预设的阈值条件;
当所述模糊图像的数据点的绝对值大于所述星载合成孔径雷达图像中对应的数据点的绝对值的预定倍数、且所述方位向时间满足所述阈值条件时,确定所述星载合成孔径雷达图像中的数据点满足所述预定抑制条件。
可选地,所述模糊图像包括负频率方位模糊图像;
所述阈值条件包括:
其中,η为所述方位向时间;fp为脉冲重复频率;fr为多普勒调频率;和/或
所述模糊图像包括正频率方位模糊图像;
所述阈值条件包括:
其中,η为所述方位向时间;Na为图像矩阵总行数;fp为脉冲重复频率;fr为多普勒调频率。
可选地,所述当所述星载合成孔径雷达图像中的数据点满足所述预定抑制条件时,利用所述模糊图像,对所述星载合成孔径雷达图像中的数据点进行模糊抑制,包括:
将所述星载合成孔径雷达图像中、满足所述预定抑制条件的数据点的绝对值,与所述模糊图像中对应的数据点的绝对值做差值,得到经模糊抑制后的图像。
本公开还提供一种星载合成孔径雷达图像的模糊抑制装置,所述装置包括:
方位模糊重建模块,用于对星载合成孔径雷达图像进行方位模糊重建,确定出所述星载合成孔径雷达图像的模糊图像;
判断模块,用于根据与所述星载合成孔径雷达图像中的数据点对应的所述模糊图像的数据点,判断所述星载合成孔径雷达图像中的数据点是否满足预定抑制条件;
模糊抑制模块,用于当所述星载合成孔径雷达图像中的数据点满足所述预定抑制条件时,利用所述模糊图像,对所述星载合成孔径雷达图像中、满足所述预定抑制条件的数据点进行模糊抑制。
可选地,所述模糊图像包括负频率方位模糊图像;
所述方位模糊重建模块用于:
对所述星载合成孔径雷达图像s0(τ,η)按列做快速傅里叶变换,得到图像S1(τ,fη),其中,j为列索引,Nr为图像矩阵总列数,fs为信号采样率,i为行索引,Na为图像矩阵总行数,fp为脉冲重复频率,τ为距离向时间,η为方位向时间,fη为方位向频率;
对所述图像S1(τ,fη)按列进行方位相位补偿和天线方向图补偿,得到图像S2(τ,fη):
S2(τ,fη)=H1·S1(τ,fη)
其中,H1为传递函数,f0为信号载频,R0为斜距,c为光速,Rref为场景中心斜距,λ为波长,Vr为等效速度,fηc为多普勒中心频率,D(fη-fp)为徙动因子,Wa(fη-fηc-fp)为卫星方位向天线方向图的频域表达式,wa(τ)为方向图的时域表达式;
对所述图像S2(τ,fη)按行做快速傅里叶变换,得到图像S3(fτ,fη);
对所述图像S3(fτ,fη)按行进行一致模糊距离徙动校正和方位向调制,得到图像S4(fτ,fη);
对所述图像S4(fτ,fη)按行做快速逆傅里叶变换,得到图像S5(τ,fη);
对所述图像S5(τ,fη)按行进行补余模糊距离徙动校正,得到图像S6(τ,fη);
对所述图像S6(τ,fη)按行做快速傅里叶变换,得到图像S7(fτ,fη);
对所述图像S7(fτ,fη)按行进行脉冲压缩和二次距离压缩误差校正,得到图像S8(fτ,fη);
对所述图像S8(fτ,fη)按行做快速逆傅里叶变换,得到图像S9(τ,fη);
对所述图像S9(τ,fη)补偿残余误差和进行方位散焦处理,得到图像S10(τ,fη);
对所述图像S10(τ,fη)按列做快速逆傅里叶变换,得到所述负频率方位模糊图像。
可选地,所述模糊图像包括正频率方位模糊图像;
所述方位模糊重建模块用于:
对所述星载合成孔径雷达图像s0(τ,η)按列做快速傅里叶变换,得到图像S11(τ,fη),其中,j为列索引,Nr为图像矩阵总列数,fs为信号采样率,i为行索引,Na为图像矩阵总行数,fp为脉冲重复频率,τ为距离向时间,η为方位向时间,fη为方位向频率;
对所述图像S11(τ,fη)按列进行方位相位补偿和天线方向图补偿,得到图像S12(τ,fη):
S12(τ,fη)=H1·S11(τ,fη)
其中,H1为传递函数,f0为信号载频,R0为斜距,c为光速,Rref为场景中心斜距,λ为波长,Vr为等效速度,fηc为多普勒中心频率,D(fη+fp)为徙动因子,Wa(fη-fηc+fp)为卫星方位向天线方向图的频域表达式,wa(τ)为方向图的时域表达式;
对所述图像S12(τ,fη)按行做快速傅里叶变换,得到图像S13(fτ,fη);
对所述图像S13(fτ,fη)按行进行一致模糊距离徙动校正和方位向调制,得到图像S14(fτ,fη);
对所述图像S14(fτ,fη)按行做快速逆傅里叶变换,得到图像S15(τ,fη);
对所述图像S15(τ,fη)按行进行补余模糊距离徙动校正,得到图像S16(τ,fη);
对所述图像S16(τ,fη)按行做快速傅里叶变换,得到图像S17(fτ,fη);
对所述图像S17(fτ,fη)按行进行脉冲压缩和二次距离压缩误差校正,得到图像S18(fτ,fη);
对所述图像S18(fτ,fη)按行做快速逆傅里叶变换,得到图像S19(τ,fη);
对所述图像S19(τ,fη)补偿残余误差和进行方位散焦处理,得到图像S20(τ,fη);
对所述图像S20(τ,fη)按列做快速逆傅里叶变换,得到所述正频率方位模糊图像。
可选地,所述判断模块包括:
第一判断子模块,用于判断所述模糊图像的数据点的绝对值是否大于所述星载合成孔径雷达图像中对应的数据点的绝对值的预定倍数,其中,所述预定倍数大于0且小于1;
第二判断子模块,用于判断所述方位向时间是否满足预设的阈值条件;
确定子模块,用于当所述模糊图像的数据点的绝对值大于所述星载合成孔径雷达图像中对应的数据点的绝对值的预定倍数、且所述方位向时间满足所述阈值条件时,确定所述星载合成孔径雷达图像中的数据点满足所述预定抑制条件。
可选地,所述模糊图像包括负频率方位模糊图像;
所述阈值条件包括:
其中,η为所述方位向时间;fp为脉冲重复频率;fr为多普勒调频率;和/或
所述模糊图像包括正频率方位模糊图像;
所述阈值条件包括:
其中,η为所述方位向时间;Na为图像矩阵总行数;fp为脉冲重复频率;fr为方位向时间。
可选地,所述模糊抑制模块用于将所述星载合成孔径雷达图像中、满足所述预定抑制条件的数据点的绝对值,与所述模糊图像中对应的数据点的绝对值做差值,得到经模糊抑制后的图像。
通过上述技术方案,对星载合成孔径雷达图像进行方位模糊重建,得到模糊图像,并利用模糊图像,对星载合成孔径雷达图像中、满足预定抑制条件的数据点进行模糊抑制。这样,能够有效去除星载合成孔径雷达图像中的方位模糊,提高星载合成孔径雷达图像的清晰度。
本公开的其他特征和优点将在随后的具体实施方式部分予以详细说明。
附图说明
附图是用来提供对本公开的进一步理解,并且构成说明书的一部分,与下面的具体实施方式一起用于解释本公开,但并不构成对本公开的限制。在附图中:
图1是根据一示例性实施例示出的一种星载合成孔径雷达图像的模糊抑制方法的流程图。
图2是根据另一示例性实施例示出的一种星载合成孔径雷达图像的模糊抑制方法的流程图。
图3A是根据一示例性实施例示出的一种星载合成孔径雷达图像的模糊抑制装置的框图。
图3B是根据另一示例性实施例示出的一种星载合成孔径雷达图像的模糊抑制装置中的判断模块的框图。
具体实施方式
以下结合附图对本公开的具体实施方式进行详细说明。应当理解的是,此处所描述的具体实施方式仅用于说明和解释本公开,并不用于限制本公开。
请参考图1,图1是一示例性实施例提供的星载合成孔径雷达图像的模糊抑制方法的流程图,如图1所示,该方法包括S11、S12以及S13。下面依次对各个步骤进行说明。
在S11中,对星载合成孔径雷达图像进行方位模糊重建,确定出星载合成孔径雷达图像的模糊图像。
本公开中,对星载合成孔径雷达图像进行方位模糊重建,重建模糊的幅度响应和相位历史,得到与星载合成孔径雷达图像对应的模糊图像。
此外,模糊图像包括负频率方位模糊图像和/或正频率方位模糊图像。
首先,对确定负频率方位模糊图像的方法进行说明:
对星载合成孔径雷达图像s0(τ,η)按列做快速傅里叶变换,得到图像S1(τ,fη):
S1(τ,fη)=FFTa[s0(τ,η)]
其中,FFTa为按列做快速傅里叶变换,j为列索引,Nr为图像矩阵总列数,fs为信号采样率,i为行索引,Na为图像矩阵总行数,fp为脉冲重复频率τ为距离向时间,η为方位向时间,fη为方位向频率;
对图像S1(τ,fη)按列进行方位相位补偿和天线方向图补偿,每一列数据乘以传递函数H1,得到图像S2(τ,fη):
S2(τ,fη)=H1·S1(τ,fη)
其中,f0为信号载频,R0为斜距,c为光速,Rref为场景中心斜距,λ为波长,Vr为等效速度,fηc为多普勒中心频率,D(fη-fp)为徙动因子,Wa(fη-fηc-fp)为卫星方位向天线方向图的频域表达式,wa(τ)为方向图的时域表达式;
对图像S2(τ,fη)按行做快速傅里叶变换,得到图像S3(fτ,fη):
S3(fτ,fη)=FFTr[S2(τ,fη)]
其中,FFTr为按行做快速傅里叶变换;
对图像S3(fτ,fη)按行进行一致模糊距离徙动校正和方位向调制,每一行数据乘以传递函数H2,得到图像S4(fτ,fη):
其中,Q1(fη)为第一散焦因子,Kr为信号调频率,D(fηc)为徙动因子,fτ为距离向频率;
对图像S4(fτ,fη)按行做快速逆傅里叶变换,得到图像S5(τ,fη):
S5(τ,fη)=IFFTr[S4(fτ,fη)]
其中,IFFTr为按行做快速逆傅里叶变换;
对图像S5(τ,fη)按行进行补余模糊距离徙动校正,每一行数据乘以校正函数H3,得到图像S6(τ,fη):
其中,Knew为信号调频率,α为线性变标因子;
对图像S6(τ,fη)按行做快速傅里叶变换,得到图像S7(fτ,fη):
S7(fτ,fη)=FFTr[S6(τ,fη)]
对图像S7(fτ,fη)按行进行脉冲压缩和二次距离压缩误差校正,每一行数据乘以校正函数H4,得到图像S8(fτ,fη):
其中,Q2(fη)为第二散焦因子,D(fη)为徙动因子,Km(fη)为距离向调频率,Km(fη-fp)为距离向调频率;
对图像S8(fτ,fη)按行做快速逆傅里叶变换,得到图像S9(τ,fη):
S9(τ,fη)=IFFTr[S8(fτ,fη)]
对图像S9(τ,fη)补偿残余误差和进行方位散焦处理,矩阵中每一个点均需要乘以对应的补偿函数H5,得到图像S10(τ,fη):
对图像S10(τ,fη)按列做快速逆傅里叶变换,得到负频率方位模糊图像:
s-1(τ,η)=IFFTa[S10(τ,fη)]
其中,s-1(τ,η)为所述负频率方位模糊图像,IFFTa为按列做快速逆傅里叶变换。
其次,对确定正频率方位模糊图像的方法进行说明:
正频率方位模糊图像也可以按照上述步骤进行,只需要将上述公式中所有的-fp替换成+fp,即可得到方位模糊重建之后的正频率方位模糊图像S+1(τ,η),本实施例中不做赘述。
在S12中,根据与星载合成孔径雷达图像中的数据点对应的模糊图像的数据点,判断星载合成孔径雷达图像中的数据点是否满足预定抑制条件。
在本公开中,根据与星载合成孔径雷达图像的数据点对应的模糊图像的数据点,判断星载合成孔径雷达图像的数据点是否满足预定抑制条件,也就是判断该数据点是否为模糊点,如果该数据点是模糊点,则对该数据点进行模糊抑制。
图2是根据另一示例性实施例示出的一种星载合成孔径雷达图像的模糊抑制方法的流程图。如图2所示,上述S12可以包括以下步骤。
在S121中,判断模糊图像的数据点的绝对值是否大于星载合成孔径雷达图像中对应的数据点的绝对值的预定倍数。
在本公开中,该预定倍数大于0且小于1,它可以例如是0.2。
示例地,当模糊图像为负频率方位模糊图像时,当模糊图像满足以下等式(1)时,可以判定模糊图像的数据点的绝对值大于星载合成孔径雷达图像中对应的数据点的绝对值的预定倍数:
|s-1(τ,η)|>m×|s0(τ,η)| (1)
其中,s-1(τ,η)为所述负频率方位模糊图像,s0(τ,η)为星载合成孔径雷达图像,m为所述预定倍数。
又示例地,当模糊图像为正频率方位模糊图像时,当模糊图像满足以下等式(2)时可以判定模糊图像的数据点的绝对值大于星载合成孔径雷达图像中对应的数据点的绝对值的预定倍数:
|s+1(τ,η)|>m×|s0(τ,η)| (2)
其中,s+1(τ,η)为所述正频率方位模糊图像,s0(τ,η)为星载合成孔径雷达图像,m为所述预定倍数。
在S122中,判断方位向时间是否满足预设的阈值条件。
在一种实施方式中,当模糊图像为负频率方位模糊图像时,该阈值条件可以是其中,η为所述方位向时间,fp为脉冲重复频率,fr为多普勒调频率。
在另一种实施方式中,当模糊图像为正频率方位模糊图像时,该阈值条件可以是其中,η为所述方位向时间,Na为图像矩阵总行数,fp为脉冲重复频率,fr为多普勒调频率。
另外,需要说明的是,上述S122可以在上述S121之前执行,也可以与S121同时执行,还可以在S121之后执行,S121和S122之间的执行顺序在此不作具体限定。
在S123中,当模糊图像的数据点的绝对值大于星载合成孔径雷达图像中对应的数据点的绝对值的预定倍数、且方位向时间满足阈值条件时,确定星载合成孔径雷达图像中的数据点满足预定抑制条件。
示例地,当模糊图像为负频率方位模糊图像时,当模糊图像满足以下等式(3)时,可以确定星载合成孔径雷达图像中的数据点满足预定抑制条件:
|s-1(τ,η)|>m×|s0(τ,η)|,且
当星载合成孔径雷达图像中的数据点满足预定抑制条件时,可以用逻辑值1来表示;当星载合成孔径雷达图像中的数据点不满足预定抑制条件时,可以用逻辑值0来表示。此外,还可以将该判定结果保存在负模糊判定矩阵中。
示例地,可以通过以下等式(4)来建立负模糊判定矩阵:
其中,B-1(τ,η)为所述负模糊判定矩阵。
又示例地,当模糊图像为正频率方位模糊图像时,当模糊图像满足以下等式(5)时,可以确定星载合成孔径雷达图像中的数据点满足预定抑制条件:
|s+1(τ,η)|>m×|s0(τ,η)|,且
同样地,可以通过以下等式(6)来建立正模糊判定矩阵:
其中,B+1(τ,η)为所述正模糊判定矩阵。
在S13中,当星载合成孔径雷达图像中的数据点满足预定抑制条件时,利用模糊图像,对星载合成孔径雷达图像中、满足预定抑制条件的数据点进行模糊抑制。
在一种实施方式中,可以通过将星载合成孔径雷达图像中、满足预定抑制条件的数据点的绝对值,与模糊图像中对应的数据点的绝对值做差值,来得到经模糊抑制后的图像。
示例地,可以通过以下等式(7)来得到经模糊抑制后的图像:
其中,sfinal(τ,η)为所述经模糊抑制后的图像。
通过上述技术方案,对星载合成孔径雷达图像进行方位模糊重建,得到模糊图像,并利用模糊图像,对星载合成孔径雷达图像中、满足预定抑制条件的数据点进行模糊抑制。这样,能够有效去除星载合成孔径雷达图像中的方位模糊,提高星载合成孔径雷达图像的清晰度。
图3A是根据一示例性实施例示出的一种星载合成孔径雷达图像的模糊抑制装置。参照图3A,该装置30可以包括:方位模糊重建模块301,用于对星载合成孔径雷达图像进行方位模糊重建,确定出所述星载合成孔径雷达图像的模糊图像;判断模块302,用于根据与所述星载合成孔径雷达图像中的数据点对应的所述模糊图像的数据点,判断所述星载合成孔径雷达图像中的数据点是否满足预定抑制条件;模糊抑制模块303,用于当所述星载合成孔径雷达图像中的数据点满足所述预定抑制条件时,利用所述模糊图像,对所述星载合成孔径雷达图像中、满足所述预定抑制条件的数据点进行模糊抑制。
可选地,所述模糊图像包括负频率方位模糊图像;所述方位模糊重建模块301用于:对所述星载合成孔径雷达图像s0(τ,η)按列做快速傅里叶变换,得到图像S1(τ,fη),其中,j为列索引,Nr为图像矩阵总列数,fs为信号采样率,i为行索引,Na为图像矩阵总行数,fp为脉冲重复频率,τ为距离向时间,η为方位向时间,fη为方位向频率;
对所述图像S1(τ,fη)按列进行方位相位补偿和天线方向图补偿,得到图像S2(τ,fη):
S2(τ,fη)=H1·S1(τ,fη)
其中,H1为传递函数,f0为信号载频,R0为斜距,c为光速,Rref为场景中心斜距,λ为波长,Vr为等效速度,fηc为多普勒中心频率,D(fη-fp)为徙动因子,Wa(fη-fηc-fp)为卫星方位向天线方向图的频域表达式,wa(τ)为方向图的时域表达式;
对所述图像S2(τ,fη)按行做快速傅里叶变换,得到图像S3(fτ,fη);
对所述图像S3(fτ,fη)按行进行一致模糊距离徙动校正和方位向调制,得到图像S4(fτ,fη);
对所述图像S4(fτ,fη)按行做快速逆傅里叶变换,得到图像S5(τ,fη);
对所述图像S5(τ,fη)按行进行补余模糊距离徙动校正,得到图像S6(τ,fη);
对所述图像S6(τ,fη)按行做快速傅里叶变换,得到图像S7(fτ,fη);
对所述图像S7(fτ,fη)按行进行脉冲压缩和二次距离压缩误差校正,得到图像S8(fτ,fη);
对所述图像S8(fτ,fη)按行做快速逆傅里叶变换,得到图像S9(τ,fη);
对所述图像S9(τ,fη)补偿残余误差和进行方位散焦处理,得到图像S10(τ,fη);
对所述图像S10(τ,fη)按列做快速逆傅里叶变换,得到所述负频率方位模糊图像。
可选地,所述模糊图像包括正频率方位模糊图像;
所述方位模糊重建模块301用于:
对所述星载合成孔径雷达图像s0(τ,η)按列做快速傅里叶变换,得到图像S11(τ,fη),其中,j为列索引,Nr为图像矩阵总列数,fs为信号采样率,i为行索引,Na为图像矩阵总行数,fp为脉冲重复频率,τ为距离向时间,η为方位向时间,fη为方位向频率;
对所述图像S11(τ,fη)按列进行方位相位补偿和天线方向图补偿,得到图像S12(τ,fη):
S12(τ,fη)=H1·S11(τ,fη)
其中,H1为传递函数,f0为信号载频,R0为斜距,c为光速,Rref为场景中心斜距,λ为波长,Vr为等效速度,fηc为多普勒中心频率,D(fη+fp)为徙动因子,Wa(fη-fηc+fp)为卫星方位向天线方向图的频域表达式,wa(τ)为方向图的时域表达式;
对所述图像S12(τ,fη)按行做快速傅里叶变换,得到图像S13(fτ,fη);
对所述图像S13(fτ,fη)按行进行一致模糊距离徙动校正和方位向调制,得到图像S14(fτ,fη);
对所述图像S14(fτ,fη)按行做快速逆傅里叶变换,得到图像S15(τ,fη);
对所述图像S15(τ,fη)按行进行补余模糊距离徙动校正,得到图像S16(τ,fη);
对所述图像S16(τ,fη)按行做快速傅里叶变换,得到图像S17(fτ,fη);
对所述图像S17(fτ,fη)按行进行脉冲压缩和二次距离压缩误差校正,得到图像S18(fτ,fη);
对所述图像S18(fτ,fη)按行做快速逆傅里叶变换,得到图像S19(τ,fη);
对所述图像S19(τ,fη)补偿残余误差和进行方位散焦处理,得到图像S20(τ,fη);
对所述图像S20(τ,fη)按列做快速逆傅里叶变换,得到所述正频率方位模糊图像。
如图3B所示,所述判断模块302可以包括:第一判断子模块3021,用于判断所述模糊图像的数据点的绝对值是否大于所述星载合成孔径雷达图像中对应的数据点的绝对值的预定倍数,其中,所述预定倍数大于0且小于1;第二判断子模块3022,用于判断所述方位向时间是否满足预设的阈值条件;确定子模块3023,用于当所述模糊图像的数据点的绝对值大于所述星载合成孔径雷达图像中对应的数据点的绝对值的预定倍数、且所述方位向时间满足所述阈值条件时,确定所述星载合成孔径雷达图像中的数据点满足所述预定抑制条件。
可选地,所述模糊图像包括负频率方位模糊图像;所述阈值条件包括:
其中,η为所述方位向时间;fp为脉冲重复频率;fr为多普勒调频率;和/或
所述模糊图像包括正频率方位模糊图像;所述阈值条件包括:
其中,η为所述方位向时间;Na为图像矩阵总行数;fp为脉冲重复频率;fr为多普勒调频率。
可选地,所述模糊抑制模块303用于将所述星载合成孔径雷达图像中、满足所述预定抑制条件的数据点的绝对值,与所述模糊图像中对应的数据点的绝对值做差值,得到经模糊抑制后的图像。
以上结合附图详细描述了本公开的优选实施方式,但是,本公开并不限于上述实施方式中的具体细节,在本公开的技术构思范围内,可以对本公开的技术方案进行多种简单变型,这些简单变型均属于本公开的保护范围。
另外需要说明的是,在上述具体实施方式中所描述的各个具体技术特征,在不矛盾的情况下,可以通过任何合适的方式进行组合。为了避免不必要的重复,本公开对各种可能的组合方式不再另行说明。
此外,本公开的各种不同的实施方式之间也可以进行任意组合,只要其不违背本公开的思想,其同样应当视为本公开所公开的内容。
Claims (12)
1.一种星载合成孔径雷达图像的模糊抑制方法,其特征在于,包括:
对星载合成孔径雷达图像进行方位模糊重建,确定出所述星载合成孔径雷达图像的模糊图像;
根据与所述星载合成孔径雷达图像中的数据点对应的所述模糊图像的数据点,判断所述星载合成孔径雷达图像中的数据点是否满足预定抑制条件;
当所述星载合成孔径雷达图像中的数据点满足所述预定抑制条件时,利用所述模糊图像,对所述星载合成孔径雷达图像中、满足所述预定抑制条件的数据点进行模糊抑制。
2.根据权利要求1所述的方法,其特征在于,所述模糊图像包括负频率方位模糊图像;
所述对星载合成孔径雷达图像进行方位模糊重建,确定出所述星载合成孔径雷达图像的模糊图像,包括:
对所述星载合成孔径雷达图像s0(τ,η)按列做快速傅里叶变换,得到图像S1(τ,fη),其中,j为列索引,Nr为图像矩阵总列数,fs为信号采样率,i为行索引,Na为图像矩阵总行数,fp为脉冲重复频率,τ为距离向时间,η为方位向时间,fη为方位向频率;
对所述图像S1(τ,fη)按列进行方位相位补偿和天线方向图补偿,得到图像S2(τ,fη):
S2(τ,fη)=H1·S1(τ,fη)
<mrow>
<msub>
<mi>H</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mi>exp</mi>
<mo>{</mo>
<mo>-</mo>
<mi>j</mi>
<mfrac>
<mrow>
<mn>4</mn>
<msub>
<mi>&pi;f</mi>
<mn>0</mn>
</msub>
<msub>
<mi>R</mi>
<mn>0</mn>
</msub>
<mi>D</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>&eta;</mi>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mi>p</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mi>c</mi>
</mfrac>
<mo>}</mo>
<mo>&CenterDot;</mo>
<mfrac>
<mrow>
<msub>
<mi>W</mi>
<mi>a</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>&eta;</mi>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>&eta;</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mi>p</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>W</mi>
<mi>a</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>&eta;</mi>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>&eta;</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>R</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>-</mo>
<mfrac>
<msub>
<mi>N</mi>
<mi>r</mi>
</msub>
<mn>2</mn>
</mfrac>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mfrac>
<mi>c</mi>
<mrow>
<mn>2</mn>
<msub>
<mi>f</mi>
<mi>s</mi>
</msub>
</mrow>
</mfrac>
</mrow>
<mrow>
<mi>D</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>&eta;</mi>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mi>p</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msqrt>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<mrow>
<msup>
<mi>&lambda;</mi>
<mn>2</mn>
</msup>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>&eta;</mi>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mi>p</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mn>4</mn>
<msubsup>
<mi>V</mi>
<mi>r</mi>
<mn>2</mn>
</msubsup>
</mrow>
</mfrac>
</mrow>
</msqrt>
</mrow>
<mrow>
<msub>
<mi>f</mi>
<mi>&eta;</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>-</mo>
<mfrac>
<msub>
<mi>N</mi>
<mi>a</mi>
</msub>
<mn>2</mn>
</mfrac>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mfrac>
<msub>
<mi>f</mi>
<mi>p</mi>
</msub>
<msub>
<mi>N</mi>
<mi>a</mi>
</msub>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>W</mi>
<mi>a</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>&eta;</mi>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>&eta;</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>w</mi>
<mi>a</mi>
</msub>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>cR</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>&eta;</mi>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>&eta;</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mn>2</mn>
<msub>
<mi>f</mi>
<mn>0</mn>
</msub>
<msubsup>
<mi>V</mi>
<mi>r</mi>
<mn>2</mn>
</msubsup>
<msqrt>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<mrow>
<msup>
<mi>c</mi>
<mn>2</mn>
</msup>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>&eta;</mi>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>&eta;</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mn>4</mn>
<msubsup>
<mi>V</mi>
<mi>r</mi>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>f</mi>
<mn>0</mn>
<mn>2</mn>
</msubsup>
</mrow>
</mfrac>
</mrow>
</msqrt>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>W</mi>
<mi>a</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>&eta;</mi>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>&eta;</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mi>p</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>w</mi>
<mi>a</mi>
</msub>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>cR</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>&eta;</mi>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>&eta;</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mi>p</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mn>2</mn>
<msub>
<mi>f</mi>
<mn>0</mn>
</msub>
<msubsup>
<mi>V</mi>
<mi>r</mi>
<mn>2</mn>
</msubsup>
<msqrt>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<mrow>
<msup>
<mi>c</mi>
<mn>2</mn>
</msup>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>&eta;</mi>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>&eta;</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mi>p</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mn>4</mn>
<msubsup>
<mi>V</mi>
<mi>r</mi>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>f</mi>
<mn>0</mn>
<mn>2</mn>
</msubsup>
</mrow>
</mfrac>
</mrow>
</msqrt>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
1
<mrow>
<msub>
<mi>w</mi>
<mi>a</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mi>&pi;</mi>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<mi>&pi;</mi>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mfrac>
</mrow>
其中,H1为传递函数,f0为信号载频,R0为斜距,c为光速,Rref为场景中心斜距,λ为波长,Vr为等效速度,fηc为多普勒中心频率,D(fη-fp)为徙动因子,Wa(fη-fηc-fp)为卫星方位向天线方向图的频域表达式,wa(τ)为方向图的时域表达式;
对所述图像S2(τ,fη)按行做快速傅里叶变换,得到图像S3(fτ,fη);
对所述图像S3(fτ,fη)按行进行一致模糊距离徙动校正和方位向调制,得到图像S4(fτ,fη);
对所述图像S4(fτ,fη)按行做快速逆傅里叶变换,得到图像S5(τ,fη);
对所述图像S5(τ,fη)按行进行补余模糊距离徙动校正,得到图像S6(τ,fη);
对所述图像S6(τ,fη)按行做快速傅里叶变换,得到图像S7(fτ,fη);
对所述图像S7(fτ,fη)按行进行脉冲压缩和二次距离压缩误差校正,得到图像S8(fτ,fη);
对所述图像S8(fτ,fη)按行做快速逆傅里叶变换,得到图像S9(τ,fη);
对所述图像S9(τ,fη)补偿残余误差和进行方位散焦处理,得到图像S10(τ,fη);
对所述图像S10(τ,fη)按列做快速逆傅里叶变换,得到所述负频率方位模糊图像。
3.根据权利要求1或2所述的方法,其特征在于,所述模糊图像包括正频率方位模糊图像;
所述对星载合成孔径雷达图像进行方位模糊重建,确定出所述星载合成孔径雷达图像的模糊图像,包括:
对所述星载合成孔径雷达图像s0(τ,η)按列做快速傅里叶变换,得到图像S11(τ,fη),其中,j为列索引,Nr为图像矩阵总列数,fs为信号采样率,i为行索引,Na为图像矩阵总行数,fp为脉冲重复频率,τ为距离向时间,η为方位向时间,fη为方位向频率;
对所述图像S11(τ,fη)按列进行方位相位补偿和天线方向图补偿,得到图像S12(τ,fη):
S12(τ,fη)=H1·S11(τ,fη)
<mrow>
<msub>
<mi>H</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mi>exp</mi>
<mo>{</mo>
<mo>-</mo>
<mi>j</mi>
<mfrac>
<mrow>
<mn>4</mn>
<msub>
<mi>&pi;f</mi>
<mn>0</mn>
</msub>
<msub>
<mi>R</mi>
<mn>0</mn>
</msub>
<mi>D</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>&eta;</mi>
</msub>
<mo>+</mo>
<msub>
<mi>f</mi>
<mi>p</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mi>c</mi>
</mfrac>
<mo>}</mo>
<mo>&CenterDot;</mo>
<mfrac>
<mrow>
<msub>
<mi>W</mi>
<mi>a</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>&eta;</mi>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>&eta;</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>f</mi>
<mi>p</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>W</mi>
<mi>a</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>&eta;</mi>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>&eta;</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>R</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>-</mo>
<mfrac>
<msub>
<mi>N</mi>
<mi>r</mi>
</msub>
<mn>2</mn>
</mfrac>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mfrac>
<mi>c</mi>
<mrow>
<mn>2</mn>
<msub>
<mi>f</mi>
<mi>s</mi>
</msub>
</mrow>
</mfrac>
</mrow>
<mrow>
<mi>D</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>&eta;</mi>
</msub>
<mo>+</mo>
<msub>
<mi>f</mi>
<mi>p</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msqrt>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<mrow>
<msup>
<mi>&lambda;</mi>
<mn>2</mn>
</msup>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>&eta;</mi>
</msub>
<mo>+</mo>
<msub>
<mi>f</mi>
<mi>p</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mn>4</mn>
<msubsup>
<mi>V</mi>
<mi>r</mi>
<mn>2</mn>
</msubsup>
</mrow>
</mfrac>
</mrow>
</msqrt>
</mrow>
<mrow>
<msub>
<mi>f</mi>
<mi>&eta;</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>-</mo>
<mfrac>
<msub>
<mi>N</mi>
<mi>a</mi>
</msub>
<mn>2</mn>
</mfrac>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mfrac>
<msub>
<mi>f</mi>
<mi>p</mi>
</msub>
<msub>
<mi>N</mi>
<mi>a</mi>
</msub>
</mfrac>
</mrow>
2
<mrow>
<msub>
<mi>W</mi>
<mi>a</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>&eta;</mi>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>&eta;</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>w</mi>
<mi>a</mi>
</msub>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>cR</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>&eta;</mi>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>&eta;</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mn>2</mn>
<msub>
<mi>f</mi>
<mn>0</mn>
</msub>
<msubsup>
<mi>V</mi>
<mi>r</mi>
<mn>2</mn>
</msubsup>
<msqrt>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<mrow>
<msup>
<mi>c</mi>
<mn>2</mn>
</msup>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>&eta;</mi>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>&eta;</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mn>4</mn>
<msubsup>
<mi>V</mi>
<mi>r</mi>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>f</mi>
<mn>0</mn>
<mn>2</mn>
</msubsup>
</mrow>
</mfrac>
</mrow>
</msqrt>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>W</mi>
<mi>a</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>&eta;</mi>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>&eta;</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>f</mi>
<mi>p</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>w</mi>
<mi>a</mi>
</msub>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>cR</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>&eta;</mi>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>&eta;</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>f</mi>
<mi>p</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mn>2</mn>
<msub>
<mi>f</mi>
<mn>0</mn>
</msub>
<msubsup>
<mi>V</mi>
<mi>r</mi>
<mn>2</mn>
</msubsup>
<msqrt>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<mrow>
<msup>
<mi>c</mi>
<mn>2</mn>
</msup>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>&eta;</mi>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>&eta;</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>f</mi>
<mi>p</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mn>4</mn>
<msubsup>
<mi>V</mi>
<mi>r</mi>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>f</mi>
<mn>0</mn>
<mn>2</mn>
</msubsup>
</mrow>
</mfrac>
</mrow>
</msqrt>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>w</mi>
<mi>a</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mi>&pi;</mi>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<mi>&pi;</mi>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mfrac>
</mrow>
其中,H1为传递函数,f0为信号载频,R0为斜距,c为光速,Rref为场景中心斜距,λ为波长,Vr为等效速度,fηc为多普勒中心频率,D(fη-fp)为徙动因子,Wa(fη-fηc+fp)为卫星方位向天线方向图的频域表达式,wa(τ)为方向图的时域表达式;
对所述图像S12(τ,fη)按行做快速傅里叶变换,得到图像S13(fτ,fη);
对所述图像S13(fτ,fη)按行进行一致模糊距离徙动校正和方位向调制,得到图像S14(fτ,fη);
对所述图像S14(fτ,fη)按行做快速逆傅里叶变换,得到图像S15(τ,fη);
对所述图像S15(τ,fη)按行进行补余模糊距离徙动校正,得到图像S16(τ,fη);
对所述图像S16(τ,fη)按行做快速傅里叶变换,得到图像S17(fτ,fη);
对所述图像S17(fτ,fη)按行进行脉冲压缩和二次距离压缩误差校正,得到图像S18(fτ,fη);
对所述图像S18(fτ,fη)按行做快速逆傅里叶变换,得到图像S19(τ,fη);
对所述图像S19(τ,fη)补偿残余误差和进行方位散焦处理,得到图像S20(τ,fη);
对所述图像S20(τ,fη)按列做快速逆傅里叶变换,得到所述正频率方位模糊图像。
4.根据权利要求1所述的方法,其特征在于,所述根据与所述星载合成孔径雷达图像中的数据点对应的所述模糊图像的数据点,判断所述星载合成孔径雷达图像中的数据点是否满足预定抑制条件,包括:
判断所述模糊图像的数据点的绝对值是否大于所述星载合成孔径雷达图像中对应的数据点的绝对值的预定倍数,其中,所述预定倍数大于0且小于1;
判断方位向时间是否满足预设的阈值条件;
当所述模糊图像的数据点的绝对值大于所述星载合成孔径雷达图像中对应的数据点的绝对值的预定倍数、且所述方位向时间满足所述阈值条件时,确定所述星载合成孔径雷达图像中的数据点满足所述预定抑制条件。
5.根据权利要求4所述的方法,其特征在于,所述模糊图像包括负频率方位模糊图像;
所述阈值条件包括:
<mrow>
<mi>&eta;</mi>
<mo>></mo>
<mfrac>
<msubsup>
<mi>f</mi>
<mi>p</mi>
<mn>2</mn>
</msubsup>
<msub>
<mi>f</mi>
<mi>r</mi>
</msub>
</mfrac>
</mrow>
其中,η为所述方位向时间;fp为脉冲重复频率;fr为多普勒调频率;和/或
所述模糊图像包括正频率方位模糊图像;
所述阈值条件包括:
<mrow>
<mi>&eta;</mi>
<mo><</mo>
<msub>
<mi>N</mi>
<mi>a</mi>
</msub>
<mo>-</mo>
<mfrac>
<msubsup>
<mi>f</mi>
<mi>p</mi>
<mn>2</mn>
</msubsup>
<msub>
<mi>f</mi>
<mi>r</mi>
</msub>
</mfrac>
</mrow>
其中,η为所述方位向时间;Na为图像矩阵总行数;fp为脉冲重复频率;fr为多普勒调频率。
6.根据权利要求1所述的方法,其特征在于,所述当所述星载合成孔径雷达图像中的数据点满足所述预定抑制条件时,利用所述模糊图像,对所述星载合成孔径雷达图像中的数据点进行模糊抑制,包括:
将所述星载合成孔径雷达图像中、满足所述预定抑制条件的数据点的绝对值,与所述模糊图像中对应的数据点的绝对值做差值,得到经模糊抑制后的图像。
7.一种星载合成孔径雷达图像的模糊抑制装置,其特征在于,包括:
方位模糊重建模块,用于对星载合成孔径雷达图像进行方位模糊重建,确定出所述星载合成孔径雷达图像的模糊图像;
判断模块,用于根据与所述星载合成孔径雷达图像中的数据点对应的所述模糊图像的数据点,判断所述星载合成孔径雷达图像中的数据点是否满足预定抑制条件;
模糊抑制模块,用于当所述星载合成孔径雷达图像中的数据点满足所述预定抑制条件时,利用所述模糊图像,对所述星载合成孔径雷达图像中、满足所述预定抑制条件的数据点进行模糊抑制。
8.根据权利要求7所述的装置,其特征在于,所述模糊图像包括负频率方位模糊图像;
所述方位模糊重建模块用于:
对所述星载合成孔径雷达图像s0(τ,η)按列做快速傅里叶变换,得到图像S1(τ,fη),其中,j为列索引,Nr为图像矩阵总列数,fs为信号采样率,i为行索引,Na为图像矩阵总行数,fp为脉冲重复频率,τ为距离向时间,η为方位向时间,fη为方位向频率;
对所述图像S1(τ,fη)按列进行方位相位补偿和天线方向图补偿,得到图像S2(τ,fη):
S2(τ,fη)=H1·S1(τ,fη)
<mrow>
<msub>
<mi>H</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mi>exp</mi>
<mo>{</mo>
<mo>-</mo>
<mi>j</mi>
<mfrac>
<mrow>
<mn>4</mn>
<msub>
<mi>&pi;f</mi>
<mn>0</mn>
</msub>
<msub>
<mi>R</mi>
<mn>0</mn>
</msub>
<mi>D</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>&eta;</mi>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mi>p</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mi>c</mi>
</mfrac>
<mo>}</mo>
<mo>&CenterDot;</mo>
<mfrac>
<mrow>
<msub>
<mi>W</mi>
<mi>a</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>&eta;</mi>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>&eta;</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mi>p</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>W</mi>
<mi>a</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>&eta;</mi>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>&eta;</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>R</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>-</mo>
<mfrac>
<msub>
<mi>N</mi>
<mi>r</mi>
</msub>
<mn>2</mn>
</mfrac>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mfrac>
<mi>c</mi>
<mrow>
<mn>2</mn>
<msub>
<mi>f</mi>
<mi>s</mi>
</msub>
</mrow>
</mfrac>
</mrow>
<mrow>
<mi>D</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>&eta;</mi>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mi>p</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msqrt>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<mrow>
<msup>
<mi>&lambda;</mi>
<mn>2</mn>
</msup>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>&eta;</mi>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mi>p</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mn>4</mn>
<msubsup>
<mi>V</mi>
<mi>r</mi>
<mn>2</mn>
</msubsup>
</mrow>
</mfrac>
</mrow>
</msqrt>
</mrow>
4
<mrow>
<msub>
<mi>f</mi>
<mi>&eta;</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>-</mo>
<mfrac>
<msub>
<mi>N</mi>
<mi>a</mi>
</msub>
<mn>2</mn>
</mfrac>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mfrac>
<msub>
<mi>f</mi>
<mi>p</mi>
</msub>
<msub>
<mi>N</mi>
<mi>a</mi>
</msub>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>W</mi>
<mi>a</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>&eta;</mi>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>&eta;</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>w</mi>
<mi>a</mi>
</msub>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>cR</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>&eta;</mi>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>&eta;</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mn>2</mn>
<msub>
<mi>f</mi>
<mn>0</mn>
</msub>
<msubsup>
<mi>V</mi>
<mi>r</mi>
<mn>2</mn>
</msubsup>
<msqrt>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<mrow>
<msup>
<mi>c</mi>
<mn>2</mn>
</msup>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>&eta;</mi>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>&eta;</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mn>4</mn>
<msubsup>
<mi>V</mi>
<mi>r</mi>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>f</mi>
<mn>0</mn>
<mn>2</mn>
</msubsup>
</mrow>
</mfrac>
</mrow>
</msqrt>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>W</mi>
<mi>a</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>&eta;</mi>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>&eta;</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mi>p</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>w</mi>
<mi>a</mi>
</msub>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>cR</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>&eta;</mi>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>&eta;</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mi>p</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mn>2</mn>
<msub>
<mi>f</mi>
<mn>0</mn>
</msub>
<msubsup>
<mi>V</mi>
<mi>r</mi>
<mn>2</mn>
</msubsup>
<msqrt>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<mrow>
<msup>
<mi>c</mi>
<mn>2</mn>
</msup>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>&eta;</mi>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>&eta;</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mi>p</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mn>4</mn>
<msubsup>
<mi>V</mi>
<mi>r</mi>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>f</mi>
<mn>0</mn>
<mn>2</mn>
</msubsup>
</mrow>
</mfrac>
</mrow>
</msqrt>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>w</mi>
<mi>a</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mi>&pi;</mi>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<mi>&pi;</mi>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mfrac>
</mrow>
其中,H1为传递函数,f0为信号载频,R0为斜距,c为光速,Rref为场景中心斜距,λ为波长,Vr为等效速度,fηc为多普勒中心频率,D(fη-fp)为徙动因子,Wa(fη-fηc-fp)为卫星方位向天线方向图的频域表达式,wa(τ)为方向图的时域表达式;
对所述图像S2(τ,fη)按行做快速傅里叶变换,得到图像S3(fτ,fη);
对所述图像S3(fτ,fη)按行进行一致模糊距离徙动校正和方位向调制,得到图像S4(fτ,fη);
对所述图像S4(fτ,fη)按行做快速逆傅里叶变换,得到图像S5(τ,fη);
对所述图像S5(τ,fη)按行进行补余模糊距离徙动校正,得到图像S6(τ,fη);
对所述图像S6(τ,fη)按行做快速傅里叶变换,得到图像S7(fτ,fη);
对所述图像S7(fτ,fη)按行进行脉冲压缩和二次距离压缩误差校正,得到图像S8(fτ,fη);
对所述图像S8(fτ,fη)按行做快速逆傅里叶变换,得到图像S9(τ,fη);
对所述图像S9(τ,fη)补偿残余误差和进行方位散焦处理,得到图像S10(τ,fη);
对所述图像S10(τ,fη)按列做快速逆傅里叶变换,得到所述负频率方位模糊图像。
9.根据权利要求7或8所述的装置,其特征在于,所述模糊图像包括正频率方位模糊图像;
所述方位模糊重建模块用于:
对所述星载合成孔径雷达图像s0(τ,η)按列做快速傅里叶变换,得到图像S11(τ,fη),其中,j为列索引,Nr为图像矩阵总列数,fs为信号采样率,i为行索引,Na为图像矩阵总行数,fp为脉冲重复频率,τ为距离向时间,η为方位向时间,fη为方位向频率;
对所述图像S11(τ,fη)按列进行方位相位补偿和天线方向图补偿,得到图像S12(τ,fη):
S12(τ,fη)=H1·S11(τ,fη)
<mrow>
<msub>
<mi>H</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mi>exp</mi>
<mo>{</mo>
<mo>-</mo>
<mi>j</mi>
<mfrac>
<mrow>
<mn>4</mn>
<msub>
<mi>&pi;f</mi>
<mn>0</mn>
</msub>
<msub>
<mi>R</mi>
<mn>0</mn>
</msub>
<mi>D</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>&eta;</mi>
</msub>
<mo>+</mo>
<msub>
<mi>f</mi>
<mi>p</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mi>c</mi>
</mfrac>
<mo>}</mo>
<mo>&CenterDot;</mo>
<mfrac>
<mrow>
<msub>
<mi>W</mi>
<mi>a</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>&eta;</mi>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>&eta;</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>f</mi>
<mi>p</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>W</mi>
<mi>a</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>&eta;</mi>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>&eta;</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>R</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>-</mo>
<mfrac>
<msub>
<mi>N</mi>
<mi>r</mi>
</msub>
<mn>2</mn>
</mfrac>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mfrac>
<mi>c</mi>
<mrow>
<mn>2</mn>
<msub>
<mi>f</mi>
<mi>s</mi>
</msub>
</mrow>
</mfrac>
</mrow>
<mrow>
<mi>D</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>&eta;</mi>
</msub>
<mo>+</mo>
<msub>
<mi>f</mi>
<mi>p</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msqrt>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<mrow>
<msup>
<mi>&lambda;</mi>
<mn>2</mn>
</msup>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>&eta;</mi>
</msub>
<mo>+</mo>
<msub>
<mi>f</mi>
<mi>p</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mn>4</mn>
<msubsup>
<mi>V</mi>
<mi>r</mi>
<mn>2</mn>
</msubsup>
</mrow>
</mfrac>
</mrow>
</msqrt>
</mrow>
<mrow>
<msub>
<mi>f</mi>
<mi>&eta;</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>-</mo>
<mfrac>
<msub>
<mi>N</mi>
<mi>a</mi>
</msub>
<mn>2</mn>
</mfrac>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mfrac>
<msub>
<mi>f</mi>
<mi>p</mi>
</msub>
<msub>
<mi>N</mi>
<mi>a</mi>
</msub>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>W</mi>
<mi>a</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>&eta;</mi>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>&eta;</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>w</mi>
<mi>a</mi>
</msub>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>cR</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>&eta;</mi>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>&eta;</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mn>2</mn>
<msub>
<mi>f</mi>
<mn>0</mn>
</msub>
<msubsup>
<mi>V</mi>
<mi>r</mi>
<mn>2</mn>
</msubsup>
<msqrt>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<mrow>
<msup>
<mi>c</mi>
<mn>2</mn>
</msup>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>&eta;</mi>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>&eta;</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mn>4</mn>
<msubsup>
<mi>V</mi>
<mi>r</mi>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>f</mi>
<mn>0</mn>
<mn>2</mn>
</msubsup>
</mrow>
</mfrac>
</mrow>
</msqrt>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>W</mi>
<mi>a</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>&eta;</mi>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>&eta;</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>f</mi>
<mi>p</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>w</mi>
<mi>a</mi>
</msub>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>cR</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>&eta;</mi>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>&eta;</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>f</mi>
<mi>p</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mn>2</mn>
<msub>
<mi>f</mi>
<mn>0</mn>
</msub>
<msubsup>
<mi>V</mi>
<mi>r</mi>
<mn>2</mn>
</msubsup>
<msqrt>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<mrow>
<msup>
<mi>c</mi>
<mn>2</mn>
</msup>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>&eta;</mi>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>&eta;</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>f</mi>
<mi>p</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mn>4</mn>
<msubsup>
<mi>V</mi>
<mi>r</mi>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>f</mi>
<mn>0</mn>
<mn>2</mn>
</msubsup>
</mrow>
</mfrac>
</mrow>
</msqrt>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>w</mi>
<mi>a</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mi>&pi;</mi>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<mi>&pi;</mi>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mfrac>
</mrow>
其中,H1为传递函数,f0为信号载频,R0为斜距,c为光速,Rref为场景中心斜距,λ为波长,Vr为等效速度,fηc为多普勒中心频率,D(fη+fp)为徙动因子,Wa(fη-fηc+fp)为卫星方位向天线方向图的频域表达式,wa(τ)为方向图的时域表达式;
对所述图像S12(τ,fη)按行做快速傅里叶变换,得到图像S13(fτ,fη);
对所述图像S13(fτ,fη)按行进行一致模糊距离徙动校正和方位向调制,得到图像S14(fτ,fη);
对所述图像S14(fτ,fη)按行做快速逆傅里叶变换,得到图像S15(τ,fη);
对所述图像S15(τ,fη)按行进行补余模糊距离徙动校正,得到图像S16(τ,fη);
对所述图像S16(τ,fη)按行做快速傅里叶变换,得到图像S17(fτ,fη);
对所述图像S17(fτ,fη)按行进行脉冲压缩和二次距离压缩误差校正,得到图像S18(fτ,fη);
对所述图像S18(fτ,fη)按行做快速逆傅里叶变换,得到图像S19(τ,fη);
对所述图像S19(τ,fη)补偿残余误差和进行方位散焦处理,得到图像S20(τ,fη);
对所述图像S20(τ,fη)按列做快速逆傅里叶变换,得到所述正频率方位模糊图像。
10.根据权利要求7所述的装置,其特征在于,所述判断模块包括:
第一判断子模块,用于判断所述模糊图像的数据点的绝对值是否大于所述星载合成孔径雷达图像中对应的数据点的绝对值的预定倍数,其中,所述预定倍数大于0且小于1;
第二判断子模块,用于判断方位向时间是否满足预设的阈值条件;
确定子模块,用于当所述模糊图像的数据点的绝对值大于所述星载合成孔径雷达图像中对应的数据点的绝对值的预定倍数、且所述方位向时间满足所述阈值条件时,确定所述星载合成孔径雷达图像中的数据点满足所述预定抑制条件。
11.根据权利要求10所述的装置,其特征在于,所述模糊图像包括负频率方位模糊图像;
所述阈值条件包括:
<mrow>
<mi>&eta;</mi>
<mo>></mo>
<mfrac>
<msubsup>
<mi>f</mi>
<mi>p</mi>
<mn>2</mn>
</msubsup>
<msub>
<mi>f</mi>
<mi>r</mi>
</msub>
</mfrac>
</mrow>
其中,η为所述方位向时间;fp为脉冲重复频率;fr为多普勒调频率;和/或
所述模糊图像包括正频率方位模糊图像;
所述阈值条件包括:
<mrow>
<mi>&eta;</mi>
<mo><</mo>
<msub>
<mi>N</mi>
<mi>a</mi>
</msub>
<mo>-</mo>
<mfrac>
<msubsup>
<mi>f</mi>
<mi>p</mi>
<mn>2</mn>
</msubsup>
<msub>
<mi>f</mi>
<mi>r</mi>
</msub>
</mfrac>
</mrow>
其中,η为所述方位向时间;Na为图像矩阵总行数;fp为脉冲重复频率;fr为多普勒调频率。
12.根据权利要求7所述的装置,其特征在于,所述模糊抑制模块用于将所述星载合成孔径雷达图像中、满足所述预定抑制条件的数据点的绝对值,与所述模糊图像中对应的数据点的绝对值做差值,得到经模糊抑制后的图像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710134700.3A CN107064895B (zh) | 2017-03-08 | 2017-03-08 | 星载合成孔径雷达图像的模糊抑制方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710134700.3A CN107064895B (zh) | 2017-03-08 | 2017-03-08 | 星载合成孔径雷达图像的模糊抑制方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107064895A true CN107064895A (zh) | 2017-08-18 |
CN107064895B CN107064895B (zh) | 2019-09-17 |
Family
ID=59621570
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710134700.3A Active CN107064895B (zh) | 2017-03-08 | 2017-03-08 | 星载合成孔径雷达图像的模糊抑制方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107064895B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107976674A (zh) * | 2017-10-31 | 2018-05-01 | 西安空间无线电技术研究所 | 一种基于最优方位模糊的双波束sar的prf选择方法及系统 |
CN108389166A (zh) * | 2017-11-21 | 2018-08-10 | 北京航空航天大学 | 图像处理方法、装置、设备及计算机可读存储介质 |
WO2022073287A1 (zh) * | 2020-10-10 | 2022-04-14 | 深圳大学 | 一种sar数据点可视化方法、智能终端及存储介质 |
-
2017
- 2017-03-08 CN CN201710134700.3A patent/CN107064895B/zh active Active
Non-Patent Citations (3)
Title |
---|
JIE CHEN ET AL.: "Accurate Reconstruction and Suppression for Azimuth Ambiguities in Spaceborne stripmap SAR Images", 《IEEE GEOSCIENCE AND REMOTE SENSING LETTERS》 * |
吴正升 等: "一种改进的SAR图像方位模糊抑制方法", 《测绘通报》 * |
邴磊 等: "星载SAR影像上船舶方位向模糊去除算法", 《中国图象图形学报》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107976674A (zh) * | 2017-10-31 | 2018-05-01 | 西安空间无线电技术研究所 | 一种基于最优方位模糊的双波束sar的prf选择方法及系统 |
CN108389166A (zh) * | 2017-11-21 | 2018-08-10 | 北京航空航天大学 | 图像处理方法、装置、设备及计算机可读存储介质 |
CN108389166B (zh) * | 2017-11-21 | 2021-08-13 | 北京航空航天大学 | 模糊覆盖区域处理方法、装置、设备及计算机可读存储介质 |
WO2022073287A1 (zh) * | 2020-10-10 | 2022-04-14 | 深圳大学 | 一种sar数据点可视化方法、智能终端及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN107064895B (zh) | 2019-09-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP3144702B1 (en) | Method and device for synthethic aperture radar imaging based on non-linear frequency modulation signal | |
Kim et al. | Correcting distortion of polarimetric SAR data induced by ionospheric scintillation | |
EP3012658B1 (en) | Method and device for implementing sar imaging | |
CN107064895A (zh) | 星载合成孔径雷达图像的模糊抑制方法和装置 | |
CN107942327A (zh) | 基于脉冲相位编码的单通道hrws‑sar成像方法 | |
CN105974376B (zh) | 一种sar射频干扰抑制方法 | |
CN104635221B (zh) | 一种基于内定标数据的子带拼接方法 | |
CN104330808B (zh) | 基于解重扩技术的多类卫星导航干扰抑制方法 | |
CN101425176A (zh) | 一种基于中值滤波的图像小波去噪方法 | |
CN103323818A (zh) | 多通道合成孔径雷达系统非均匀采样奇异点的方法和装置 | |
CN103901416A (zh) | 一种基于稳健主成分分析法的多通道杂波抑制方法 | |
US7460059B1 (en) | Removing interfering clutter associated with radar pulses that an airborne radar receives from a radar transponder | |
CN102520396A (zh) | 基于复数经验模态分解的时变窄带干扰抑制方法 | |
CN105223572A (zh) | 一种基于pfa算法的正前视双基sar成像处理方法 | |
CN103760534B (zh) | 一种星载sar数据的电离层色散效应校正方法 | |
CN103885061A (zh) | 基于改进相位梯度自聚焦的双基sar运动补偿方法 | |
CN106443721A (zh) | 一种快速收敛的卫星导航天线阵抗干扰方法 | |
CN104393883A (zh) | 基于改进的功率倒置法的ads-b压制式干扰抑制方法 | |
CN107102328B (zh) | 基于fpga的实时成像信号处理方法及fpga | |
Rabaste et al. | Rejection of doppler shifted multipaths in airborne passive radar | |
CN103576153A (zh) | 一种方位向多波束合成孔径雷达及其实现方法和装置 | |
CN103630911B (zh) | 一种导航信号的处理方法及装置 | |
Li et al. | An autofocus scheme of bistatic SAR considering cross-cell residual range migration | |
CN103323820B (zh) | 一种抑制频率步进sar一维距离图像栅瓣的方法 | |
CN115616628B (zh) | 基于角跟踪环路的gnss天线阵接收机盲波束形成方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |