CN109343056A - 一种针对非线性调频sar的rd成像方法和装置 - Google Patents

一种针对非线性调频sar的rd成像方法和装置 Download PDF

Info

Publication number
CN109343056A
CN109343056A CN201811290311.0A CN201811290311A CN109343056A CN 109343056 A CN109343056 A CN 109343056A CN 201811290311 A CN201811290311 A CN 201811290311A CN 109343056 A CN109343056 A CN 109343056A
Authority
CN
China
Prior art keywords
data
pulse
distance
apart
frequency domain
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
CN201811290311.0A
Other languages
English (en)
Other versions
CN109343056B (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.)
Institute of Electronics of CAS
Original Assignee
Institute of Electronics of CAS
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 Institute of Electronics of CAS filed Critical Institute of Electronics of CAS
Priority to CN201811290311.0A priority Critical patent/CN109343056B/zh
Publication of CN109343056A publication Critical patent/CN109343056A/zh
Application granted granted Critical
Publication of CN109343056B publication Critical patent/CN109343056B/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/9004—SAR image acquisition techniques
    • 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)
  • Electromagnetism (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明实施例公开了一种针对非线性调频SAR的RD成像方法和装置，该方法包括：在发射非线性调频信号后，对回波数据进行距离向傅里叶变换，得到距离频域的回波数据；根据所述非线性调频信号，对所述距离频域的回波数据进行距离向脉冲压缩处理，得到距离向脉冲压缩的距离频域数据；将所述距离向脉冲压缩的距离频域数据变换至距离多普勒域，得到距离多普勒域的数据；对所述距离多普勒域的数据依次进行距离徙动校正RCMC和方位向脉冲压缩处理，得到方位向脉冲压缩处理后的数据；根据所述方位向脉冲压缩处理后的数据，得到SAR图像数据；如此，可以对非线性调频信号为发射波形的SAR，提供改进的RD算法的成像方法，并获得聚焦良好的图像。

Description

一种针对非线性调频SAR的RD成像方法和装置
技术领域
本发明涉及合成孔径雷达(Synthetic Aperture Radar，SAR)成像技术，尤其涉及一种针对非线性调频SAR的距离多普勒(Range Doppler，RD)成像方法和装置，可以应用于以非线性调频信号为发射波形的SAR在小斜视角下的成像。
背景技术
SAR是一种高分辨率成像雷达，可以在能见度极低的气象条件下得到类似光学照相的高分辨雷达图像，SAR通常装在飞机或卫星上，按平台的运动航迹来测距和二维成像。线性调频(Linear Frequency Modulation，LFM)信号是SAR中最常用的发射波形，这种波形通过匹配滤波后，生成的响应函数，经归一化后的峰值旁瓣比(PSLR，Peak Side LobeRatio)为-13dB。为了抑制旁瓣的高度，通常采用加权窗函数、自适应滤波和优化算法，但这些抑制旁瓣的高度的方法会使得匹配滤波器失配，进而降低了输出的信噪比(Signal toNoise Ratio，SNR)。
相比于LFM信号，非线性调频(Non-Linear Frequency Modulation，NLFM)信号的脉冲压缩结果可以获得很低的PSLR，并且，并没有损失输出SNR。通过相关实验，NLFM信号可以避免SNR的1-2dB损失，相当于节省25％的天线发射功率。对于能量紧缺的雷达系统中，采用NLFM信号作为发射波形可以提升系统性能。
RD算法是SAR成像处理中最直观、最基本的经典方法，可以在两个一维操作之间使用距离徙动校正，对距离和方位进行了分离处理。该算法能兼顾成熟、简便和高效等因素，至今仍是使用最广泛的成像处理算法之一。但是RD算法要求回波信号具有线性调频信号特性，对于非线性调频信号不适用。
发明内容
本发明实施例期望提供一种针对非线性调频SAR的RD成像方法和装置，能够对非线性调频信号为发射波形的SAR，提供基于RD算法的成像方法。
本发明实施例提供了一种针对非线性调频SAR的RD成像方法，所述方法包括：
在发射非线性调频信号后，对回波数据进行距离向傅里叶变换，得到距离频域的回波数据；
根据所述非线性调频信号，对所述距离频域的回波数据进行距离向脉冲压缩处理，得到距离向脉冲压缩的距离频域数据；
将所述距离向脉冲压缩的距离频域数据变换至距离多普勒域，得到距离多普勒域的数据；对所述距离多普勒域的数据依次进行距离徙动校正(Range Cell MigrationCorrection，RCMC)和方位向脉冲压缩处理，得到方位向脉冲压缩处理后的数据；根据所述方位向脉冲压缩处理后的数据，得到SAR图像数据。
上述方案中，所述根据所述非线性调频信号，对所述距离频域的回波数据进行距离向脉冲压缩处理，得到距离向脉冲压缩的数据，包括：
对所述非线性调频信号进行补零，使补零后的非线性调频信号的长度等于所述非线性调频信号对应的原始回波数据的长度；对所述补零后的非线性调频信号进行傅里叶变换和取共轭操作，得到处理后的非线性调频信号；对所述处理后的非线性调频信号与所述距离频域的回波数据进行相乘，得到距离向脉冲压缩的距离频域数据。
上述方案中，所述将所述距离向脉冲压缩的距离频域数据变换至距离多普勒域，得到距离多普勒域的数据，包括：
对所述距离向脉冲压缩的距离频域数据进行运动误差补偿，得到运动误差补偿后的数据，对所述运动误差补偿后的数据进行距离向傅里叶逆变换和方位向傅里叶变换，得到距离多普勒域的数据。
上述方案中，所述对所述距离向脉冲压缩的距离频域数据进行运动误差补偿，得到运动误差补偿后的数据，包括：
确定SAR对应的斜距误差，根据所述斜距误差对所述距离向脉冲压缩的距离频域数据进行运动误差补偿，得到运动误差补偿后的数据。
上述方案中，所述根据所述方位向脉冲压缩处理后的数据，得到SAR图像数据，包括：
对所述方位向脉冲压缩处理后的数据依次进行残余运动误差补偿和方位向傅里叶逆变换，得到SAR图像数据。
本发明实施例还提供了一种针对非线性调频SAR的RD成像装置，所述装置包括：获取模块、第一处理模块和第二处理模块；其中，
获取模块，用于在SAR发射非线性调频信号后，对回波数据进行距离向傅里叶变换，得到距离频域的回波数据；
第一处理模块，用于根据所述非线性调频信号，对所述距离频域的回波数据进行距离向脉冲压缩处理，得到距离向脉冲压缩的距离频域数据；
第二处理模块，用于将所述距离向脉冲压缩的距离频域数据变换至距离多普勒域，得到距离多普勒域的数据；对所述距离多普勒域的数据依次进行RCMC和方位向脉冲压缩处理，得到方位向脉冲压缩处理后的数据；根据所述方位向脉冲压缩处理后的数据，得到SAR图像数据。
上述方案中，所述第一处理模块，具体用于对所述非线性调频信号进行补零，使补零后的非线性调频信号的长度等于所述非线性调频信号对应的原始回波数据的长度；对所述补零后的非线性调频信号进行傅里叶变换和取共轭操作，得到处理后的非线性调频信号；对所述处理后的非线性调频信号与所述距离频域的回波数据进行相乘，得到距离向脉冲压缩的距离频域数据。
上述方案中，所述第二处理模块，具体用于对所述距离向脉冲压缩的距离频域数据进行运动误差补偿，得到运动误差补偿后的数据，对所述运动误差补偿后的数据进行距离向傅里叶逆变换和方位向傅里叶变换，得到距离多普勒域的数据。
本发明实施例还提供了一种针对非线性调频SAR的RD成像设备，所述设备包括：处理器和用于存储能够在处理器上运行的计算机程序的存储器；其中，
所述处理器用于运行所述计算机程序时，执行上述任意一种针对非线性调频SAR的RD成像方法的步骤。
本发明实施例还提供了一种计算机存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现上述任意一种针对非线性调频SAR的RD成像方法的步骤。
本发明实施例中，在发射非线性调频信号后，对回波数据进行距离向傅里叶变换，得到距离频域的回波数据；根据所述非线性调频信号，对所述距离频域的回波数据进行距离向脉冲压缩处理，得到距离向脉冲压缩的距离频域数据；将所述距离向脉冲压缩的距离频域数据变换至距离多普勒域，得到距离多普勒域的数据；对所述距离多普勒域的数据依次进行RCMC和方位向脉冲压缩处理，得到方位向脉冲压缩处理后的数据；根据所述方位向脉冲压缩处理后的数据，得到SAR图像数据。可以看出，在SAR成像处理过程中，对距离向脉冲压缩处理进行改进，即，将非线性调频信号作为距离向脉冲压缩处理的依据，从而，可以使基于RD算法的成像方法适用于非线性调频信号为发射波形的SAR中；也就是说，对于非线性调频信号为发射波形的SAR，提供基于RD算法的成像方法，获得聚焦良好的图像，从而提升SAR的性能。
附图说明
图1为本发明实施例的针对非线性调频SAR的RD成像方法的流程图；
图2为本发明实施例中设计的一个成像场景的示意图；
图3为本发明实施例中对于图2所示的成像场景得出的聚焦后的点目标成像结果示意图；
图4a表示本发明实施例的图2所示的成像场景中P1点的幅度轮廓图；
图4b表示本发明实施例的图2所示的成像场景中P1点的距离向剖面示意图；
图4c表示本发明实施例的图2所示的成像场景中P1点的方位向剖面示意图；
图4d表示本发明实施例的图2所示的成像场景中P5点的幅度轮廓图；
图4e表示本发明实施例的图2所示的成像场景中P5点的距离向剖面示意图；
图4f表示本发明实施例的图2所示的成像场景中P5点的方位向剖面示意图；
图5为本发明实施例的针对非线性调频SAR的RD成像装置的组成结构示意图；
图6为本发明实施例的针对非线性调频SAR的RD成像设备的硬件结构示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述。
本发明实施例中，在发射非线性调频信号后，对回波数据进行距离向傅里叶变换，得到距离频域的回波数据；根据所述非线性调频信号，对所述距离频域的回波数据进行距离向脉冲压缩处理，得到距离向脉冲压缩的距离频域数据；将所述距离向脉冲压缩的距离频域数据变换至距离多普勒域，得到距离多普勒域的数据；对所述距离多普勒域的数据依次进行距离徙动校正RCMC和方位向脉冲压缩处理，得到方位向脉冲压缩处理后的数据；根据所述方位向脉冲压缩处理后的数据，得到SAR图像数据。
下面结合实施例对本发明作进一步详细的说明。
图1为本发明实施例的针对非线性调频SAR的RD成像方法的流程图，如图1所示，该流程可以包括：
步骤101：在发射非线性调频信号后，对回波数据进行距离向傅里叶变换，得到距离频域的回波数据；
实际应用中，可以控制SAR发射非线性调频信号，然后利用SAR接收回波数据；然后，对接收的回波数据进行解调以及距离向傅里叶变换，得到距离频域的回波数据。
在一个具体的示例中，SAR发射信号在遇到障碍物后形成反射，对于一个反射点其回波模型经过解调后，可以用式(1)表示。
其中，r(t,η)表示解调后的回波模型，A0是复常数，表示信号的幅度；t表示距离向时间，η表示方位时间，ηc表示波束中心偏移时间，f0表示载频，R(η)表示瞬时斜距，wr和wa分别表示距离包络和方位包络，φ(t-2R(η)/c)表示时移后的非线性调频信号相位，c表示光速。
在理想情况下，SAR被安装在飞行器上，飞行器的飞行轨迹为直线，那么，SAR的飞行轨迹同样为直线，此时点目标的瞬时斜距可以表示为其中，R0表示最近斜距，Vr表示SAR的速度。
对解调后的回波数据进行距离向傅里叶变换后，可以将式(1)表示的回波模型变换到距离频域，得到距离频域的回波数据R(ft,η)，R(ft,η)可以用式(2)表示。
其中，A1表示频谱的幅度，Wr表示距离向的频率包络，ft表示距离向频率，f0表示载频，其中Θ(ft)表示非线性调频信号基带频谱的相位项。
步骤102：根据所述非线性调频信号，对所述距离频域的回波数据进行距离向脉冲压缩处理，得到距离向脉冲压缩的距离频域数据；
在本步骤的一种实现方式中，可以对所述非线性调频信号进行两端补零，使补零后的非线性调频信号的长度等于所述非线性调频信号对应的原始回波数据的长度；对所述补零后的非线性调频信号进行傅里叶变换和取共轭操作，得到处理后的非线性调频信号；对所述处理后的非线性调频信号与所述距离频域的回波数据进行相乘，得到距离向脉冲压缩的距离频域数据。
具体地说，SAR发射的非线性调频信号可以表示为s(t)，通过对s(t)两端补零，使补零后的信号达到回波数据的长度，补零后的信号表示为s′(t)，对s′(t)进行傅里叶变换，并取共轭，得到处理后的非线性调频信号Hr(ft)，Hr(ft)可以用式(3)表示。
Hr(ft)＝exp(-jΘ(ft))＝[FT[s′(t)]]* (3)
其中，表示傅里叶变换。
可以看出，将处理后的非线性调频信号Hr(ft)与距离频域的回波数据R(ft,η)相乘，可以消除Θ(ft)项。
步骤103：将所述距离向脉冲压缩的距离频域数据变换至距离多普勒域，得到距离多普勒域的数据；对所述距离多普勒域的数据依次进行RCMC和方位向脉冲压缩处理，得到方位向脉冲压缩处理后的数据；根据所述方位向脉冲压缩处理后的数据，得到SAR图像数据。
对于将所述距离向脉冲压缩的距离频域数据变换至距离多普勒域，得到距离多普勒域的数据的实现方式，示例性地，可以对所述距离向脉冲压缩的距离频域数据进行运动误差补偿，得到运动误差补偿后的数据，对所述运动误差补偿后的数据进行距离向傅里叶逆变换和方位向傅里叶变换，得到距离多普勒域的数据。
这里，对所述距离向脉冲压缩的距离频域数据进行运动误差补偿的一种实现方式可以是：确定SAR对应的斜距误差，根据所述斜距误差对所述距离向脉冲压缩的距离频域数据进行运动误差补偿，得到运动误差补偿后的数据。
可以理解的是，由于SAR的载机飞行不平稳，实际飞行轨迹会偏离理想航迹，实际实施时，可以记录载机在SAR的每个发射脉冲时的东、北、天向速度，之后，通过计算，可以拟合出载机的真实运动轨迹。从得出的真实运动轨迹中，计算载机与SAR发射信号的覆盖范围中心(即景中心目标)沿视线方向的真实斜距，将计算得出的真实斜距与理想斜距相减，得到真实的斜距误差。这个斜距误差在时域相当于延时误差，但是这个延时通常带有小数延迟，因此补偿的办法可以是：通过将所述经过距离向脉冲压缩后的距离频域数据，乘以一个对应的相位，补偿掉这个斜距误差。通过上述步骤，可以补偿参考距离处的运动误差，得到运动误差补偿后的数据。本发明实施例中，可以采用载机搭载的惯性导航系统，记录载机在每个发射脉冲时的东、北、天向速度。
在一实施方式中，对所述运动误差补偿后的数据进行距离向傅里叶逆变换后，得到的时域信号可以用R(t,η)表示，R(t,η)可以用式(4)表示。
其中，pr为脉冲压缩后的距离向包络。
然后对其进行方位向傅里叶变换得距离多普勒域的数据
其中Rrd(fη)为距离多普勒域中的斜距表达，为
RCMC的目的就是使得Rrd(fη)＝R0，具体的，在对所述距离多普勒域的数据进行RCMC的一种实现方式中，可以对在距离多普勒域进行距离插值运算完成RCMC操作，进行RCMC时需要校正的徙动量ΔR(fη)为：
其中，λ表示非线性调频信号的波长，fη表示方位频率。
在对所述距离多普勒域的数据进行RCMC，可以得到RCMC后的距离多普勒域的信号R(t,fη)，R(t,fη)的表达式为：
其中，A2代表频谱幅度，表示多普勒中心频率，Ka表示方位向调频斜率，Wa表示方位向的频率包络
在一个示例中，
在得到RCMC后的距离多普勒域的信号，可以对RCMC后的距离多普勒域的信号进行方位向脉冲压缩处理，得到方位向脉冲压缩处理后的数据。
在对RCMC后的距离多普勒域的信号进行方位向脉冲压缩处理的一种实现方式中，可以在距离多普勒域进行方位向匹配滤波实现方位脉冲压缩操作，方位向脉冲压缩处理需要使用的匹配滤波器可以用式(7)表示。
其中，Ha(fη)表示方位向脉冲压缩处理需要使用的匹配滤波器。
方位向脉冲压缩处理后的数据可以用Rac(t,fη)表示。
对于根据所述方位向脉冲压缩处理后的数据，得到SAR图像数据的实现方式，在一个示例中，可以对所述方位向脉冲压缩处理后的数据依次进行残余运动误差补偿和方位向傅里叶逆变换，得到SAR图像数据。
在一实施方式中，可以采用空变的运动补偿方法补偿残余运动误差，得到补偿残余运动误差后的数据，具体地说，可以对方位向脉冲压缩处理后的数据，采用空变的运动补偿方法补偿残余运动误差，得到补偿残余运动误差后的数据；本发明实施例中，空变的运动补偿方法可以是相位梯度自聚焦(Phase Gradient Autofocus，PGA)方法。
最后，通过进行方位向傅里叶逆变换，即可以完成二维图像聚焦，得到聚焦良好的SAR图像数据，再进行方位向傅里叶逆变换，表明已完成了非线性调频信号的成像。
需要说明的是，SAR图像数据中，数据存在除参考距离处的残余运动误差，并且还存在一些未知的误差，比如系统本身的误差，还有非运动误差。
下面通过仿真对本发明实施例的效果进行说明。
在采用非线性调频信号作为发射波形(其自相关函数旁瓣为-40.5dB)，仿真的系统参数如表1所示。
工作模式 条带
雷达速度 125m/s
载频 9.6GHz
脉冲重复频率(PRF) 1500Hz
发射信号带宽 100MHz
脉冲宽度 10us
表1
图2为本发明实施例中设计的一个成像场景的示意图；图2中，横坐标表示距离向，纵坐标表示方位向，图2所示的成像场景共有九个点目标，分别表示为P1、P2、P3、P4、P5、P6、P7、P8和P9。
对于图2所示的成像场景，可以采用前述记载的针对非线性调频SAR的RD成像方法进行成像，得出图3所示的聚焦后的点目标成像结果示意图，图3中，横坐标表示距离向，单位为km，纵坐标表示方位向，单位为km，从图3中可以看出，聚焦后的点目标包括P1和P5。
这里，在进行仿真时，没有引入运动误差，因此，在对图2所示的成像场景采用前述记载的SAR成像方法进行成像的过程中，可以省略前述记载的运动误差补偿和残余运动误差补偿。
对图2所示的点P1、P5分别进行插值，可以得到相应的点目标幅度轮廓图、距离向剖面示意图和方位向剖面示意图。图4a表示本发明实施例的图2所示的成像场景中P1点的幅度轮廓图，图4b表示本发明实施例的图2所示的成像场景中P1点的距离向剖面示意图，图4c表示本发明实施例的图2所示的成像场景中P1点的方位向剖面示意图，图4d表示本发明实施例的图2所示的成像场景中P5点的幅度轮廓图，图4e表示本发明实施例的图2所示的成像场景中P5点的距离向剖面示意图，图4f表示本发明实施例的图2所示的成像场景中P5点的方位向剖面示意图；图4a和图4d中，横坐标表示距离向采样点，纵坐标表示方位向采样点，图4b和图4e中，横坐标表示距离向采样点，纵坐标表示幅度，单位为dB，图4c和图4f中，横坐标表示方位向采样点，纵坐标表示幅度，单位为dB。
对于图3聚焦后的点目标P1和P5，质量评估参数如表2所示。
表2
其中，ISLR表示积分旁瓣比(Integrated Side Lobe Ratio)。
基于前述记载的针对非线性调频SAR的RD成像方法，本发明实施例还提出了一种针对非线性调频SAR的RD成像装置。
图5为本发明实施例的针对非线性调频SAR的RD成像装置的组成结构示意图，如图5所示，该装置包括获取模块501、第一处理模块502和第二处理模块503；其中，
获取模块501，用于在SAR发射非线性调频信号后，对回波数据进行距离向傅里叶变换，得到距离频域的回波数据；
第一处理模块502，用于根据所述非线性调频信号，对所述距离频域的回波数据进行距离向脉冲压缩处理，得到距离向脉冲压缩的距离频域数据；
第二处理模块503，用于将所述距离向脉冲压缩的距离频域数据变换至距离多普勒域，得到距离多普勒域的数据；对所述距离多普勒域的数据依次进行距离徙动校正RCMC和方位向脉冲压缩处理，得到方位向脉冲压缩处理后的数据；根据所述方位向脉冲压缩处理后的数据，得到SAR图像数据。
在一实施方式中，所述第一处理模块502，具体用于对所述非线性调频信号进行补零，使补零后的非线性调频信号的长度等于所述非线性调频信号对应的原始回波数据的长度；对所述补零后的非线性调频信号进行傅里叶变换和取共轭操作，得到处理后的非线性调频信号；对所述处理后的非线性调频信号与所述距离频域的回波数据进行相乘，得到距离向脉冲压缩的距离频域数据。
在一实施方式中，所述第二处理模块503，具体用于对所述距离向脉冲压缩的距离频域数据进行运动误差补偿，得到运动误差补偿后的数据，对所述运动误差补偿后的数据进行距离向傅里叶逆变换和方位向傅里叶变换，得到距离多普勒域的数据。
在一实施方式中，所述第二处理模块503，具体用于确定SAR对应的斜距误差，根据所述斜距误差对所述距离向脉冲压缩的距离频域数据进行运动误差补偿，得到运动误差补偿后的数据。
在一实施方式中，所述第二处理模块503，具体用于对所述方位向脉冲压缩处理后的数据依次进行残余运动误差补偿和方位向傅里叶逆变换，得到SAR图像数据。
实际应用中，上述获取模块501、第一处理模块502和第二处理模块503均可由中央处理器(Central Processing Unit，CPU)、微处理器(Micro Processor Unit，MPU)、数字信号处理器(Digital Signal Processor，DSP)、或现场可编程门阵列(Field ProgrammableGate Array，FPGA)等实现。
另外，在本实施例中的各功能模块可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现，也可以采用软件功能模块的形式实现。
所述集成的单元如果以软件功能模块的形式实现并非作为独立的产品进行销售或使用时，可以存储在一个计算机可读取存储介质中，基于这样的理解，本实施例的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的全部或部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)或processor(处理器)执行本实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(Read Only Memory，ROM)、随机存取存储器(Random Access Memory，RAM)、磁碟或者光盘等各种可以存储程序代码的介质。
具体来讲，本实施例中的一种针对非线性调频SAR的RD成像方法对应的计算机程序指令可以被存储在光盘，硬盘，U盘等存储介质上，当存储介质中的与一种针对非线性调频SAR的RD成像方法对应的计算机程序指令被一电子设备读取或被执行时，实现前述实施例的任意一种针对非线性调频SAR的RD成像方法的步骤。
基于前述实施例相同的技术构思，参见图6，其示出了本发明实施例提供的一种针对非线性调频SAR的RD成像设备60，可以包括：存储器61、处理器62和总线63；其中，
所述总线63用于连接所述存储器61、处理器62和这些器件之间的相互通信；
所述存储器61，用于存储计算机程序和数据；
所述处理器62，用于执行所述存储器中存储的计算机程序，以实现前述实施例的任意一种针对非线性调频SAR的RD成像方法的步骤。
在实际应用中，上述存储器61可以是易失性存储器(volatile memory)，例如RAM；或者非易失性存储器(non-volatile memory)，例如ROM，快闪存储器(flash memory)，硬盘(Hard Disk Drive，HDD)或固态硬盘(Solid-State Drive，SSD)；或者上述种类的存储器的组合，并向处理器62提供指令和数据。
上述处理器62可以为特定用途集成电路(Application Specific IntegratedCircuit，ASIC)、DSP、数字信号处理装置(Digital Signal Processing Device，DSPD)、可编程逻辑装置(Programmable Logic Device，PLD)、FPGA、CPU、控制器、微控制器、微处理器中的至少一种。可以理解地，对于不同的设备，用于实现上述处理器功能的电子器件还可以为其它，本发明实施例不作具体限定。
上述本发明实施例序号仅仅为了描述，不代表实施例的优劣。
通过以上的实施方式的描述，本领域的技术人员可以清楚地了解到上述实施例方法可借助软件加必需的通用硬件平台的方式来实现，当然也可以通过硬件，但很多情况下前者是更佳的实施方式。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质(如ROM/RAM、磁碟、光盘)中，包括若干指令用以使得一台终端(可以是手机，计算机，服务器，空调器，或者网络设备等)执行本发明各个实施例所述的方法。
上面结合附图对本发明的实施例进行了描述，但是本发明并不局限于上述的具体实施方式，上述的具体实施方式仅仅是示意性的，而不是限制性的，本领域的普通技术人员在本发明的启示下，在不脱离本发明宗旨和权利要求所保护的范围情况下，还可做出很多形式，这些均属于本发明的保护之内。

