CN114624707A - 基于柱面扫描体制下的mimo阵列三维成像方法及成像装置 - Google Patents

基于柱面扫描体制下的mimo阵列三维成像方法及成像装置 Download PDF

Info

Publication number
CN114624707A
CN114624707A CN202210187239.9A CN202210187239A CN114624707A CN 114624707 A CN114624707 A CN 114624707A CN 202210187239 A CN202210187239 A CN 202210187239A CN 114624707 A CN114624707 A CN 114624707A
Authority
CN
China
Prior art keywords
target
wave number
radius
cylindrical
echo signal
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.)
Pending
Application number
CN202210187239.9A
Other languages
English (en)
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.)
National Space Science Center of CAS
Original Assignee
National Space Science Center 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 National Space Science Center of CAS filed Critical National Space Science Center of CAS
Priority to CN202210187239.9A priority Critical patent/CN114624707A/zh
Publication of CN114624707A publication Critical patent/CN114624707A/zh
Pending legal-status Critical Current

Links

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/28—Details of pulse systems
    • G01S7/285—Receivers
    • G01S7/295—Means for transforming co-ordinates or for evaluating data, e.g. using computers
    • 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/35—Details of non-pulse systems
    • G01S7/352—Receivers
    • 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

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

本发明属于MIMO雷达信号处理和成像技术领域，具体涉及一种基于柱面扫描体制下的MIMO阵列三维成像方法，包括：根据柱面扫描几何下，MIMO阵列中的发射天线阵元与目标之间的距离，以及MIMO阵列中的接收天线阵元与目标之间的距离，得到距离衰减因子，构建考虑距离衰减的回波模型；对所述回波模型相对于发射天线阵元位置、接收天线阵元位置及方位角度变量分别作一维傅里叶变换，得到波数域回波信号；根据频谱回波信号，建立与重建目标半径相关的聚焦函数表达式；对所述聚焦函数表达式求逆后执行一维傅里叶变换，得到波数域聚焦函数；根据得到的波数域回波信号和波数域聚焦函数，将二者相乘后，再进行降维重排、累积和逆傅里叶变换处理，得到目标反射率图像。

Description

