CN109343056A - 一种针对非线性调频sar的rd成像方法和装置 - Google Patents
一种针对非线性调频sar的rd成像方法和装置 Download PDFInfo
- 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
- distance
- range
- pulse
- 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
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 54
- 238000007906 compression Methods 0.000 claims abstract description 66
- 238000000034 method Methods 0.000 claims abstract description 24
- 238000013508 migration Methods 0.000 claims abstract description 7
- 230000005012 migration Effects 0.000 claims abstract description 7
- 238000012937 correction Methods 0.000 claims abstract description 6
- 230000009466 transformation Effects 0.000 claims abstract description 3
- 230000006835 compression Effects 0.000 claims description 61
- 238000012545 processing Methods 0.000 claims description 59
- 230000015654 memory Effects 0.000 claims description 14
- 238000004590 computer program Methods 0.000 claims description 12
- 230000001131 transforming effect Effects 0.000 claims description 10
- 238000004422 calculation algorithm Methods 0.000 abstract description 9
- 230000008569 process Effects 0.000 abstract description 5
- 238000010586 diagram Methods 0.000 description 13
- 230000003287 optical effect Effects 0.000 description 4
- 238000005070 sampling Methods 0.000 description 4
- 230000005540 biological transmission Effects 0.000 description 3
- 238000001228 spectrum Methods 0.000 description 3
- 230000006870 function Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 101000977638 Homo sapiens Immunoglobulin superfamily containing leucine-rich repeat protein Proteins 0.000 description 1
- 102100023538 Immunoglobulin superfamily containing leucine-rich repeat protein Human genes 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 238000005311 autocorrelation function Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000002595 magnetic resonance imaging Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000001303 quality assessment method Methods 0.000 description 1
- 238000005316 response function Methods 0.000 description 1
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
技术领域
本发明涉及合成孔径雷达(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任一项所述方法的步骤。
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 (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112363143A (zh) * | 2020-11-24 | 2021-02-12 | 北京小米移动软件有限公司 | 一种空调设备基于毫米波进行空间识别的方法和系统 |
| CN113433543A (zh) * | 2020-03-05 | 2021-09-24 | 上海新微技术研发中心有限公司 | 一种检测人体久坐状态的智能座位牌及检测方法 |
| CN116148856A (zh) * | 2023-04-17 | 2023-05-23 | 南京邮电大学 | 一种sar动目标二维自聚焦成像处理方法 |
Citations (6)
| 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雷达系统的冰层截面成像方法 |
-
2018
- 2018-10-31 CN CN201811290311.0A patent/CN109343056B/zh active Active
Patent Citations (6)
| 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 (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113433543A (zh) * | 2020-03-05 | 2021-09-24 | 上海新微技术研发中心有限公司 | 一种检测人体久坐状态的智能座位牌及检测方法 |
| 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) | 一种基于非线性调频信号的合成孔径雷达成像方法和装置 | |
| CN111142105B (zh) | 复杂运动目标isar成像方法 | |
| US8044846B1 (en) | Method for deblurring radar range-doppler images | |
| CN108051809B (zh) | 基于Radon变换的运动目标成像方法、装置及电子设备 | |
| US9291711B2 (en) | Compressive radar imaging technology | |
| CN109343056B (zh) | 一种针对非线性调频sar的rd成像方法和装置 | |
| CN108051812B (zh) | 基于二维速度搜索的星载sar运动目标检测方法 | |
| EP2998763A1 (en) | Phase calibration of a stepped-chirp signal for a synthetic aperture radar | |
| CN103454632B (zh) | 一站固定式调频连续波双基地sar成像方法 | |
| CN106932778B (zh) | 方位多通道fmcw sar滑动聚束成像方法 | |
| CN102393518A (zh) | 一种适用于大斜视角的机载sar成像方法 | |
| US8760340B2 (en) | Processing radar return signals to detect targets | |
| CN108562884A (zh) | 一种基于最大后验概率的机载前视海面目标角超分辨方法 | |
| CN106338731A (zh) | 一种多子带并发的mimo‑sar雷达成像方法 | |
| CN113702974B (zh) | 一种机载/弹载合成孔径雷达图像快速优化方法 | |
| CN112882031A (zh) | 一种离散频率编码波形的isar成像方法、装置及存储介质 | |
| CN113640807B (zh) | 多子阵合成孔径声纳脉内多普勒频移补偿逐线成像方法 | |
| CN109343057B (zh) | 一种针对非线性调频sar的cs成像方法和装置 | |
| CN117826155A (zh) | 一种车载雷达成像信号处理方法、设备及介质 | |
| CN112415512B (zh) | 基于进退法和黄金分割法的sar运动目标聚焦方法 | |
| CN116953684B (zh) | 频率捷变fmcw雷达目标速度估计方法及系统 | |
| CN115755046A (zh) | 一种用于条带sar大方位幅宽成像的扩展pfa算法 | |
| CN111007512B (zh) | 车载雷达成像方法及装置、电子设备 | |
| CN104076361B (zh) | 无人机载战场监视雷达超分辨广域成像方法 | |
| CN114488054B (zh) | 一种计算高效的合成孔径雷达地面运动目标聚焦方法 |
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 |