Claims (10)

1.一种针对非线性调频合成孔径雷达SAR的距离多普勒RD成像方法，其特征在于，所述方法包括：
在发射非线性调频信号后，对回波数据进行距离向傅里叶变换，得到距离频域的回波数据；
根据所述非线性调频信号，对所述距离频域的回波数据进行距离向脉冲压缩处理，得到距离向脉冲压缩的距离频域数据；
将所述距离向脉冲压缩的距离频域数据变换至距离多普勒域，得到距离多普勒域的数据；对所述距离多普勒域的数据依次进行距离徙动校正RCMC和方位向脉冲压缩处理，得到方位向脉冲压缩处理后的数据；根据所述方位向脉冲压缩处理后的数据，得到SAR图像数据。
2.根据权利要求1所述的方法，其特征在于，所述根据所述非线性调频信号，对所述距离频域的回波数据进行距离向脉冲压缩处理，得到距离向脉冲压缩的距离频域数据，包括：
对所述非线性调频信号进行补零，使补零后的非线性调频信号的长度等于所述非线性调频信号对应的原始回波数据的长度；对所述补零后的非线性调频信号进行傅里叶变换和取共轭操作，得到处理后的非线性调频信号；对所述处理后的非线性调频信号与所述距离频域的回波数据进行相乘，得到距离向脉冲压缩的距离频域数据。
3.根据权利要求1所述的方法，其特征在于，所述将所述距离向脉冲压缩的距离频域数据变换至距离多普勒域，得到距离多普勒域的数据，包括：
对所述距离向脉冲压缩的距离频域数据进行运动误差补偿，得到运动误差补偿后的数据，对所述运动误差补偿后的数据进行距离向傅里叶逆变换和方位向傅里叶变换，得到距离多普勒域的数据。
4.根据权利要求3所述的方法，其特征在于，所述对所述距离向脉冲压缩的数据进行运动误差补偿，得到运动误差补偿后的数据，包括：
确定SAR对应的斜距误差，根据所述斜距误差对所述距离向脉冲压缩的距离频域数据进行运动误差补偿，得到运动误差补偿后的数据。
5.根据权利要求1所述的方法，其特征在于，所述根据所述方位向脉冲压缩处理后的数据，得到SAR图像数据，包括：
对所述方位向脉冲压缩处理后的数据依次进行残余运动误差补偿和方位向傅里叶逆变换，得到SAR图像数据。
6.一种针对非线性调频合成孔径雷达SAR的距离多普勒RD成像装置，其特征在于，所述装置包括：获取模块、第一处理模块和第二处理模块；其中，
获取模块，用于在SAR发射非线性调频信号后，对回波数据进行距离向傅里叶变换，得到距离频域的回波数据；
第一处理模块，用于根据所述非线性调频信号，对所述距离频域的回波数据进行距离向脉冲压缩处理，得出距离向脉冲压缩的距离频域数据；
第二处理模块，用于将所述距离向脉冲压缩的距离频域数据变换至距离多普勒域，得到距离多普勒域的数据；对所述距离多普勒域的数据依次进行距离徙动校正RCMC和方位向脉冲压缩处理，得到方位向脉冲压缩处理后的数据；根据所述方位向脉冲压缩处理后的数据，得到SAR图像数据。
7.根据权利要求6所述的装置，其特征在于，所述第一处理模块，具体用于对所述非线性调频信号进行补零，使补零后的非线性调频信号的长度等于所述非线性调频信号对应的原始回波数据的长度；对所述补零后的非线性调频信号进行傅里叶变换和取共轭操作，得到处理后的非线性调频信号；对所述处理后的非线性调频信号与所述距离频域的回波数据进行相乘，得到距离向脉冲压缩的距离频域数据。
8.根据权利要求6所述的装置，其特征在于，所述第二处理模块，具体用于对所述距离向脉冲压缩的距离频域数据进行运动误差补偿，得到运动误差补偿后的数据，对所述运动误差补偿后的数据进行距离向傅里叶逆变换和方位向傅里叶变换，得到距离多普勒域的数据。
9.一种针对非线性调频合成孔径雷达SAR的距离多普勒RD成像设备，其特征在于，所述设备包括：处理器和用于存储能够在处理器上运行的计算机程序的存储器；其中，
所述处理器用于运行所述计算机程序时，执行权利要求1至5任一项所述方法的步骤。
10.一种计算机存储介质，其上存储有计算机程序，其特征在于，该计算机程序被处理器执行时实现权利要求1至5任一项所述方法的步骤。
CN201811290311.0A 2018-10-31 2018-10-31 一种针对非线性调频sar的rd成像方法和装置 Active CN109343056B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811290311.0A CN109343056B (zh) 2018-10-31 2018-10-31 一种针对非线性调频sar的rd成像方法和装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811290311.0A CN109343056B (zh) 2018-10-31 2018-10-31 一种针对非线性调频sar的rd成像方法和装置

