CN107064895B - 星载合成孔径雷达图像的模糊抑制方法和装置 - Google Patents

星载合成孔径雷达图像的模糊抑制方法和装置 Download PDF

Info

Publication number
CN107064895B
CN107064895B CN201710134700.3A CN201710134700A CN107064895B CN 107064895 B CN107064895 B CN 107064895B CN 201710134700 A CN201710134700 A CN 201710134700A CN 107064895 B CN107064895 B CN 107064895B
Authority
CN
China
Prior art keywords
image
satellite
aperture radar
synthetic aperture
borne synthetic
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.)
Active
Application number
CN201710134700.3A
Other languages
English (en)
Other versions
CN107064895A (zh
Inventor
陈杰
王凯
匡辉
曾虹程
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beihang University
Original Assignee
Beihang University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beihang University filed Critical Beihang University
Priority to CN201710134700.3A priority Critical patent/CN107064895B/zh
Publication of CN107064895A publication Critical patent/CN107064895A/zh
Application granted granted Critical
Publication of CN107064895B publication Critical patent/CN107064895B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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 (10)

1.一种星载合成孔径雷达图像的模糊抑制方法，其特征在于，包括：
对星载合成孔径雷达图像进行方位模糊重建，确定出所述星载合成孔径雷达图像的模糊图像；
判断所述模糊图像的数据点的绝对值是否大于所述星载合成孔径雷达图像中对应的数据点的绝对值的预定倍数，其中，所述预定倍数大于0且小于1；
判断方位向时间是否满足预设的阈值条件；
当所述模糊图像的数据点的绝对值大于所述星载合成孔径雷达图像中对应的数据点的绝对值的预定倍数、且所述方位向时间满足所述阈值条件时，确定所述星载合成孔径雷达图像中的数据点满足预定抑制条件；
当所述星载合成孔径雷达图像中的数据点满足所述预定抑制条件时，利用所述模糊图像，对所述星载合成孔径雷达图像中、满足所述预定抑制条件的数据点进行模糊抑制。
2.根据权利要求1所述的方法，其特征在于，所述模糊图像包括负频率方位模糊图像；
所述对星载合成孔径雷达图像进行方位模糊重建，确定出所述星载合成孔径雷达图像的模糊图像，包括：
对所述星载合成孔径雷达图像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η)按列做快速逆傅里叶变换，得到所述负频率方位模糊图像。
3.根据权利要求1或2所述的方法，其特征在于，所述模糊图像包括正频率方位模糊图像；
所述对星载合成孔径雷达图像进行方位模糊重建，确定出所述星载合成孔径雷达图像的模糊图像，包括：
对所述星载合成孔径雷达图像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η)按列做快速逆傅里叶变换，得到所述正频率方位模糊图像。
4.根据权利要求1所述的方法，其特征在于，所述模糊图像包括负频率方位模糊图像；
所述阈值条件包括：
其中，η为所述方位向时间；fp为脉冲重复频率；fr为多普勒调频率；和/或
所述模糊图像包括正频率方位模糊图像；
所述阈值条件包括：
其中，η为所述方位向时间；Na为图像矩阵总行数；fp为脉冲重复频率；fr为多普勒调频率。
5.根据权利要求1所述的方法，其特征在于，所述当所述星载合成孔径雷达图像中的数据点满足所述预定抑制条件时，利用所述模糊图像，对所述星载合成孔径雷达图像中的数据点进行模糊抑制，包括：
将所述星载合成孔径雷达图像中、满足所述预定抑制条件的数据点的绝对值，与所述模糊图像中对应的数据点的绝对值做差值，得到经模糊抑制后的图像。
6.一种星载合成孔径雷达图像的模糊抑制装置，其特征在于，包括：
方位模糊重建模块，用于对星载合成孔径雷达图像进行方位模糊重建，确定出所述星载合成孔径雷达图像的模糊图像；
判断模块，用于根据与所述星载合成孔径雷达图像中的数据点对应的所述模糊图像的数据点，判断所述星载合成孔径雷达图像中的数据点是否满足预定抑制条件，其中，所述判断模块包括：
第一判断子模块，用于判断所述模糊图像的数据点的绝对值是否大于所述星载合成孔径雷达图像中对应的数据点的绝对值的预定倍数，其中，所述预定倍数大于0且小于1；
第二判断子模块，用于判断方位向时间是否满足预设的阈值条件；
确定子模块，用于当所述模糊图像的数据点的绝对值大于所述星载合成孔径雷达图像中对应的数据点的绝对值的预定倍数、且所述方位向时间满足所述阈值条件时，确定所述星载合成孔径雷达图像中的数据点满足所述预定抑制条件；
模糊抑制模块，用于当所述星载合成孔径雷达图像中的数据点满足所述预定抑制条件时，利用所述模糊图像，对所述星载合成孔径雷达图像中、满足所述预定抑制条件的数据点进行模糊抑制。
7.根据权利要求6所述的装置，其特征在于，所述模糊图像包括负频率方位模糊图像；
所述方位模糊重建模块用于：
对所述星载合成孔径雷达图像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η)按列做快速逆傅里叶变换，得到所述负频率方位模糊图像。
8.根据权利要求6或7所述的装置，其特征在于，所述模糊图像包括正频率方位模糊图像；
所述方位模糊重建模块用于：
对所述星载合成孔径雷达图像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η)按列做快速逆傅里叶变换，得到所述正频率方位模糊图像。
9.根据权利要求6所述的装置，其特征在于，所述模糊图像包括负频率方位模糊图像；
所述阈值条件包括：
其中，η为所述方位向时间；fp为脉冲重复频率；fr为多普勒调频率；和/或
所述模糊图像包括正频率方位模糊图像；
所述阈值条件包括：
其中，η为所述方位向时间；Na为图像矩阵总行数；fp为脉冲重复频率；fr为多普勒调频率。
10.根据权利要求6所述的装置，其特征在于，所述模糊抑制模块用于将所述星载合成孔径雷达图像中、满足所述预定抑制条件的数据点的绝对值，与所述模糊图像中对应的数据点的绝对值做差值，得到经模糊抑制后的图像。
CN201710134700.3A 2017-03-08 2017-03-08 星载合成孔径雷达图像的模糊抑制方法和装置 Active CN107064895B (zh)

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 CN107064895A (zh) 2017-08-18
CN107064895B true 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)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107976674A (zh) * 2017-10-31 2018-05-01 西安空间无线电技术研究所 一种基于最优方位模糊的双波束sar的prf选择方法及系统
CN108389166B (zh) * 2017-11-21 2021-08-13 北京航空航天大学 模糊覆盖区域处理方法、装置、设备及计算机可读存储介质
CN112462366B (zh) * 2020-10-10 2023-08-08 深圳大学 一种sar数据点可视化方法、智能终端及存储介质

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Accurate Reconstruction and Suppression for Azimuth Ambiguities in Spaceborne stripmap SAR Images;Jie Chen et al.;《IEEE GEOSCIENCE AND REMOTE SENSING LETTERS》;20170131;第14卷(第1期);102-106
一种改进的SAR图像方位模糊抑制方法;吴正升 等;《测绘通报》;20141231;149-152,162
星载SAR影像上船舶方位向模糊去除算法;邴磊 等;《中国图象图形学报》;20160731;第21卷(第7期);951-958

