CN107037434B - 基于和差波束单脉冲测角的成像方法 - Google Patents
基于和差波束单脉冲测角的成像方法 Download PDFInfo
- Publication number
- CN107037434B CN107037434B CN201710269899.0A CN201710269899A CN107037434B CN 107037434 B CN107037434 B CN 107037434B CN 201710269899 A CN201710269899 A CN 201710269899A CN 107037434 B CN107037434 B CN 107037434B
- Authority
- CN
- China
- Prior art keywords
- imaging
- sum
- azimuth
- difference
- difference beam
- 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
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
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
- G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Acoustics & Sound (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明公开了基于和差波束单脉冲测角的成像方法,属于信息传输和图像处理的技术领域。本发明对通道内的原始数据进行加权处理以形成和波束及差波束,对成像波束扫描获取的原始图像进行方位向上的阈值处理以提取强散射点,根据和差波束比估计强散射点方位位置的偏移,结合强散射点方位位置的偏移在成像波束方位向上对强散射点进行重定位,对重定位得到的距离方位图进行扇形视图转换得到二维声学图像,能够明显提高特征目标的分辨效果,从而改善整个图像的清晰度,解决了声呐对探测区域方位成像分辨率不佳的技术问题。
Description
技术领域
本发明公开了基于和差波束单脉冲测角的成像方法,属于信息传输和图像处理的技术领域。
背景技术
针对多波束前视声呐的成像特点,一般采用实波束扫描的成像方法,但是这种方法的方位向分辨率受限于成像波束的主瓣宽度,难以获得更为清晰的水下声学图像。在运动平台声呐对水下成像的过程中,常常采用合成孔径声呐(Synthetic Aperture Sonar,SAS)或多普勒波束锐化(Doppler Beam Sharpening,DBS)技术来实现方位向分辨率的提高,但是仅针对成像区域在声呐系统运动轨迹侧向的情况,常用的声呐系统成像方法并不适用于前视成像。现有的前视成像技术对具有强散射点目标的成像不是太理想,方位向分辨率较低,本发明旨在将单脉冲测角引入到成像过程以解决这个问题。
发明内容
本发明的发明目的是针对上述背景技术的不足,提供了基于和差波束单脉冲测角的成像方法,利用和差波束的比值对特征目标进行精确定位,能够明显提高特征目标的分辨效果,根据方位角估计值对特征目标成像结果的方位向位置进行重新分配,从而改善整个图像的清晰度,解决了声呐对探测区域方位成像分辨率不佳的技术问题。
本发明为实现上述发明目的采用如下技术方案:
基于和差波束单脉冲测角的成像方法,对通道内的原始数据进行加权处理以形成和波束及差波束,对成像波束扫描获取的原始图像进行方位向上的阈值处理以提取强散射点,根据和差波束比估计强散射点方位位置的偏移,结合强散射点方位位置的偏移在成像波束方位向上对强散射点进行重定位,对重定位得到的距离方位图进行扇形视图转换得到二维声学图像。
作为基于和差波束单脉冲测角的成像方法的进一步优化方案,对通道内的原始数据进行加权处理以形成和波束及差波束,具体方法为:以扫描在目标方向形成主瓣的成像波束为和波束,对接收到的回波信号进行差波束加权处理以使回波信号在波束扫描的方向形成零陷,回波信号在波束扫描方向上形成零陷的成像波束即为差波束。
作为基于和差波束单脉冲测角的成像方法的进一步优化方案,采用分水岭阈值算法对成像波束扫描获取的原始图像进行方位向上的阈值处理以提取强散射点。
再进一步的,基于和差波束单脉冲测角的成像方法,根据和差波束比估计强散射点每个距离单元上方位位置的偏移,结合强散射点每个距离单元上方位位置的偏移在成像波束方位向上对强散射点进行重定位。
再进一步的,基于和差波束单脉冲测角的成像方法,根据和差波束比估计强散射点每个距离单元上方位位置的偏移的具体方法为:
首先,构建第m成像波束在第n距离单元上的和差波束比K(n,m):D(n,m)为第m成像波束在第n距离单元上的差波束加权值,S(n,m)为第m成像波束在第n距离单元上的和波束值,Re[·]表示取实部数据,
然后,根据和差波束比与方位角的线性关系确定第n距离单元相对于第m成像波束中心的方位角Δθ(n,m):γ为表征和差波束比与方位角线性关系的斜率常数。
再进一步的,基于和差波束单脉冲测角的成像方法,结合强散射点每个距离单元上方位位置的偏移在成像波束方位向上对强散射点进行重定位的具体方法为:根据第n距离单元相对于第m成像波束中心的方位角Δθ(n,m)以及第m成像波束的指向方位θs(m)确定第m成像波束在第n距离单元上的实际方位位置θ(n,m):θ(n,m)=Δθ(n,m)+θs(m), 为成像波束的方位向间隔。
再进一步的,基于和差波束单脉冲测角的成像方法,重定位得到的距离方位图是通过逐波束非相干累加的方式获得的,逐波束非相干累加的递推公式为:
其中,为第m成像波束扫描第n距离单元的第方位向像素单元形成的距离方位图,为第m-1成像波束扫描第n距离单元的第方位向像素单元形成的距离方位图,为第0成像波束扫描第n距离单元的第方位向像素单元形成的距离方位图,S(n,0)为第0成像波束在第n距离单元上的和波束值,|·|为取绝对值计算,[·]为取整计算,θ(n,0)为第0成像波束在第n距离单元上的实际方位位置,θ0为第0方位向像素单元对应的方位位置,ρ为方位向像素单元的大小。
更进一步的,基于和差波束单脉冲测角的成像方法,方位向像素单元的大小ρ为:θML为和波束的3dB主瓣宽度,N为子孔径阵元数。
本发明采用上述技术方案,具有以下有益效果:
(1)将单脉冲测角引入到成像过程,可以改善水下特征目标的分辨效果,显著提高声呐对探测区域的成像清晰度;
(2)单脉冲成像过程中并不涉及复杂的信号处理运算,有利于工程实现,具有较高的实用价值。
附图说明
图1为本发明成像方法的信号处理流程图。
图2(a)为实波束成像的仿真结果,图2(b)为单脉冲前视成像的仿真结果。
图3(a)为实波束成像的测验结果,图3(b)为单脉冲前视成像的测验结果。
具体实施方式
下面结合附图对发明的技术方案进行详细说明。
本发明的发明构思是在多波束扫描成像方法的基础上,利用和差波束单脉冲测角对扫描波束内强散射点的精确定位性能,将其引入声呐系统对探测区域的成像过程中,可以显著提高声呐图像的质量,使某些具有特征的目标位置更加精确,提出了一种基于和差波束单脉冲测角的前视成像方法。
扫描在目标方向形成主瓣的成像波束为和波束。单脉冲前视成像在波束扫描成像的基础上,对接收到的180路回波信号进行差波束加权以便使回波信号在波束扫描方向形成零陷,得到的差波束加权值表示为D(n,m),回波信号在波束扫描方向上形成零陷的成像波束即为差波束。利用和差波束的比值对目标进行精确定位,这样可以明显提高特征目标(即,强散射点)的分辨效果,从而实现对整个图像清晰度的改善,其信号处理流程如图1所示。
首先,进行方位向上的阈值处理,利用分水岭阈值算法寻找背景噪声区域并提取强散射点信息,假定最大峰值30dB以下的像素对应背景噪声,最大峰值10dB以上对应强散射点。
然后,在强散射点每个距离单元上构造和差波束比,从而对强散射点每个距离单元的精确方位位置进行估计,K(n,m)表示第m成像波束在第n距离单元上的单脉冲和差波束比,则,
其中,S(n,m)表示第m成像波束在第n距离单元上的和波束值,D(n,m)表示第m成像波束在第n距离单元上的差波束加权值,Re[·]表示取实部数据。
接着,根据和差波束比与方位角的关系可得第n距离单元相对于第m成像波束中心的方位角Δθ(n,m)为:
其中,γ为表征和差波束比与方位角线性关系的斜率常数,可以通过实验确定具体数值,
读入成像波束序号,第m成像波束在第n距离单元上的实际方位位置θ(n,m)可表示为:
θ(n,m)=Δθ(n,m)+θs(m),
式中:θs(m)为第m成像波束的指向方位, 为成像波束的方位向间隔,
最后,定义重新定位后的距离方位图像中的每个像素单元为I(n,k),k为方位向像素单元序号,为了实现图像方位向分辨率的提高,方位向像素单元的大小ρ选择为其中,θML为和波束的3dB主瓣宽度,N为子孔径阵元数,N=0,1,……,91,设0号方位向像素单元对应的方位位置为θ0,则,
上述处理过程仅针对方位向上的成像数据,而对应的距离位置不变。完成上式处理后,成像波束的回波信号重新定位于距离方位图像中的I(n,k)像素单元内,重定位得到的距离方位图通过逐波束非相干累加的方式生成,即如下递推公式:
递推公式中,m=1,2,……,540,I0表示第0成像波束扫描距离单元形成的距离方位图,Im表示第m成像波束扫描距离单元形成的距离方位图,Im-1表示第m-1成像波束扫描距离单元形成的距离方位图,θ(n,0)第0成像波束在第n距离单元上的实际方位位置,S(n,0)为第0成像波束在第n距离单元上的和波束值。
对强散射点每个距离单元都按照上述步骤进行精确的方位估计并更新距离方位图,通过逐波束非相干累加的方式生成成像波束扫描每个距离单元的距离方位图像。
对成像波束扫描每个距离单元的距离方位图像做扇形视图转换,把成像数据从极坐标转换到直角坐标,即,将距离方位图像转换为最终的二维声学视图。
下面通过仿真实验和实际测验来阐述本发明技术方案的优越性。
按照表1设置的仿真参数对实波束成像结果以及单脉冲前视成像结果进行仿真,仿真结果如图2(a)、图2(b)所示。
参数名称 | 参数设置 |
和波束主瓣宽度 | 1.2° |
扫描波束个数 | 540 |
相邻波束间隔 | 0.168° |
脉冲重复周期 | 100ms |
脉冲宽度 | 0.05ms |
中心频率 | 450kHz |
表1 仿真参数的设置值
由仿真结果对比可以看出:由于实波束成像过程中,目标场景的方位向分辨率完全依赖于扫描波束的主瓣宽度,因此每个点目标的成像数据均在方位向上占据了若干个位置,使得图像难以反映其精确的方位;单脉冲前视成像方法可以明显改善这样的情况,经过单脉冲成像处理之后,图像中原来实波束成像得到的方位向上的“线”压缩成了“点”,并且点的位置精确反映了目标的真实位置。
为了验证单脉冲成像算法的实用性,对声呐系统在消声水池录制的两组180通道原始回波数据同时进行实波束成像和单脉冲成像处理,三角架和圆环成像结果的矩形视图如图3(a)、图3(b)所示。
Claims (6)
1.基于和差波束单脉冲测角的成像方法,其特征在于,以扫描在目标方向形成主瓣的成像波束为和波束,对接收到的多通道回波信号进行加权处理以形成差波束加权处理以使回波信号在波束扫描的方向形成零陷,对成像波束扫描获取的原始图像进行方位向上的阈值处理以提取强散射点,根据和差波束比估计强散射点每个距离单元上方位位置的偏移,结合强散射点每个距离单元上方位位置的偏移在成像波束方位向上对强散射点进行重定位,对重定位得到的距离方位图进行扇形视图转换得到二维声学图像。
2.根据权利要求1所述基于和差波束单脉冲测角的成像方法,其特征在于,采用分水岭阈值算法对成像波束扫描获取的原始图像进行方位向上的阈值处理以提取强散射点。
3.根据权利要求1所述基于和差波束单脉冲测角的成像方法,其特征在于,根据和差波束比估计强散射点每个距离单元上方位位置的偏移的具体方法为:
首先,构建第m成像波束在第n距离单元上的和差波束比K(n,m):D(n,m)为第m成像波束在第n距离单元上的差波束加权值,S(n,m)为第m成像波束在第n距离单元上的和波束值,Re[·]表示取实部数据,
然后,根据和差波束比与方位角的线性关系确定第n距离单元相对于第m成像波束中心的方位角Δθ(n,m):γ为表征和差波束比与方位角线性关系的斜率常数。
4.根据权利要求3所述基于和差波束单脉冲测角的成像方法,其特征在于,结合强散射点每个距离单元上方位位置的偏移在成像波束方位向上对强散射点进行重定位的具体方法为:根据第n距离单元相对于第m成像波束中心的方位角Δθ(n,m)以及第m成像波束的指向方位θs(m)确定第m成像波束在第n距离单元上的实际方位位置θ(n,m):θ(n,m)=Δθ(n,m)+θs(m), 为成像波束的方位向间隔。
5.根据权利要求4所述基于和差波束单脉冲测角的成像方法,其特征在于,重定位得到的距离方位图是通过逐波束非相干累加的方式获得的,逐波束非相干累加的递推公式为:
其中,为第m成像波束扫描第n距离单元的第方位向像素单元形成的距离方位图,为第m-1成像波束扫描第n距离单元的第方位向像素单元形成的距离方位图,为第0成像波束扫描第n距离单元的第方位向像素单元形成的距离方位图,S(n,0)为第0成像波束在第n距离单元上的和波束值,|·|为取绝对值计算,[·]为取整计算,θ(n,0)为第0成像波束在第n距离单元上的实际方位位置,θ0为第0方位向像素单元对应的方位位置,ρ为方位向像素单元的大小。
6.根据权利要求5所述基于和差波束单脉冲测角的成像方法,其特征在于,所述方位向像素单元的大小ρ为:θML为和波束的3dB主瓣宽度,N为子孔径阵元数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710269899.0A CN107037434B (zh) | 2017-04-24 | 2017-04-24 | 基于和差波束单脉冲测角的成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710269899.0A CN107037434B (zh) | 2017-04-24 | 2017-04-24 | 基于和差波束单脉冲测角的成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107037434A CN107037434A (zh) | 2017-08-11 |
CN107037434B true CN107037434B (zh) | 2019-12-24 |
Family
ID=59535542
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710269899.0A Expired - Fee Related CN107037434B (zh) | 2017-04-24 | 2017-04-24 | 基于和差波束单脉冲测角的成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107037434B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109471093A (zh) * | 2018-11-07 | 2019-03-15 | 中国人民解放军国防科技大学 | 单脉冲雷达和差相关目标检测方法及系统 |
CN115201821B (zh) * | 2022-07-25 | 2023-09-05 | 苏州威陌电子信息科技有限公司 | 基于强目标成像对消的小目标检测方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS6050415A (ja) * | 1983-08-31 | 1985-03-20 | Unyusho Kowan Gijutsu Kenkyusho | 超音波海底沈埋異常物探査装置 |
US7667638B1 (en) * | 2009-02-26 | 2010-02-23 | Northrop Grumman Systems Corporation | Detection and resolution of closely spaced targets in a monopulse system |
CN101937085B (zh) * | 2010-09-01 | 2012-02-22 | 华南理工大学 | 一种基于超分辨率图像重建的合成孔径声纳成像方法 |
CN103630887B (zh) * | 2013-11-11 | 2015-12-09 | 南京航空航天大学 | 一种多波束形成方法及使用该方法的多波束声纳 |
CN105093229B (zh) * | 2014-05-19 | 2018-03-27 | 日本电气株式会社 | 多个目标的定位方法及装置 |
CN105607055B (zh) * | 2015-12-24 | 2017-12-19 | 西安电子科技大学 | 一种基于天线方向图的机载雷达单脉冲前视成像方法 |
-
2017
- 2017-04-24 CN CN201710269899.0A patent/CN107037434B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN107037434A (zh) | 2017-08-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101887120B (zh) | 一种用于单波束机械扫描雷达测量目标方位角的方法 | |
CN111352107B (zh) | 基于多通道数字和差的单脉冲跟踪与成像方法 | |
CN112444811A (zh) | 一种综合mimo雷达和isar的目标检测和成像方法 | |
CN102608597A (zh) | 基于非完全数据解卷积的实孔径前视成像方法 | |
CN1769925A (zh) | 一种合成孔径雷达运动目标成像方法 | |
CN110146889B (zh) | 基于最优回波子区域选择的大转角isar成像方法 | |
CN109655819B (zh) | 一种基于实孔径多普勒波束锐化的杂波抑制三维成像方法 | |
CN111352105B (zh) | 用于目标跟踪的测角方法 | |
CN109061638B (zh) | 相控阵近距离数字成像方法 | |
CN114002673A (zh) | 星载无源sar的非合作信号感知系统与多维参数估计方法 | |
CN110907938B (zh) | 一种近场的快速下视合成孔径三维成像方法 | |
CN107037434B (zh) | 基于和差波束单脉冲测角的成像方法 | |
CN112034459A (zh) | 基于天线方向图补偿的直线合成孔径雷达成像方法和系统 | |
CN108107427A (zh) | 基于超分辨技术的机载/弹载阵列雷达前视成像方法 | |
CN112433210A (zh) | 一种双站前视探地雷达快速时域成像方法 | |
CN113156444B (zh) | 一种基于运动补偿的多波束声呐高精度成像方法 | |
CN111273290A (zh) | 基于预成像曲线轨迹的三维sar成像方法 | |
CN111929682B (zh) | 基于距离向分块的机载dbf-tops sar成像方法 | |
CN114252878A (zh) | 一种基于逆合成孔径雷达对运动目标进行成像及横向定标的方法 | |
CN111707996B (zh) | 基于改进grft-stap的geo星机sar动目标检测方法 | |
CN113534077A (zh) | 一种雷达辐射源威力反演方法、装置及电子设备 | |
CN115865142B (zh) | 一种雷达天线波束指向校准方法、系统、装置及存储介质 | |
Al-Ibadi et al. | DEM extraction of the basal topography of the Canadian archipelago ICE caps via 2D automated layer-tracker | |
CN113359196B (zh) | 基于子空间法和dbf的多目标生命体征探测方法 | |
CN115390070A (zh) | 一种机动平台单通道短孔径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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20191224 |