基于柱面扫描体制下的MIMO阵列三维成像方法及成像装置
技术领域
本发明属于MIMO雷达信号处理和成像技术领域，具体地说，涉及一种基于柱面扫描体制下的MIMO阵列三维成像方法及成像装置。
背景技术
柱面扫描体制可实现对目标360度的观测，在安检及医学等领域可提供全方位的成像结果。
近年来，多输入多输出(Multiple-Input MultipleOutput,MIMO)雷达技术由于其诸多优点得到快速发展。其中，MIMO阵列相比于单站阵列，可大幅度的减少阵元数量，具有高分辨率成像的潜力。柱面扫描体制下的MIMO阵列是将MIMO阵列与柱面扫描方式相结合，可实现对目标全方位高分辨的重建。MIMO体制下的数据录取方式相对于单站更为复杂，因此，对该几何在近距离场景下的三维成像方法需要进一步的研究。
传统的距离徙动算法(Range Migration Algorithm,RMA)是基于单站数据获取模式下的成像算法，不能直接应用于MIMO阵列信号的处理。在近场条件下，由于波前是球面而不是平面的，因此，在发射和接收天线之间假设虚拟相位中心的成像算法对于近场MIMO重建是不准确的。
全息成像系统的多数重建算法为了简化推导过程忽略了距离衰减幅度项，从而推导过程引入了一定的近似，而对于近距离成像场景，忽略距离衰减的成像算法将在一定程度上影响成像质量。现有的成像算法在构建回波模型时大多忽略距离衰减幅度项，而忽略该项会一定程度上影响成像质量。
发明内容
为了解决现有技术存在的问题，本发明提出了一种基于柱面扫描体制下的MIMO阵列三维成像方法，该方法包括：
根据柱面扫描几何下，MIMO阵列中的发射天线阵元与目标之间的距离，以及MIMO阵列中的接收天线阵元与目标之间的距离，得到距离衰减因子，构建考虑距离衰减的回波模型；
对所述回波模型相对于发射天线阵元位置、接收天线阵元位置及方位角度变量分别作一维傅里叶变换，得到波数域回波信号；根据频谱回波信号，建立与重建目标半径相关的聚焦函数表达式；对所述聚焦函数表达式求逆后执行一维傅里叶变换，得到波数域聚焦函数；
根据得到的波数域回波信号和波数域聚焦函数，将二者相乘后，再进行降维重排、累积和逆傅里叶变换处理得到目标反射率图像。
作为上述技术方案的改进之一，所述MIMO阵列是由多个发射天线阵元和多个接收天线阵元沿某一方向直线排列而成，并以半径R0作圆周运动，以实现柱面扫描，对柱面内目标区域内的目标实现观测。
作为上述技术方案的改进之一，所述根据柱面扫描几何下，MIMO阵列中的发射天线阵元与目标之间的距离，以及MIMO阵列中的接收天线阵元与目标之间的距离，得到距离衰减因子，构建考虑距离衰减的回波模型；其具体过程包括：
根据柱面扫描几何下，发射天线阵元之间的距离RT和接收天线阵元与目标之间的距离RR，得到距离衰减因子，记为1/(4πRTRR)；
其中，
Figure BDA0003523214120000021
其中，(r,θ,z)为预先建立的柱坐标系O-rθz区域中的任一点；
Figure BDA0003523214120000022
Figure BDA0003523214120000023
分别为MIMO阵列中发射天线阵元和接收天线阵元的柱坐标系下的坐标，其中，R0为MIMO阵列柱面扫描的半径；
Figure BDA0003523214120000024
为柱面扫描过程中MIMO阵列同X轴正向的夹角；zT、zR分别代表发射天线阵元在MIMO阵列方向的坐标和接收天线阵元在MIMO阵列方向的坐标；
根据上述得到的距离衰减因子1/(4πRTRR)，构建考虑距离衰减的回波模型：
Figure BDA0003523214120000025
其中，
Figure BDA0003523214120000026
为快时间波数域回波信号；其中，kR为快时间波数；σ(r,θ,z)表示待重建的目标反射率函数；j表示虚部单元；(r,θ,z)为柱坐标系O-rθz区域中的任一点。
作为上述技术方案的改进之一，所述对回波模型相对于发射天线阵元位置、接收天线阵元位置及方位角度变量分别作一维傅里叶变换，得到波数域回波信号；其具体过程包括：
根据回波模型，获取快时间波数域回波信号
Figure BDA0003523214120000031
对快时间波数域回波信号
Figure BDA0003523214120000032
两侧相对于发射天线阵元位置变量zT、接收天线阵元位置变量zR分别执行一维傅里叶变换，得到变换后的频谱回波信号
Figure BDA0003523214120000033
Figure BDA0003523214120000034
其中，kzT、kzR分别为代表变量zT、zR傅里叶变化后对应的波数域变量；σ(r,θ,z)表示待重建的目标反射率函数；
Figure BDA0003523214120000035
为相对于发射天线阵元位置变量zT执行一维傅里叶变换操作；
Figure BDA0003523214120000036
为相对于接收天线阵元位置变量zR执行一维傅里叶变换操作；
利用驻定相位原理求解上式
Figure BDA0003523214120000037
中
Figure BDA0003523214120000038
Figure BDA0003523214120000039
Figure BDA00035232141200000310
Figure BDA00035232141200000311
根据
Figure BDA0003523214120000041
的结果得到另一表示形式下的频谱回波信号
Figure BDA0003523214120000042
Figure BDA0003523214120000043
其中，R0为MIMO阵列柱面扫描的半径；σ(r,θ,z)为待重建的目标反射率函数；
对
Figure BDA0003523214120000044
中的角度变量
Figure BDA0003523214120000045
执行一维傅里叶变换的操作，得到波数域回波信号SSS(kR,kθ,kzT,kzR)；
其中，kθ为角度变量
Figure BDA0003523214120000046
对应的波数域变量。
作为上述技术方案的改进之一，根据所述频谱回波信号，建立与重建目标半径相关的聚焦函数表达式；对所述聚焦函数表达式求逆后执行一维傅里叶变换，得到波数域聚焦函数；包括：
根据所述频谱回波信号
Figure BDA0003523214120000047
建立半径r处聚焦函数表达式为：
Figure BDA0003523214120000048
对所述半径r处聚焦函数表达式求逆后相对于角度变量
Figure BDA0003523214120000049
执行一维傅里叶变换，得到半径r处的波数域聚焦函数
Figure BDA00035232141200000410
作为上述技术方案的改进之一，所述根据得到的波数域回波信号和波数域聚焦函数，将二者相乘后，再进行降维重排、累积和逆傅里叶变换处理得到目标反射率图像；其具体过程包括：
将所述波数域回波信号与某半径处波数域聚焦函数表达式相乘得到聚焦后的波数域回波信号；对所述聚焦后的波数域回波信号执行降维重排操作，得到降维处理后的信号；
对所述降维处理后的信号在快时间波数上进行累积操作得到半径r处的柱面目标的二维频谱；
对所述半径r处的柱面目标的二维频谱执行角度维和阵列维的逆傅里叶变换，得到所述半径r处的柱面目标的重建图像；
对每个半径处的柱面目标重复上述过程得到每个半径位置处的柱面目标图像，即可得到三维目标区域的反射率图像。
作为上述技术方案的改进之一，所述将所述波数域回波信号与某半径处波数域聚焦函数表达式相乘得到聚焦后的波数域回波信号；对所述聚焦后的波数域回波信号执行降维重排操作，得到降维处理后的信号；其具体过程包括：
将所述波数域回波信号SSS(kR,kθ,kzT,kzR)与半径r处波数域聚焦函数
Figure BDA0003523214120000051
相乘，得到聚焦后的波数域回波信号表达式：
Figure BDA0003523214120000052
对所述聚焦后的波数域回波信号SSSr(kR,kθ,kzT,kzR)执行降维重排操作：按照设定的色散关系kz＝kzT+kzR对该信号进行降维重排操作，得到降维处理后的信号
Figure BDA0003523214120000053
其中，kz为定义的阵列向波数；
Figure BDA0003523214120000054
作为上述技术方案的改进之一，所述对所述降维处理后的信号在快时间波数上进行累积操作得到半径r处的柱面目标的二维频谱；对所述半径r处的柱面目标的二维频谱执行角度维和阵列维的逆傅里叶变换，得到所述半径r处的柱面目标的重建图像；其具体过程包括：
对所述降维处理后的信号
Figure BDA0003523214120000055
在快时间波数kR上进行累积操作得到半径处的柱面目标的二维频谱表达式：
Figure BDA0003523214120000056
对所述半径r处的柱面目标的二维频谱Σr(kθ,kz)执行角度维和阵列维的逆傅里叶变换，得到半径r处的柱面目标的重建图像
Figure BDA0003523214120000061
Figure BDA0003523214120000062
本发明还提供了一种基于柱面扫描体制下的MIMO阵列三维成像装置，该装置包括：回波模型构建模块、信号获取模块和图像三维重建模块；
所述回波模型构建模块，用于根据柱面扫描几何下，MIMO阵列中的发射天线阵元与目标之间的距离，以及MIMO阵列中的接收天线阵元与目标之间的距离，得到距离衰减因子，构建考虑距离衰减的回波模型；
所述信号获取模块，所述对回波模型相对于发射天线阵元位置、接收天线阵元位置及方位角度变量分别作一维傅里叶变换，得到波数域回波信号；根据频谱回波信号，建立与重建目标半径相关的聚焦函数表达式；对所述聚焦函数表达式求逆后执行一维傅里叶变换，得到波数域聚焦函数；
所述图像三维重建模块，用于根据得到的波数域回波信号和波数域聚焦函数，将二者相乘后，再进行降维重排、累积和逆傅里叶变换处理，得到目标反射率图像。
本发明与现有技术相比的有益效果是：
1、本发明的方法可有效解决对近场场景下目标的三维成像问题，与传统的时域成像算法相比大大提高的成像的效率；
2、本发明的方法补偿了波方程中的固有因子——距离衰减项，可提高重建目标图像的动态范围。
附图说明
图1为本发明的一种用于柱面扫描体制下的MIMO阵列三维成像方法的流程图；
图2是本发明的柱面扫描体制下的MIMO阵列成像几何的示意图；
图3a是本发明的方法在仿真实验中的MIMO阵列的排布示意图；
图3b是本发明的方法在仿真实验中的仿真示意图；
图4是本发明MIMO阵列柱面扫描系统下对点目标的三维成像结果；
图5a是三维成像结果的二维切面图像中的xy平面的切片图；
图5b是三维成像结果的二维切面图像中的xz方向的切片图；
图6是本发明的一种基于柱面扫描体制下的MIMO阵列三维成像系统的结构示意图。
具体实施方式
现结合附图对本发明作进一步的描述。
如图1所示，本发明提供了一种基于柱面扫描体制下的MIMO阵列三维成像方法，该方法适用于MIMO阵列在近场场景下对柱面孔径内目标进行高分辨三维成像，采用逐个二维柱面成像堆叠的方法实现了目标的三维成像，算法过程中补偿了距离衰减项，提高了图像的质量；在建立回波模型时考虑了距离衰减幅度项，因此本发明的成像方法会一定程度上提高算法的准确度。该方法包括：
在MIMO阵列柱面扫描几何下，构建了考虑距离衰减的回波模型；
对所述回波模型相对于发射天线阵元位置、接收天线阵元位置及方位角度变量分别作一维傅里叶变换，得到波数域回波信号；
根据频谱回波信号，建立与重建目标半径相关的聚焦函数表达式；对所述聚焦函数表达式求逆后执行一维傅里叶变换，得到波数域聚焦函数；
重复上述操作，遍历成像区域中所有半径即可得到完整的三维反射率图像。
如图1所示，该方法具体包括：
根据柱面扫描几何下，MIMO阵列中的发射天线阵元与目标之间的距离，以及MIMO阵列中的接收天线阵元与目标之间的距离，得到距离衰减因子，构建考虑距离衰减的回波模型；
具体地，根据柱面扫描几何下，发射天线阵元之间的距离RT和接收天线阵元与目标之间的距离RR，得到距离衰减因子，记为1/(4πRTRR)；
其中，
Figure BDA0003523214120000071
其中，(r,θ,z)为预先建立的柱坐标系O-rθz区域中的任一点；
Figure BDA0003523214120000072
Figure BDA0003523214120000073
分别为MIMO阵列中发射天线阵元和接收天线阵元的柱坐标系下的坐标，其中，R0为MIMO阵列柱面扫描的半径；
Figure BDA0003523214120000074
为柱面扫描过程中MIMO阵列同X轴正向的夹角；zT、zR分别代表发射天线阵元在MIMO阵列方向的坐标和接收天线阵元在MIMO阵列方向的坐标；
根据上述得到的距离衰减因子1/(4πRTRR)，构建考虑距离衰减的回波模型：
Figure BDA0003523214120000081
其中，
Figure BDA0003523214120000082
为快时间波数域回波信号；其中，kR为快时间波数；σ(r,θ,z)表示待重建的目标反射率函数；j表示虚部单元；(r,θ,z)为柱坐标系O-rθz区域中的任一点。
对所述回波模型相对于发射天线阵元位置、接收天线阵元位置及方位角度变量分别作一维傅里叶变换，得到波数域回波信号；
具体地，根据回波模型，获取快时间波数域回波信号
Figure BDA0003523214120000083
对快时间波数域回波信号
Figure BDA0003523214120000084
两侧相对于发射天线阵元位置变量zT、接收天线阵元位置变量zR分别执行一维傅里叶变换，得到变换后的频谱回波信号
Figure BDA0003523214120000085
Figure BDA0003523214120000086
其中，kzT、kzR分别为代表变量zT、zR傅里叶变化后对应的波数域变量；σ(r,θ,z)表示待重建的目标反射率函数；FTzT为相对于发射天线阵元位置变量zT执行一维傅里叶变换操作；FTzR为相对于接收天线阵元位置变量zR执行一维傅里叶变换操作；
利用驻定相位原理求解上式
Figure BDA0003523214120000087
中
Figure BDA0003523214120000088
Figure BDA0003523214120000089
Figure BDA0003523214120000091
Figure BDA0003523214120000092
根据
Figure BDA0003523214120000093
的结果得到另一表示形式下的频谱回波信号
Figure BDA0003523214120000094
Figure BDA0003523214120000095
其中，R0为MIMO阵列柱面扫描的半径；σ(r,θ,z)为待重建的目标反射率函数；
对
Figure BDA0003523214120000096
中的角度变量
Figure BDA0003523214120000097
执行一维傅里叶变换的操作，得到波数域回波信号
Figure BDA0003523214120000098
其中，kθ为角度变量
Figure BDA0003523214120000099
对应的波数域变量。
根据频谱回波信号，建立与重建目标半径相关的聚焦函数表达式；对所述聚焦函数表达式求逆后执行一维傅里叶变换，得到波数域聚焦函数；
具体地，根据所述频谱回波信号
Figure BDA00035232141200000910
建立半径r处聚焦函数表达式为：
Figure BDA0003523214120000101
对所述半径r处聚焦函数表达式求逆后相对于角度变量
Figure BDA0003523214120000102
执行一维傅里叶变换，得到半径r处的波数域聚焦函数
Figure BDA0003523214120000103
根据得到的波数域回波信号和波数域聚焦函数，将二者相乘后，再进行降维重排、累积和逆傅里叶变换处理，得到目标反射率图像。
具体地，将所述波数域回波信号与某半径处波数域聚焦函数表达式相乘得到聚焦后的波数域回波信号；对所述聚焦后的波数域回波信号执行降维重排操作，得到降维处理后的信号；
具体地，将所述波数域回波信号SSS(kR,kθ,kzT,kzR)与半径r处波数域聚焦函数
Figure BDA0003523214120000104
相乘，得到聚焦后的波数域回波信号表达式：
Figure BDA0003523214120000105
对所述聚焦后的波数域回波信号SSSr(kR,kθ,kzT,kzR)执行降维重排操作：按照设定的色散关系kz＝kzT+kzR进对该信号进行降维重排操作得到降维处理后的信号
Figure BDA0003523214120000106
其中，kz为定义的阵列向波数；
Figure BDA0003523214120000107
对所述降维处理后的信号在快时间波数上进行累积操作得到半径r处的柱面目标的二维频谱；
具体地，对所述降维处理后的信号
Figure BDA0003523214120000108
在快时间波数kR上进行累积操作得到半径处的柱面目标的二维频谱表达式：
Figure BDA0003523214120000109
对所述半径r处的柱面目标的二维频谱执行角度维和阵列维的逆傅里叶变换，得到所述半径r处的柱面目标的重建图像；
对每个半径处的柱面目标重复上述相乘、降维重排、累积、逆傅里叶变换等过程得到每个半径位置处的柱面目标图像，即可得到三维目标区域的反射率图像。
具体地，对所述半径r处的柱面目标的二维频谱Σr(kθ,kz)执行角度维和阵列维的逆傅里叶变换，得到半径r处的柱面目标的重建图像
Figure BDA0003523214120000111
Figure BDA0003523214120000112
遍历目标区域，划分N个柱面单元，对每个半径位置重复上述相乘、降维重排、累积、逆傅里叶变换等过程得到每个半径位置处的柱面目标图像，得到三维目标区域的反射率图像
Figure BDA0003523214120000113
其中，r＝rn，n＝1,2,...N。
如图2所示，所述MIMO阵列是由多个发射天线阵元和多个接收天线阵元沿某一方向直线排列而成，并以半径R0作圆周运动，以实现柱面扫描，对柱面内目标区域内的目标实现观测。
仿真实验
1、仿真参数
为了验证本发明方法的有效性，此处给出了表1中的仿真数据参数。
表1仿真数据参数
Figure BDA0003523214120000114
2、实验内容
在本实施例中，所述MIMO阵列是由76个均匀分布的接收天线阵元和12个非均匀分布的发射天线阵元组成，如图3a所示。应用的MIMO阵列和柱面扫描系统的几何示意如图3b所示。当MIMO阵列位于某一方位角度位置时，发射天线阵元依次发射探测信号，对每一个发射天线阵元有全部的接收天线阵元对目标的回波信号进行接收。当所有发射天线阵元完成信号发射后，阵列运动到下一方位位置处，重复上述发射接收过程。最终整个过程会形成360度的柱面孔径。
利用本发明的方法所获得的实验处理结果如图4、图5所示。图4为原点处，点目标的三维成像效果，图5a为散射点在xy平面的切面，图5b为散射点在xz平面的切面。从图5a和5b中可以看出本发明方法能准确重建目标位置的成像效果。
综上所述，实验验证了本发明的正确性、有效性和可靠性。
在另一个实施例中，如图6所示，本发明提供了一种基于柱面扫描体制下的MIMO阵列三维成像装置，该装置包括：回波模型构建模块1102，信号获取模块1104和图像三维重建模块1106；
所述回波模型构建模块1102，用于根据柱面扫描几何下，MIMO阵列中的发射天线阵元与目标之间的距离，以及MIMO阵列中的接收天线阵元与目标之间的距离，得到距离衰减因子，构建考虑距离衰减的回波模型；
具体地，所述回波模型构建模块1102包括：距离衰减因子获取单元和模型构建单元；
所述距离衰减因子获取单元，用于根据柱面扫描几何下，发射天线阵元之间的距离RT和接收天线阵元与目标之间的距离RR，得到距离衰减因子，记为1/(4πRTRR)；
其中，
Figure BDA0003523214120000121
其中，(r,θ,z)为预先建立的柱坐标系O-rθz区域中的任一点；
Figure BDA0003523214120000122
Figure BDA0003523214120000123
分别为MIMO阵列中发射天线阵元和接收天线阵元的柱坐标系下的坐标，其中，R0为MIMO阵列柱面扫描的半径；
Figure BDA0003523214120000124
为柱面扫描过程中MIMO阵列同X轴正向的夹角；zT、zR分别代表发射天线阵元在MIMO阵列方向的坐标和接收天线阵元在MIMO阵列方向的坐标；
所述模型构建单元，用于根据上述得到的距离衰减因子1/(4πRTRR)，构建考虑距离衰减的回波模型：
Figure BDA0003523214120000125
其中，
Figure BDA0003523214120000126
为快时间波数域回波信号；其中，kR为快时间波数；σ(r,θ,z)表示待重建的目标反射率函数；j表示虚部单元；(r,θ,z)为柱坐标系O-rθz区域中的任一点。
所述信号获取模块1104，用于对回波模型相对于发射天线阵元位置、接收天线阵元位置及方位角度变量分别作一维傅里叶变换，得到波数域回波信号；根据频谱回波信号，建立与重建目标半径相关的聚焦函数表达式；对所述聚焦函数表达式求逆后执行一维傅里叶变换，得到波数域聚焦函数；
具体地，所述信号获取模块1104包括：波数域回波信号获取单元和波数域聚焦函数获取单元；
所述波数域回波信号获取单元，用于根据回波模型，获取快时间波数域回波信号
Figure BDA0003523214120000131
对快时间波数域回波信号
Figure BDA0003523214120000132
两侧相对于发射天线阵元位置变量zT、接收天线阵元位置变量zR分别执行一维傅里叶变换，得到变换后的频谱回波信号
Figure BDA0003523214120000133
Figure BDA0003523214120000134
其中，kzT、kzR分别为代表变量zT、zR傅里叶变化后对应的波数域变量；σ(r,θ,z)表示待重建的目标反射率函数；
Figure BDA0003523214120000135
为相对于发射天线阵元位置变量zT执行一维傅里叶变换操作；
Figure BDA0003523214120000136
为相对于接收天线阵元位置变量zR执行一维傅里叶变换操作；
利用驻定相位原理求解上式
Figure BDA0003523214120000137
中
Figure BDA0003523214120000138
Figure BDA0003523214120000139
Figure BDA00035232141200001310
Figure BDA0003523214120000141
根据
Figure BDA0003523214120000142
的结果得到另一表示形式下的频谱回波信号
Figure BDA0003523214120000143
Figure BDA0003523214120000144
其中，R0为MIMO阵列柱面扫描的半径；σ(r,θ,z)为待重建的目标反射率函数；
对
Figure BDA0003523214120000145
中的角度变量
Figure BDA0003523214120000146
执行一维傅里叶变换的操作，得到波数域回波信号SSS(kR,kθ,kzT,kzR)；
其中，kθ为角度变量
Figure BDA0003523214120000147
对应的波数域变量。
所述波数域聚焦函数获取单元，用于根据所述频谱回波信号
Figure BDA0003523214120000148
建立半径r处聚焦函数表达式为：
Figure BDA0003523214120000149
对所述半径r处聚焦函数表达式求逆后相对于角度变量
Figure BDA00035232141200001411
执行一维傅里叶变换，得到半径r处的波数域聚焦函数
Figure BDA00035232141200001410
所述图像三维重建模块1106，用于根据得到的波数域回波信号和波数域聚焦函数，将二者相乘后，再进行降维重排、累积和逆傅里叶变换处理，得到目标反射率图像。
具体地，所述图像三维重建模块1106包括：降维处理单元、累积处理单元和图像重建单元；
所述降维处理单元，用于将所述波数域回波信号与某半径处波数域聚焦函数表达式相乘得到聚焦后的波数域回波信号；对所述聚焦后的波数域回波信号执行降维重排操作，得到降维处理后的信号；
将所述波数域回波信号SSS(kR,kθ,kzT,kzR)与半径r处波数域聚焦函数
Figure BDA0003523214120000151
相乘，得到聚焦后的波数域回波信号表达式：
Figure BDA0003523214120000152
对所述聚焦后的波数域回波信号SSSr(kR,kθ,kzT,kzR)执行降维重排操作：按照设定的色散关系kz＝kzT+kzR进对该信号进行降维重排操作得到降维处理后的信号
Figure BDA0003523214120000153
其中，kz为定义的阵列向波数；
Figure BDA0003523214120000154
所述累积处理单元，用于对所述降维处理后的信号在快时间波数上进行累积操作得到半径r处的柱面目标的二维频谱；
具体地，对所述降维处理后的信号
Figure BDA0003523214120000155
在快时间波数kR上进行累积操作得到半径处的柱面目标的二维频谱表达式：
Figure BDA0003523214120000156
所述图像重建单元，用于对所述半径r处的柱面目标的二维频谱执行角度维和阵列维的逆傅里叶变换，得到所述半径r处的柱面目标的重建图像；
对每个半径处的柱面目标重复上述过程得到每个半径位置处的柱面目标图像，即可得到三维目标区域的反射率图像。
具体地，对所述半径r处的柱面目标的二维频谱Σr(kθ,kz)执行角度维和阵列维的逆傅里叶变换，得到半径r处的柱面目标的重建图像
Figure BDA0003523214120000157
Figure BDA0003523214120000158
遍历目标区域，划分N个柱面单元，对每个半径位置重复上述过程得到每个半径位置处的柱面目标图像，得到三维目标区域的反射率图像
Figure BDA0003523214120000161
其中，r＝rn，n＝1,2,...N。
最后所应说明的是，以上实施例仅用以说明本发明的技术方案而非限制。尽管参照实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，对本发明的技术方案进行修改或者等同替换，都不脱离本发明技术方案的精神和范围，其均应涵盖在本发明的权利要求范围当中。