Publications (2)

Publication Number Publication Date
CN109343056A true CN109343056A (zh) 2019-02-15
CN109343056B CN109343056B (zh) 2020-05-08

Family

ID=65312806

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811290311.0A Active CN109343056B (zh) 2018-10-31 2018-10-31 一种针对非线性调频sar的rd成像方法和装置

Country Status (1)

Country Link
CN (1) CN109343056B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112363143A (zh) * 2020-11-24 2021-02-12 北京小米移动软件有限公司 一种空调设备基于毫米波进行空间识别的方法和系统
CN116148856A (zh) * 2023-04-17 2023-05-23 南京邮电大学 一种sar动目标二维自聚焦成像处理方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6072419A (en) * 1997-05-13 2000-06-06 Thomson-Csf Method for the processing of the reception signal of a deramp type synthetic aperture radar
CN102129066A (zh) * 2010-01-13 2011-07-20 中国科学院电子学研究所 一种宽幅星载sar快速地理编码方法
CN102147469A (zh) * 2010-12-29 2011-08-10 电子科技大学 一种双基地前视合成孔径雷达成像方法
CN105259552A (zh) * 2015-09-17 2016-01-20 中国科学院电子学研究所 一种基于非线性调频信号的合成孔径雷达成像方法和装置
CN105301588A (zh) * 2015-10-10 2016-02-03 中国测绘科学研究院 一种联合StereoSAR和InSAR的DEM提取方法
CN106093941A (zh) * 2016-06-07 2016-11-09 电子科技大学 一种基于fmcw雷达系统的冰层截面成像方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6072419A (en) * 1997-05-13 2000-06-06 Thomson-Csf Method for the processing of the reception signal of a deramp type synthetic aperture radar
CN102129066A (zh) * 2010-01-13 2011-07-20 中国科学院电子学研究所 一种宽幅星载sar快速地理编码方法
CN102147469A (zh) * 2010-12-29 2011-08-10 电子科技大学 一种双基地前视合成孔径雷达成像方法
CN105259552A (zh) * 2015-09-17 2016-01-20 中国科学院电子学研究所 一种基于非线性调频信号的合成孔径雷达成像方法和装置
CN105301588A (zh) * 2015-10-10 2016-02-03 中国测绘科学研究院 一种联合StereoSAR和InSAR的DEM提取方法
CN106093941A (zh) * 2016-06-07 2016-11-09 电子科技大学 一种基于fmcw雷达系统的冰层截面成像方法

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112363143A (zh) * 2020-11-24 2021-02-12 北京小米移动软件有限公司 一种空调设备基于毫米波进行空间识别的方法和系统
CN112363143B (zh) * 2020-11-24 2023-11-07 北京小米移动软件有限公司 一种空调设备基于毫米波进行空间识别的方法和系统
CN116148856A (zh) * 2023-04-17 2023-05-23 南京邮电大学 一种sar动目标二维自聚焦成像处理方法

