CN117289276A - Sar图像的处理方法、装置、电子设备和存储介质 - Google Patents

Sar图像的处理方法、装置、电子设备和存储介质 Download PDF

Info

Publication number
CN117289276A
CN117289276A CN202311576620.5A CN202311576620A CN117289276A CN 117289276 A CN117289276 A CN 117289276A CN 202311576620 A CN202311576620 A CN 202311576620A CN 117289276 A CN117289276 A CN 117289276A
Authority
CN
China
Prior art keywords
envelope
part data
imaginary
distance
data
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
Application number
CN202311576620.5A
Other languages
English (en)
Other versions
CN117289276B (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.)
Aerospace Hongtu Information Technology Co Ltd
Original Assignee
Aerospace Hongtu Information Technology Co Ltd
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 Aerospace Hongtu Information Technology Co Ltd filed Critical Aerospace Hongtu Information Technology Co Ltd
Priority to CN202311576620.5A priority Critical patent/CN117289276B/zh
Publication of CN117289276A publication Critical patent/CN117289276A/zh
Application granted granted Critical
Publication of CN117289276B publication Critical patent/CN117289276B/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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00—Geometric image transformations in the plane of the image
    • G06T3/40—Scaling of whole images or parts thereof, e.g. expanding or contracting
    • G06T3/4053—Scaling of whole images or parts thereof, e.g. expanding or contracting based on super-resolution, i.e. the output image resolution being higher than the sensor resolution

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Theoretical Computer Science (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明提供了一种SAR图像的处理方法、装置、电子设备和存储介质，属于图像处理的技术领域，该处理方法中，是先对原始虚部数据矩阵和原始实部数据矩阵进行包络对消处理后，再对包络对消后的虚部数据矩阵和包络对消后的实部数据矩阵进行空间变迹处理，最后再对空间变迹后的虚部数据矩阵和空间变迹后的实部数据矩阵进行合并，得到旁瓣抑制和图像增强后的单视复SAR图像，由于包络对消处理能够消除原始虚部数据矩阵和原始实部数据矩阵中的低频分量，进而消除这些低频分量所导致的回波不在0值附近的值，后续再对包络对消后的虚部数据矩阵和包络对消后的实部数据矩阵进行逐行逐列空间变迹处理时，准确性好，单视复SAR图像的旁瓣抑制效果好，分辨率高。

Description

SAR图像的处理方法、装置、电子设备和存储介质
技术领域
本发明涉及图像处理的技术领域，尤其是涉及一种SAR图像的处理方法、装置、电子设备和存储介质。
背景技术
空间变迹算法（Spatially Variant Apodization，SVA）是一种不需要先验知识，基于余弦类频域加权的非线性加权方法，广泛应用于合成孔径雷达（SAR，SyntheticAperture Radar）的图像滤波、旁瓣抑制及高分辨率成像。目前，针对空间变迹法的改进主要集中在选取滤波器窗函数类型及滤波器的点数。
传统空间变迹算法在对SAR数据进行实际处理时，均是直接在时域进行，对单个距离向和方位向数据的回波数值分别进行处理，因此，对整个波形的具体取值要求比较高，具体来说要求点函数脉冲响应函数的旁瓣都应该在0值上下周期循环。然而，当原始SAR数据中存在偏差及干扰时，利用传统的空间变迹算法处理时，往往难以取得理想的效果。
综上，现有的SAR图像的处理方法存在旁瓣抑制效果差、处理后的SAR图像分辨率不高的技术问题。
发明内容
有鉴于此，本发明的目的在于提供一种SAR图像的处理方法、装置、电子设备和存储介质，以缓解现有的SAR图像的处理方法旁瓣抑制效果差、处理后的SAR图像分辨率不高的技术问题。
第一方面，本发明实施例提供了一种SAR图像的处理方法，包括：
获取单视复SAR图像，并对所述单视复SAR图像进行虚部数据和实部数据的提取，得到所述单视复SAR图像对应的原始虚部数据矩阵和原始实部数据矩阵；
对所述原始虚部数据矩阵进行逐行逐列包络对消处理，以消除所述原始虚部数据矩阵中的低频分量，进而得到包络对消后的虚部数据矩阵；
对所述原始实部数据矩阵进行逐行逐列包络对消处理，以消除所述原始实部数据矩阵中的低频分量，进而得到包络对消后的实部数据矩阵；
对所述包络对消后的虚部数据矩阵进行逐行逐列空间变迹处理，得到空间变迹后的虚部数据矩阵；
对所述包络对消后的实部数据矩阵进行逐行逐列空间变迹处理，得到空间变迹后的实部数据矩阵；
对所述空间变迹后的虚部数据矩阵和所述空间变迹后的实部数据矩阵进行合并，得到旁瓣抑制和图像增强后的单视复SAR图像。
进一步的，对所述原始虚部数据矩阵进行逐行逐列包络对消处理，包括：
对所述原始虚部数据矩阵的每行原始虚部数据中的每个原始虚部数据分别进行上包络提取和下包络提取，并根据提取得到的每个所述原始虚部数据的上包络和下包络对对应的原始虚部数据进行距离向包络对消处理，得到距离向包络对消后的虚部数据矩阵；
对所述距离向包络对消后的虚部数据矩阵的每列距离向包络对消后的虚部数据中的每个距离向包络对消后的虚部数据分别进行上包络提取和下包络提取，并根据提取得到的每个所述距离向包络对消后的虚部数据的上包络和下包络对对应的距离向包络对消后的虚部数据进行方位向包络对消处理，得到包络对消后的虚部数据矩阵；
对所述原始实部数据矩阵进行逐行逐列包络对消处理，包括：
对所述原始实部数据矩阵的每行原始实部数据中的每个原始实部数据分别进行上包络提取和下包络提取，并根据提取得到的每个所述原始实部数据的上包络和下包络对对应的原始实部数据进行距离向包络对消处理，得到距离向包络对消后的实部数据矩阵；
对所述距离向包络对消后的实部数据矩阵的每列距离向包络对消后的实部数据中的每个距离向包络对消后的实部数据分别进行上包络提取和下包络提取，并根据提取得到的每个所述距离向包络对消后的实部数据的上包络和下包络对对应的距离向包络对消后的实部数据进行方位向包络对消处理，得到包络对消后的实部数据矩阵。
进一步的，根据提取得到的每个所述原始虚部数据的上包络和下包络对对应的原始虚部数据进行距离向包络对消处理，包括：
若提取得到的第一原始虚部数据的上包络和下包络的符号相同，则将所述第一原始虚部数据与对应的上包络作差，并将得到的差值与对应的上包络和下包络的平均值作和，进而实现对所述第一原始虚部数据的距离向包络对消处理；
若提取得到的第二原始虚部数据的上包络和下包络的符号不同，则不对所述第二原始虚部数据进行距离向包络对消处理；
根据提取得到的每个所述距离向包络对消后的虚部数据的上包络和下包络对对应的距离向包络对消后的虚部数据进行方位向包络对消处理，包括：
若提取得到的第一距离向包络对消后的虚部数据的上包络和下包络的符号相同，则将所述第一距离向包络对消后的虚部数据与对应的上包络作差，并将得到的差值与对应的上包络和下包络的平均值作和，进而实现对所述第一距离向包络对消后的虚部数据的方位向包络对消处理；
若提取得到的第二距离向包络对消后的虚部数据的上包络和下包络的符号不同，则不对所述第二距离向包络对消后的虚部数据进行方位向包络对消处理；
根据提取得到的每个所述原始实部数据的上包络和下包络对对应的原始实部数据进行距离向包络对消处理，包括：
若提取得到的第一原始实部数据的上包络和下包络的符号相同，则将所述第一原始实部数据与对应的上包络作差，并将得到的差值与对应的上包络和下包络的平均值作和，进而实现对所述第一原始实部数据的距离向包络对消处理；
若提取得到的第二原始实部数据的上包络和下包络的符号不同，则不对所述第二原始实部数据进行距离向包络对消处理；
根据提取得到的每个所述距离向包络对消后的实部数据的上包络和下包络对对应的距离向包络对消后的实部数据进行方位向包络对消处理，包括：
若提取得到的第一距离向包络对消后的实部数据的上包络和下包络的符号相同，则将所述第一距离向包络对消后的实部数据与对应的上包络作差，并将得到的差值与对应的上包络和下包络的平均值作和，进而实现对所述第一距离向包络对消后的实部数据的方位向包络对消处理；
若提取得到的第二距离向包络对消后的实部数据的上包络和下包络的符号不同，则不对所述第二距离向包络对消后的实部数据进行方位向包络对消处理。
进一步的，对所述包络对消后的虚部数据矩阵进行逐行逐列空间变迹处理，包括：
对所述包络对消后的虚部数据矩阵的每行包络对消后的虚部数据，遍历其中的每个包络对消后的虚部数据，计算第一权重下的当前包络对消后的虚部数据的第一加权值和第二权重下的所述当前包络对消后的虚部数据的第二加权值；
若所述当前包络对消后的虚部数据的第一加权值与所述当前包络对消后的虚部数据的第二加权值的符号不同，则将所述当前包络对消后的虚部数据置为0；
若所述当前包络对消后的虚部数据的第一加权值与所述当前包络对消后的虚部数据的第二加权值的符号相同，则将所述当前包络对消后的虚部数据置为第一目标值，进而得到距离向空间变迹后的虚部数据矩阵，其中，所述第一目标值为基于所述当前包络对消后的虚部数据的第一加权值与所述当前包络对消后的虚部数据的第二加权值计算得到的；
对所述距离向空间变迹后的虚部数据矩阵的每列距离向空间变迹后的虚部数据，遍历其中的每个距离向空间变迹后的虚部数据，计算第一权重下的当前距离向空间变迹后的虚部数据的第一加权值和第二权重下的所述当前距离向空间变迹后的虚部数据的第二加权值；
若所述当前距离向空间变迹后的虚部数据的第一加权值与所述当前距离向空间变迹后的虚部数据的第二加权值的符号不同，则将所述当前距离向空间变迹后的虚部数据置为0；
若所述当前距离向空间变迹后的虚部数据的第一加权值与所述当前距离向空间变迹后的虚部数据的第二加权值的符号相同，则将所述当前距离向空间变迹后的虚部数据置为第二目标值，进而得到所述空间变迹后的虚部数据矩阵，其中，所述第二目标值为基于所述当前距离向空间变迹后的虚部数据的第一加权值与所述当前距离向空间变迹后的虚部数据的第二加权值计算得到的；
对所述包络对消后的实部数据矩阵进行逐行逐列空间变迹处理，包括：
对所述包络对消后的实部数据矩阵的每行包络对消后的实部数据，遍历其中的每个包络对消后的实部数据，计算第一权重下的当前包络对消后的实部数据的第一加权值和第二权重下的所述当前包络对消后的实部数据的第二加权值；
若所述当前包络对消后的实部数据的第一加权值与所述当前包络对消后的实部数据的第二加权值的符号不同，则将所述当前包络对消后的实部数据置为0；
若所述当前包络对消后的实部数据的第一加权值与所述当前包络对消后的实部数据的第二加权值的符号相同，则将所述当前包络对消后的实部数据置为第三目标值，进而得到距离向空间变迹后的实部数据矩阵，其中，所述第三目标值为基于所述当前包络对消后的实部数据的第一加权值与所述当前包络对消后的实部数据的第二加权值计算得到的；
对所述距离向空间变迹后的实部数据矩阵的每列距离向空间变迹后的实部数据，遍历其中的每个距离向空间变迹后的实部数据，计算第一权重下的当前距离向空间变迹后的实部数据的第一加权值和第二权重下的所述当前距离向空间变迹后的实部数据的第二加权值；
若所述当前距离向空间变迹后的实部数据的第一加权值与所述当前距离向空间变迹后的实部数据的第二加权值的符号不同，则将所述当前距离向空间变迹后的实部数据置为0；
若所述当前距离向空间变迹后的实部数据的第一加权值与所述当前距离向空间变迹后的实部数据的第二加权值的符号相同，则将所述当前距离向空间变迹后的实部数据置为第四目标值，进而得到所述空间变迹后的实部数据矩阵，其中，所述第四目标值为基于所述当前距离向空间变迹后的实部数据的第一加权值与所述当前距离向空间变迹后的实部数据的第二加权值计算得到的。
进一步的，计算第一权重下的当前包络对消后的虚部数据的第一加权值和第二权重下的所述当前包络对消后的虚部数据的第二加权值，包括：
根据加权值计算算式计算第一权重下的所述当前包络对消后的虚部数据的第一加权值/>和第二权重下的所述当前包络对消后的虚部数据的第二加权值/>，其中，/>表示所述当前包络对消后的虚部数据的加权值，/>表示所述当前包络对消后的虚部数据，/>表示所述第一加权值0或所述第二加权值/>，/>，/>，表示合成孔径雷达工作时的采样率，/>表示合成孔径雷达工作时的处理带宽，/>表示所述当前包络对消后的虚部数据之前的第k个包络对消后的虚部数据，/>表示所述当前包络对消后的虚部数据之后的第k个包络对消后的虚部数据，，/>为向下取整。
进一步的，所述第一目标值的计算算式如下：
其中，表示所述第一目标值，/>表示所述当前包络对消后的虚部数据的第一加权值，/>表示所述当前包络对消后的虚部数据的第二加权值，sign()函数表示取数值符号，min()函数表示取最小值，abs()函数表示取绝对值。
进一步的，对所述空间变迹后的虚部数据矩阵和所述空间变迹后的实部数据矩阵进行合并，包括：
采用实部和虚部加和的方式对所述空间变迹后的虚部数据矩阵和所述空间变迹后的实部数据矩阵进行合并，得到处理后的单视复数据；
将所述处理后的单视复数据可视化，得到所述旁瓣抑制和图像增强后的单视复SAR图像。
第二方面，本发明实施例还提供了一种SAR图像的处理装置，包括：
提取单元，用于获取单视复SAR图像，并对所述单视复SAR图像进行虚部数据和实部数据的提取，得到所述单视复SAR图像对应的原始虚部数据矩阵和原始实部数据矩阵；
虚部包络对消处理单元，用于对所述原始虚部数据矩阵进行逐行逐列包络对消处理，以消除所述原始虚部数据矩阵中的低频分量，进而得到包络对消后的虚部数据矩阵；
实部包络对消处理单元，用于对所述原始实部数据矩阵进行逐行逐列包络对消处理，以消除所述原始实部数据矩阵中的低频分量，进而得到包络对消后的实部数据矩阵；
虚部空间变迹处理单元，用于对所述包络对消后的虚部数据矩阵进行逐行逐列空间变迹处理，得到空间变迹后的虚部数据矩阵；
实部空间变迹处理单元，用于对所述包络对消后的实部数据矩阵进行逐行逐列空间变迹处理，得到空间变迹后的实部数据矩阵；
合并单元，用于对所述空间变迹后的虚部数据矩阵和所述空间变迹后的实部数据矩阵进行合并，得到旁瓣抑制和图像增强后的单视复SAR图像。
第三方面，本发明实施例还提供了一种电子设备，包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述第一方面任一项所述的方法的步骤。
第四方面，本发明实施例还提供了一种计算机可读存储介质，所述计算机可读存储介质存储有机器可运行指令，所述机器可运行指令在被处理器调用和运行时，所述机器可运行指令促使所述处理器运行上述第一方面任一项所述的方法。
在本发明实施例中，提供了一种SAR图像的处理方法，包括：获取单视复SAR图像，并对单视复SAR图像进行虚部数据和实部数据的提取，得到单视复SAR图像对应的原始虚部数据矩阵和原始实部数据矩阵；对原始虚部数据矩阵进行逐行逐列包络对消处理，以消除原始虚部数据矩阵中的低频分量，进而得到包络对消后的虚部数据矩阵；对原始实部数据矩阵进行逐行逐列包络对消处理，以消除原始实部数据矩阵中的低频分量，进而得到包络对消后的实部数据矩阵；对包络对消后的虚部数据矩阵进行逐行逐列空间变迹处理，得到空间变迹后的虚部数据矩阵；对包络对消后的实部数据矩阵进行逐行逐列空间变迹处理，得到空间变迹后的实部数据矩阵；对空间变迹后的虚部数据矩阵和空间变迹后的实部数据矩阵进行合并，得到旁瓣抑制和图像增强后的单视复SAR图像。通过上述描述可知，本发明的SAR图像的处理方法中，是先对原始虚部数据矩阵和原始实部数据矩阵进行逐行逐列包络对消处理后，再对包络对消后的虚部数据矩阵和包络对消后的实部数据矩阵进行逐行逐列空间变迹处理，最后再对空间变迹后的虚部数据矩阵和空间变迹后的实部数据矩阵进行合并，得到旁瓣抑制和图像增强后的单视复SAR图像，由于包络对消处理能够消除原始虚部数据矩阵和原始实部数据矩阵中的低频分量，进而消除这些低频分量所导致的回波不在0值附近的值，这样，后续再对包络对消后的虚部数据矩阵和包络对消后的实部数据矩阵进行逐行逐列空间变迹处理时，准确性好，单视复SAR图像的旁瓣抑制效果好，并且最终得到的旁瓣抑制和图像增强后的单视复SAR图像的分辨率高，缓解了现有的SAR图像的处理方法旁瓣抑制效果差、处理后的SAR图像分辨率不高的技术问题。
附图说明
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案，下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施方式，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的一种SAR图像的处理方法的流程图；
图2为本发明实施例提供的单道数据包络对消处理的示意图；
图3为本发明实施例提供的单道数据包络对消处理前后对比的示意图；
图4为本发明实施例提供的单视复SAR图像的处理结果对比示意图；
图5为本发明实施例提供的一种SAR图像的处理装置的示意图；
图6为本发明实施例提供的一种电子设备的示意图。
具体实施方式
下面将结合实施例对本发明的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
现有的SAR图像的处理方法旁瓣抑制效果差、处理后的SAR图像分辨率不高。
基于此，本发明的SAR图像的处理方法中，是先对原始虚部数据矩阵和原始实部数据矩阵进行逐行逐列包络对消处理后，再对包络对消后的虚部数据矩阵和包络对消后的实部数据矩阵进行逐行逐列空间变迹处理，最后再对空间变迹后的虚部数据矩阵和空间变迹后的实部数据矩阵进行合并，得到旁瓣抑制和图像增强后的单视复SAR图像，由于包络对消处理能够消除原始虚部数据矩阵和原始实部数据矩阵中的低频分量，进而消除这些低频分量所导致的回波不在0值附近的值，这样，后续再对包络对消后的虚部数据矩阵和包络对消后的实部数据矩阵进行逐行逐列空间变迹处理时，准确性好，单视复SAR图像的旁瓣抑制效果好，并且最终得到的旁瓣抑制和图像增强后的单视复SAR图像的分辨率高。
为便于对本实施例进行理解，首先对本发明实施例所公开的一种SAR图像的处理方法进行详细介绍。
实施例一：
根据本发明实施例，提供了一种SAR图像的处理方法的实施例，需要说明的是，在附图的流程图示出的步骤可以在诸如一组计算机可执行指令的计算机系统中执行，并且，虽然在流程图中示出了逻辑顺序，但是在某些情况下，可以以不同于此处的顺序执行所示出或描述的步骤。
图1是根据本发明实施例的一种SAR图像的处理方法的流程图，如图1所示，该方法包括如下步骤：
步骤S102，获取单视复SAR图像，并对单视复SAR图像进行虚部数据和实部数据的提取，得到单视复SAR图像对应的原始虚部数据矩阵和原始实部数据矩阵；
具体的，单视复SAR图像是SAR产品级别的通用说法。单视SAR图像是指只用一段合成孔径长度所成的SAR图像，没有和其它SAR图像进行叠加，通常所用的合成孔径长度比较长，这样方位分辨率比较高；多视是指将整个有效合成孔径长度分成多段分别对同一场景进行成像，然后将所得的图像求和叠加得到一幅SAR图像，好处是：可以提高SAR图像的信噪比，有效抑制斑点噪声，代价是降低方位分辨率。
单视复SAR图像为单视复图像数据（Single Look Complex，SLC，里面包含实部和虚部），即在单视复SAR图像中，每个像素点为复数，包含实部和虚部。在单视复SAR图像中，提取数据的实部矩阵和虚部矩阵，也就是对这个复数矩阵，取实部形成原始实部数据矩阵，取虚部形成原始虚部数据矩阵。
单视复SAR图像通常为tiff格式，获取得到对应的原始SAR数据矩阵data，其维度为[row,col,2]。row为SAR图像方位向宽度；col为SAR图像的距离向宽度，对原始SAR数据矩阵data，分别取虚部和实部部分，得到原始虚部数据矩阵data_img和原始实部数据矩阵data_real，其大小均为[row，col]。
步骤S104，对原始虚部数据矩阵进行逐行逐列包络对消处理，以消除原始虚部数据矩阵中的低频分量，进而得到包络对消后的虚部数据矩阵；
具体的，包络对消处理的目的在于消除原始虚部数据矩阵中引入的低频分量导致回波不在0值附近。
步骤S106，对原始实部数据矩阵进行逐行逐列包络对消处理，以消除原始实部数据矩阵中的低频分量，进而得到包络对消后的实部数据矩阵；
具体的，包络对消处理的目的在于消除原始实部数据矩阵中引入的低频分量导致回波不在0值附近。
步骤S108，对包络对消后的虚部数据矩阵进行逐行逐列空间变迹处理，得到空间变迹后的虚部数据矩阵；
具体的，本发明的空间变迹处理不仅能够实现旁瓣抑制，还能在一定程度上起到图像增强的效果。
步骤S110，对包络对消后的实部数据矩阵进行逐行逐列空间变迹处理，得到空间变迹后的实部数据矩阵；
步骤S112，对空间变迹后的虚部数据矩阵和空间变迹后的实部数据矩阵进行合并，得到旁瓣抑制和图像增强后的单视复SAR图像。
在本发明实施例中，提供了一种SAR图像的处理方法、装置、电子设备和存储介质，包括：获取单视复SAR图像，并对单视复SAR图像进行虚部数据和实部数据的提取，得到单视复SAR图像对应的原始虚部数据矩阵和原始实部数据矩阵；对原始虚部数据矩阵进行逐行逐列包络对消处理，以消除原始虚部数据矩阵中的低频分量，进而得到包络对消后的虚部数据矩阵；对原始实部数据矩阵进行逐行逐列包络对消处理，以消除原始实部数据矩阵中的低频分量，进而得到包络对消后的实部数据矩阵；对包络对消后的虚部数据矩阵进行逐行逐列空间变迹处理，得到空间变迹后的虚部数据矩阵；对包络对消后的实部数据矩阵进行逐行逐列空间变迹处理，得到空间变迹后的实部数据矩阵；对空间变迹后的虚部数据矩阵和空间变迹后的实部数据矩阵进行合并，得到旁瓣抑制和图像增强后的单视复SAR图像。通过上述描述可知，本发明的SAR图像的处理方法中，是先对原始虚部数据矩阵和原始实部数据矩阵进行逐行逐列包络对消处理后，再对包络对消后的虚部数据矩阵和包络对消后的实部数据矩阵进行逐行逐列空间变迹处理，最后再对空间变迹后的虚部数据矩阵和空间变迹后的实部数据矩阵进行合并，得到旁瓣抑制和图像增强后的单视复SAR图像，由于包络对消处理能够消除原始虚部数据矩阵和原始实部数据矩阵中的低频分量，进而消除这些低频分量所导致的回波不在0值附近的值，这样，后续再对包络对消后的虚部数据矩阵和包络对消后的实部数据矩阵进行逐行逐列空间变迹处理时，准确性好，单视复SAR图像的旁瓣抑制效果好，并且最终得到的旁瓣抑制和图像增强后的单视复SAR图像的分辨率高，缓解了现有的SAR图像的处理方法旁瓣抑制效果差、处理后的SAR图像分辨率不高的技术问题。
上述内容对本发明的SAR图像的处理方法进行了简要介绍，下面对其中涉及到的具体内容进行详细处理。
在本发明的一个可选实施例中，对原始虚部数据矩阵进行逐行逐列包络对消处理，具体包括如下步骤：
（1）对原始虚部数据矩阵的每行原始虚部数据中的每个原始虚部数据分别进行上包络提取和下包络提取，并根据提取得到的每个原始虚部数据的上包络和下包络对对应的原始虚部数据进行距离向包络对消处理，得到距离向包络对消后的虚部数据矩阵；
具体的，在原始虚部数据矩阵data_img中，先提取第一行原始虚部数据data_img_one，如图2中的实线。
对第一行（行即为距离向）原始虚部数据data_img_one中的每个原始虚部数据分别进行上包络提取（如图2中的点线）和下包络提取（如图2中的虚线），该过程可以是并行提取各个原始虚部数据的上包络和下包络，本发明对其实现方式不进行具体限制。
进而根据提取得到的每个原始虚部数据的上包络和下包络对对应的原始虚部数据进行距离向包络对消处理，具体过程包括：
（11）若提取得到的第一原始虚部数据的上包络和下包络的符号相同，则将第一原始虚部数据与对应的上包络作差，并将得到的差值与对应的上包络和下包络的平均值作和，进而实现对第一原始虚部数据的距离向包络对消处理；
上述符号相同是指上包络和下包络均为正数或负数。
（12）若提取得到的第二原始虚部数据的上包络和下包络的符号不同，则不对第二原始虚部数据进行距离向包络对消处理；
上述符号不同是指上包络和下包络一个为正数，一个为负数，其对应的物理含义为：一个标准的点脉冲响应回波，其波峰和波谷应分别位于0值的两侧，正常应该一正一负交替出现。如果同时为正或同时为负，说明该点脉冲响应回波存在低频分量、偏置或干扰，需要进行抑制。
遍历完第一行原始虚部数据data_img_one中的每个原始虚部数据，得到第一行距离向包络对消后的虚部数据data_img_one_pro，如图3中的虚线。
然后遍历所有距离向（即所有行），得到距离向包络对消后的虚部数据矩阵data_img_row。
（2）对距离向包络对消后的虚部数据矩阵的每列距离向包络对消后的虚部数据中的每个距离向包络对消后的虚部数据分别进行上包络提取和下包络提取，并根据提取得到的每个距离向包络对消后的虚部数据的上包络和下包络对对应的距离向包络对消后的虚部数据进行方位向包络对消处理，得到包络对消后的虚部数据矩阵；
具体的，因为单视复SAR图像分为距离向和方位向，距离向和方位向上都会存在干扰，对应数据矩阵就是行和列。包络对消处理是分别对每一行进行一遍处理，然后对每一列再进行一遍处理。
在距离向包络对消后的虚部数据矩阵data_img_row上，在依次提取每一个方位向（列）数据data_img_row_one，即每列距离向包络对消后的虚部数据，对其中的每个距离向包络对消后的虚部数据分别进行上包络提取和下包络提取，并根据提取得到的每个距离向包络对消后的虚部数据的上包络和下包络对对应的距离向包络对消后的虚部数据进行方位向包络对消处理（与上述（1）的过程相似），得到包络对消后的虚部数据矩阵data_img_row_col。
其中，根据提取得到的每个距离向包络对消后的虚部数据的上包络和下包络对对应的距离向包络对消后的虚部数据进行方位向包络对消处理，具体包括如下步骤：
（21）若提取得到的第一距离向包络对消后的虚部数据的上包络和下包络的符号相同，则将第一距离向包络对消后的虚部数据与对应的上包络作差，并将得到的差值与对应的上包络和下包络的平均值作和，进而实现对第一距离向包络对消后的虚部数据的方位向包络对消处理；
（22）若提取得到的第二距离向包络对消后的虚部数据的上包络和下包络的符号不同，则不对第二距离向包络对消后的虚部数据进行方位向包络对消处理；
对原始实部数据矩阵进行逐行逐列包络对消处理，具体包括如下步骤：
（3）对原始实部数据矩阵的每行原始实部数据中的每个原始实部数据分别进行上包络提取和下包络提取，并根据提取得到的每个原始实部数据的上包络和下包络对对应的原始实部数据进行距离向包络对消处理，得到距离向包络对消后的实部数据矩阵；
其中，根据提取得到的每个原始实部数据的上包络和下包络对对应的原始实部数据进行距离向包络对消处理，具体包括如下步骤：
（31）若提取得到的第一原始实部数据的上包络和下包络的符号相同，则将第一原始实部数据与对应的上包络作差，并将得到的差值与对应的上包络和下包络的平均值作和，进而实现对第一原始实部数据的距离向包络对消处理；
（32）若提取得到的第二原始实部数据的上包络和下包络的符号不同，则不对第二原始实部数据进行距离向包络对消处理；
（4）对距离向包络对消后的实部数据矩阵的每列距离向包络对消后的实部数据中的每个距离向包络对消后的实部数据分别进行上包络提取和下包络提取，并根据提取得到的每个距离向包络对消后的实部数据的上包络和下包络对对应的距离向包络对消后的实部数据进行方位向包络对消处理，得到包络对消后的实部数据矩阵。
其中，根据提取得到的每个距离向包络对消后的实部数据的上包络和下包络对对应的距离向包络对消后的实部数据进行方位向包络对消处理，具体包括如下步骤：
（41）若提取得到的第一距离向包络对消后的实部数据的上包络和下包络的符号相同，则将第一距离向包络对消后的实部数据与对应的上包络作差，并将得到的差值与对应的上包络和下包络的平均值作和，进而实现对第一距离向包络对消后的实部数据的方位向包络对消处理；
（42）若提取得到的第二距离向包络对消后的实部数据的上包络和下包络的符号不同，则不对第二距离向包络对消后的实部数据进行方位向包络对消处理。
具体的，对原始实部数据矩阵进行逐行逐列包络对消处理的过程与对原始虚部数据矩阵进行逐行逐列包络对消处理的过程相同，因此不再展开描述。
本发明中，对原始实部数据矩阵data_real进行包络对消处理后，得到包络对消后的实部数据矩阵data_real_row_col。
在本发明的一个可选实施例中，对包络对消后的虚部数据矩阵进行逐行逐列空间变迹处理，具体包括如下步骤：
（1）对包络对消后的虚部数据矩阵的每行包络对消后的虚部数据，遍历其中的每个包络对消后的虚部数据，计算第一权重下的当前包络对消后的虚部数据的第一加权值和第二权重下的当前包络对消后的虚部数据的第二加权值；
具体的，计算第一权重下的当前包络对消后的虚部数据的第一加权值和第二权重下的当前包络对消后的虚部数据的第二加权值，包括：
根据加权值计算算式计算第一权重下的当前包络对消后的虚部数据的第一加权值/>和第二权重下的当前包络对消后的虚部数据的第二加权值/>，其中，/>表示当前包络对消后的虚部数据的加权值，/>表示当前包络对消后的虚部数据，/>表示第一加权值0或第二加权值/>，/>，/>，表示合成孔径雷达工作时的采样率，/>表示合成孔径雷达工作时的处理带宽，/>表示当前包络对消后的虚部数据之前的第k个包络对消后的虚部数据，/>表示当前包络对消后的虚部数据之后的第k个包络对消后的虚部数据，/>，为向下取整。
（2）若当前包络对消后的虚部数据的第一加权值与当前包络对消后的虚部数据的第二加权值的符号不同，则将当前包络对消后的虚部数据置为0；
具体的，与/>的数值符号不同（如一个为正数，另一个为负数），则/>=0。也就是符号不同时，确定当前包络对消后的虚部数据位于旁瓣处，需要置为0。
（3）若当前包络对消后的虚部数据的第一加权值与当前包络对消后的虚部数据的第二加权值的符号相同，则将当前包络对消后的虚部数据置为第一目标值，进而得到距离向空间变迹后的虚部数据矩阵，其中，第一目标值为基于当前包络对消后的虚部数据的第一加权值与当前包络对消后的虚部数据的第二加权值计算得到的；
具体的，与/>的数值符号相同（两个都为正数或负数），将当前包络对消后的虚部数据置为第一目标值。
其中，第一目标值的计算算式如下：
，表示第一目标值，/>表示当前包络对消后的虚部数据的第一加权值，表示当前包络对消后的虚部数据的第二加权值，sign()函数表示取数值符号，min()函数表示取最小值，abs()函数表示取绝对值。
（4）对距离向空间变迹后的虚部数据矩阵的每列距离向空间变迹后的虚部数据，遍历其中的每个距离向空间变迹后的虚部数据，计算第一权重下的当前距离向空间变迹后的虚部数据的第一加权值和第二权重下的当前距离向空间变迹后的虚部数据的第二加权值；
该过程与上述计算第一权重下的当前包络对消后的虚部数据的第一加权值和第二权重下的当前包络对消后的虚部数据的第二加权值的过程相同，在此不再赘述。
（5）若当前距离向空间变迹后的虚部数据的第一加权值与当前距离向空间变迹后的虚部数据的第二加权值的符号不同，则将当前距离向空间变迹后的虚部数据置为0；
（6）若当前距离向空间变迹后的虚部数据的第一加权值与当前距离向空间变迹后的虚部数据的第二加权值的符号相同，则将当前距离向空间变迹后的虚部数据置为第二目标值，进而得到空间变迹后的虚部数据矩阵，其中，第二目标值为基于当前距离向空间变迹后的虚部数据的第一加权值与当前距离向空间变迹后的虚部数据的第二加权值计算得到的；
该过程与上述计算第一目标值的过程相同，在此不再赘述。
上述对包络对消后的虚部数据矩阵data_img_row_col进行逐行逐列空间变迹处理后，得到data_img_row_col_sva。
对包络对消后的实部数据矩阵进行逐行逐列空间变迹处理，包括：
（1）对包络对消后的实部数据矩阵的每行包络对消后的实部数据，遍历其中的每个包络对消后的实部数据，计算第一权重下的当前包络对消后的实部数据的第一加权值和第二权重下的当前包络对消后的实部数据的第二加权值；
（2）若当前包络对消后的实部数据的第一加权值与当前包络对消后的实部数据的第二加权值的符号不同，则将当前包络对消后的实部数据置为0；
（3）若当前包络对消后的实部数据的第一加权值与当前包络对消后的实部数据的第二加权值的符号相同，则将当前包络对消后的实部数据置为第三目标值，进而得到距离向空间变迹后的实部数据矩阵，其中，第三目标值为基于当前包络对消后的实部数据的第一加权值与当前包络对消后的实部数据的第二加权值计算得到的；
（4）对距离向空间变迹后的实部数据矩阵的每列距离向空间变迹后的实部数据，遍历其中的每个距离向空间变迹后的实部数据，计算第一权重下的当前距离向空间变迹后的实部数据的第一加权值和第二权重下的当前距离向空间变迹后的实部数据的第二加权值；
（5）若当前距离向空间变迹后的实部数据的第一加权值与当前距离向空间变迹后的实部数据的第二加权值的符号不同，则将当前距离向空间变迹后的实部数据置为0；
（6）若当前距离向空间变迹后的实部数据的第一加权值与当前距离向空间变迹后的实部数据的第二加权值的符号相同，则将当前距离向空间变迹后的实部数据置为第四目标值，进而得到空间变迹后的实部数据矩阵，其中，第四目标值为基于当前距离向空间变迹后的实部数据的第一加权值与当前距离向空间变迹后的实部数据的第二加权值计算得到的。
该对包络对消后的实部数据矩阵进行逐行逐列空间变迹处理的过程与上述对包络对消后的虚部数据矩阵进行逐行逐列空间变迹处理的过程相同，也不再展开描述。
该对包络对消后的实部数据矩阵data_real_row_col进行逐行逐列空间变迹处理后，得到data_real_row_col_sva。
在本发明的一个可选实施例中，对空间变迹后的虚部数据矩阵和空间变迹后的实部数据矩阵进行合并，具体包括如下步骤：
（1）采用实部和虚部加和的方式对空间变迹后的虚部数据矩阵和空间变迹后的实部数据矩阵进行合并，得到处理后的单视复数据；
具体的，data_final = data_img_row_col_sva×j+data_real_row_col_sva，其中，j为虚数，等于，data_final表示处理后的单视复数据矩阵，data_img_row_col_sva表示空间变迹后的虚部数据矩阵，data_real_row_col_sva表示空间变迹后的实部数据矩阵。
（2）将处理后的单视复数据可视化，得到旁瓣抑制和图像增强后的单视复SAR图像。
本发明的SAR图像的处理方法在空间变迹处理之前引入了包络对消处理，解决了原始虚部数据矩阵和原始实部数据矩阵存在低频分量导致直接利用空间变迹处理效果不理想的问题。
图4为对比的处理结果图，（a）为原始的单视复SAR图像，（b）为传统的空间变迹法的处理结果，（c）为加入包络对消后的空间变迹处理结果。结果表明，加入包络对消处理后，对存在旁瓣的单视复SAR图像有更加明显的抑制效果，且图像分辨率也有所提高。
实施例二：
本发明实施例还提供了一种SAR图像的处理装置，该SAR图像的处理装置主要用于执行本发明实施例一中所提供的SAR图像的处理方法，以下对本发明实施例提供的SAR图像的处理装置做具体介绍。
图5是根据本发明实施例的一种SAR图像的处理装置的示意图，如图5所示，该装置主要包括：提取单元10、虚部包络对消处理单元20、实部包络对消处理单元30、虚部空间变迹处理单元40、实部空间变迹处理单元50和合并单元60，其中：
提取单元，用于获取单视复SAR图像，并对单视复SAR图像进行虚部数据和实部数据的提取，得到单视复SAR图像对应的原始虚部数据矩阵和原始实部数据矩阵；
虚部包络对消处理单元，用于对原始虚部数据矩阵进行逐行逐列包络对消处理，以消除原始虚部数据矩阵中的低频分量，进而得到包络对消后的虚部数据矩阵；
实部包络对消处理单元，用于对原始实部数据矩阵进行逐行逐列包络对消处理，以消除原始实部数据矩阵中的低频分量，进而得到包络对消后的实部数据矩阵；
虚部空间变迹处理单元，用于对包络对消后的虚部数据矩阵进行逐行逐列空间变迹处理，得到空间变迹后的虚部数据矩阵；
实部空间变迹处理单元，用于对包络对消后的实部数据矩阵进行逐行逐列空间变迹处理，得到空间变迹后的实部数据矩阵；
合并单元，用于对空间变迹后的虚部数据矩阵和空间变迹后的实部数据矩阵进行合并，得到旁瓣抑制和图像增强后的单视复SAR图像。
在本发明实施例中，提供了一种SAR图像的处理装置，包括：获取单视复SAR图像，并对单视复SAR图像进行虚部数据和实部数据的提取，得到单视复SAR图像对应的原始虚部数据矩阵和原始实部数据矩阵；对原始虚部数据矩阵进行逐行逐列包络对消处理，以消除原始虚部数据矩阵中的低频分量，进而得到包络对消后的虚部数据矩阵；对原始实部数据矩阵进行逐行逐列包络对消处理，以消除原始实部数据矩阵中的低频分量，进而得到包络对消后的实部数据矩阵；对包络对消后的虚部数据矩阵进行逐行逐列空间变迹处理，得到空间变迹后的虚部数据矩阵；对包络对消后的实部数据矩阵进行逐行逐列空间变迹处理，得到空间变迹后的实部数据矩阵；对空间变迹后的虚部数据矩阵和空间变迹后的实部数据矩阵进行合并，得到旁瓣抑制和图像增强后的单视复SAR图像。通过上述描述可知，本发明的SAR图像的处理装置中，是先对原始虚部数据矩阵和原始实部数据矩阵进行逐行逐列包络对消处理后，再对包络对消后的虚部数据矩阵和包络对消后的实部数据矩阵进行逐行逐列空间变迹处理，最后再对空间变迹后的虚部数据矩阵和空间变迹后的实部数据矩阵进行合并，得到旁瓣抑制和图像增强后的单视复SAR图像，由于包络对消处理能够消除原始虚部数据矩阵和原始实部数据矩阵中的低频分量，进而消除这些低频分量所导致的回波不在0值附近的值，这样，后续再对包络对消后的虚部数据矩阵和包络对消后的实部数据矩阵进行逐行逐列空间变迹处理时，准确性好，单视复SAR图像的旁瓣抑制效果好，并且最终得到的旁瓣抑制和图像增强后的单视复SAR图像的分辨率高，缓解了现有的SAR图像的处理方法旁瓣抑制效果差、处理后的SAR图像分辨率不高的技术问题。
可选地，虚部包络对消处理单元还用于：对原始虚部数据矩阵的每行原始虚部数据中的每个原始虚部数据分别进行上包络提取和下包络提取，并根据提取得到的每个原始虚部数据的上包络和下包络对对应的原始虚部数据进行距离向包络对消处理，得到距离向包络对消后的虚部数据矩阵；对距离向包络对消后的虚部数据矩阵的每列距离向包络对消后的虚部数据中的每个距离向包络对消后的虚部数据分别进行上包络提取和下包络提取，并根据提取得到的每个距离向包络对消后的虚部数据的上包络和下包络对对应的距离向包络对消后的虚部数据进行方位向包络对消处理，得到包络对消后的虚部数据矩阵；
实部包络对消处理单元还用于：对原始实部数据矩阵的每行原始实部数据中的每个原始实部数据分别进行上包络提取和下包络提取，并根据提取得到的每个原始实部数据的上包络和下包络对对应的原始实部数据进行距离向包络对消处理，得到距离向包络对消后的实部数据矩阵；对距离向包络对消后的实部数据矩阵的每列距离向包络对消后的实部数据中的每个距离向包络对消后的实部数据分别进行上包络提取和下包络提取，并根据提取得到的每个距离向包络对消后的实部数据的上包络和下包络对对应的距离向包络对消后的实部数据进行方位向包络对消处理，得到包络对消后的实部数据矩阵。
可选地，虚部包络对消处理单元还用于：若提取得到的第一原始虚部数据的上包络和下包络的符号相同，则将第一原始虚部数据与对应的上包络作差，并将得到的差值与对应的上包络和下包络的平均值作和，进而实现对第一原始虚部数据的距离向包络对消处理；若提取得到的第二原始虚部数据的上包络和下包络的符号不同，则不对第二原始虚部数据进行距离向包络对消处理；若提取得到的第一距离向包络对消后的虚部数据的上包络和下包络的符号相同，则将第一距离向包络对消后的虚部数据与对应的上包络作差，并将得到的差值与对应的上包络和下包络的平均值作和，进而实现对第一距离向包络对消后的虚部数据的方位向包络对消处理；若提取得到的第二距离向包络对消后的虚部数据的上包络和下包络的符号不同，则不对第二距离向包络对消后的虚部数据进行方位向包络对消处理；
实部包络对消处理单元还用于：若提取得到的第一原始实部数据的上包络和下包络的符号相同，则将第一原始实部数据与对应的上包络作差，并将得到的差值与对应的上包络和下包络的平均值作和，进而实现对第一原始实部数据的距离向包络对消处理；若提取得到的第二原始实部数据的上包络和下包络的符号不同，则不对第二原始实部数据进行距离向包络对消处理；若提取得到的第一距离向包络对消后的实部数据的上包络和下包络的符号相同，则将第一距离向包络对消后的实部数据与对应的上包络作差，并将得到的差值与对应的上包络和下包络的平均值作和，进而实现对第一距离向包络对消后的实部数据的方位向包络对消处理；若提取得到的第二距离向包络对消后的实部数据的上包络和下包络的符号不同，则不对第二距离向包络对消后的实部数据进行方位向包络对消处理。
可选地，虚部空间变迹处理单元还用于：对包络对消后的虚部数据矩阵的每行包络对消后的虚部数据，遍历其中的每个包络对消后的虚部数据，计算第一权重下的当前包络对消后的虚部数据的第一加权值和第二权重下的当前包络对消后的虚部数据的第二加权值；若当前包络对消后的虚部数据的第一加权值与当前包络对消后的虚部数据的第二加权值的符号不同，则将当前包络对消后的虚部数据置为0；若当前包络对消后的虚部数据的第一加权值与当前包络对消后的虚部数据的第二加权值的符号相同，则将当前包络对消后的虚部数据置为第一目标值，进而得到距离向空间变迹后的虚部数据矩阵，其中，第一目标值为基于当前包络对消后的虚部数据的第一加权值与当前包络对消后的虚部数据的第二加权值计算得到的；对距离向空间变迹后的虚部数据矩阵的每列距离向空间变迹后的虚部数据，遍历其中的每个距离向空间变迹后的虚部数据，计算第一权重下的当前距离向空间变迹后的虚部数据的第一加权值和第二权重下的当前距离向空间变迹后的虚部数据的第二加权值；若当前距离向空间变迹后的虚部数据的第一加权值与当前距离向空间变迹后的虚部数据的第二加权值的符号不同，则将当前距离向空间变迹后的虚部数据置为0；若当前距离向空间变迹后的虚部数据的第一加权值与当前距离向空间变迹后的虚部数据的第二加权值的符号相同，则将当前距离向空间变迹后的虚部数据置为第二目标值，进而得到空间变迹后的虚部数据矩阵，其中，第二目标值为基于当前距离向空间变迹后的虚部数据的第一加权值与当前距离向空间变迹后的虚部数据的第二加权值计算得到的；
实部空间变迹处理单元还用于：对包络对消后的实部数据矩阵的每行包络对消后的实部数据，遍历其中的每个包络对消后的实部数据，计算第一权重下的当前包络对消后的实部数据的第一加权值和第二权重下的当前包络对消后的实部数据的第二加权值；若当前包络对消后的实部数据的第一加权值与当前包络对消后的实部数据的第二加权值的符号不同，则将当前包络对消后的实部数据置为0；若当前包络对消后的实部数据的第一加权值与当前包络对消后的实部数据的第二加权值的符号相同，则将当前包络对消后的实部数据置为第三目标值，进而得到距离向空间变迹后的实部数据矩阵，其中，第三目标值为基于当前包络对消后的实部数据的第一加权值与当前包络对消后的实部数据的第二加权值计算得到的；对距离向空间变迹后的实部数据矩阵的每列距离向空间变迹后的实部数据，遍历其中的每个距离向空间变迹后的实部数据，计算第一权重下的当前距离向空间变迹后的实部数据的第一加权值和第二权重下的当前距离向空间变迹后的实部数据的第二加权值；若当前距离向空间变迹后的实部数据的第一加权值与当前距离向空间变迹后的实部数据的第二加权值的符号不同，则将当前距离向空间变迹后的实部数据置为0；若当前距离向空间变迹后的实部数据的第一加权值与当前距离向空间变迹后的实部数据的第二加权值的符号相同，则将当前距离向空间变迹后的实部数据置为第四目标值，进而得到空间变迹后的实部数据矩阵，其中，第四目标值为基于当前距离向空间变迹后的实部数据的第一加权值与当前距离向空间变迹后的实部数据的第二加权值计算得到的。
可选地，虚部空间变迹处理单元还用于：根据加权值计算算式计算第一权重下的当前包络对消后的虚部数据的第一加权值/>和第二权重下的当前包络对消后的虚部数据的第二加权值/>，其中，/>表示当前包络对消后的虚部数据的加权值，/>表示当前包络对消后的虚部数据，/>表示第一加权值0或第二加权值/>，，/>，/>表示合成孔径雷达工作时的采样率，/>表示合成孔径雷达工作时的处理带宽，/>表示当前包络对消后的虚部数据之前的第k个包络对消后的虚部数据，/>表示当前包络对消后的虚部数据之后的第k个包络对消后的虚部数据，/>，为向下取整。
可选地，第一目标值的计算算式如下：，其中，表示第一目标值，/>表示当前包络对消后的虚部数据的第一加权值，表示当前包络对消后的虚部数据的第二加权值，sign()函数表示取数值符号，min()函数表示取最小值，abs()函数表示取绝对值。
可选地，合并单元还用于：采用实部和虚部加和的方式对空间变迹后的虚部数据矩阵和空间变迹后的实部数据矩阵进行合并，得到处理后的单视复数据；将处理后的单视复数据可视化，得到旁瓣抑制和图像增强后的单视复SAR图像。
本发明实施例所提供的装置，其实现原理及产生的技术效果和前述方法实施例相同，为简要描述，装置实施例部分未提及之处，可参考前述方法实施例中相应内容。
如图6所示，本申请实施例提供的一种电子设备600，包括：处理器601、存储器602和总线，所述存储器602存储有所述处理器601可执行的机器可读指令，当电子设备运行时，所述处理器601与所述存储器602之间通过总线通信，所述处理器601执行所述机器可读指令，以执行如上述SAR图像的处理方法的步骤。
具体地，上述存储器602和处理器601能够为通用的存储器和处理器，这里不做具体限定，当处理器601运行存储器602存储的计算机程序时，能够执行上述SAR图像的处理方法。
处理器601可能是一种集成电路芯片，具有信号的处理能力。在实现过程中，上述方法的各步骤可以通过处理器601中的硬件的集成逻辑电路或者软件形式的指令完成。上述的处理器601可以是通用处理器，包括中央处理器(Central Processing Unit，简称CPU)、网络处理器(Network Processor，简称NP)等；还可以是数字信号处理器(DigitalSignal Processing，简称DSP)、专用集成电路(Application Specific IntegratedCircuit，简称ASIC)、现成可编程门阵列(Field-Programmable Gate Array，简称FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。可以实现或者执行本申请实施例中的公开的各方法、步骤及逻辑框图。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。结合本申请实施例所公开的方法的步骤可以直接体现为硬件译码处理器执行完成，或者用译码处理器中的硬件及软件模块组合执行完成。软件模块可以位于随机存储器，闪存、只读存储器，可编程只读存储器或者电可擦写可编程存储器、寄存器等本领域成熟的存储介质中。该存储介质位于存储器602，处理器601读取存储器602中的信息，结合其硬件完成上述方法的步骤。
对应于上述SAR图像的处理方法，本申请实施例还提供了一种计算机可读存储介质，所述计算机可读存储介质存储有机器可运行指令，所述计算机可运行指令在被处理器调用和运行时，所述计算机可运行指令促使所述处理器运行上述SAR图像的处理方法的步骤。
本申请实施例所提供的SAR图像的处理装置可以为设备上的特定硬件或者安装于设备上的软件或固件等。本申请实施例所提供的装置，其实现原理及产生的技术效果和前述方法实施例相同，为简要描述，装置实施例部分未提及之处，可参考前述方法实施例中相应内容。所属领域的技术人员可以清楚地了解到，为描述的方便和简洁，前述描述的系统、装置和单元的具体工作过程，均可以参考上述方法实施例中的对应过程，在此不再赘述。
在本申请所提供的实施例中，应该理解到，所揭露装置和方法，可以通过其它的方式实现。以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，又例如，多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些通信接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。
再例如，附图中的流程图和框图显示了根据本申请的多个实施例的装置、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上，流程图或框图中的每个方框可以代表一个模块、程序段或代码的一部分，所述模块、程序段或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。也应当注意，在有些作为替换的实现方式中，方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如，两个连续的方框实际上可以基本并行地执行，它们有时也可以按相反的顺序执行，这依所涉及的功能而定。也要注意的是，框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合，可以用执行规定的功能或动作的专用的基于硬件的系统来实现，或者可以用专用硬件与计算机指令的组合来实现。
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
另外，在本申请提供的实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。
所述功能如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本申请的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台电子设备（可以是个人计算机，服务器，或者网络设备等）执行本申请各个实施例所述SAR图像的处理方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器（Read-Only Memory，简称ROM）、随机存取存储器（Random Access Memory，简称RAM）、磁碟或者光盘等各种可以存储程序代码的介质。
应注意到：相似的标号和字母在下面的附图中表示类似项，因此，一旦某一项在一个附图中被定义，则在随后的附图中不需要对其进行进一步定义和解释，此外，术语“第一”、“第二”、“第三”等仅用于区分描述，而不能理解为指示或暗示相对重要性。
最后应说明的是：以上所述实施例，仅为本申请的具体实施方式，用以说明本申请的技术方案，而非对其限制，本申请的保护范围并不局限于此，尽管参照前述实施例对本申请进行了详细的说明，本领域的普通技术人员应当理解：任何熟悉本技术领域的技术人员在本申请揭露的技术范围内，其依然可以对前述实施例所记载的技术方案进行修改或可轻易想到变化，或者对其中部分技术特征进行等同替换；而这些修改、变化或者替换，并不使相应技术方案的本质脱离本申请实施例技术方案的范围。都应涵盖在本申请的保护范围之内。因此，本申请的保护范围应以权利要求的保护范围为准。

Claims (10)

1.一种SAR图像的处理方法，其特征在于，包括：
获取单视复SAR图像，并对所述单视复SAR图像进行虚部数据和实部数据的提取，得到所述单视复SAR图像对应的原始虚部数据矩阵和原始实部数据矩阵；
对所述原始虚部数据矩阵进行逐行逐列包络对消处理，以消除所述原始虚部数据矩阵中的低频分量，进而得到包络对消后的虚部数据矩阵；
对所述原始实部数据矩阵进行逐行逐列包络对消处理，以消除所述原始实部数据矩阵中的低频分量，进而得到包络对消后的实部数据矩阵；
对所述包络对消后的虚部数据矩阵进行逐行逐列空间变迹处理，得到空间变迹后的虚部数据矩阵；
对所述包络对消后的实部数据矩阵进行逐行逐列空间变迹处理，得到空间变迹后的实部数据矩阵；
对所述空间变迹后的虚部数据矩阵和所述空间变迹后的实部数据矩阵进行合并，得到旁瓣抑制和图像增强后的单视复SAR图像。
2.根据权利要求1所述的方法，其特征在于，对所述原始虚部数据矩阵进行逐行逐列包络对消处理，包括：
对所述原始虚部数据矩阵的每行原始虚部数据中的每个原始虚部数据分别进行上包络提取和下包络提取，并根据提取得到的每个所述原始虚部数据的上包络和下包络对对应的原始虚部数据进行距离向包络对消处理，得到距离向包络对消后的虚部数据矩阵；
对所述距离向包络对消后的虚部数据矩阵的每列距离向包络对消后的虚部数据中的每个距离向包络对消后的虚部数据分别进行上包络提取和下包络提取，并根据提取得到的每个所述距离向包络对消后的虚部数据的上包络和下包络对对应的距离向包络对消后的虚部数据进行方位向包络对消处理，得到包络对消后的虚部数据矩阵；
对所述原始实部数据矩阵进行逐行逐列包络对消处理，包括：
对所述原始实部数据矩阵的每行原始实部数据中的每个原始实部数据分别进行上包络提取和下包络提取，并根据提取得到的每个所述原始实部数据的上包络和下包络对对应的原始实部数据进行距离向包络对消处理，得到距离向包络对消后的实部数据矩阵；
对所述距离向包络对消后的实部数据矩阵的每列距离向包络对消后的实部数据中的每个距离向包络对消后的实部数据分别进行上包络提取和下包络提取，并根据提取得到的每个所述距离向包络对消后的实部数据的上包络和下包络对对应的距离向包络对消后的实部数据进行方位向包络对消处理，得到包络对消后的实部数据矩阵。
3.根据权利要求2所述的方法，其特征在于，根据提取得到的每个所述原始虚部数据的上包络和下包络对对应的原始虚部数据进行距离向包络对消处理，包括：
若提取得到的第一原始虚部数据的上包络和下包络的符号相同，则将所述第一原始虚部数据与对应的上包络作差，并将得到的差值与对应的上包络和下包络的平均值作和，进而实现对所述第一原始虚部数据的距离向包络对消处理；
若提取得到的第二原始虚部数据的上包络和下包络的符号不同，则不对所述第二原始虚部数据进行距离向包络对消处理；
根据提取得到的每个所述距离向包络对消后的虚部数据的上包络和下包络对对应的距离向包络对消后的虚部数据进行方位向包络对消处理，包括：
若提取得到的第一距离向包络对消后的虚部数据的上包络和下包络的符号相同，则将所述第一距离向包络对消后的虚部数据与对应的上包络作差，并将得到的差值与对应的上包络和下包络的平均值作和，进而实现对所述第一距离向包络对消后的虚部数据的方位向包络对消处理；
若提取得到的第二距离向包络对消后的虚部数据的上包络和下包络的符号不同，则不对所述第二距离向包络对消后的虚部数据进行方位向包络对消处理；
根据提取得到的每个所述原始实部数据的上包络和下包络对对应的原始实部数据进行距离向包络对消处理，包括：
若提取得到的第一原始实部数据的上包络和下包络的符号相同，则将所述第一原始实部数据与对应的上包络作差，并将得到的差值与对应的上包络和下包络的平均值作和，进而实现对所述第一原始实部数据的距离向包络对消处理；
若提取得到的第二原始实部数据的上包络和下包络的符号不同，则不对所述第二原始实部数据进行距离向包络对消处理；
根据提取得到的每个所述距离向包络对消后的实部数据的上包络和下包络对对应的距离向包络对消后的实部数据进行方位向包络对消处理，包括：
若提取得到的第一距离向包络对消后的实部数据的上包络和下包络的符号相同，则将所述第一距离向包络对消后的实部数据与对应的上包络作差，并将得到的差值与对应的上包络和下包络的平均值作和，进而实现对所述第一距离向包络对消后的实部数据的方位向包络对消处理；
若提取得到的第二距离向包络对消后的实部数据的上包络和下包络的符号不同，则不对所述第二距离向包络对消后的实部数据进行方位向包络对消处理。
4.根据权利要求1所述的方法，其特征在于，对所述包络对消后的虚部数据矩阵进行逐行逐列空间变迹处理，包括：
对所述包络对消后的虚部数据矩阵的每行包络对消后的虚部数据，遍历其中的每个包络对消后的虚部数据，计算第一权重下的当前包络对消后的虚部数据的第一加权值和第二权重下的所述当前包络对消后的虚部数据的第二加权值；
若所述当前包络对消后的虚部数据的第一加权值与所述当前包络对消后的虚部数据的第二加权值的符号不同，则将所述当前包络对消后的虚部数据置为0；
若所述当前包络对消后的虚部数据的第一加权值与所述当前包络对消后的虚部数据的第二加权值的符号相同，则将所述当前包络对消后的虚部数据置为第一目标值，进而得到距离向空间变迹后的虚部数据矩阵，其中，所述第一目标值为基于所述当前包络对消后的虚部数据的第一加权值与所述当前包络对消后的虚部数据的第二加权值计算得到的；
对所述距离向空间变迹后的虚部数据矩阵的每列距离向空间变迹后的虚部数据，遍历其中的每个距离向空间变迹后的虚部数据，计算第一权重下的当前距离向空间变迹后的虚部数据的第一加权值和第二权重下的所述当前距离向空间变迹后的虚部数据的第二加权值；
若所述当前距离向空间变迹后的虚部数据的第一加权值与所述当前距离向空间变迹后的虚部数据的第二加权值的符号不同，则将所述当前距离向空间变迹后的虚部数据置为0；
若所述当前距离向空间变迹后的虚部数据的第一加权值与所述当前距离向空间变迹后的虚部数据的第二加权值的符号相同，则将所述当前距离向空间变迹后的虚部数据置为第二目标值，进而得到所述空间变迹后的虚部数据矩阵，其中，所述第二目标值为基于所述当前距离向空间变迹后的虚部数据的第一加权值与所述当前距离向空间变迹后的虚部数据的第二加权值计算得到的；
对所述包络对消后的实部数据矩阵进行逐行逐列空间变迹处理，包括：
对所述包络对消后的实部数据矩阵的每行包络对消后的实部数据，遍历其中的每个包络对消后的实部数据，计算第一权重下的当前包络对消后的实部数据的第一加权值和第二权重下的所述当前包络对消后的实部数据的第二加权值；
若所述当前包络对消后的实部数据的第一加权值与所述当前包络对消后的实部数据的第二加权值的符号不同，则将所述当前包络对消后的实部数据置为0；
若所述当前包络对消后的实部数据的第一加权值与所述当前包络对消后的实部数据的第二加权值的符号相同，则将所述当前包络对消后的实部数据置为第三目标值，进而得到距离向空间变迹后的实部数据矩阵，其中，所述第三目标值为基于所述当前包络对消后的实部数据的第一加权值与所述当前包络对消后的实部数据的第二加权值计算得到的；
对所述距离向空间变迹后的实部数据矩阵的每列距离向空间变迹后的实部数据，遍历其中的每个距离向空间变迹后的实部数据，计算第一权重下的当前距离向空间变迹后的实部数据的第一加权值和第二权重下的所述当前距离向空间变迹后的实部数据的第二加权值；
若所述当前距离向空间变迹后的实部数据的第一加权值与所述当前距离向空间变迹后的实部数据的第二加权值的符号不同，则将所述当前距离向空间变迹后的实部数据置为0；
若所述当前距离向空间变迹后的实部数据的第一加权值与所述当前距离向空间变迹后的实部数据的第二加权值的符号相同，则将所述当前距离向空间变迹后的实部数据置为第四目标值，进而得到所述空间变迹后的实部数据矩阵，其中，所述第四目标值为基于所述当前距离向空间变迹后的实部数据的第一加权值与所述当前距离向空间变迹后的实部数据的第二加权值计算得到的。
5.根据权利要求4所述的方法，其特征在于，计算第一权重下的当前包络对消后的虚部数据的第一加权值和第二权重下的所述当前包络对消后的虚部数据的第二加权值，包括：
根据加权值计算算式计算第一权重下的所述当前包络对消后的虚部数据的第一加权值/>和第二权重下的所述当前包络对消后的虚部数据的第二加权值/>，其中，/>表示所述当前包络对消后的虚部数据的加权值，/>表示所述当前包络对消后的虚部数据，/>表示所述第一加权值0或所述第二加权值/>，/>，，/>表示合成孔径雷达工作时的采样率，表示合成孔径雷达工作时的处理带宽，/>表示所述当前包络对消后的虚部数据之前的第k个包络对消后的虚部数据，/>表示所述当前包络对消后的虚部数据之后的第k个包络对消后的虚部数据，/>，/>为向下取整。
6.根据权利要求4所述的方法，其特征在于，所述第一目标值的计算算式如下：
其中，表示所述第一目标值，/>表示所述当前包络对消后的虚部数据的第一加权值，/>表示所述当前包络对消后的虚部数据的第二加权值，sign()函数表示取数值符号，min()函数表示取最小值，abs()函数表示取绝对值。
7.根据权利要求1所述的方法，其特征在于，对所述空间变迹后的虚部数据矩阵和所述空间变迹后的实部数据矩阵进行合并，包括：
采用实部和虚部加和的方式对所述空间变迹后的虚部数据矩阵和所述空间变迹后的实部数据矩阵进行合并，得到处理后的单视复数据；
将所述处理后的单视复数据可视化，得到所述旁瓣抑制和图像增强后的单视复SAR图像。
8.一种SAR图像的处理装置，其特征在于，包括：
提取单元，用于获取单视复SAR图像，并对所述单视复SAR图像进行虚部数据和实部数据的提取，得到所述单视复SAR图像对应的原始虚部数据矩阵和原始实部数据矩阵；
虚部包络对消处理单元，用于对所述原始虚部数据矩阵进行逐行逐列包络对消处理，以消除所述原始虚部数据矩阵中的低频分量，进而得到包络对消后的虚部数据矩阵；
实部包络对消处理单元，用于对所述原始实部数据矩阵进行逐行逐列包络对消处理，以消除所述原始实部数据矩阵中的低频分量，进而得到包络对消后的实部数据矩阵；
虚部空间变迹处理单元，用于对所述包络对消后的虚部数据矩阵进行逐行逐列空间变迹处理，得到空间变迹后的虚部数据矩阵；
实部空间变迹处理单元，用于对所述包络对消后的实部数据矩阵进行逐行逐列空间变迹处理，得到空间变迹后的实部数据矩阵；
合并单元，用于对所述空间变迹后的虚部数据矩阵和所述空间变迹后的实部数据矩阵进行合并，得到旁瓣抑制和图像增强后的单视复SAR图像。
9.一种电子设备，包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现上述权利要求1至7中任一项所述的方法的步骤。
10.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有机器可运行指令，所述机器可运行指令在被处理器调用和运行时，所述机器可运行指令促使所述处理器运行上述权利要求1至7中任一项所述的方法。
CN202311576620.5A 2023-11-24 2023-11-24 Sar图像的处理方法、装置、电子设备和存储介质 Active CN117289276B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202311576620.5A CN117289276B (zh) 2023-11-24 2023-11-24 Sar图像的处理方法、装置、电子设备和存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202311576620.5A CN117289276B (zh) 2023-11-24 2023-11-24 Sar图像的处理方法、装置、电子设备和存储介质

Publications (2)

Publication Number Publication Date
CN117289276A true CN117289276A (zh) 2023-12-26
CN117289276B CN117289276B (zh) 2024-02-27

Family

ID=89252032

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202311576620.5A Active CN117289276B (zh) 2023-11-24 2023-11-24 Sar图像的处理方法、装置、电子设备和存储介质

Country Status (1)

Country Link
CN (1) CN117289276B (zh)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5349359A (en) * 1993-05-10 1994-09-20 Environmental Research Institute Of Michigan Spatially variant apodization
CN101957458A (zh) * 2009-07-21 2011-01-26 北京大学 一种处理超低频数据的方法
CN106483506A (zh) * 2016-10-26 2017-03-08 河海大学 一种认知型自适应干扰抑制方法
CN106780329A (zh) * 2016-12-07 2017-05-31 华中科技大学 一种基于反透视平面变换的超声平面波成像方法
CN110632596A (zh) * 2019-10-09 2019-12-31 上海无线电设备研究所 一种太赫兹sar多频振动误差补偿方法
CN114706075A (zh) * 2022-03-08 2022-07-05 重庆邮电大学 一种毫米波近场sar图像旁瓣抑制方法、设备及存储介质
CN115079171A (zh) * 2022-06-22 2022-09-20 青岛大学 基于混合范数的自适应稀疏约束逆合成孔径雷达成像方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5349359A (en) * 1993-05-10 1994-09-20 Environmental Research Institute Of Michigan Spatially variant apodization
CN101957458A (zh) * 2009-07-21 2011-01-26 北京大学 一种处理超低频数据的方法
CN106483506A (zh) * 2016-10-26 2017-03-08 河海大学 一种认知型自适应干扰抑制方法
CN106780329A (zh) * 2016-12-07 2017-05-31 华中科技大学 一种基于反透视平面变换的超声平面波成像方法
CN110632596A (zh) * 2019-10-09 2019-12-31 上海无线电设备研究所 一种太赫兹sar多频振动误差补偿方法
CN114706075A (zh) * 2022-03-08 2022-07-05 重庆邮电大学 一种毫米波近场sar图像旁瓣抑制方法、设备及存储介质
CN115079171A (zh) * 2022-06-22 2022-09-20 青岛大学 基于混合范数的自适应稀疏约束逆合成孔径雷达成像方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
张平;杨汝良;: "合成孔径雷达空间变迹带宽外推超分辨算法", 测绘学报, no. 01 *
徐光耀;樊霞良;: "基于变迹滤波的SAR图像旁瓣抑制技术", 指挥控制与仿真, no. 02 *
杨科;廖桂生;徐青;王伟伟;: "改进的合成孔径雷达旁瓣抑制空间变迹算法", 电波科学学报, no. 06 *

Also Published As

Publication number Publication date
CN117289276B (zh) 2024-02-27

Similar Documents

Publication Publication Date Title
CN111242137B (zh) 一种基于形态成分分析的椒盐噪声滤波方法、装置
US6061100A (en) Noise reduction for video signals
EP1411471B1 (en) Removing color aliasing artifacts from color digital images
Zhu et al. A new algorithm for SAR image despeckling using an enhanced Lee filter and median filter
CN108109123B (zh) 一种图像去噪方法
US9319602B2 (en) Method and apparatus for processing image
CN111179186A (zh) 一种保护图像细节的图像去噪系统
US8055093B2 (en) Method and system for noise level detection in image data
CN117289276B (zh) Sar图像的处理方法、装置、电子设备和存储介质
US8090215B2 (en) Method and system for noise level detection in image data
US20040046894A1 (en) Reducing color aliasing artifacts from color digital images
CN111723753B (zh) 卫星遥感图像的去条带方法、装置和电子设备
US8208753B2 (en) Method and system for noise level detection in image data
US20110254728A1 (en) Three quarter spatially variant apodization
Choppala et al. Isolated vector median filtering for noise reduction in digital color images
KR101412076B1 (ko) 전방 관측 영상 레이더의 영상 품질 향상 방법
CN111652811B (zh) 一种基于刃边函数和最优窗维纳滤波的运动模糊图像复原方法
CN110889878B (zh) 结合图像平滑处理的ct图像环形伪影抑制方法
CN114820376A (zh) 条带噪声的融合校正方法、装置、电子设备及存储介质
CN109615597B (zh) 一种基于中心仿射滤波的图像斑点噪声抑制方法
US8180169B2 (en) System and method for multi-scale sigma filtering using quadrature mirror filters
CN110766005B (zh) 目标特征提取方法、装置和终端设备
CN114280571A (zh) 一种雨杂波信号的处理方法、装置及设备
CN107154019B (zh) 一种基于lr-几何插值的电磁成像超分辨方法
Rajesh et al. An effective denoising and enhancement technique for removal of random impulse noise in images

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