CN108761457B - 基于mimo阵列合成孔径的高精度三维快速成像方法及装置 - Google Patents
基于mimo阵列合成孔径的高精度三维快速成像方法及装置 Download PDFInfo
- Publication number
- CN108761457B CN108761457B CN201810577125.9A CN201810577125A CN108761457B CN 108761457 B CN108761457 B CN 108761457B CN 201810577125 A CN201810577125 A CN 201810577125A CN 108761457 B CN108761457 B CN 108761457B
- Authority
- CN
- China
- Prior art keywords
- spectrum
- synthetic aperture
- mimo array
- total
- original
- 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
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
- 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
-
- 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/904—SAR modes
-
- 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阵列合成孔径的高精度三维快速成像方法,其包括以下步骤:S1、向目标发送宽带信号,并接收所述宽带信号经过目标散射后得到的原始回波信号;S2、将所述原始回波信号变换到MIMO阵列方向以及合成孔径方向所对应的空间频率域,确定原始空间谱;S3、将所述原始空间谱结合相位偏移因子,确定包括距离向的总空间谱;以及S4、根据所述总空间谱确定目标的成像函数。本发明还涉及一种基于MIMO阵列合成孔径的高精度三维快速成像装置。
Description
技术领域
本发明涉及信号处理技术领域,尤其涉及一种基于MIMO阵列合成孔径的高精度三维快速成像方法及装置。
背景技术
雷达三维实时成像装置最常用的是通过二维阵列实现方位向快速聚焦,以及通过宽带信号实现距离向的聚焦。然而在频段比较高的应用场景下,例如毫米波太赫兹频段,二维阵列由于阵元个数较多,在现有条件下,给系统增加了诸多成本,尤其是目前毫米波太赫兹器件本身成本就相对较高。虽然只采用单个阵元进行二维扫描,结合宽带信号也可获取三维聚焦图像,但是数据的获取时间太长,做不到实时成像。
结合MIMO阵列与合成孔径技术的多输入多输出合成孔径雷达(Multiple InputMultiple Output Synthetic Aperture Radar,MIMO-SAR),可减少阵元个数,降低系统成本。然而目前基于此装置的成像算法,不能同时满足高精度成像与快速成像的要求。
发明内容
本发明的目的在于提供一种基于MIMO阵列合成孔径的高精度三维快速成像方法及装置,以解决上述的至少一项技术问题。
本发明的一方面,提供了一种基于MIMO阵列合成孔径的高精度三维快速成像方法,包括如下步骤:
S1、向目标发送宽带信号,并接收所述宽带信号经过目标散射后得到的原始回波信号;
S2、将所述原始回波信号变换到MIMO阵列方向以及合成孔径方向所对应的空间频率域,确定原始空间谱;
S3、将所述原始空间谱结合相位偏移因子,确定总空间谱;以及
S4、根据所述总空间谱确定目标的成像函数。
在一些实施例中,在步骤S1中,构建三维坐标系,以x方向表示MIMO阵列方向,y方向表示合成孔径方向,z方向表示距离向,所述原始回波信号为:
S(xt,xr,y,0,k)
其中,MIMO阵列的发射天线位于(xt,y,z)处,接收天线位于(xr,y,z)处,k为所述宽带信号不同发射频率所对应的波数,MIMO阵列合成孔径所在的平面的距离向位于z=0处。
在一些实施例中,在步骤S2中,通过对所述原始回波信号进行三维傅里叶变换以将所述原始回波信号变换到MIMO阵列方向以及合成孔径方向所对应的空间频率域,得到的所述原始空间谱的公式为:
其中,kxt,kxr,ky分别表示xt,xr,y对应的空间频率域坐标。
在一些实施例中,步骤S3包括:
S31、根据空间频率域的坐标关系将所述原始空间谱结合相位偏移因子,得到初始总空间谱;
S32、对初始总空间谱进行数据重排,并将重复的空间频率域坐标对应的空间谱能量进行叠加,得到总空间谱。
在一些实施例中,在步骤S31中:
所述相位偏移因子为exp(jkzz),
所述空间频率域的坐标关系包括kx=kxt+kxr;
所述初始总空间谱的公式为:s(kxt,kxr,ky,z,k)=s(kxt,kxr,ky,0,k)·exp(jkzz),
其中,z=zN,zN为距离向z的N个实数取值。
在一些实施例中,步骤S32中:
对所述初始总空间谱进行数据重排后得到的总空间谱的公式为:
S(kx,ky,z,k)=S(kxt,kxr,ky,z,k)rearrange。
在一些实施例中,步骤S4包括:
S41、对所述总空间谱进行波数域积分,确定平均波数域总空间谱;以及
S42、对所述平均波数域总空间谱进行逆傅里叶变换,确定目标的成像函数。
在一些实施例中:
所述目标的成像函数为:
本发明的另一方面,还提供了一种基于MIMO阵列合成孔径的高精度三维快速成像装置,包括:
存储器,用于存储指令;以及
处理器,用于根据所述指令,执行前述的基于MIMO阵列合成孔径的高精度三维快速成像方法。
本发明的基于MIMO阵列合成孔径的高精度三维快速成像方法及装置,相较于现有技术,至少具有以下有益效果其中之一或其中的一部分:
1、相较于经典的可用于任意阵列形式的反向传播(BP)算法,成像速度得到了提升。
2、相较于成像速度较快的经过一定条件近似的快速傅氏变换(FFT)算法,成像质量得到改善。
附图说明
图1为本发明实施例的基于MIMO阵列合成孔径的高精度三维快速成像方法的步骤的流程图;
图2为本发明实施例的基于MIMO阵列合成孔径的高精度三维快速成像装置的结构示意图;
图3为MIMO线性阵列实例示意图;
图4为本发明实施例的步骤S3的具体步骤的流程图;
图5为本发明实施例的步骤S4的具体步骤的流程图;
图6为本发明的实施例的目标三维成像结果示意图;
图7为本发明的实施例的距离为R处的目标二维成像结果示意图;
图8为本发明实施例的用基于MIMO-SAR的高精度三维快速成像装置的结构示意图。
具体实施方式
为使本公开的目的、技术方案和优点更加清楚明白,以下结合具体实施例,并参照附图,对本公开进一步详细说明。
本发明的一方面,提供了一种基于MIMO阵列合成孔径的高精度三维快速成像方法,图1为本发明实施例的基于MIMO阵列合成孔径的高精度三维快速成像方法的步骤示意图,如图1所示,该方法包括以下步骤:
S1、向目标发送宽带信号,并接收宽带信号经过目标散射后得到的原始回波信号;
S2、将原始回波信号变换到MIMO阵列方向以及合成孔径方向所对应的空间频率域,确定原始空间谱;
S3、将原始空间谱结合相位偏移因子,确定总空间谱;
S4、根据总空间谱确定目标的成像函数。
其中,在步骤S1中,宽带信号由MIMO阵列的发射阵元发射,MIMO阵列的接收阵元接收原始回波信号。
在本发明实施例中,MIMO阵列选择发射在接收阵列两端形式的线性阵列,可以理解的是,在其他实施例中也可以选择其他形式的MIMO阵列。图2为本发明实施例的MIMO线性阵列合成孔径仿真场景示意图,图3为本发明实施例的MIMO线性阵列实例示意图,如图2和图3所示,在该MIMO线性阵列中,接收阵元以间距为dR排列,在第一个和最后一个接收阵元的两端,各有发射阵元以间距为dT排列。接收阵列长度为LR,接收阵列和发射阵列总长度为LT。有两组目标点组成的线性目标以交叉形式排列,目标到MIMO阵列合成孔径平面的距离为R。
在一些实施例中,参照图2,构建三维坐标系,以x方向表示MIMO阵列方向,y方向表示合成孔径方向,z方向表示距离向,原始回波信号为:
S(xt,xr,y,O,k)
其中,MIMO阵列的发射天线位于(xt,y,z)处,接收天线位于(xr,y,z)处,k为宽带信号不同发射频率所对应的波数,MIMO阵列合成孔径所在的平面的距离向位于z=0处。
根据一些实施例,在步骤S2中,通过对原始回波信号进行三维傅里叶变换以将原始回波信号变换到MIMO阵列方向以及合成孔径方向所对应的空间频率域,以得到原始空间谱,原始空间谱的公式为:
其中,kxt,kxr,ky分别表示xt,xr,y对应的空间频率域坐标。
图4为本发明实施例的步骤S3的具体步骤示意图,如图4所示,根据一些实施例,步骤S3可以包括以下子步骤:
S31、根据空间频率域的坐标关系将原始空间谱结合相位偏移因子,得到初始总空间谱;
例如,相位偏移因子为exp(jkzz),
空间频率域的坐标关系包括kx=kxt+kxr,
初始总空间谱的公式为:S(kxt,kxr,ky,z,k)=S(kxt,kxr,ky,0,k)·exp(jkzz),
其中,z=zN,zN为距离向z的N个实数取值。
S32、将初始总空间谱进行数据重排,并将重复的空间频率域坐标对应的空间谱能量进行叠加,得到总空间谱。
例如,重排后得到的总空间谱的公式为:
S(kx,ky,z,k)=S(kxt,kxr,ky,z,k)rearrange
图5为本发明实施例的步骤S4的具体步骤示意图,如图5所示,该步骤S4可以包括以下子步骤:
S41、对总空间谱进行波数域积分,确定平均波数域总空间谱。
S42、对平均波数域总空间谱进行逆傅里叶变换,最终确定目标的成像函数。
例如,目标的成像函数公式为:
下面结合附图介绍本发明的一个具体实施例。参照图2,目标区域的图像如图2所示,按照以下步骤以实现基于MIMO阵列合成孔径的高精度三维快速成像。
S1、向目标发送宽带信号,并接收宽带信号经过目标散射后得到的原始回波信号。
构建三维坐标系,以x方向表示MIMO阵列方向,y方向表示合成孔径方向,z方向表示距离向,如图2所示,MIMO阵列的发射天线位于(xt,y,z)处,接收天线位于(xr,y,z)处,MIMO-SAR所在的平面的距离向位于z=0处,以k来表示所述宽带信号不同发射频率所对应的波数,
则所述原始回波信号可表示为:S(xt,xr,y,0,k)。
S2、将原始回波信号变换到MIMO阵列方向以及合成孔径方向所对应的空间频率域,以得到原始空间谱。
可以理解的是,在本实施例中,将所述原始回波信号变换到空间频率域是可以通过对所述简化公式进行三维傅里叶变换实现的,则原始空间谱的公式为:
其中,kxt,kxr,ky分别表示xt,xr,y对应的空间频率域坐标。
S3、将原始空间谱结合相位偏移因子,确定总空间谱。
本实施例中,步骤S3包括以下子步骤:
S31、根据空间频率域的坐标关系将所述原始空间谱结合相位偏移因子,得到初始总空间谱,相位偏移因子为exp(jkzz),
空间频率域的坐标关系包括kx=kxt+kxr;
初始总空间谱的公式为:S(kxt,kxr,ky,z,k)=S(kxt,kxr,ky,0,k)·exp(jkzz),
其中,z=zN,zN为距离向z的N个实数取值;
S32、对初始总空间谱进行数据重排,并将重复的空间频率域坐标对应的空间谱能量进行叠加,得到总空间谱,排后总空间谱的公式为:
S(kx,ky,z,k)=s(kxt,kxr,ky,z,k)rearrange。
S4、根据总空间谱确定目标的成像函数。
本实施例中,步骤S4包括以下子步骤:
S42、对平均波数域总空间谱进行逆傅里叶变换,确定目标的成像函数,目标的成像函数为:
根据上述方法,参照图6和图7。图6为本实施例中得到的目标三维成像结果示意图,图7为本实施例中距离为R处的目标二维成像结果示意图。
可见,本发明的基于MIMO阵列合成孔径的高精度三维快速成像方法缩短了成像时间,提高了成像质量。
本发明的另一方面,还提供了一种基于MIMO阵列合成孔径的高精度三维快速成像装置,图8为本发明实施例的用基于MIMO阵列合成孔径的高精度三维快速成像装置的结构示意图,如图8所示,该装置包括:
存储器81,用于存储指令;以及
处理器82,用于根据存储器81中的指令,执行前述的基于MIMO阵列合成孔径的高精度三维快速成像方法。
综上,本发明的基于MIMO阵列合成孔径的高精度三维快速成像方法及装置,相较于经典的可用于任意阵列形式的BP算法,成像速度得到了提升,相较于成像速度较快的经过一定条件近似的FFT快速算法,成像质量得到了改善。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (7)
1.一种基于MIMO阵列合成孔径的高精度三维快速成像方法,包括:
S1、向目标发送宽带信号,并接收所述宽带信号经过目标散射后得到的原始回波信号,其中,构建三维坐标系,以x方向表示MIMO阵列方向,y方向表示合成孔径方向,z方向表示距离向,所述原始回波信号为:
S(xt,xr,y,0,k)
其中,MIMO阵列的发射天线位于(xt,y,z)处,接收天线位于(xr,y,z)处,k为所述宽带信号不同发射频率所对应的波数,MIMO阵列合成孔径所在的平面的距离向位于z=0处;
S2、通过对所述原始回波信号进行三维傅里叶变换以将所述原始回波信号变换到MIMO阵列方向以及合成孔径方向所对应的空间频率域,得到的原始空间谱的公式为:
其中,kxt,kxr,ky分别表示xt,xr,y对应的空间频率域坐标;
S3、将所述原始空间谱结合相位偏移因子,确定总空间谱;以及
S4、根据所述总空间谱确定目标的成像函数。
2.根据权利要求1所述的方法,其中,步骤S3包括:
S31、根据空间频率域的坐标关系将所述原始空间谱结合所述相位偏移因子,得到初始总空间谱;
S32、对所述初始总空间谱进行数据重排,并将重复的空间频率域坐标对应的空间谱能量进行叠加,得到总空间谱。
4.根据权利要求3所述的方法,其中,在步骤S32中:
对所述初始总空间谱进行数据重排后得到的总空间谱的公式为:
S(kx,ky,z,k)=S(kxt,kxr,ky,z,k)rearrange。
5.根据权利要求4所述的方法,其中,步骤S4包括:
S41、对所述总空间谱进行波数域积分,确定平均波数域总空间谱;以及
S42、对所述平均波数域总空间谱进行逆傅里叶变换,确定目标的成像函数。
7.一种基于MIMO阵列合成孔径的高精度三维快速成像装置,包括:
存储器,用于存储指令;以及
处理器,用于根据所述指令,执行如权利要求1至6中的任一项所述的基于MIMO阵列合成孔径的高精度三维快速成像方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810577125.9A CN108761457B (zh) | 2018-06-05 | 2018-06-05 | 基于mimo阵列合成孔径的高精度三维快速成像方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810577125.9A CN108761457B (zh) | 2018-06-05 | 2018-06-05 | 基于mimo阵列合成孔径的高精度三维快速成像方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108761457A CN108761457A (zh) | 2018-11-06 |
CN108761457B true CN108761457B (zh) | 2021-04-16 |
Family
ID=64000340
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810577125.9A Active CN108761457B (zh) | 2018-06-05 | 2018-06-05 | 基于mimo阵列合成孔径的高精度三维快速成像方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108761457B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110850401B (zh) * | 2019-08-27 | 2022-06-28 | 天津大学 | 一种基于运动模型与合成孔径的rfid标签定位方法 |
CN111650585B (zh) * | 2020-08-04 | 2020-11-03 | 中国人民解放军国防科技大学 | 近场毫米波稀疏mimo扫描阵列全聚焦成像方法和装置 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101190731B1 (ko) * | 2010-06-28 | 2012-10-16 | 한국과학기술원 | 광역 고해상도 영상을 위한 다중 입력 다중 출력 영상 레이더 이용방법 및 이를 이용한 시스템 |
CN104237888B (zh) * | 2014-10-20 | 2016-08-24 | 内蒙古工业大学 | 弧形阵列mimo-sar的成像方法 |
CN104267401B (zh) * | 2014-10-21 | 2016-03-02 | 内蒙古工业大学 | 线性阵列天线同时mimo-sar成像系统和方法 |
CN108107428B (zh) * | 2017-12-07 | 2021-09-10 | 中国科学院电子学研究所 | 用于mimo阵列的相移偏移成像方法及装置 |
-
2018
- 2018-06-05 CN CN201810577125.9A patent/CN108761457B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN108761457A (zh) | 2018-11-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108828593B (zh) | 一种随机辐射雷达关联成像方法 | |
CN107976673B (zh) | 提高大场景目标成像质量的mimo雷达成像方法 | |
CN108107428B (zh) | 用于mimo阵列的相移偏移成像方法及装置 | |
KR101603630B1 (ko) | 송신신호 파형의 직교성을 이용한 재머 신호 억제 및 원신호 도래각 추정 방법 및 장치 | |
CN108919260B (zh) | 用于mimo阵列的相移偏移成像方法及装置 | |
CN110596646A (zh) | 一种基于mimo体制的提高雷达角度分辨率的布局及方法 | |
CN108957449A (zh) | 基于mimo阵列合成孔径的三维快速成像方法及装置 | |
CN110346793A (zh) | 一种分布式阵列雷达迭代自适应高分辨成像方法 | |
CN111537997B (zh) | 一种基于mimo和压缩感知技术的三维雷达成像方法 | |
CN109884627B (zh) | 任意线阵构型的近程毫米波快速三维成像方法 | |
JP2016130726A (ja) | 3次元(3d)シーン画像を生成するシステム及び方法 | |
JP5489813B2 (ja) | レーダ画像処理装置 | |
CN108761457B (zh) | 基于mimo阵列合成孔径的高精度三维快速成像方法及装置 | |
Skouroliakou et al. | Fourier-based radar processing for multistatic millimetre-wave imaging with sparse apertures | |
CN110988869B (zh) | 基于mimo阵列的成像方法及装置 | |
JP7228567B2 (ja) | センサアレイ撮像装置 | |
CN108181624B (zh) | 一种差分计算成像装置及方法 | |
CN110554383B (zh) | 用于微波频段的mimo环形阵列方位向成像方法及装置 | |
CN108693529B (zh) | 基于mimo-sar的高精度三维快速成像方法及装置 | |
CN108761433B (zh) | 一种使用mimo声纳差合阵处理的高分辨成像方法 | |
CN111090102B (zh) | 一种超分辨反射式太赫兹三维目标重建成像方法 | |
CN110658502A (zh) | 一种幅相误差校正方法 | |
CN113671477B (zh) | 一种基于图信号处理的雷达目标距离估计方法 | |
CN110554387B (zh) | 综合孔径干涉近场主动源成像方法和装置 | |
CN113176566B (zh) | 一种用于mimo-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 |