CN103698761B - 一种基于椭球切点包络面的雷达三维成像方法及系统 - Google Patents
一种基于椭球切点包络面的雷达三维成像方法及系统 Download PDFInfo
- Publication number
- CN103698761B CN103698761B CN201310651858.XA CN201310651858A CN103698761B CN 103698761 B CN103698761 B CN 103698761B CN 201310651858 A CN201310651858 A CN 201310651858A CN 103698761 B CN103698761 B CN 103698761B
- Authority
- CN
- China
- Prior art keywords
- ellipsoid
- antenna
- enveloping surface
- axis
- dual
- 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.)
- Expired - Fee Related
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 27
- 238000000034 method Methods 0.000 claims abstract description 17
- 238000002592 echocardiography Methods 0.000 claims abstract description 11
- 230000036039 immunity Effects 0.000 abstract description 2
- 230000003321 amplification Effects 0.000 description 1
- 230000037237 body shape Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000003199 nucleic acid amplification method 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
-
- 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)
- Aerials With Secondary Devices (AREA)
Abstract
本发明公开了一种基于椭球切点包络面的雷达三维成像方法及系统。所述方法包括:将收发天线分置;于Δt时间间隔在m×n天线坐标矩阵扫描域内间隔移动收发天线,得到m×n组电磁回波数据;分别以收发天线坐标为椭圆焦点、目标体反射回波路径长度为椭圆长轴,构造出m×n个椭圆并分别求出各个椭圆半短轴;以椭圆半长轴为X轴和Y轴的椭球赤道半径,以椭圆半短轴为Z轴的椭球极半径,构造出m×n个椭球组成椭球集;构造与椭球集的各椭球相切的包络面,选取包络最大Z轴方向切点的包络面,在所述包络最大Z轴方向切点的包络面中,以所有的任意4个相邻切点的几何中心组成的包络面作为目标体边界实现成像过程。本发明具有强抗噪性,能提高成像速度和精确度,满足实时性应用要求。
Description
技术领域
本发明涉及雷达成像技术领域,具体涉及一种基于椭球切点包络面的雷达三维成像方法及系统。
背景技术
目前,雷达二维成像技术理论框架已趋于成熟,然而二维成像仅能还原目标体横截面形状,无法对其整体进行成像与重建,导致雷达系统对物体辨识度较低。因此,雷达三维成像技术成为了近几年国内外学者研究的热点与难点。与二维成像不同,三维成像需要考虑目标体的距离-方位-俯仰三个参数。在数据处理方面,增加了表征高度的坐标轴来进行目标体重建,必须将其扩展至面阵天线或合成孔径雷达等效阵面,以采集多个水平面上的回波信息。
已有的雷达三维成像技术主要是通过回波信号的能量叠加原理,对整个成像域进行目标位置还原。但是,此类方法由于计算量过大,不利于实时成像的应用。同时,对于目标体边界形状的还原精确度不高,因此成像辨识度也受到了一定的影响。
发明内容
针对现有技术的不足,本发明提供一种基于椭球切点包络面的雷达三维成像方法及系统。本发明根据回波路径通过构造目标体边界点与雷达天线之间的几何关系,重建目标体边界,其中的数学算法过程简单、计算复杂度低、无任何放大噪声的数学式,具有较强的抗噪性,从而实现实时性应用。另外,结合收发分置天线尤其是一发多收天线进行发射与回波接收过程,能提高成像的速度以及成像的精确度,更利于应用在实时成像中。
本发明的技术方案如下所述。
一种基于椭球切点包络面的雷达三维成像方法,所述方法包括:
将收发天线分置;
于Δt时间间隔在m×n天线坐标矩阵扫描域内间隔移动收发天线,得到m×n组电磁回波数据;
分别以收发天线坐标为椭圆焦点、从发射天线至目标体后反射回接收天线的电磁波路径长度的1/2为椭圆半长轴,在目标体与收发天线构成的二维平面上的构造出m×n个椭圆并分别求出各个椭圆半短轴;
分别将收发天线坐标视为X轴方向的椭球焦点,以所述椭圆半长轴为X轴方向的椭球赤道半径,以所述椭圆半短轴为Y轴方向的椭球赤道半径,同时以所述椭圆半短轴为Z轴方向的椭球极半径,构造出m×n个椭球组成椭球集;
构造与椭球集的各椭球相切的包络面,选取包络最大Z轴方向切点的包络面,在所述包络最大Z轴方向切点的包络面中,以所有的任意4个相邻切点的几何中心组成的包络面作为目标体边界,实现成像过程;
其中:m=1,2,3,…;n=1,2,3,…。
所述收发天线可以是一发多收天线,也可以是阵列面天线;阵列面天线按Δt时间间隔顺序扫描获取目标体电磁回波数据或同时获取目标体电磁回波数据。
所述电磁回波数据中的载波频率是500MHz~7.5GHz。
一种基于椭球切点包络面的雷达三维成像系统,所述系统包括:
分置的收发天线;
于Δt时间间隔在m×n天线坐标矩阵扫描域内间隔移动收发天线,得到m×n组电磁回波数据的装置;
分别以收发天线坐标为椭圆焦点、从发射天线至目标体后反射回接收天线的电磁波路径长度的1/2为椭圆半长轴,在目标体与收发天线构成的二维平面上的构造出m×n个椭圆并分别求出各个椭圆半短轴的装置;
分别将收发天线坐标视为X轴方向的椭球焦点,以所述椭圆半长轴为X轴方向的椭球赤道半径,以所述椭圆半短轴为Y轴方向的椭球赤道半径,同时以所述椭圆半短轴为Z轴方向的椭球极半径,构造出m×n个椭球组成椭球集的装置;
构造与椭球集的各椭球相切的包络面,选取包络最大Z轴方向切点的包络面,在所述包络最大Z轴方向切点的包络面中,以所有的任意4个相邻切点的几何中心组成的包络面作为目标体边界,实现成像过程的装置。
附图说明
图1描述了收发天线与目标体的电磁波路径的几何关系。
图2A、图2B描述了本发明所构造椭球与目标体的几何关系。
图3、图4是本发明的FDTD仿真三维成像效果图。
具体实施方式
下面结合附图和实施例对本发明做进一步的详细描述。
发射与接收天线坐标分别为FT(XT,YT,0)和FR(XR,YR,0),电磁波经过发射天线传播至目标体后反射回到接收天线处。设电磁波在目标体边界点C(x,y,z)反射,那么电磁波路径是FT(XT,YT,0)→C(x,y,z)→FR(XR,YR,0),即从发射天线至目标体后反射回接收天线的电磁波路径长度可表示成R=|FTC|+|CFR|,如图1所示。将发射天线与接收天线坐标FT(XT,YT,0)和FR(XR,YR,0)视为椭圆的两个焦点坐标,其焦距为|XT-XR|。以FTCFR为参考平面,构造出焦点坐标为FT(XT,YT,0)和FR(XR,YR,0)、长轴为R=|FTC|+|CFR|的椭圆,根据椭圆方程式求出椭圆的短轴为
以上述椭圆的半长轴R/2为X轴方向的椭球赤道半径,以上述椭圆半短轴为Y轴方向的椭球赤道半径,同时也以上述椭圆半短轴为Z轴方向的椭球极半径,构造出的椭球表达式为:
通过上述的构造方法,在每个Δt时间间隔都可构造出一个椭球。设在Δtij所构造的椭球为那么经过总时间T={Δtij,i=1,2,…,m;j=1,2,…,n}后,得到
的椭球集为:
其中:i=1,2,…,m;j=1,2,…,n,m和n分别为天线移动位置的行列数。
构造与椭球集的各椭球相切的包络面,选取包络最大Z轴方向切点的包络面,在所述包络最大Z轴方向切点的包络面中,以所有的任意4个相邻切点的几何中心组成的包络面作为目标体边界,实现成像过程。
Claims (6)
1.一种基于椭球切点包络面的雷达三维成像方法,所述方法包括:
将收发天线分置;
于△t时间间隔在m×n天线坐标矩阵扫描域内间隔移动收发天线,得到
m×n组电磁回波数据;
分别以收发天线坐标为椭圆焦点、从发射天线至目标体后反射回接收天线的
电磁波路径长度的1/2为椭圆半长轴,在目标体与收发天线构成的二维平面上构
造出m×n个椭圆并分别求出各个椭圆半短轴;
分别将收发天线坐标视为X轴方向的椭球焦点,以所述椭圆半长轴为X轴方
向的椭球赤道半径,以所述椭圆半短轴为Y轴方向的椭球赤道半径,同时以所述
椭圆半短轴为z轴方向的椭球极半径,构造出m×n个椭球组成椭球集;
构造与椭球集的各椭球相切的包络面,选取包络最大z轴方向切点的包络面,
在所述包络最大z轴方向切点的包络面中,以所有的任意4个相邻切点的几何中
心组成的包络面作为目标体边界,实现成像过程;
其中:m=1,2,3,....;n=1,2,3,...。
2.根据权利要求1所述的方法,其中所述收发天线是一发多收天线。
3.根据权利要求1所述的方法,其中所述收发天线是阵列面天线。
4.根据权利要求3所述的方法,其中所述阵列面天线按Δt时间间隔顺序扫
描获取目标体电磁回波数据或同时获取目标体电磁回波数据。
5.根据权利要求1,2,3或4所述的方法,其中所述电磁回波数据中的载波
频率是500MHz一7.5GHz。
6.一种基于椭球切点包络面的雷达三维成像系统,所述系统包括:
于Δt时间间隔在m×n天线坐标矩阵扫描域内间隔移动收发天线,得到
m×n组电磁回波数据的装置;
分别以收发天线坐标为椭圆焦点、从发射天线至目标体后反射回接收天线的
电磁波路径长度的1/2为椭圆半长轴,在目标体与收发天线构成的二维平面上构
造出m×n个椭圆并分别求出各个椭圆半短轴的装置;
分别将收发天线坐标视为X轴方向的椭球焦点,以所述椭圆半长轴为X轴
方向的椭球赤道半径,以所述椭圆半短轴为Y轴方向的椭球赤道半径,同时以所
述椭圆半短轴为z轴方向的椭球极半径,构造出m×n个椭球组成椭球集的装置;
构造与椭球集的各椭球相切的包络面,选取包络最大Z轴方向切点的包络
面,在所述包络最大Z轴方向切点的包络面中,以所有的任意4个相邻切点的几
何中心组成的包络面作为目标体边界,实现成像过程的装置;
其中:m=1,2,3,…;n=1,2,3,…。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310651858.XA CN103698761B (zh) | 2013-12-06 | 2013-12-06 | 一种基于椭球切点包络面的雷达三维成像方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310651858.XA CN103698761B (zh) | 2013-12-06 | 2013-12-06 | 一种基于椭球切点包络面的雷达三维成像方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103698761A CN103698761A (zh) | 2014-04-02 |
CN103698761B true CN103698761B (zh) | 2016-01-20 |
Family
ID=50360355
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310651858.XA Expired - Fee Related CN103698761B (zh) | 2013-12-06 | 2013-12-06 | 一种基于椭球切点包络面的雷达三维成像方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103698761B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3422955B1 (en) * | 2016-02-29 | 2023-10-18 | Packsize International, LLC | System and method for assisted 3d scanning |
CN109194338A (zh) * | 2018-11-07 | 2019-01-11 | 桂林电子科技大学 | 一种混合节点多比特处理的极化码译码方法 |
CN116148862B (zh) * | 2023-01-16 | 2024-04-02 | 无锡市雷华科技有限公司 | 一种探鸟雷达飞鸟综合预警评估方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6340946B1 (en) * | 2000-08-03 | 2002-01-22 | Massachusetts Institue Of Technology | Method for determining storm predictability |
CN101738614A (zh) * | 2008-11-17 | 2010-06-16 | 清华大学 | 基于时空图像序列的逆合成孔径雷达的目标转动估计方法 |
CN102253386A (zh) * | 2011-03-29 | 2011-11-23 | 西安电子科技大学 | 基于发射波束扫描的机载下视三维合成孔径雷达成像系统 |
CN102866401A (zh) * | 2012-08-06 | 2013-01-09 | 西北工业大学 | 一种基于mimo技术的三维成像方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7609198B2 (en) * | 2007-05-21 | 2009-10-27 | Spatial Digital Systems, Inc. | Apparatus and method for radar imaging by measuring spatial frequency components |
-
2013
- 2013-12-06 CN CN201310651858.XA patent/CN103698761B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6340946B1 (en) * | 2000-08-03 | 2002-01-22 | Massachusetts Institue Of Technology | Method for determining storm predictability |
CN101738614A (zh) * | 2008-11-17 | 2010-06-16 | 清华大学 | 基于时空图像序列的逆合成孔径雷达的目标转动估计方法 |
CN102253386A (zh) * | 2011-03-29 | 2011-11-23 | 西安电子科技大学 | 基于发射波束扫描的机载下视三维合成孔径雷达成像系统 |
CN102866401A (zh) * | 2012-08-06 | 2013-01-09 | 西北工业大学 | 一种基于mimo技术的三维成像方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103698761A (zh) | 2014-04-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102640014B (zh) | 图像生成装置和图像生成方法 | |
CN105572652B (zh) | 一种利用外推获得具有多次散射目标远场rcs的方法 | |
CN106405548A (zh) | 基于多任务贝叶斯压缩感知的逆合成孔径雷达成像方法 | |
CN112099011B (zh) | 基于focuss算法的全息sar子孔径三维重建方法 | |
CN103698761B (zh) | 一种基于椭球切点包络面的雷达三维成像方法及系统 | |
CN110426707A (zh) | 基于轨道角动量的涡旋sar成像方法及成像系统 | |
CN102426360A (zh) | 空中微动旋转目标的二维isar成像方法 | |
CN102914773B (zh) | 一种多航过圆周sar三维成像方法 | |
CN102393518A (zh) | 一种适用于大斜视角的机载sar成像方法 | |
CN103793059A (zh) | 一种基于时域多普勒手势恢复识别方法 | |
CN103630900B (zh) | 3-d sar波数域快速成像的方法 | |
JP2013174584A (ja) | シーン内の反射物を再構成する方法 | |
CN104101864A (zh) | 基于eof分解的导航x波段雷达海浪参数反演算法 | |
US20210275141A1 (en) | Ultrasound method and apparatus | |
CN103913733B (zh) | 极地冰川厚度探测方法 | |
CN103300826A (zh) | 一种用于早期乳腺肿瘤超宽带微波检测的成像方法 | |
CN106680796A (zh) | 基于频率干涉的平面全息阵列目标三维表面重构方法 | |
CN102596011A (zh) | 光声测量装置 | |
CN103675817B (zh) | 一种基于变换域稀疏的合成孔径雷达侧视三维成像方法 | |
CN109541594A (zh) | 基于涡旋电磁波的条带sar三维成像方法 | |
CN103300825B (zh) | 一种用于超宽带微波检测的乳腺肿瘤特征信号提取方法 | |
CN105677937A (zh) | 一种电磁逆散射重构介质目标的方法 | |
Peng et al. | Convolution back-projection imaging algorithm for downward-looking sparse linear array three dimensional synthetic aperture radar | |
CN103605130B (zh) | 超宽带穿墙雷达隐藏运动目标检测与成像方法 | |
CN101694525A (zh) | 非相干合成孔径激光成像雷达及其聚束模式成像方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160120 |