Also Published As

Publication number Publication date
CN107064895A (zh) 2017-08-18

Similar Documents

Publication Publication Date Title
Yatawatta et al. Initial deep LOFAR observations of epoch of reionization windows-I. The north celestial pole
Romeiser et al. First analysis of TerraSAR-X along-track InSAR-derived current fields
CN108051809B (zh) 基于Radon变换的运动目标成像方法、装置及电子设备
Guarnieri Adaptive removal of azimuth ambiguities in SAR images
CN107064895B (zh) 星载合成孔径雷达图像的模糊抑制方法和装置
EP3012658B1 (en) Method and device for implementing sar imaging
Kapahi et al. The Molonglo reference catalog 1 Jy radio source survey. II. Radio structures of galaxy identifications
US20080297405A1 (en) Synthetic Aperture focusing techniques
CN107942327A (zh) 基于脉冲相位编码的单通道hrws‑sar成像方法
Fan et al. Estimation of wind direction in tropical cyclones using C-band dual-polarization synthetic aperture radar
Cotton Fourier Plane Image Combination by Feathering
CN106646529A (zh) 一种基于多波束优选的gnss天线阵抗干扰方法
DE102015122420A1 (de) Sendeanordnung zum Erzeugen eines für eine Lokalisierung geeigneten Signalmusters und Empfangsanordnung zum Durchführen einer Lokalisierung
CN105676190A (zh) 一种校正合成孔径雷达回波数据的方法和装置
CN103576153A (zh) 一种方位向多波束合成孔径雷达及其实现方法和装置
Camps et al. An imaging algorithm for synthetic aperture interferometric radiometers with built-in RFI mitigation
CN108375779B (zh) 一种基于Kalman滤波器的导航接收机干扰抑制方法
Corbella et al. Reduction of secondary lobes in aperture synthesis radiometry
Zheng Adaptive antenna array processing for GPS receivers.
Campbell et al. Arecibo radar maps of Venus from 1988 to 2020
Raout et al. Passive bistatic noise radar using DVB-T signals
CN103323820A (zh) 一种抑制频率步进sar一维距离图像栅瓣的方法
Zhang et al. A new SAR–GMTI high-accuracy focusing and relocation method using instantaneous interferometry
Zhu et al. Super-resolution for 4-D SAR tomography via compressive sensing
Hellbourg et al. Spatial filtering experiment with the askap beta array

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