Claims (9)

1.一种基于柱面扫描体制下的MIMO阵列三维成像方法，该方法包括：
根据柱面扫描几何下，MIMO阵列中的发射天线阵元与目标之间的距离，以及MIMO阵列中的接收天线阵元与目标之间的距离，得到距离衰减因子，构建考虑距离衰减的回波模型；
对所述回波模型相对于发射天线阵元位置、接收天线阵元位置及方位角度变量分别作一维傅里叶变换，得到波数域回波信号；根据频谱回波信号，建立与重建目标半径相关的聚焦函数表达式；对所述聚焦函数表达式求逆后执行一维傅里叶变换，得到波数域聚焦函数；
根据得到的波数域回波信号和波数域聚焦函数，将二者相乘后，再进行降维重排、累积和逆傅里叶变换处理，得到目标反射率图像。
2.根据权利要求1所述的基于柱面扫描体制下的MIMO阵列三维成像方法，其特征在于，所述MIMO阵列是由多个发射天线阵元和多个接收天线阵元沿某一方向直线排列而成，并以半径R0作圆周运动，以实现柱面扫描，对柱面内目标区域内的目标实现观测。
3.根据权利要求1所述的基于柱面扫描体制下的MIMO阵列三维成像方法，其特征在于，所述根据柱面扫描几何下，MIMO阵列中的发射天线阵元与目标之间的距离，以及MIMO阵列中的接收天线阵元与目标之间的距离，得到距离衰减因子，构建考虑距离衰减的回波模型；其具体过程包括：
根据柱面扫描几何下，发射天线阵元之间的距离RT和接收天线阵元与目标之间的距离RR，得到距离衰减因子，记为1/(4πRTRR)；
其中，
Figure FDA0003523214110000011
其中，(r,θ,z)为预先建立的柱坐标系O-rθz区域中的任一点；
Figure FDA0003523214110000012
Figure FDA0003523214110000013
分别为MIMO阵列中发射天线阵元和接收天线阵元的柱坐标系下的坐标，其中，R0为MIMO阵列柱面扫描的半径；
Figure FDA0003523214110000014
为柱面扫描过程中MIMO阵列同X轴正向的夹角；zT、zR分别代表发射天线阵元在MIMO阵列方向的坐标和接收天线阵元在MIMO阵列方向的坐标；
根据上述得到的距离衰减因子1/(4πRTRR)，构建考虑距离衰减的回波模型：
Figure FDA0003523214110000021
其中，
Figure FDA0003523214110000022
为快时间波数域回波信号；其中，kR为快时间波数；σ(r,θ,z)表示待重建的目标反射率函数；j表示虚部单元；(r,θ,z)为柱坐标系O-rθz区域中的任一点。
4.根据权利要求1所述的基于柱面扫描体制下的MIMO阵列三维成像方法，其特征在于，所述对回波模型相对于发射天线阵元位置、接收天线阵元位置及方位角度变量分别作一维傅里叶变换，得到波数域回波信号；其具体过程包括：
根据回波模型，获取快时间波数域回波信号
Figure FDA0003523214110000023
对快时间波数域回波信号
Figure FDA0003523214110000024
两侧相对于发射天线阵元位置变量zT、接收天线阵元位置变量zR分别执行一维傅里叶变换，得到变换后的频谱回波信号
Figure FDA0003523214110000025
Figure FDA0003523214110000026
其中，kzT、kzR分别为代表变量zT、zR傅里叶变化后对应的波数域变量；σ(r,θ,z)表示待重建的目标反射率函数；
Figure FDA0003523214110000027
为相对于发射天线阵元位置变量zT执行一维傅里叶变换操作；
Figure FDA0003523214110000028
为相对于接收天线阵元位置变量zR执行一维傅里叶变换操作；
利用驻定相位原理求解上式
Figure FDA0003523214110000029
中
Figure FDA00035232141100000210
Figure FDA00035232141100000211
Figure FDA0003523214110000031
Figure FDA0003523214110000032
根据
Figure FDA0003523214110000033
的结果得到另一表示形式下的频谱回波信号
Figure FDA0003523214110000034
Figure FDA0003523214110000035
其中，R0为MIMO阵列柱面扫描的半径；σ(r,θ,z)为待重建的目标反射率函数；
对
Figure FDA0003523214110000036
中的角度变量
Figure FDA0003523214110000037
执行一维傅里叶变换的操作，得到波数域回波信号SSS(kR,kθ,kzT,kzR)；
其中，kθ为角度变量
Figure FDA0003523214110000038
对应的波数域变量。
5.根据权利要求4所述的基于柱面扫描体制下的MIMO阵列三维成像方法，其特征在于，根据所述频谱回波信号，建立与重建目标半径相关的聚焦函数表达式；对所述聚焦函数表达式求逆后执行一维傅里叶变换，得到波数域聚焦函数；包括：
根据所述频谱回波信号
Figure FDA0003523214110000039
建立半径r处聚焦函数表达式为：
Figure FDA0003523214110000041
对所述半径r处聚焦函数表达式求逆后相对于角度变量
Figure FDA0003523214110000042
执行一维傅里叶变换，得到半径r处的波数域聚焦函数
Figure FDA0003523214110000043
6.根据权利要求1所述的基于柱面扫描体制下的MIMO阵列三维成像方法，其特征在于，所述根据得到的波数域回波信号和波数域聚焦函数，将二者相乘后，再进行降维重排、累积和逆傅里叶变换处理得到目标反射率图像；其具体过程包括：
将所述波数域回波信号与某半径处波数域聚焦函数表达式相乘得到聚焦后的波数域回波信号；对所述聚焦后的波数域回波信号执行降维重排操作，得到降维处理后的信号；
对所述降维处理后的信号在快时间波数上进行累积操作得到半径r处的柱面目标的二维频谱；对所述半径r处的柱面目标的二维频谱执行角度维和阵列维的逆傅里叶变换，得到所述半径r处的柱面目标的重建图像；
对每个半径处的柱面目标重复上述过程得到每个半径位置处的柱面目标图像，得到三维目标区域的反射率图像。
7.根据权利要求6所述的基于柱面扫描体制下的MIMO阵列三维成像方法，其特征在于，所述将所述波数域回波信号与某半径处波数域聚焦函数表达式相乘得到聚焦后的波数域回波信号；对所述聚焦后的波数域回波信号执行降维重排操作，得到降维处理后的信号；其具体过程包括：
将所述波数域回波信号SSS(kR,kθ,kzT,kzR)与半径r处波数域聚焦函数
Figure FDA0003523214110000044
相乘，得到聚焦后的波数域回波信号表达式：
Figure FDA0003523214110000045
对所述聚焦后的波数域回波信号SSSr(kR,kθ,kzT,kzR)执行降维重排操作：按照设定的色散关系kz＝kzT+kzR对该信号进行降维重排操作，得到降维处理后的信号
Figure FDA0003523214110000046
其中，kz为定义的阵列向波数；
Figure FDA0003523214110000051
8.根据权利要求6所述的基于柱面扫描体制下的MIMO阵列三维成像方法，其特征在于，所述对所述降维处理后的信号在快时间波数上进行累积操作得到半径r处的柱面目标的二维频谱；对所述半径r处的柱面目标的二维频谱执行角度维和阵列维的逆傅里叶变换，得到所述半径r处的柱面目标的重建图像；其具体过程包括：
对所述降维处理后的信号
Figure FDA0003523214110000052
在快时间波数kR上进行累积操作得到半径处的柱面目标的二维频谱表达式：
Figure FDA0003523214110000053
对所述半径r处的柱面目标的二维频谱Σr(kθ,kz)执行角度维和阵列维的逆傅里叶变换，得到半径r处的柱面目标的重建图像
Figure FDA0003523214110000054
Figure FDA0003523214110000055
9.一种基于柱面扫描体制下的MIMO阵列三维成像装置，其特征在于，该装置包括：回波模型构建模块、信号获取模块和图像三维重建模块；
所述回波模型构建模块，用于根据柱面扫描几何下，MIMO阵列中的发射天线阵元与目标之间的距离，以及MIMO阵列中的接收天线阵元与目标之间的距离，得到距离衰减因子，构建考虑距离衰减的回波模型；
所述信号获取模块，所述对回波模型相对于发射天线阵元位置、接收天线阵元位置及方位角度变量分别作一维傅里叶变换，得到波数域回波信号；根据频谱回波信号，建立与重建目标半径相关的聚焦函数表达式；对所述聚焦函数表达式求逆后执行一维傅里叶变换，得到波数域聚焦函数；
所述图像三维重建模块，用于根据得到的波数域回波信号和波数域聚焦函数，将二者相乘后，再进行降维重排、累积和逆傅里叶变换处理，得到目标反射率图像。
CN202210187239.9A 2022-02-28 2022-02-28 基于柱面扫描体制下的mimo阵列三维成像方法及成像装置 Pending CN114624707A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210187239.9A CN114624707A (zh) 2022-02-28 2022-02-28 基于柱面扫描体制下的mimo阵列三维成像方法及成像装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210187239.9A CN114624707A (zh) 2022-02-28 2022-02-28 基于柱面扫描体制下的mimo阵列三维成像方法及成像装置