Also Published As

Publication number Publication date
CN109343056B (zh) 2020-05-08

Similar Documents

Publication Publication Date Title
CN105259552B (zh) 一种基于非线性调频信号的合成孔径雷达成像方法和装置
CN108051809B (zh) 基于Radon变换的运动目标成像方法、装置及电子设备
EP2998763B1 (en) Phase calibration of a stepped-chirp signal for a synthetic aperture radar
CN107462887B (zh) 基于压缩感知的宽幅星载合成孔径雷达成像方法
US9291711B2 (en) Compressive radar imaging technology
US9329264B2 (en) SAR image formation
CN104316924B (zh) 机载超高分辨率sar反投影图像的自聚焦运动补偿方法
CN106338731B (zh) 一种多子带并发的mimo-sar雷达成像方法
CN103454632B (zh) 一站固定式调频连续波双基地sar成像方法
CN106772253B (zh) 一种非均匀杂波环境下的雷达杂波抑制方法
CN104977582A (zh) 一种实现扫描雷达方位超分辨成像的解卷积方法
CN103605116B (zh) 基于稀疏分析的成像雷达通道参数在线补偿方法
CN109597072A (zh) 一种双基合成孔径雷达sar系统的成像处理方法及装置
CN106950565A (zh) 星载sar成像抖动补偿方法、成像方法
CN109343056A (zh) 一种针对非线性调频sar的rd成像方法和装置
CN111273292B (zh) 合成孔径雷达高频振动补偿方法、装置、电子设备及介质
CN110879391A (zh) 基于电磁仿真和弹载回波仿真的雷达图像数据集制作方法
CN108845318A (zh) 基于Relax算法的星载高分宽幅成像方法
CN109343057A (zh) 一种针对非线性调频sar的cs成像方法和装置
Saeedi et al. Improved navigation-based motion compensation for LFMCW synthetic aperture radar imaging
US9103918B2 (en) Enhanced radar range resolution
CN115755046A (zh) 一种用于条带sar大方位幅宽成像的扩展pfa算法
CN111007512B (zh) 车载雷达成像方法及装置、电子设备
KR20150055279A (ko) 방위각 고분해능 신호처리 알고리즘을 이용하는 차량용 레이더 및 그 운영 방법
CN111638516A (zh) 基于双频共轭处理技术的太赫兹频段sar运动补偿算法

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