CN110244285A - 一种合成孔径成像声呐斜距向非均匀分段补偿方法及系统 - Google Patents
一种合成孔径成像声呐斜距向非均匀分段补偿方法及系统 Download PDFInfo
- Publication number
- CN110244285A CN110244285A CN201910553855.XA CN201910553855A CN110244285A CN 110244285 A CN110244285 A CN 110244285A CN 201910553855 A CN201910553855 A CN 201910553855A CN 110244285 A CN110244285 A CN 110244285A
- Authority
- CN
- China
- Prior art keywords
- distance
- subsegment
- echo data
- homogeneous
- oblique distance
- 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.)
- Granted
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
-
- 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/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52004—Means for monitoring or calibrating
-
- 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/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/539—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明属于成像声呐信号处理技术领域,特别涉及一种合成孔径成像声呐斜距向非均匀分段补偿方法,该方法包括:获取合成孔径声呐在整个大测绘带区域内接收地距向的回波数据;将获取的回波数据沿地距向均匀分成N段,通过地距向斜距映射,获得每个斜距向的非均匀子段的回波数据;根据横荡误差和升沉误差,计算每个斜距向的非均匀子段的声程差;基于获得的每个斜距向的非均匀子段的声程差和斜距向的非均匀子段的回波数据,对N个斜距向的非均匀子段的回波数据进行运动补偿,从而实现整个测绘带区域内的接收的地距向的回波数据的运动补偿。
Description
技术领域
本发明属于成像声呐信号处理技术领域,特别涉及一种合成孔径成像声呐斜距向非均匀分段补偿方法及系统。
背景技术
合成孔径声呐利用小孔径基阵沿方位向的移动合成虚拟大孔径,然后根据空间位置和相位关系对不同位置处的接收回波进行相干叠加处理,获得高分辨率的图像。合成孔径声呐对运动状态要求比较苛刻,要想实现高分辨率成像,需要做理想的匀速直线运动。
但是,在实际外场环境中,受海浪、暗涌和母船等因素影响,合成孔径声呐不可避免地会产生运动误差,运动误差将会引起声呐回波的幅度和相位发生变化,从而导致目标出现散焦、重影等现象,严重时甚至不能成像。因此,需要对声呐回波进行运动补偿。
目前,现有的运动补偿方法包括运动误差估计和运动误差补偿两个步骤。针对运动误差补偿,当测绘带较窄时,可认为测绘带内的目标具有相同的声程差,选取测绘带中心为参考距离,计算声程差,对整个测绘带的回波数据补偿。但是,对于大测绘带合成孔径声呐,运动误差在近距和远距处引起的声程差有所不同,具有明显的斜距向空变性,窄带假设不再成立,无法对不同斜距处的回波数据区别处理,进而无法实现对声呐回波准确的运动补偿。
发明内容
本发明的目的在于,为解决现有的运动补偿方法存在的上述缺陷,本发明提出了一种合成孔径成像声呐斜距向非均匀分段补偿方法,该方法基于大测绘带;克服了大测绘带合成孔径声呐运动误差引起的声程差在斜距向存在空变性的问题,通过沿地距向对回波数据均匀分段,借助距离压缩效应,得到多个斜距向的非均匀子段,顺应运动误差引起的声程差“先快后慢”的空变规律,实现更准确的运动补偿。该方法适用于大测绘带合成孔径成像声呐的运动补偿,对空变性较大的近距目标补偿效果尤其明显,同时计算量不大,便于实时操作。
为了实现上述目的,本发明提出了一种合成孔径成像声呐斜距向非均匀分段补偿方法,其特征在于,该方法包括:
获取合成孔径声呐在整个大测绘带区域内接收地距向的回波数据;
将获取的回波数据沿地距向均匀分成N段,通过地距向斜距映射,获得每个斜距向的非均匀子段的回波数据;
根据横荡误差和升沉误差,计算每个斜距向的非均匀子段的声程差;
基于获得的每个斜距向的非均匀子段的声程差和斜距向的非均匀子段的回波数据,对N个斜距向的非均匀子段的回波数据进行运动补偿,从而实现整个测绘带区域内的接收的地距向的回波数据的运动补偿。
作为上述技术方案的改进之一,所述将获取的回波数据沿地距向均匀分成N段,获得每个斜距向的非均匀子段的回波数据;具体包括:
利用合成孔径声呐斜距成像的距离压缩效应,通过将获取的地距向的回波数据沿地距向均匀分成N段,则第i个斜距向的非均匀子段的回波数据
其中,i表示第i个子段,1≤i≤N;为合成孔径声呐在整个测绘带区域内接收地距向的回波数据;τ为时间变量;xb表示第i个子段的起始点,xe表示第i个子段的终点,Rmin为声呐测绘带最近端与声呐之间的距离,称为最短距离;h为声呐距底高度;u(τ)为阶跃函数;c为声速,c=1500m/s;δx为每个地距向均匀分段的子段长度;
Rmax为声呐测绘带最远端与声呐之间的距离,称为最长距离。
作为上述技术方案的改进之一,所述根据横荡误差和升沉误差,计算每个斜距向的非均匀子段的声程差;具体包括:
以第i个斜距向的非均匀子段中心为参考距离,计算第i个斜距向的非均匀子段的声程差Δri为:
其中,Δx为声呐运动产生的横荡误差;Δh为声呐运动产生的升沉误差;
作为上述技术方案的改进之一,所述基于获得的每个斜距向的非均匀子段的声程差和斜距向的非均匀子段的回波数据,对N个斜距向的非均匀子段的回波数据进行运动补偿,从而实现整个测绘带区域内的接收的地距向的回波数据的运动补偿;具体包括:
通过时域移位的方法,基于获得的第i个斜距向的非均匀子段的声程差和第i个斜距向的非均匀子段的回波数据,对N个斜距向的非均匀子段的回波数据进行运动补偿:
其中,s′(τ)为N个斜距向的非均匀子段的回波数据进行运动补偿后的回波数据,fc为载波频率,j为虚数单位;
基于对N个斜距向的非均匀子段的回波数据进行运动补偿,从而实现整个测绘带区域内的接收的地距向的回波数据的运动补偿。
作为上述技术方案的改进之一,所述基于获得的每个斜距向的非均匀子段的声程差和斜距向的非均匀子段的回波数据,对N个斜距向的非均匀子段的回波数据进行运动补偿,从而实现整个测绘带区域内的接收的地距向的回波数据的运动补偿;具体包括:
通过频域补相的方法,基于获得的第i个斜距向的非均匀子段的声程差和第i个斜距向的非均匀子段的回波数据,对N个斜距向的非均匀子段的回波数据进行运动补偿:
其中,s′(τ)为N个斜距向的非均匀子段的回波数据进行运动补偿后的回波数据;ifft()为傅里叶反变换;f为斜距向频率变量;
基于对N个斜距向的非均匀子段的回波数据进行运动补偿,从而实现整个测绘带区域内的接收的地距向的回波数据的运动补偿。
一种合成孔径成像声呐斜距向非均匀分段补偿系统,该系统包括:数据获取模块、分段模块、处理模块和运动补偿模块;
所述数据获取模块,用于获取合成孔径声呐在整个大测绘带区域内接收地距向的回波数据;
所述分段模块,用于将获取的回波数据沿地距向均匀分成N段,通过地距向斜距映射,获得每个斜距向的非均匀子段的回波数据;
所述处理模块,用于根据横荡误差和升沉误差,计算每个斜距向的非均匀子段的声程差;
所述运动补偿模块,用于基于获得的每个斜距向的非均匀子段的声程差和斜距向的非均匀子段的回波数据,对N个斜距向的非均匀子段的回波数据进行运动补偿,从而实现整个测绘带区域内的接收的地距向的回波数据的运动补偿。
作为上述技术方案的改进之一,所述分段模块具体包括:
数据接收单元,用于接收地距向的回波数据;和
处理单元,用于根据合成孔径声呐斜距成像的距离压缩效应,通过将接收的地距向的回波数据沿地距向均匀分成N段,则第i个斜距向的非均匀子段的回波数据
其中,i表示第i个子段,1≤i≤N;为合成孔径声呐在整个测绘带区域内接收地距向的回波数据;τ为时间变量;xb表示第i个子段的起始点,xe表示第i个子段的终点,Rmin为声呐测绘带最近端与声呐之间的距离,称为最短距离;h为声呐距底高度;u(τ)为阶跃函数;c为声速,c=1500m/s;δx为每个地距向均匀分段的子段长度;
Rmax为声呐测绘带最远端与声呐之间的距离,称为最长距离。
作为上述技术方案的改进之一,所述处理模块具体包括:
根据公式(6),以第i个斜距向的非均匀子段中心为参考距离,计算第i个斜距向的非均匀子段的声程差Δri:
其中,Δx为声呐运动产生的横荡误差;Δh为声呐运动产生的升沉误差;
作为上述技术方案的改进之一,所述运动补偿模块具体包括:
根据时域移位的方法,对N个斜距向的非均匀子段的回波数据进行运动补偿:
其中,s′(τ)为N个斜距向的非均匀子段的回波数据进行运动补偿后的回波数据,fc为载波频率,j为虚数单位;
基于对N个斜距向的非均匀子段的回波数据进行运动补偿,从而实现整个测绘带区域内的接收的地距向的回波数据的运动补偿。
作为上述技术方案的改进之一,所述运动补偿模块具体包括:
根据频域补相的方法,对N个斜距向的非均匀子段的回波数据进行运动补偿:
其中,s′(τ)为N个斜距向的非均匀子段的回波数据进行运动补偿后的回波数据;ifft()为傅里叶反变换;f为斜距向频率变量;
基于对N个斜距向的非均匀子段的回波数据进行运动补偿,从而实现整个测绘带区域内的接收的地距向的回波数据的运动补偿。
本发明相比于现有技术的有益效果在于:
本发明的方法可以有效消除大测绘带合成孔径声呐运动误差在斜距向引起的空变性影响,提高运动补偿的准确性,尤其对近处目标成像结果改善明显;同时计算量不大,便于实时操作。
附图说明
图1是本发明的一种合成孔径成像声呐斜距向非均匀分段补偿方法考虑存在横荡误差和升沉误差时成像声呐几何端射图;
图2是本发明的一种合成孔径成像声呐斜距向非均匀分段补偿方法采用的分段策略的示意图。
具体实施方式
现结合附图对本发明作进一步的描述。
如图1所示,本发明提出了一种合成孔径成像声呐斜距向非均匀分段补偿方法,该方法具体包括:
获取合成孔径声呐在整个大测绘带区域内接收地距向的回波数据;其中,所述整个测绘带区域具体包括:最短距离Rmin与最长距离Rmax之间的成像区域;其中,最短距离Rmin为声呐测绘带最近端与声呐之间的距离,称为最短距离;最长距离Rmax为声呐测绘带最远端与声呐之间的距离,称为最长距离。
将获取的回波数据沿地距向均匀分成N段,通过地距向斜距的映射,获得每个斜距向的非均匀子段的回波数据;具体包括:
利用合成孔径声呐斜距成像的距离压缩效应,通过将获取的地距向的回波数据沿地距向均匀分成N段,则第i个斜距向的非均匀子段的回波数据
其中,i表示第i个子段,1≤i≤N;为合成孔径声呐在整个测绘带区域内接收地距向的回波数据;τ为时间变量;xb表示第i个子段的起始点,xe表示第i个子段的终点,Rmin为声呐测绘带最近端与声呐之间的距离,称为最短距离;h为声呐距底高度;u(τ)为阶跃函数;c为声速,c=1500m/s;δx为每个地距向均匀分段的子段长度;
Rmax为声呐测绘带最远端与声呐之间的距离,称为最远距离。
根据横荡误差和升沉误差,计算每个斜距向的非均匀子段的声程差;具体地,
如图1所示,图1是合成孔径声呐成像几何端射图,在xoz坐标系下,假设某时刻声呐理想位置处于点A,坐标为(0,h),由于横荡误差Δx和升沉误差Δh的影响,阵元实际处于点B,坐标为(-Δx,h+Δh)。在声呐波束照射范围内,存在某个目标M,坐标为(x,0)。在理想位置点A处的声呐与目标M的斜距为
而实际处于点B时的斜距为
将实际斜距与理想斜距之间的差值定义为声程差,则目标M与声呐之间的声程差可近似表示为
首先考虑存在横荡误差和升沉误差的情况下,运动误差在斜距向引起的空变效应。基于公式(3)可以发现,声程差不仅仅是关于横荡误差和升沉误差的函数,同时也是关于斜距的函数,不同斜距处的目标,声程差各不相同,这就是运动误差在斜距向引起的空变效应。对于窄测绘带合成孔径声呐,可以忽略声程差随目标斜距变化的空变特性,直接利用测绘带中心处的声程差对整个测绘带回波数据进行运动补偿。但对于大测绘带合成孔径声呐,近距和远距处的声程差相差较大,空变效应明显,不能简单地以测绘带中心处的声程差一概等效,需要区别处理不同斜距处的回波数据。
因此,对于大测绘带合成孔径声呐,以第i个斜距向的非均匀子段中心为参考距离,计算第i个斜距向的非均匀子段的声程差Δri为:
其中,Δx为声呐运动产生的横荡误差;Δh为声呐运动产生的升沉误差;
基于获得的每个斜距向的非均匀子段的声程差和斜距向的非均匀子段的回波数据,对N个斜距向的非均匀子段的回波数据进行运动补偿,从而实现整个测绘带区域内的接收的地距向的回波数据的运动补偿。具体地,
通过时域移位的方法,基于获得的第i个斜距向的非均匀子段的声程差和第i个斜距向的非均匀子段的回波数据,对N个斜距向的非均匀子段的回波数据进行运动补偿:
其中,s′(τ)为N个斜距向的非均匀子段的回波数据进行运动补偿后的回波数据,fc为载波频率,j为虚数单位;τ为时间变量;
基于对N个斜距向的非均匀子段的回波数据进行运动补偿,从而实现整个测绘带区域内的接收的地距向的回波数据的运动补偿;
或者通过频域补相的方法,基于获得的第i个斜距向的非均匀子段的声程差和第i个斜距向的非均匀子段的回波数据,对N个斜距向的非均匀子段的回波数据进行运动补偿:
其中,ifft()为傅里叶反变换;f为斜距向频率变量;
基于对N个斜距向的非均匀子段的回波数据进行运动补偿,从而实现整个测绘带区域内的接收的地距向的回波数据的运动补偿。
基于补偿准确度和运算量的综合考虑,通常将最短距离Rmin与最长距离Rmax之间的成像区域沿斜距向均匀划分为多个斜距向的非均匀子段,当第i个斜距向的非均匀子段足够窄时,该子段对应的回波数据的声程差相差不大,可按照第i个斜距向的非均匀子段内中心点作为参考距离,计算第i个斜距向的非均匀子段的回波数据的声程差对该段回波数据一致补偿,这样大大节省了运算量,也非常方便。但是,按照斜距向均匀分段的策略,各子段的长度相同,空变性却差异较大。斜距向均匀分段对回波数据做一致的处理,势必会造成近距处的子段残留相对较大的补偿误差,因此有必要区别处理最短距离Rmin与最长距离Rmax的斜距向的非均匀子段。
由于运动误差引起的空变呈现“先快后慢”的规律,最短距离Rmin的斜距向的非均匀子段空变程度大于最长距离Rmax的斜距向的非均匀子段,为了顺应运动误差空变规律,实现更精确的运动补偿,应当缩短最短距离Rmin处空变性较大的子段长度,适当延长最长距离Rmax处空变性较小的子段长度,即采用非均匀分段策略,如图2所示,图2给出了分段策略的示意图,在地距向均匀分段的子段,映射到斜距向为非均匀的子段,其中,斜距向最短距离处的子段长度较短,记为Rmin,斜距向最大距离处的子段长度较长,记为Rmax。利用合成孔径声呐斜距成像的距离压缩效应,通过对地距向的回波数据沿地距向均匀分成N段,获得第i个斜距向的非均匀子段的回波数据。
利用合成孔径声呐斜距成像的距离压缩效应,通过对地距向的回波数据沿斜距向均匀分成N段,获取N个斜距向的非均匀子段。距离合成孔径声呐越近的斜距向的非均匀子段,在斜距向上的投影被压缩越严重,即缩短了最短距离Rmin处空变性较大的子段长度;距离合成孔径声呐较远的子段,由于其空变性较小,适当延长并不会带来太大影响。
本发明还提供了一种合成孔径成像声呐斜距向非均匀分段补偿系统,其特征在于,该系统包括:数据获取模块、分段模块、处理模块和运动补偿模块;
所述数据获取模块,用于获取合成孔径声呐在整个大测绘带区域内接收地距向的回波数据;
所述分段模块,用于将获取的回波数据沿地距向均匀分成N段,通过地距向斜距映射,获得每个斜距向的非均匀子段的回波数据;
所述分段模块具体包括:
数据接收单元,用于接收地距向的回波数据;和
处理单元,用于根据合成孔径声呐斜距成像的距离压缩效应,通过将接收的地距向的回波数据沿地距向均匀分成N段,则第i个斜距向的非均匀子段的回波数据
其中,i表示第i个子段,1≤i≤N;为合成孔径声呐在整个测绘带区域内接收地距向的回波数据;τ为时间变量;xb表示第i个子段的起始点,xe表示第i个子段的终点,Rmin为声呐测绘带最近端与声呐之间的距离,称为最短距离;h为声呐距底高度;u(τ)为阶跃函数;c为声速,c=1500m/s;δx为每个地距向均匀分段的子段长度;
Rmax为声呐测绘带最远端与声呐之间的距离,称为最长距离。
所述处理模块,用于根据横荡误差和升沉误差,计算每个斜距向的非均匀子段的声程差;
所述处理模块具体包括:
根据公式(6),以第i个斜距向的非均匀子段中心为参考距离,计算第i个斜距向的非均匀子段的声程差Δri:
其中,Δx为声呐运动产生的横荡误差;Δh为声呐运动产生的升沉误差。
所述运动补偿模块,用于基于获得的每个斜距向的非均匀子段的声程差和斜距向的非均匀子段的回波数据,对N个斜距向的非均匀子段的回波数据进行运动补偿,从而实现整个测绘带区域内的接收的地距向的回波数据的运动补偿。
所述运动补偿模块具体包括:
根据时域移位的方法,基于获得的第i个斜距向的非均匀子段的声程差和第i个斜距向的非均匀子段的回波数据,对N个斜距向的非均匀子段的回波数据进行运动补偿:
其中,s′(τ)为N个斜距向的非均匀子段的回波数据进行运动补偿后的回波数据,fc为载波频率,j为虚数单位;
基于对N个斜距向的非均匀子段的回波数据进行运动补偿,从而实现整个测绘带区域内的接收的地距向的回波数据的运动补偿。
所述运动补偿模块具体还包括:
根据频域补相的方法,基于获得的第i个斜距向的非均匀子段的声程差和第i个斜距向的非均匀子段的回波数据,对N个斜距向的非均匀子段的回波数据进行运动补偿:
其中,s′(τ)为N个斜距向的非均匀子段的回波数据进行运动补偿后的回波数据;ifft()为傅里叶反变换;f为斜距向频率变量;
基于对N个斜距向的非均匀子段的回波数据进行运动补偿,从而实现整个测绘带区域内的接收的地距向的回波数据的运动补偿。
最后所应说明的是,以上实施例仅用以说明本发明的技术方案而非限制。尽管参照实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,对本发明的技术方案进行修改或者等同替换,都不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (10)
1.一种合成孔径成像声呐斜距向非均匀分段补偿方法,其特征在于,该方法包括:
获取合成孔径声呐在整个大测绘带区域内接收地距向的回波数据;
将回波数据沿地距向均匀分成N段,通过地距向斜距映射,获得每个斜距向的非均匀子段的回波数据;
根据横荡误差和升沉误差,计算每个斜距向的非均匀子段的声程差;
基于获得的每个斜距向的非均匀子段的声程差和斜距向的非均匀子段的回波数据,对N个斜距向的非均匀子段的回波数据进行运动补偿,从而实现整个测绘带区域内的接收的地距向的回波数据的运动补偿。
2.根据权利要求1所述的方法,其特征在于,所述将获取的回波数据沿地距向均匀分成N段,获得每个斜距向的非均匀子段的回波数据;具体包括:
利用合成孔径声呐斜距成像的距离压缩效应,通过将获取的地距向的回波数据沿地距向均匀分成N段,则第i个斜距向的非均匀子段的回波数据
其中,i表示第i个子段,1≤i≤N;为合成孔径声呐在整个测绘带区域内接收地距向的回波数据;τ为时间变量;xb表示第i个子段的起始点,xe表示第i个子段的终点,Rmin为声呐测绘带最近端与声呐之间的距离,称为最短距离;h为声呐距底高度;u(τ)为阶跃函数;c为声速,c=1500m/s;δx为每个地距向均匀分段的子段长度;
Rmax为声呐测绘带最远端与声呐之间的距离,称为最长距离。
3.根据权利要求2所述的方法,其特征在于,所述根据横荡误差和升沉误差,计算每个斜距向的非均匀子段的声程差;具体包括:
以第i个斜距向的非均匀子段中心为参考距离,计算第i个斜距向的非均匀子段的声程差Δri为:
其中,Δx为声呐运动产生的横荡误差;Δh为声呐运动产生的升沉误差。
4.根据权利要求2或3所述的方法,其特征在于,所述基于获得的每个斜距向的非均匀子段的声程差和斜距向的非均匀子段的回波数据,对N个斜距向的非均匀子段的回波数据进行运动补偿,从而实现整个测绘带区域内的接收的地距向的回波数据的运动补偿;具体包括:
通过时域移位的方法,基于获得的第i个斜距向的非均匀子段的声程差和第i个斜距向的非均匀子段的回波数据,对N个斜距向的非均匀子段的回波数据进行运动补偿:
其中,s′(τ)为N个斜距向的非均匀子段的回波数据进行运动补偿后的回波数据,fc为载波频率,j为虚数单位;
基于对N个斜距向的非均匀子段的回波数据进行运动补偿,从而实现整个测绘带区域内的接收的地距向的回波数据的运动补偿。
5.根据权利要求2或3所述的方法,其特征在于,所述基于获得的每个斜距向的非均匀子段的声程差和斜距向的非均匀子段的回波数据,对N个斜距向的非均匀子段的回波数据进行运动补偿,从而实现整个测绘带区域内的接收的地距向的回波数据的运动补偿;具体包括:
通过频域补相的方法,基于获得的第i个斜距向的非均匀子段的声程差和第i个斜距向的非均匀子段的回波数据,对N个斜距向的非均匀子段的回波数据进行运动补偿:
其中,s′(τ)为N个斜距向的非均匀子段的回波数据进行运动补偿后的回波数据;ifft()为傅里叶反变换;f为斜距向频率变量,fc为载波频率;
基于对N个斜距向的非均匀子段的回波数据进行运动补偿,从而实现整个测绘带区域内的接收的地距向的回波数据的运动补偿。
6.一种合成孔径成像声呐斜距向非均匀分段补偿系统,其特征在于,该系统包括:数据获取模块、分段模块、处理模块和运动补偿模块;
所述数据获取模块,用于获取合成孔径声呐在整个大测绘带区域内接收地距向的回波数据;
所述分段模块,用于将获取的回波数据沿地距向均匀分成N段,通过地距向斜距映射,获得每个斜距向的非均匀子段的回波数据;
所述声程差计算模块,用于根据横荡误差和升沉误差,计算每个斜距向的非均匀子段的声程差;
所述运动补偿模块,用于基于获得的每个斜距向的非均匀子段的声程差和斜距向的非均匀子段的回波数据,对N个斜距向的非均匀子段的回波数据进行运动补偿,从而实现整个测绘带区域内的接收的地距向的回波数据的运动补偿。
7.根据权利要求6所述的系统,其特征在于,所述分段模块具体包括:
数据接收单元,用于接收地距向的回波数据;和
处理单元,用于根据合成孔径声呐斜距成像的距离压缩效应,通过将接收的地距向的回波数据沿地距向均匀分成N段,则第i个斜距向的非均匀子段的回波数据
其中,i表示第i个子段,1≤i≤N;为合成孔径声呐在整个测绘带区域内接收地距向的回波数据;τ为时间变量;xb表示第i个子段的起始点,xe表示第i个子段的终点,Rmin为声呐测绘带最近端与声呐之间的距离,称为最短距离;h为声呐距底高度;u(τ)为阶跃函数;c为声速,c=1500m/s;δx为每个地距向均匀分段的子段长度;
Rmax为声呐测绘带最远端与声呐之间的距离,称为最长距离。
8.根据权利要求7所述的系统,其特征在于,所述声程差计算模块具体包括:
根据公式(6),以第i个斜距向的非均匀子段中心为参考距离,计算第i个斜距向的非均匀子段的声程差Δri:
其中,Δx为声呐运动产生的横荡误差;Δh为声呐运动产生的升沉误差。
9.根据权利要求7或8所述的系统,其特征在于,所述运动补偿模块具体包括:
根据时域移位的方法,对N个斜距向的非均匀子段的回波数据进行运动补偿:
其中,s′(τ)为N个斜距向的非均匀子段的回波数据进行运动补偿后的回波数据,fc为载波频率,j为虚数单位;
基于对N个斜距向的非均匀子段的回波数据进行运动补偿,从而实现整个测绘带区域内的接收的地距向的回波数据的运动补偿。
10.根据权利要求7或8所述的系统,其特征在于,所述运动补偿模块具体包括:
根据频域补相的方法,对N个斜距向的非均匀子段的回波数据进行运动补偿:
其中,s′(τ)为N个斜距向的非均匀子段的回波数据进行运动补偿后的回波数据;ifft()为傅里叶反变换;f为斜距向频率变量;
基于对N个斜距向的非均匀子段的回波数据进行运动补偿,从而实现整个测绘带区域内的接收的地距向的回波数据的运动补偿。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910553855.XA CN110244285B (zh) | 2019-06-25 | 2019-06-25 | 一种合成孔径成像声呐斜距向非均匀分段补偿方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910553855.XA CN110244285B (zh) | 2019-06-25 | 2019-06-25 | 一种合成孔径成像声呐斜距向非均匀分段补偿方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110244285A true CN110244285A (zh) | 2019-09-17 |
CN110244285B CN110244285B (zh) | 2021-02-09 |
Family
ID=67889327
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910553855.XA Active CN110244285B (zh) | 2019-06-25 | 2019-06-25 | 一种合成孔径成像声呐斜距向非均匀分段补偿方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110244285B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113156444A (zh) * | 2021-06-02 | 2021-07-23 | 杭州电子科技大学 | 一种基于运动补偿的多波束声呐高精度成像方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101551452A (zh) * | 2008-04-01 | 2009-10-07 | 中国科学院声学研究所 | 一种合成孔径声纳联合运动补偿方法及系统 |
CN101762810A (zh) * | 2009-12-15 | 2010-06-30 | 中国科学院声学研究所 | 一种大测绘带下合成孔径声纳运动补偿方法 |
CN101799537A (zh) * | 2009-12-15 | 2010-08-11 | 中国科学院声学研究所 | 一种基于双侧回波数据的合成孔径声纳运动补偿方法 |
CN103744081A (zh) * | 2014-01-27 | 2014-04-23 | 中国科学院电子学研究所 | 一种机载圆迹合成孔径雷达高精度三维成像补偿方法 |
WO2017062553A1 (en) * | 2015-10-08 | 2017-04-13 | Mayo Foundation For Medical Education And Research | Systems and methods for ultrasound elastography with continuous transducer vibration |
CN107367722A (zh) * | 2016-05-13 | 2017-11-21 | 中国科学院声学研究所 | 一种减小dpc法累积误差的合成孔径声呐运动补偿方法 |
CN109444901A (zh) * | 2018-11-14 | 2019-03-08 | 杭州电子科技大学 | 一种异构环境下多子阵sas子孔径成像方法 |
US20190129026A1 (en) * | 2015-06-04 | 2019-05-02 | Chikayoshi Sumi | Measurement and imaging instruments and beamforming method |
-
2019
- 2019-06-25 CN CN201910553855.XA patent/CN110244285B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101551452A (zh) * | 2008-04-01 | 2009-10-07 | 中国科学院声学研究所 | 一种合成孔径声纳联合运动补偿方法及系统 |
CN101762810A (zh) * | 2009-12-15 | 2010-06-30 | 中国科学院声学研究所 | 一种大测绘带下合成孔径声纳运动补偿方法 |
CN101799537A (zh) * | 2009-12-15 | 2010-08-11 | 中国科学院声学研究所 | 一种基于双侧回波数据的合成孔径声纳运动补偿方法 |
CN103744081A (zh) * | 2014-01-27 | 2014-04-23 | 中国科学院电子学研究所 | 一种机载圆迹合成孔径雷达高精度三维成像补偿方法 |
US20190129026A1 (en) * | 2015-06-04 | 2019-05-02 | Chikayoshi Sumi | Measurement and imaging instruments and beamforming method |
WO2017062553A1 (en) * | 2015-10-08 | 2017-04-13 | Mayo Foundation For Medical Education And Research | Systems and methods for ultrasound elastography with continuous transducer vibration |
CN107367722A (zh) * | 2016-05-13 | 2017-11-21 | 中国科学院声学研究所 | 一种减小dpc法累积误差的合成孔径声呐运动补偿方法 |
CN109444901A (zh) * | 2018-11-14 | 2019-03-08 | 杭州电子科技大学 | 一种异构环境下多子阵sas子孔径成像方法 |
Non-Patent Citations (6)
Title |
---|
HAYDEN J. CALLOW 等: ""Motion-Compensation Improvement for Widebeam, Multiple-Receiver SAS Systems"", 《IEEE JOURNAL OF OCEANIC ENGINEERING》 * |
MICHAEL P.HAYES 等: ""Synthetic aperture sonar:a review of current status"", 《IEEE JOURNAL OF OCEANIC ENGINEERING》 * |
SANWEN ZHU 等: ""Segment DPCA Motion Compensation for SAS"", 《2009 2ND INTERNATIONAL CONGRESS ON IMAGE AND SIGNAL PROCESSING》 * |
刘纪元: ""合成孔径声呐技术研究进展"", 《中国科学院院刊》 * |
张学波 等: ""多子阵合成孔径声纳系统中的侧摆运动误差补偿"", 《电讯技术》 * |
张鹏飞 等: ""一种适用于大测绘带合成孔径声呐的运动补偿方法"", 《声学学报》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113156444A (zh) * | 2021-06-02 | 2021-07-23 | 杭州电子科技大学 | 一种基于运动补偿的多波束声呐高精度成像方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110244285B (zh) | 2021-02-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108490443B (zh) | 基于解析解及NUFFT的多子阵合成孔径声纳ωk成像算法 | |
US9075145B2 (en) | Combining data from multiple radar signals on a single plan position indicator (PPI) display | |
CN110146889B (zh) | 基于最优回波子区域选择的大转角isar成像方法 | |
CN104898119B (zh) | 一种基于相关函数的动目标参数估计方法 | |
CN109633642B (zh) | 一种太赫兹高速目标雷达成像方法 | |
CN106125078B (zh) | 一种水下多维声学成像系统和方法 | |
CN111856461A (zh) | 基于改进pfa的聚束sar成像方法及其dsp实现 | |
CN109655819B (zh) | 一种基于实孔径多普勒波束锐化的杂波抑制三维成像方法 | |
CN116500626B (zh) | 一种多接收阵元数据收发合置转换方法 | |
CN111090094B (zh) | 脉冲多普勒雷达的双波束角度测量方法、系统及存储介质 | |
CN108845301B (zh) | 一种用于双基地isar的目标等效旋转中心估计方法 | |
US20090099455A1 (en) | Ultrasound diagnosis method and apparatus | |
CN115061130A (zh) | 基于车载非平稳平台的搜跟一体化设计方法及系统 | |
CN110244285A (zh) | 一种合成孔径成像声呐斜距向非均匀分段补偿方法及系统 | |
US9035820B2 (en) | Measurement device, measurement system, measurement method, and program | |
CN107367722B (zh) | 一种减小dpc法累积误差的合成孔径声呐运动补偿方法 | |
CN106908789A (zh) | 一种基于空间域支撑区融合的sar成像方法 | |
CN110308447A (zh) | 合成孔径雷达方位向参数控制方法、装置及存储介质 | |
CN107329128B (zh) | 一种基于快速存储技术的超分辨高精度低空测角方法 | |
CN109613474A (zh) | 一种适用于短距离车载雷达的测角补偿方法 | |
CN112711018B (zh) | 一种近场空间网状结构聚焦波束形成方法 | |
JP2730296B2 (ja) | グランドマッピングレーダ信号処理方法とその装置 | |
JP7030629B2 (ja) | レーダ信号処理装置 | |
CN113075633A (zh) | 一种基于距离-俯仰联合字典的目标定位方法 | |
CN113960599B (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 |