Publications (1)

Publication Number Publication Date
CN114624707A true CN114624707A (zh) 2022-06-14

Family

ID=81900603

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210187239.9A Pending CN114624707A (zh) 2022-02-28 2022-02-28 基于柱面扫描体制下的mimo阵列三维成像方法及成像装置

Country Status (1)

Country Link
CN (1) CN114624707A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2024154212A1 (ja) * 2023-01-17 2024-07-25 株式会社三井E&S データ処理方法、計測システム、及び、プログラム
WO2024154211A1 (ja) * 2023-01-17 2024-07-25 株式会社三井E&S データ処理方法、計測システム、及び、プログラム

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2024154212A1 (ja) * 2023-01-17 2024-07-25 株式会社三井E&S データ処理方法、計測システム、及び、プログラム
WO2024154211A1 (ja) * 2023-01-17 2024-07-25 株式会社三井E&S データ処理方法、計測システム、及び、プログラム

Similar Documents

Publication Publication Date Title
CN110275166B (zh) 基于admm的快速稀疏孔径isar自聚焦与成像方法
CN112180329B (zh) 一种基于阵元随机均匀分布球阵反卷积波束形成的汽车噪声源声成像方法
CN113820711B (zh) 阵列快速成像方法及其应用
CN114624707A (zh) 基于柱面扫描体制下的mimo阵列三维成像方法及成像装置
WO2018196254A1 (zh) 一种微波成像系统的幅相校正方法及系统
CN111896951B (zh) 一种毫米波柱面全息成像系统的三维成像与重构方法
CN106772380A (zh) 一种圆周合成孔径雷达成像方法
CN106569191A (zh) 一种利用高分辨率成像获取目标rcs的方法
CN107390215A (zh) 一种高速超分辨率mimo阵列成像方法
CN110609282B (zh) 基于后向投影的太赫兹孔径编码三维成像方法及装置
CN112462363B (zh) 非均匀稀疏极化阵列相干目标参数估计方法
CN113933834A (zh) 圆柱扫描微波成像方法
CN109884627A (zh) 任意线阵构型的近程毫米波快速三维成像方法
CN108107428B (zh) 用于mimo阵列的相移偏移成像方法及装置
CN105699970A (zh) 一种用于合成孔径雷达大场景欺骗干扰的逆距离徙动方法
CN110764089A (zh) 一种超分辨率毫米波mimo阵列实时成像方法
Wang et al. Near-field synthetic spectrum imaging algorithm based on single-frequency sparse MIMO array
CN111812644B (zh) 基于稀疏估计的mimo雷达成像方法
CN113835222A (zh) 曲面阵列快速成像方法
CN103630907B (zh) 近距离主动式毫米波圆柱扫描成像系统的免插值重构方法
CN111090102B (zh) 一种超分辨反射式太赫兹三维目标重建成像方法
CN102183755A (zh) 一种基于柯西-高斯模型的新型高分辨方位估计方法
Cao et al. Microwave coincidence imaging based on attributed scattering model
Zhang et al. Mid-range 2-D fast imaging in pseudo-polar coordinate format for wideband array
CN114994668A (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