CN110967692B - 一种成像方法 - Google Patents
一种成像方法 Download PDFInfo
- Publication number
- CN110967692B CN110967692B CN201911347403.2A CN201911347403A CN110967692B CN 110967692 B CN110967692 B CN 110967692B CN 201911347403 A CN201911347403 A CN 201911347403A CN 110967692 B CN110967692 B CN 110967692B
- Authority
- CN
- China
- Prior art keywords
- electromagnetic wave
- vortex electromagnetic
- scattering
- scattering point
- detection target
- 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
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 40
- 238000001514 detection method Methods 0.000 claims abstract description 53
- 238000013507 mapping Methods 0.000 claims abstract description 7
- 238000007781 pre-processing Methods 0.000 claims abstract description 7
- 239000011159 matrix material Substances 0.000 claims description 21
- 238000006243 chemical reaction Methods 0.000 claims description 20
- 238000012545 processing Methods 0.000 claims description 12
- 238000001228 spectrum Methods 0.000 claims description 12
- 238000005070 sampling Methods 0.000 claims description 5
- 230000009466 transformation Effects 0.000 claims description 4
- 230000007547 defect Effects 0.000 description 7
- 238000000034 method Methods 0.000 description 7
- 230000005855 radiation Effects 0.000 description 7
- 238000011160 research Methods 0.000 description 5
- 238000004891 communication Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 230000033001 locomotion Effects 0.000 description 4
- 230000005540 biological transmission Effects 0.000 description 3
- 238000013461 design Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000000007 visual effect Effects 0.000 description 3
- 238000002592 echocardiography Methods 0.000 description 2
- 230000010287 polarization Effects 0.000 description 2
- BSYNRYMUTXBXSQ-UHFFFAOYSA-N Aspirin Chemical compound CC(=O)OC1=CC=CC=C1C(O)=O BSYNRYMUTXBXSQ-UHFFFAOYSA-N 0.000 description 1
- LFQSCWFLJHTTHZ-UHFFFAOYSA-N Ethanol Chemical compound CCO LFQSCWFLJHTTHZ-UHFFFAOYSA-N 0.000 description 1
- 239000000969 carrier Substances 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000002860 competitive effect Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000007123 defense Effects 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
- 230000008569 process Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
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
-
- 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
- G01S7/418—Theoretical aspects
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
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)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种成像方法，包含：步骤S1、获取雷达探测目标的涡旋电磁波回波信号；步骤S2、对所述涡旋电磁波回波信号进行预处理；步骤S3、根据经所述预处理的涡旋电磁波回波信号，获取所述雷达探测目标的位置信息；步骤S4、根据经所述预处理的涡旋电磁波回波信号，获取所述雷达探测目标的若干个散射点的方位角信息；步骤S5、根据所述位置信息和所述若干个散射点的方位角信息，获得与每一所述散射点相对应的散射点俯仰角信息；步骤S6、将所有所述散射点俯仰角信息映射至三维坐标中，以重构出所述雷达探测目标的立体图像。本发明提出了一种适应不同距离探测的涡旋电磁波三维成像方法。
Description
技术领域
本发明涉及成像技术领域，尤其是涉及一种雷达成像方法。
背景技术
雷达成像技术具备全天时、全天候的优势，不受自然条件的影响，高分辨率成像雷达在战场侦察、防空反导、精确打击等军事领域中的作用逐步突出，其已成为诸多国家未来发展竞相争夺的制高点。
雷达高分辨成像理论与技术历经数十年的发展已取得长足进步，部分成果也在雷达成像系统中得到了应用，统筹来看大多数方法通常都是依据距离-多普勒原理展开，其高分辨率通过发射宽带信号、雷达与目标的相对运动形成大的虚拟合成孔径来获取。其中方位向分辨率的获得取决于雷达与目标的横向相对运动，然而在雷达与目标间的横向相对运动较小的一些应用场景中，传统的SAR/ISAR成像所需的大虚拟孔径无法保证，雷达探测能力受到了限制，因此，探索新的信息获取方式是提升雷达成像分辨率的一条重要途径。
雷达领域的信息调制主要利用远场平面波近似，当电磁波的相位波前具有非平面结构时，可通过在其上调制所需信息来提高信息传递、获取能力。相比于传统电磁波，涡旋电磁波则有着更高维度的信息调制自由度，因其轨道角动量拥有不同于线性动量已有自由度的新自由度，有望在不增加采样点的情况下提升雷达分辨率，并具备抗干扰能力，对于雷达特性测量及目标识别能力的提高具有诱人的应用前景。
从电磁涡旋近况可以看出，它是量子信息技术不可分割的组成部分。轨道角动量具有不同的模式，利用不同的OAM状态可作为不同信息的载体。由于其波前是连续变化的，相当于能在短时间内实现二维角度空间连续采样，比现行的MIMO雷达在空间分集方面具有明显的优势，可望实现更为理想的目标多维成像能力。
目前已有的涡旋电磁波研究主要涉及涡旋电磁波的产生、轨道角动量天线设计、涡旋电磁波模式复用等方面。雷达和通信主要集中的微波波段电磁波也具有涡旋的性质，且目前将涡旋电磁波应用主要集中在通信领域，其在雷达领域的研究尚不多见，且存在诸多弊端。由此，需要提出一种应用在雷达领域的较为完善的基于涡旋电磁波的新型成像方法。
发明内容
本发明的目的在于提供一种成像方法，用以解决现有技术在雷达成像应用领域中没有适应不同探测距离的涡旋电磁波三维成像方法的问题。
为了解决上述问题，本发明通过以下技术方案实现：
一种成像方法，包含：步骤S1、获取雷达探测目标的涡旋电磁波回波信号；步骤S2、对所述涡旋电磁波回波信号进行预处理；步骤S3、根据经所述预处理的涡旋电磁波回波信号，获取所述雷达探测目标的位置信息；步骤S4、根据经所述预处理的涡旋电磁波回波信号，获取所述雷达探测目标的若干个散射点的方位角信息；步骤S5、根据所述位置信息和所述若干个散射点的方位角信息，获得与每一所述散射点相对应的散射点俯仰角信息；步骤S6、将所有所述散射点俯仰角信息映射至三维坐标中，以重构出所述雷达探测目标的立体图像。
进一步的，所述步骤S2包括：根据不同的涡旋电磁波回波信号特性，对所述涡旋电磁波回波信号进行模拟正交解调后直接进行模数变换，得到回波频谱，并将所述回波频谱搬移至预设的处理频带内；
或者，采用中频采样后数字下变频的方式以对所述涡旋电磁波回波信号进行模数转换，得到回波频谱，并将所述回波频谱搬移至预设的处理频带内。
进一步的，所述步骤S2还包括：
将用于接收所述涡旋电磁波回波信号的接收系统设置为基模模式，以线性调频前端去斜方式对经所述模数转换的涡旋电磁波回波信号进行快时间高阶相位项补偿和包络对齐处理。
进一步的，所述步骤S3包括：所述雷达探测目标包括P个散射点，每个散射点后向散射强度表示为σp，忽略初相，则所述涡旋电磁波回波信号srm(t)为：
式中，fp表示散射点的频率，t表示时间脉宽，Lm(.)表示Lagurre多项式，ρ表示散射点到雷达视轴的距离，ω3dB表示发射波束r处照射的3dB波束宽度所对应的束腰波束大小，m表示模式的阶数，
表示散射点的相位；
对所述涡旋电磁波回波信号srm(t)进行快时间维FFT变换，获得不同散射点回波距离对应的频率值；
根据频率与距离关系计算出所述雷达探测目标的位置信息。
进一步的，所述步骤S4包括：采用慢时间维FFT变换可获得每个所述散射点的空间相位分布信息，所述雷达探测目标的第p散射点方位信息Apm采用如下公式表示：
式中，σp表示第P散射点后向散射强度；Lm(.)表示Lagurre多项式，ρ表示散射点到雷达视轴的距离，ω3dB表示发射波束r处照射的3dB波束宽度所对应的束腰波束大小，m表示涡旋电磁波的基模模式的阶数，
表示散射点的相位，α表示涡旋电磁波的轨道角动量的本征值，用于描述涡旋电磁波的轨道角动量的状态；根据三角函数特性，当m≠α时，Apm为零，则可确定该第p散射点的方位角信息。
进一步的，所述步骤S5包括：根据基模模式分布的特征构建同一散射点下的散射点俯仰角和模式幅度的转换矩阵D，可利用Laguerre多项式直接求得；模式数为M，俯仰向可分辨数为N，则所述转换矩阵D为M×N的矩阵，表示为：
根据所述雷达探测目标的第p个散射点方位信息Apm，结合所述转换矩阵D获得对应的散射点俯仰角信息。其中，所述转换矩阵D中的矩阵[A0…Am…AM]表示不同模式散射点对应的回波强度，矩阵[θ0…θn…θN]表示散射点对应模式的俯仰角度，Dmn表示某种俯仰角为θn对应幅度为Am时的系数。
本发明与现有技术相比具有以下优点：
本发明提供了一种成像方法，包含：获取雷达探测目标的涡旋电磁波回波信号；对所述涡旋电磁波回波信号进行预处理；根据经所述预处理的涡旋电磁波回波信号，获取所述雷达探测目标的位置信息；根据经所述预处理的涡旋电磁波回波信号，获取所述雷达探测目标的若干个散射点的方位角信息；根据所述位置信息和所述若干个散射点的方位角信息，获得与每一所述散射点相对应的散射点俯仰角信息；将所有所述散射点俯仰角信息映射至三维坐标中，以重构出所述雷达探测目标的立体图像。由此可知，从已有的发明来开，传统的雷达信息调制主要在时域、频域或极化域上进行，雷达图像距离向分辨率取决于发射信号带宽，方位向分辨率依赖于雷达观测角度范围，而本发明基于涡旋电磁波相位螺旋分布的辐射特性探索雷达信息获取新方式，不但可以突破传统的雷达成像方式主要依托于距离-多普勒而存在的探测场景应用约束，还可以实现成像方位高分辨，本发明可为新体制的雷达设计、目标识别技术的发展提供参考和借鉴。
对于现有已公开的涡旋电磁波相关的成像方法，由于目标俯仰角靠近高阶模态涡旋电磁波主瓣区域，从涡旋电磁波辐射强度方向图可知，在此俯仰角附近较大的范围内，高阶模态回波包络均为正值，应用贝塞尔函数进行相位补偿误差的可能性较小，只能实现近距离、小范围探测。在高模态情况下，无法实现对于多个散射点的高分辨，本发明提出了克服此弊端的处理方式，避免了贝塞尔函数调制导致的螺旋波强度中空的弊端，可以达到扩大探测范围、提高方位向分辨力的目的。
附图说明
图1为本发明一实施例提供的一种成像方法的流程图；
图2为基于不同模式的归一化强度与偏离波束中心角度关系曲线；
图3为基于不同模式的空间二维相位调制分布图。
具体实施方式
承如背景技术所述，目前已有的涡旋电磁波研究主要涉及涡旋电磁波的产生、轨道角动量天线设计、涡旋电磁波模式复用等方面。雷达和通信主要集中的微波波段电磁波也具有涡旋的性质，且目前将涡旋电磁波应用主要集中在通信领域，其在雷达领域的研究尚不多见。例如公开号为CN106199600B，名称为：基于多普勒估计的方位多通道合成孔径雷达成像方法的中国专利提出了一种基于多普勒估计的方位多通道合成孔径雷达成像方法，其属于传统的依据距离-多普勒原理展开的雷达成像方法，雷达探测能力受到了限制。名称为《涡旋电磁波在雷达成像中的应用研究》毕业论文文献中对涡旋电磁波在雷达成像中的应用进行了初步探索，根据涡旋电磁波产生原理，提出了一种用螺旋形阵列产生涡旋电磁波的新方法，通过改变频率控制产生的涡旋电磁波模态。但由于辐射能量弱和模态混叠的原因，大俯仰角并不适合实际应用，无法实现高分辨成像。由此可知，在雷达成像应用领域，至今还没有较为完善的应用涡旋电磁波的三维成像方法，需要提出适应不同距离探测的涡旋电磁波三维成像方法。
下面将结合示意图对本发明的具体实施方式进行更详细的描述。根据下列描述和权利要求书，本发明的优点和特征将更清楚。需说明的是，附图均采用非常简化的形式且均使用非精准的比例，仅用以方便、明晰地辅助说明本发明实施例的目的。
如图1所示，本实施例一种成像方法包括以下过程：
步骤S1、获取雷达探测目标的涡旋电磁波回波信号；
步骤S2、对所述涡旋电磁波回波信号进行预处理；
步骤S3、根据经所述预处理的涡旋电磁波回波信号，获取所述雷达探测目标的位置信息；
步骤S4、根据经所述预处理的涡旋电磁波回波信号，获取所述雷达探测目标的若干个散射点的方位角信息；
步骤S5、根据所述位置信息和所述若干个散射点的方位角信息，获得与每一所述散射点相对应的散射点俯仰角信息；
步骤S6、将所有所述散射点俯仰角信息映射至三维坐标中，以重构出所述雷达探测目标的立体图像。
具体的，在本实施例中，以理想的散射点为例对所述成像方法进行进一步说明:基于幅度模态变换和空间相位调制的工作方式能够获得雷达探测目标的真实三维图像，通过信号处理可获得雷达探测目标所包括的不同散射点的距离、俯仰角度信息和方位角度信息。距离维采用传统的宽带线性调频或者步进频等工作体制能够对距离方向或者雷达视轴方向实现距离高分辨，理论分辨率为c/2B(c为光速，B为信号带宽)；方位维采用多模态空间相位调制方式能够对目标方位向散射点进行分辨，理论分辨率为2π/Δp(Δp为空间相位调制最高模态特征值和最低模态特征值之差)；俯仰维采用同一散射点不同模态回波幅度信息进行解算，理论分辨率为θ3dB/Δp(θ3dB为辐射口径所对应的3dB波束宽度)。
雷达系统发射信号形式采用如下公式进行表示：
式中，stp(t)为雷达系统发射信号形式，t表示时间脉宽，u(t)为信号复包络(可为线性调频或者步进频信号包络)，fc为信号载波频率，A(r)为辐射场在远场r处任一点
的信号，m表示涡旋电磁波的基模模式的阶数，其通常取整数，θ为俯仰角，
为方位角，i表示虚数单位，j表示虚数单位。
式中，Lm(·)为Lagurre多项式，ρ为散射点到雷达视轴的距离，ω3dB为发射波束r处照射的3dB波束宽度所对应的束腰波束大小，k为波数，
为初相。
进一步的，所述步骤S2包括：根据不同的涡旋电磁波回波信号特性，对所述涡旋电磁波回波信号进行模拟正交解调后直接进行模数变换，得到回波频谱，并将所述回波频谱搬移至预设的处理频带内；
或者，采用中频采样后数字下变频的方式以对所述涡旋电磁波回波信号进行模数转换，得到回波频谱，并将所述回波频谱搬移至预设的处理频带内。
进一步的，所述步骤S2还包括：将用于接收所述涡旋电磁波回波信号的接收系统设置为基模模式，由于雷达系统发射波束照射强度分布近似为高斯性，散射点空间相位(空间二维相位)无调制。以线性调频前端去斜方式对经所述模数转换的涡旋电磁波回波信号进行快时间高阶相位项补偿和包络对齐等的预处理。
进一步的，上述涡旋电磁波回波信号经过上述步骤S2的预处理之后，其消除了所述雷达探测目标相对雷达运动带来的空间相位变化的影响。所述步骤S3包括：所述雷达探测目标包括P个散射点，每个散射点后向散射强度表示为σp，忽略初相
则所述涡旋电磁波回波信号srm(t)采用如下公式进行表示：
式中，fp表示第p散射点的频率，t表示时间脉宽，Lm(.)表示Lagurre多项式，ρ表示散射点到雷达视轴的距离，ω3dB表示发射波束r处照射的3dB波束宽度所对应的束腰波束大小，m表示涡旋电磁波的基模模式的阶数，
表示第p散射点的相位。Lagurre多项式表征涡旋电磁波幅度分布。此时如图2所示，其给出了m为0～3阶模时，涡旋电磁波幅度分布情况，表明了该种类型涡旋电磁波幅度无中空现象，由此，能够适应不同距离的目标探测；
表征涡旋电磁波空间相位分布，如图3所示，其给出了m为0～3阶模时的涡旋电磁波的空间二维相位调制分布并表明该种类型涡旋电磁波具备方位向角度分析能力。对式(3)进行快时间维FFT变换，即快速傅里叶变换，获得不同散射点回波距离对应的频率值；根据频率与距离关系计算出所述雷达探测目标的位置信息。
进一步的，所述步骤S4包括：由于空间相位分布具有以2π为周期的周期性，则对于涡旋电磁波回波信号即式(3)用求和形式。由于式(3)中前三项与第四项相互独立，之间不存在干扰或者耦合，则本实施例采用慢时间维FFT变换可获得每个所述散射点的空间相位分布信息，所述雷达探测目标的第p散射点方位信息Apm采用如下公式表示：
式中，σp表示第p散射点后向散射强度，
表示散射点的相位，Lm(.)表示Lagurre多项式，ρ表示散射点到雷达视轴的距离，ω3dB表示发射波束r处照射的3dB波束宽度所对应的束腰波束大小，m表示涡旋电磁波基模模式的阶数(模式的阶数)；α表示涡旋电磁波的轨道角动量的本征值，用于描述涡旋电磁波的轨道角动量的状态。根据三角函数特性，当m≠α时，Apm为零，则可确定该第p散射点的方位角信息。
进一步的，所述步骤S5包括：确定雷达探测目标的距离(即位置信息)和方位角信息后，取出涡旋电磁波回波信号的强度信息。
根据涡旋电磁波的基模模式分布的特征构建同一散射点下的散射点俯仰角和模式幅度的转换矩阵D，可利用Laguerre多项式直接求得；模式数为M，俯仰向可分辨数为N，则所述转换矩阵D为M×N的矩阵，表示为：
转换矩阵D中的矩阵[A0…Am…AM]表示不同模式散射点对应的回波强度，矩阵[θ0…θn…θN]表示散射点对应模式的俯仰角度，Dmn表示某种俯仰角为θn对应幅度为Am时的系数，根据所述雷达探测目标的第p散射点方位信息Apm，结合所述转换矩阵D获得对应的散射点俯仰角信息。
由此可知，本实施例提出了克服此弊端的处理方式，避免了贝塞尔函数调制导致的涡旋电磁波强度中空的弊端，可以达到扩大探测范围、提高方位向分辨力的目的。
然后，将获取的不同散射点的
信息映射至三维坐标中，即可重构出目标的立体图像。由此可知本实施例解决了现有技术在雷达成像应用领域中没有适应不同探测距离的涡旋电磁波三维成像方法的问题。
综上所述，本发明提供了一种成像方法，包含：获取雷达探测目标的涡旋电磁波回波信号；对所述涡旋电磁波回波信号进行预处理；根据经所述预处理的涡旋电磁波回波信号，获取所述雷达探测目标的位置信息；根据经所述预处理的涡旋电磁波回波信号，获取所述雷达探测目标的若干个散射点的方位角信息；根据所述位置信息和所述若干个散射点的方位角信息，获得与每一所述散射点相对应的散射点俯仰角信息；将所有所述散射点俯仰角信息映射至三维坐标中，以重构出所述雷达探测目标的立体图像。由此可知，从已有的发明来开，传统的雷达信息调制主要在时域、频域或极化域上进行，雷达图像距离向分辨率取决于发射信号带宽，方位向分辨率依赖于雷达观测角度范围，而本发明基于涡旋电磁波相位螺旋分布的辐射特性探索雷达信息获取新方式，不但可以突破传统的雷达成像方式主要依托于距离-多普勒而存在的探测场景应用约束，还可以实现成像方位高分辨，本发明可为新体制的雷达设计、目标识别技术的发展提供参考和借鉴。
对于现有已公开的涡旋电磁波相关的成像方法，由于目标俯仰角靠近高阶模态涡旋电磁波主瓣区域，从涡旋电磁波辐射强度方向图可知，在此俯仰角附近较大的范围内，高阶模态回波包络均为正值，应用贝塞尔函数进行相位补偿误差的可能性较小，只能实现近距离、小范围探测。在高模态情况下，无法实现对于多个散射点的高分辨，本发明提出了克服此弊端的处理方式，避免了贝塞尔函数调制导致的涡旋电磁波强度中空的弊端，可以达到扩大探测范围、提高方位向分辨力的目的。
上述仅为本发明的优选实施例而已，并不对本发明起到任何限制作用。任何所属技术领域的技术人员，在不脱离本发明的技术方案的范围内，对本发明揭露的技术方案和技术内容做任何形式的等同替换或修改等变动，均属未脱离本发明的技术方案的内容，仍属于本发明的保护范围之内。
Claims (4)
1.一种成像方法，其特征在于，包含：
步骤S1、获取雷达探测目标的涡旋电磁波回波信号；
步骤S2、对所述涡旋电磁波回波信号进行预处理；
所述步骤S2还包括：将用于接收所述涡旋电磁波回波信号的接收系统设置为基模模式，以前端去斜方式，并对经模数转换的涡旋电磁波回波信号进行快时间高阶相位项补偿和包络对齐处理；
步骤S3、根据经所述预处理的涡旋电磁波回波信号，获取所述雷达探测目标的位置信息；
步骤S4、根据经所述预处理的涡旋电磁波回波信号，获取所述雷达探测目标的若干个散射点的方位角信息；
步骤S5、根据所述位置信息和所述若干个散射点的方位角信息，获得与每一所述散射点相对应的散射点俯仰角信息；
所述步骤S5包括：根据基模模式分布的特征构建同一散射点下的散射点俯仰角和模式幅度的转换矩阵D，可利用Laguerre多项式直接求得；模式数为M，俯仰向可分辨数为N，则所述转换矩阵D为M×N的矩阵，表示为：
根据所述雷达探测目标的第p散射点方位信息Apm，结合所述转换矩阵D获得对应的散射点俯仰角信息；其中，所述转换矩阵D中的矩阵[A0···Am···AM]表示不同模式散射点对应的回波强度，矩阵[θ0···θn···θN]表示散射点对应模式的俯仰角度，Dmn表示某种俯仰角为θn对应幅度为Am时的系数；
步骤S6、将所有所述散射点俯仰角信息映射至三维坐标中，以重构出所述雷达探测目标的立体图像。
2.如权利要求1所述的成像方法，其特征在于，所述步骤S2包括：根据不同的涡旋电磁波回波信号特性，对所述涡旋电磁波回波信号进行模拟正交解调后直接进行模数变换，得到回波频谱，并将所述回波频谱搬移至预设的处理频带内；
或者，采用中频采样后数字下变频的方式以对所述涡旋电磁波回波信号进行模数转换，得到回波频谱，并将所述回波频谱搬移至预设的处理频带内。
3.如权利要求2所述的成像方法，其特征在于，所述步骤S3包括：所述雷达探测目标包括P个散射点，每个散射点后向散射强度表示为σp，忽略初相，则所述涡旋电磁波回波信号srm(t)为：
式中，fp表示散射点的频率，t表示时间脉宽，Lm(.)表示Lagurre多项式，ρ表示散射点到雷达视轴的距离，ω3dB表示发射波束r处照射的3dB波束宽度所对应的束腰波束大小，m表示涡旋电磁波的基模模式的阶数，
表示散射点的相位；
对所述涡旋电磁波回波信号srm(t)进行快时间维FFT变换，获得不同散射点回波距离对应的频率值；
根据频率与距离关系计算出所述雷达探测目标的位置信息。
4.如权利要求3所述的成像方法，其特征在于，所述步骤S4包括：
采用慢时间维FFT变换可获得每个所述散射点的空间相位分布信息，所述雷达探测目标的第p个散射点方位信息Apm采用如下公式表示：
Lm(.)表示Lagurre多项式，ρ表示散射点到雷达视轴的距离，ω3dB表示发射波束r处照射的3dB波束宽度所对应的束腰波束大小，m表示涡旋电磁波的基模模式的阶数，
表示散射点的相位，α表示涡旋电磁波的轨道角动量的本征值，用于描述涡旋电磁波的轨道角动量的状态；根据三角函数特性，当m≠α时，Apm为零，则可确定该第p散射点的方位角信息。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911347403.2A CN110967692B (zh) | 2019-12-24 | 2019-12-24 | 一种成像方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911347403.2A CN110967692B (zh) | 2019-12-24 | 2019-12-24 | 一种成像方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110967692A CN110967692A (zh) | 2020-04-07 |
| CN110967692B true CN110967692B (zh) | 2023-04-11 |
Family
ID=70036228
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911347403.2A Active CN110967692B (zh) | 2019-12-24 | 2019-12-24 | 一种成像方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110967692B (zh) |
Families Citing this family (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111474543B (zh) * | 2020-04-13 | 2022-08-09 | 中国人民解放军国防科技大学 | 一种涡旋电磁波干涉的目标三维成像方法及装置 |
| CN114076958A (zh) * | 2020-08-17 | 2022-02-22 | 华为技术有限公司 | 电磁波成像方法、装置及系统 |
| CN112083432B (zh) * | 2020-09-10 | 2024-01-09 | 天津水聿方舟海洋工程技术有限公司 | 基于声学轨道角动量的超精细三维成像方法 |
| CN112526512B (zh) * | 2020-11-23 | 2022-07-22 | 电子科技大学 | 大功率大口径宽带毫米波空馈相控阵雷达系统及成像方法 |
| CN112327299B (zh) * | 2021-01-04 | 2021-03-16 | 中国人民解放军国防科技大学 | 电磁涡旋前视目标三维成像方法 |
| CN112904347B (zh) * | 2021-01-19 | 2023-05-26 | 鹏城实验室 | 成像系统和方法 |
| CN112764041B (zh) * | 2021-01-19 | 2023-05-26 | 鹏城实验室 | 成像系统和方法 |
| CN113608191B (zh) * | 2021-07-30 | 2024-01-02 | 上海无线电设备研究所 | 一种近场目标方位识别方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106526589A (zh) * | 2016-12-26 | 2017-03-22 | 浙江大学 | 一种基于涡旋电磁波的雷达目标二维成像方法 |
| CN106886020A (zh) * | 2017-03-06 | 2017-06-23 | 中国人民解放军国防科学技术大学 | 一种单天线接收条件下的电磁涡旋成像方法 |
| CN110221293A (zh) * | 2019-05-10 | 2019-09-10 | 电子科技大学 | 一种基于涡旋电磁波的俯仰角成像方法 |
| CN110412571A (zh) * | 2019-07-19 | 2019-11-05 | 西安电子科技大学 | 基于电磁涡旋波的合成孔径雷达三维成像方法 |
| CN110426707A (zh) * | 2019-08-12 | 2019-11-08 | 中国人民解放军国防科技大学 | 基于轨道角动量的涡旋sar成像方法及成像系统 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2017149526A2 (en) * | 2016-03-04 | 2017-09-08 | May Patents Ltd. | A method and apparatus for cooperative usage of multiple distance meters |
-
2019
- 2019-12-24 CN CN201911347403.2A patent/CN110967692B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106526589A (zh) * | 2016-12-26 | 2017-03-22 | 浙江大学 | 一种基于涡旋电磁波的雷达目标二维成像方法 |
| CN106886020A (zh) * | 2017-03-06 | 2017-06-23 | 中国人民解放军国防科学技术大学 | 一种单天线接收条件下的电磁涡旋成像方法 |
| CN110221293A (zh) * | 2019-05-10 | 2019-09-10 | 电子科技大学 | 一种基于涡旋电磁波的俯仰角成像方法 |
| CN110412571A (zh) * | 2019-07-19 | 2019-11-05 | 西安电子科技大学 | 基于电磁涡旋波的合成孔径雷达三维成像方法 |
| CN110426707A (zh) * | 2019-08-12 | 2019-11-08 | 中国人民解放军国防科技大学 | 基于轨道角动量的涡旋sar成像方法及成像系统 |
Non-Patent Citations (4)
| Title |
|---|
| Structure Safety Evaluation of Transmission Corridor Considering Micro-landform Under Super Typhoon Environment;Zebin Wang et al.;《IEEE》;20181231;全文 * |
| 基于单脉冲测角技术的雷达高分辨三维成像方法研究;雷璐;《科技创新导报》;20100321(第09期);全文 * |
| 涡旋电磁波及其在雷达中应用研究进展;刘康等;《电子学报》;20180915(第09期);全文 * |
| 涡旋电磁波的理论与应用研究进展;李龙等;《微波学报》;20180419(第02期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110967692A (zh) | 2020-04-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110967692B (zh) | 一种成像方法 | |
| CN112444811B (zh) | 一种综合mimo雷达和isar的目标检测和成像方法 | |
| CN109581352B (zh) | 一种基于毫米波雷达的超分辨测角系统 | |
| Feng et al. | An extended fast factorized back projection algorithm for missile-borne bistatic forward-looking SAR imaging | |
| CN107390181B (zh) | 一种基于多波束扫描的雷达高分辨成像方法 | |
| CN102156279A (zh) | 基于mimo的双基地雷达地面动目标检测方法 | |
| Sakamoto et al. | Fast imaging method for security systems using ultrawideband radar | |
| CN109597073A (zh) | 一种无人机载微型近距极化干涉合成孔径雷达系统 | |
| CN111239721A (zh) | 车载mimo雷达求熵解速度模糊的方法 | |
| CN102565772B (zh) | 基于sar子孔径序列图像的海洋动态信息提取方法 | |
| CN112834980A (zh) | 一种基于传播算子的涡旋电磁波超分辨测向方法 | |
| CN113917424B (zh) | 一种基于ddma与盲源分离的地基mimo雷达地杂波抑制方法 | |
| Chen et al. | Very High-Resolution Synthetic Aperture Radar Systems and Imaging: A Review | |
| CN212845925U (zh) | 一种船载sar系统 | |
| Kari et al. | Evolutionary developments of today’s remote sensing radar technology—right from the telemobiloscope: a review | |
| CN113359196A (zh) | 基于子空间法和dbf的多目标生命体征探测方法 | |
| Zhang et al. | Design and preliminary experiment of China imaging altimeter | |
| Zhao et al. | A novel method of ship detection by combining space-borne SAR and GNSS-R | |
| Ma et al. | High-resolution imaging using a narrowband MIMO radar system | |
| Hellsten | CARABAS-an UWB low frequency SAR | |
| CN116500550B (zh) | 一种星载sar距离模糊抑制方法 | |
| CN118671763B (zh) | 一种实现回波压缩的频扫sar信号生成方法 | |
| Zhao et al. | A Novel Beam-Steering Method for Spaceborne Interferometric Synthetic Aperture Radar Altimeter | |
| Chen et al. | Multi-target super-resolution technology of remote sensing based on L regularization | |
| JP2019158356A (ja) | 地形計測方法および地形計測装置 |
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 |