CN113608225A - 一种基于合成孔径声纳的水底目标埋深深度计算方法 - Google Patents
一种基于合成孔径声纳的水底目标埋深深度计算方法 Download PDFInfo
- Publication number
- CN113608225A CN113608225A CN202110896058.9A CN202110896058A CN113608225A CN 113608225 A CN113608225 A CN 113608225A CN 202110896058 A CN202110896058 A CN 202110896058A CN 113608225 A CN113608225 A CN 113608225A
- Authority
- CN
- China
- Prior art keywords
- buried
- depth
- target
- towed body
- coordinate
- 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
- 238000000034 method Methods 0.000 title claims abstract description 50
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims description 28
- 238000001514 detection method Methods 0.000 claims abstract description 87
- 238000003384 imaging method Methods 0.000 claims abstract description 21
- 230000008569 process Effects 0.000 claims description 31
- 238000004364 calculation method Methods 0.000 claims description 28
- 238000009933 burial Methods 0.000 claims description 19
- 238000012937 correction Methods 0.000 claims description 11
- 239000000523 sample Substances 0.000 claims description 5
- 230000009466 transformation Effects 0.000 claims description 5
- 230000001131 transforming effect Effects 0.000 claims 1
- 238000010586 diagram Methods 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 244000239635 ulla Species 0.000 description 2
- 101000939348 Homo sapiens NEDD8-activating enzyme E1 regulatory subunit Proteins 0.000 description 1
- 102100029781 NEDD8-activating enzyme E1 regulatory subunit Human genes 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 238000004873 anchoring Methods 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012552 review Methods 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
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
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Mathematical Physics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Algebra (AREA)
- Operations Research (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本申请公开了一种基于合成孔径声纳的水底目标埋深深度计算方法,该方法包括:步骤1:获取至少两次合成孔径声纳在探测区域的探测数据,根据探测数据,确定掩埋目标在成像图中的图像坐标,并根据图像坐标确定探测信息;步骤2,根据探测信息,采用遍历掩埋深度变量的方式,计算使掩埋目标定位坐标的位置偏差最小的掩埋预计深度,将掩埋预计深度记作掩埋目标的掩埋深度。通过本申请中的技术方案,基于双频合成孔径声纳的成像特点和目标定位特点,确定掩埋目标相应的定位信息,以有效计算掩埋目标的掩埋深度信息。
Description
技术领域
本申请涉及目标定位的技术领域,具体而言,涉及一种基于合成孔径声纳的水底目标埋深深度计算方法。
背景技术
合成孔径声纳是一种可以对水底进行实时声成像的装备,可以让人们清晰直观地感知水底大范围场景。目前主要用于水下搜救、水下考古、海底管缆探查、海洋工程作业前期探查等领域。
双频合成孔径声纳是将工作频率在低频(10kHz左右)和高频(100kHz左右)的两套声纳安装在一起,同时对水下成像。由于声波在低频段具有一定的穿透性,能够对处于海底一定深度范围内的目标成像,而高频段的声波穿透能力相对较弱,结合两个频段合成孔径声纳所成图像,就能将只在低频声图中存在,而在高频声图中同一位置没有显现的目标判断为掩埋目标。
但是,在掩埋目标定位的实际应用场景中,不仅要知道掩埋目标的掩埋状态,还需要明确掩埋的深度。比如海底管缆,为防止海底管缆被船只下锚刮断,工程上都会将缆线深埋于海底平面以下,埋缆作业过程先是挖沟,然后在沟中排缆,最后将挖沟挖出来的土回填。为检查管缆埋深是否符合设计要求,需要使用如浅地层剖面仪之类的设备进行掩埋深度复查。但浅地层剖面仪每次只能对一个切片内的目标深度进行判读,效率比较低。
经文献检索发现,有以下文献对目标掩埋深度问题进行了研究:
文献1:《合成孔径成像中掩埋目标的深度误差分析》,该方案中在面阵合成孔径声纳三维成像算法仿真的基础上,利用下视面阵合成孔径声纳在声纳照射范围内有规律地进行多次数据采集,而后进行相干累加获取点目标的高度信息,这种方式不仅与浅地层剖面仪相比工作效率没有明显提升,而且其基于的是仿真的面阵合成孔径声纳,采样位置必须严格指定,而实际海洋作业中这都是难以控制的,导致在实际工程作业中难以实现。
文献2:《掩埋目标探测方法研究》,该方案中提出了一种基于声学参量阵的水下掩埋目标成像方法,针对低频声纳对掩埋物的探测与成像,基于水下掩埋物探测及成像中发射声源特性、水下掩埋物的声散射机理,对掩埋物探测深度进行分析。但是,实际工程作业中,由于水底底质复杂,导致散射回波信号会受到干扰,仍旧无法从接收阵接收到的目标后向散射信号反推出目标掩埋深度。
发明内容
本申请的目的在于:解决双频合成孔径声纳水下掩埋目标探测过程中无法确定掩埋深度的问题。
本申请的技术方案是:提供了一种基于合成孔径声纳的水底目标埋深深度计算方法,该计算方法适用于根据拖体对探测区域的探测数据对掩埋目标深度的计算,拖体由母船拖动驶过探测区域,拖体上安装有合成孔径声纳,该方法包括:步骤1:获取至少两次合成孔径声纳在探测区域的探测数据,根据探测数据,确定掩埋目标在成像图中的图像坐标,并根据图像坐标确定探测信息;步骤2,根据探测信息,采用遍历掩埋深度变量的方式,计算使掩埋目标定位坐标的位置偏差最小的掩埋预计深度,将掩埋预计深度记作掩埋目标的掩埋深度。
上述任一项技术方案中,进一步地,探测信息至少包括:掩埋目标与拖体阵面的斜距、探测到掩埋目标时的拖体UTM坐标和航向角,步骤2具体包括:
步骤21,根据斜距,计算第i次探测时掩埋目标与拖体阵面之间的平距Hd(i),对应的计算公式为:
式中,tgtDpt(i)为第i次探测时对应的掩埋深度变量,Ri为第i次探测时的斜距,tawHgt(i)为第i次探测时对应的拖体距水底高度;
步骤22,根据平距Hd(i)以及对应的航向角θ(i),计算掩埋目标定位坐标,对应的计算公式为:
UlatT(i)=UlatA(i)+δy′(i)
UlonT(i)=UlonA(i)+δx′(i)
δy′(i)=Hd(i)×sin(θ(i))
δx′(i)=Hd(i)×cos(θ(i))
式中,(UlatA(i),UlonA(i))为第i次探测过程的拖体UTM坐标,(δx′(i),δy′(i))为第i次探测过程的掩埋目标UTM坐标修正值,(UlatT(i),UlonT(i))为第i次探测过程的掩埋目标定位坐标;
步骤23,计算掩埋目标定位坐标的位置偏差,采用遍历掩埋深度变量的方式,计算使掩埋目标定位坐标的位置偏差最小的掩埋预计深度,将掩埋预计深度记作掩埋目标的掩埋深度。
上述任一项技术方案中,进一步地,探测次数为2次,掩埋预计深度的计算公式为:
式中,tgtDpt为掩埋深度变量,tgtDpt′为掩埋预计深度。
上述任一项技术方案中,进一步地,步骤23具体包括:步骤A,在设定的掩埋深度范围内,选取第一掩埋深度变量;步骤B,根据第一掩埋深度变量,计算第一次和第二次探测信息中斜距R1和R2对应的掩埋目标与拖体阵面之间的平距Hd(1)和Hd(2);步骤C,根据计算出的平距Hd(1)和Hd(2),结合对应的航向角θ(1)和θ(2),计算掩埋目标定位坐标(UlatT(1),UlonT(1))和(UlatT(2),UlonT(2));步骤D,计算掩埋目标定位坐标(UlatT(1),UlonT(1))和(UlatT(2),UlonT(2))中对应坐标值的坐标差值,并将计算出的坐标差值求和,得到第一掩埋深度变量对应的第一位置偏差;步骤E,在设定的掩埋深度范围内,选取第二掩埋深度变量,重复执行步骤A至步骤D,计算第二掩埋深度变量对应的第二位置偏差,比较第一位置偏差和第二位置偏差之间的数值大小,保留较小值;步骤F,重新选取第三掩埋深度变量执行上述过程,直至位置偏差的取值为0或者得到最小的位置偏差,此时,该位置偏差对应的掩埋深度变量即为掩埋预计深度,将该掩埋预计深度记作掩埋目标的掩埋深度。
上述任一项技术方案中,进一步地,探测信息至少包括:掩埋目标与拖体阵面的斜距、拖体深度、拖体距水底高度、探测到掩埋目标时的拖体UTM坐标和航向角。
上述任一项技术方案中,进一步地,图像坐标包括和图像横坐标和图像纵坐标,步骤1中,根据图像坐标确定探测信息,具体包括:根据图像纵坐标和成像图的分辨率,确定斜距;根据图像横坐标和航行速度,计算拖体航行至此处的通过时间;根据通过时间,查询拖体航行至此处母船的经度、母船的纬度、航向角、拖体深度以及拖体距水底高度;通过坐标变换的方式,将母船的经度和母船的纬度,转换为拖体UTM坐标。
上述任一项技术方案中,进一步地,拖体UTM坐标的计算过程具体包括:对母船的经度和母船的纬度进行坐标变换,计算母船UTM坐标;计算拖体与定位装置沿航行方向的距离,并根据距离和航向角,计算拖体UTM坐标修正量;根据拖体UTM坐标修正量和母船UTM坐标,计算拖体UTM坐标。
本申请的有益效果是:
本申请中的技术方案,根据从不同方向至少两次扫测到同一个掩埋目标的探测信息,通过结合航向角和母船与拖体间的距离运算,采用遍历掩埋深度变量的方式,计算出使掩埋目标定位坐标的位置偏差最小的掩埋预计深度,以作为掩埋目标的埋深。
本申请中的掩埋目标埋深的计算过程也是掩埋目标精确定位的过程。之前合成孔径声纳在对水下掩埋目标定位的计算过程中是不考虑掩埋深度的,主要是因为掩埋深度很难测量,导致多次扫测同一掩埋目标定位位置之间距离相差较大,而本申请正是利用这种定位位置上的差异,通过两次定位位置距离差最小化的方法反过来计算掩埋深度,并同时得到掩埋目标精确的定位位置。
本申请的优点具体还表现在:
1、实际工程作业可操作性强,可使用现有设备,无需改造已有双频合成孔径声纳设备;
2、实现方法简单,节约探测成本,仅需对同一掩埋目标进行两次不同路径的扫测,就能根据两次得到的探测信息,获取掩埋目标的掩埋深度;
3、利用对两次扫测中掩埋目标定位位置坐标差的最小化,获得掩埋目标的掩埋深度,而其中双频合成孔径声纳对水下沉底目标定位的方法已经成熟,本申请是水下目标定位方法针对掩埋目标定位的进一步优化。
在本申请的优选实现方式中,首先根据航向角对拖体UTM坐标进行修订;之后,通过利用掩埋目标与拖体之间的位置关系,计算出掩埋目标与拖体阵面之间的平距;进而结合航向角和拖体坐标,计算出本次探测过程对应的准确的掩埋目标定位坐标,保证了掩埋深度计算的准确性。
附图说明
本申请的上述和/或附加方面的优点在结合下面附图对实施例的描述中将变得明显和容易理解,其中:
图1是根据本申请的一个实施例的水下掩埋目标探测过程的示意图;
图2是根据本申请的一个实施例的基于合成孔径声纳的水底目标埋深深度计算方法的示意流程图;
图3是根据本申请的一个实施例的双次扫测目标过程俯视示意图;
图4是根据本申请的一个实施例的掩埋目标定位过程的位置示意图;
图5是根据本申请的一个实施例的拖体UTM坐标修订的位置示意图;
图6是根据本申请的一个实施例的拖体UTM坐标的修正量计算的示意图;
图7是根据本申请的一个实施例的目标UTM坐标的计算示意图;
图8是根据本申请的一个实施例的仿真实例结果的示意图。
具体实施方式
为了能够更清楚地理解本申请的上述目的、特征和优点,下面结合附图和具体实施方式对本申请进行进一步的详细描述。需要说明的是,在不冲突的情况下,本申请的实施例及实施例中的特征可以相互结合。
在下面的描述中,阐述了很多具体细节以便于充分理解本申请,但是,本申请还可以采用其他不同于在此描述的其他方式来实施,因此,本申请的保护范围并不受下面公开的具体实施例的限制。
如图1所示,本实施例中的水底目标埋深深度计算方法适用于拖曳式双频合成孔径声纳对水下掩埋目标的深度探测,其中,拖曳式双频合成孔径声纳安装在拖体上,拖体通过拖缆连接于母船,以便在母船的拖动下对作业区域内的掩埋目标进行探测,拖曳式双频合成孔径声纳包括安装在拖体两侧的收发装置的基阵,拖体内的发射机和接收机,以及干端的信号处理机和显控机。
需要说明的是,为了利用合成孔径声纳对水下掩埋目标进行定位,在拖体内部一般都安装有测高计和深度计,以便测量拖体的距水底高度和距水面距离(拖体深度),同时,在母船甲板上一般都会放置GPS或者北斗定位装置用于测量母船航速、航向角、经纬度,这些信息都会被信号处理机汇总用于成像和目标定位。
如图2所示,本实施例提供了一种基于合成孔径声纳的水底目标埋深深度计算方法,该方法包括:
步骤1:获取至少两次合成孔径声纳在探测区域的探测数据,根据探测数据,确定掩埋目标在成像图中的图像坐标,并根据图像坐标确定探测信息,其中,探测信息至少包括掩埋目标与拖体阵面的斜距、探测到掩埋目标时的拖体UTM坐标、母船的航向角、拖体深度以及拖体距水底高度,拖体UTM坐标至少包括拖体第一UTM坐标和拖体第二UTM坐标。
具体的,设定x方向是成像图中与航向平行的方向,即拖体航行方向,而y方向是垂直于航向的方向,为距离向。成像图的分辨率为δ。
如图3所示,当母船拖动拖体从不同方向至少两次扫测同一掩埋目标时,根据拖体上合成孔径声纳在探测区域的探测数据,通过统计合成孔径声纳扫测过程中所成的成像区中掩埋目标处x方向和y方向像素点的个数,或者读取双频合成孔径声纳显控中的声图坐标,确定掩埋目标在成像图中xy方向的图像坐标(x0,y0)。
如图4所示,由于阵面法线方向与水平面呈一定倾角,因此,将掩埋目标与拖体阵面的距离记作斜距,该值表示目标距离拖体阵面的距离,由掩埋目标在成像图中的图像纵坐标y0和成像图的分辨率δ确定,对应的计算公式为:
Ri=y0i×δ
式中,δ为分辨率,Ri为斜距,i=1,2,i为探测次数的编号,即R1为第一次探测过程中的斜距,R2为第二次探测过程中的斜距,y0i为第i次探测时的图像纵坐标。
需要说明的是,探测次数的编号与获取合成孔径声纳在探测区域的探测数据的次数相同。
由于图像横坐标x0表示的是掩埋目标出现在拖体阵面正横方向时拖体的位置,设定拖体是处于匀速直线运动的状态,已知初始UTC时间t0,结合拖体航行速度v,因此,通过对应的第i次探测时的图像横坐标x0i可以计算出拖体航行至此处的通过时间t,对应的计算公式为:
再利用该通过时间t,通过查询的方式,可以确定该时刻母船的经度lon、纬度lat、航向角θ、拖体深度tawDpt和拖体距水底高度tawHgt。
进而通过坐标转换的方式,将该时刻母船的经度lon、纬度lat转换至水下拖体的UTM(Universal Transverse Mercatar Projection)坐标,以获取扫过同一掩埋目标时拖体合成孔径声纳的经纬度数据。
需要说明的是,在精度符合要求的前提下,可以直接将母船经纬度,通过坐标变换的方式,转换为UTM坐标系下的结果直接作为拖体UTM坐标。
本实施例中,对UTM坐标转换的方式并不限定。
进一步的,该步骤1中还包括:根据母船的航向角,对拖体UTM坐标进行修订。
具体的,如图5所示,对于任一次探测而言,设定坐标修订过程中能够确定以下参数:
1)定位设备接收天线与缆线导缆轮出缆口之间的水平距离dgps;
2)导缆轮离水面高度wheelHgt;
3)从拖体的深度计获取的拖体深度值tawDpt,即拖体距水底高度;
4)从绞车处读取的缆长值cableLen。
结合图5中的三角函数关系,可以计算得到拖体阵面(拖体)至定位装置天线沿航行方向的距离distGPS2taw为:
设定母船的经度lon、纬度lat转换后的UTM坐标系中的母船的UTM坐标为(UlatB,UlonB),如图6所示,设定本实施例中的航向角为航行方向与UTM坐标系中的y轴方向(正北方向)之间的夹角。
根据第i次探测时对应的航向角θ(i)以及计算出的拖体阵面至定位装置天线沿航行方向的距离distGPS2taw,计算拖体UTM坐标修正量,对应的计算公式为:
δx(i)=-distGPS2taw×sin(θ(i))
δy(i)=-distGPS2taw×cos(θ(i))
需要说明的是,由于拖体位于母船的后方,所以拖体UTM坐标修正量前取负号。
式中,δx为UTM横坐标修正量,δy为UTM纵坐标修正量。
之后根据母船的UTM坐标为(UlatB(i),UlonB(i)),计算修正后的拖体UTM坐标,对应的计算公式为:
UlatA(i)=UlatB(i)+δy(i)
UlonA(i)=UlonB(i)+δx(i)
式中,UlatA(i)为第i次探测时对应的拖体UTM坐标中的横坐标,UlonA(i)为第i次探测时对应的拖体UTM坐标中的纵坐标。
本实施例中,设定第一次探测扫过掩埋目标时得到拖体UTM坐标,即拖体第一UTM坐标为(ULatA1,ULonA1),第二次探测扫过掩埋目标时得到拖体UTM坐标,即拖体第二UTM坐标为(UlatA2,UlonA2)。
步骤2,根据探测信息,采用遍历掩埋深度变量的方式,计算使掩埋目标定位坐标的位置偏差最小的掩埋预计深度,将掩埋预计深度记作掩埋目标的掩埋深度,同时获取掩埋目标的准确定位。
如图4所示,在阵面位置1处,通过上述过程,可计算出掩埋目标距离阵面的斜距为R1,如果不考虑掩埋目标的掩埋深度,即认为其是个沉底目标,则推测出其位置将处于“虚目标1”的位置处。同样的,对于阵面位置2,计算出的斜距为R2,如果认为其是沉底目标,则推测出其位置会在“虚目标2”处。
因此,在不考虑掩埋深度的情况下,拖体两次从不同方向下对掩埋目标的扫测,同一个掩埋目标将会出现在两个位置相距较远的虚目标处,即两个阵面对掩埋目标定位坐标之间的位置偏差较大。
通过采用遍历掩埋深度变量的方式,改变掩埋目标定位坐标推算过程中的掩埋深度变量,最小化两次掩埋目标定位坐标之间的位置偏差,即可得到最优化的掩埋目标的掩埋深度。
设定相同的掩埋深度变量tgtDpt下,两次扫测到同一个掩埋目标后,得到的两组对应的掩埋目标UTM坐标,作为掩埋目标定位坐标,即掩埋目标定位坐标分别为:(UlatT1,UlonT1)和(UlatT2,UlonT2),因此,掩埋预计深度的计算公式为:
式中,tgtDpt为掩埋深度变量,tgtDpt′为掩埋预计深度。
本实施例以任一组掩埋目标定位坐标为例,对掩埋目标定位坐标的计算过程进行说明。
在计算掩埋目标定位坐标过程中需要用到的参数包括:航向角θ、掩埋目标与拖体阵面之间的斜距Ri、拖体距水底高度tawHgt。
在计算过程中,设定掩埋目标所在区域是个大平层,即拖体正下方的高程与掩埋目标所在高程相同。
因此,对于任一次探测而言,该步骤2中具体包括:
步骤21,在距离向上,根据所述斜距,计算第i次探测时掩埋目标与拖体阵面之间的平距Hd(i),对应的计算公式为:
式中,tgtDpt(i)为第i次探测时对应的掩埋深度变量,Ri为第i次探测时的斜距,tawHgt(i)为第i次探测时对应的拖体距水底高度,i为探测次数,i=1,2,…。
步骤22,根据所述平距Hd(i)以及对应的航向角θ(i),计算所述掩埋目标定位坐标。
如图7所示,在通过上述过程得到的第i次探测过程的拖体UTM坐标(UlatA(i),UlonA(i))后,基于三角函数原理,根据平距Hd(i)和对应的航向角θ(i),即可计算掩埋目标定位坐标(UlatT(i),UlonT(i)),对应的计算公式为:
UlatT(i)=UlatA(i)+δy′(i)
UlonT(i)=UlonA(i)+δx′(i)
δy′(i)=Hd(i)×sin(θ(i))
δx′(i)=Hd(i)×cos(θ(i))
式中,δy′(i)为掩埋目标UTM横坐标修订值,δx′(i)为掩埋目标UTM纵坐标修订值。
因此,通过上述过程可计算出两次扫测对应的掩埋目标定位坐标:(UlatT1,UlonT1)和(UlatT2,UlonT2)。
步骤23,计算所述掩埋目标定位坐标的位置偏差,采用遍历掩埋深度变量的方式,计算使所述掩埋目标定位坐标的位置偏差最小的掩埋预计深度,将所述掩埋预计深度记作所述掩埋目标的掩埋深度,其中,掩埋深度的计算公式:
采用遍历掩埋深度变量tgtDpt的方式,计算使掩埋目标定位坐标的位置偏差最小的掩埋预计深度,即掩埋目标的掩埋深度。
优选的,对第i次探测得到的探测信息,采用遍历掩埋深度变量tgtDpt的方式,计算使掩埋目标定位坐标的位置偏差最小的掩埋预计深度,具体包括:
步骤A,在设定的掩埋深度范围内,选取第一掩埋深度变量;
步骤B,根据所述第一掩埋深度变量,计算第一次和第二次探测信息中斜距R1和R2对应的掩埋目标与拖体阵面之间的平距Hd(1)和Hd(2);
步骤C,根据计算出的平距Hd(1)和Hd(2),结合对应的航向角θ(1)和θ(2),计算掩埋目标定位坐标(UlatT(1),UlonT(1))和(UlatT(2),UlonT(2));
步骤D,计算掩埋目标定位坐标(UlatT(1),UlonT(1))和(UlatT(2),UlonT(2))中对应坐标值的坐标差值,并将计算出的坐标差值求和,得到第一掩埋深度变量对应的第一位置偏差;
步骤E,在设定的掩埋深度范围内,选取第二掩埋深度变量,重复执行步骤A至步骤D,计算第二掩埋深度变量对应的第二位置偏差,比较第一位置偏差和第二位置偏差之间的数值大小,保留较小值;
步骤F,重新选取第三掩埋深度变量执行上述过程,直至位置偏差的取值为0或者得到最小的位置偏差,此时,该位置偏差对应的掩埋深度变量即为掩埋预计深度,将该掩埋预计深度记作掩埋目标的掩埋深度。
需要说明的是,掩埋深度变量的选取可以是先从掩埋深度范围内的最大值开始选取,即第一掩埋深度变量的取值小于第二掩埋深度变量的取值,第二掩埋深度变量的取值小于第三掩埋深度变量的取值,依次类推。也可以是先从掩埋深度范围内的最小值开始选取,即第一掩埋深度变量的取值大于第二掩埋深度变量的取值,第二掩埋深度变量的取值大于第三掩埋深度变量的取值,依次类推。还可以是随机选取。
根据得到该掩埋预计深度对应的掩埋目标定位坐标和掩埋深度,即可对掩埋目标进行准确定位。
为了验证本实施例中的掩埋深度计算方法,设置验证条件为:第1次扫测到掩埋目标时拖体的坐标为(400,300),第2次扫测到掩埋目标时拖体的坐标为(401,301),假设掩埋目标实际位置在(500,550)处,掩埋深度5m,拖体距水底高度为50m。
由上述假设条件可以计算出斜距R1为274.8m,斜距R2为273.5m。
从0至10m遍历掩埋深度变量tgtDpt,可以得到两组掩埋目标UTM坐标之间的定位坐标差曲线,如图8所示,从图8中可以看出通过上述过程计算出的掩埋深度最优值,即使得两次扫测对掩埋目标定位距离差最小的掩埋深度5m,与设定的掩埋深度一致。
以上结合附图详细说明了本申请的技术方案,本申请提出了一种基于合成孔径声纳的水底目标埋深深度计算方法,该方法包括:步骤1:获取至少两次合成孔径声纳在探测区域的探测数据,根据探测数据,确定掩埋目标在成像图中的图像坐标,并根据图像坐标确定探测信息;步骤2,根据探测信息,采用遍历掩埋深度变量的方式,计算使掩埋目标定位坐标的位置偏差最小的掩埋预计深度,将掩埋预计深度记作掩埋目标的掩埋深度。通过本申请中的技术方案,基于双频合成孔径声纳的成像特点和目标定位特点,确定掩埋目标相应的定位信息,以有效计算目标的掩埋深度信息。
本申请中的步骤可根据实际需求进行顺序调整、合并和删减。
本申请装置中的单元可根据实际需求进行合并、划分和删减。
尽管参考附图详地公开了本申请,但应理解的是,这些描述仅仅是示例性的,并非用来限制本申请的应用。本申请的保护范围由附加权利要求限定,并可包括在不脱离本申请保护范围和精神的情况下针对发明所作的各种变型、改型及等效方案。
Claims (7)
1.一种基于合成孔径声纳的水底目标埋深深度计算方法,所述计算方法适用于根据拖体对探测区域的探测数据对掩埋目标深度的计算,所述拖体由母船拖动驶过所述探测区域,所述拖体上安装有合成孔径声纳,其特征在于,所述方法包括:
步骤1:获取至少两次所述合成孔径声纳在所述探测区域的所述探测数据,根据所述探测数据,确定所述掩埋目标在成像图中的图像坐标,并根据所述图像坐标确定探测信息;
步骤2,根据所述探测信息,采用遍历掩埋深度变量的方式,计算使所述掩埋目标定位坐标的位置偏差最小的掩埋预计深度,将所述掩埋预计深度记作所述掩埋目标的掩埋深度。
2.如权利要求1所述的基于合成孔径声纳的水底目标埋深深度计算方法,其特征在于,所述探测信息至少包括:所述掩埋目标与拖体阵面的斜距、探测到所述掩埋目标时的拖体UTM坐标和航向角,所述步骤2具体包括:
步骤21,根据所述斜距,计算第i次探测时掩埋目标与拖体阵面之间的平距Hd(i),对应的计算公式为:
式中,tgtDpt(i)为第i次探测时对应的掩埋深度变量,Ri为第i次探测时的斜距,tawHgt(i)为第i次探测时对应的拖体距水底高度;
步骤22,根据所述平距Hd(i)以及对应的航向角θ(i),计算所述掩埋目标定位坐标,对应的计算公式为:
UlatT(i)=UlatA(i)+δy′(i)
UlonT(i)=UlonA(i)+δx′(i)
δy′(i)=Hd(i)×sin(θ(i))
δx′(i)=Hd(i)×cos(θ(i))
式中,(UlatA(i),UlonA(i))为第i次探测过程的拖体UTM坐标,(δx′(i),δy′(i))为第i次探测过程的掩埋目标UTM坐标修正值,(UlatT(i),UlonT(i))为第i次探测过程的掩埋目标定位坐标;
步骤23,计算所述掩埋目标定位坐标的位置偏差,采用遍历掩埋深度变量的方式,计算使所述掩埋目标定位坐标的位置偏差最小的掩埋预计深度,将所述掩埋预计深度记作所述掩埋目标的掩埋深度。
4.如权利要求2所述的基于合成孔径声纳的水底目标埋深深度计算方法,其特征在于,所述步骤23具体包括:
步骤A,在所述掩埋深度范围内选取第一掩埋深度变量;
步骤B,根据所述第一掩埋深度变量,计算第一次和第二次探测信息中斜距R1和R2对应的掩埋目标与拖体阵面之间的平距Hd(1)和Hd(2);
步骤C,根据计算出的所述平距Hd(1)和Hd(2),结合对应的航向角θ(1)和θ(2),计算所述掩埋目标定位坐标(UlatT(1),UlonT(1))和(UlatT(2),UlonT(2));
步骤D,计算所述掩埋目标定位坐标(UlatT(1),UlonT(1))和(UlatT(2),UlonT(2))中对应坐标值的坐标差值,并将计算出的坐标差值求和,得到第一掩埋深度变量对应的第一位置偏差;
步骤E,在所述掩埋深度范围内选取第二掩埋深度变量,重复执行步骤A至步骤D,计算所述第二掩埋深度变量对应的第二位置偏差,比较所述第一位置偏差和所述第二位置偏差之间的数值大小,保留较小值;
步骤F,重新选取第三掩埋深度变量执行上述过程,直至位置偏差的取值为0或者得到最小的位置偏差,将所述位置偏差对应的掩埋深度变量记作所述掩埋预计深度。
5.如权利要求1所述的基于合成孔径声纳的水底目标埋深深度计算方法,其特征在于,所述探测信息至少包括:所述掩埋目标与拖体阵面的斜距、拖体深度、拖体距水底高度、探测到所述掩埋目标时的拖体UTM坐标和航向角。
6.如权利要求5所述的基于合成孔径声纳的水底目标埋深深度计算方法,其特征在于,所述图像坐标包括和图像横坐标x0图像纵坐标y0,所述步骤1中,根据所述图像坐标确定探测信息,具体包括:
根据所述图像纵坐标y0和所述成像图的分辨率δ,确定所述斜距;
根据所述图像横坐标x0和航行速度,计算所述拖体航行至此处的通过时间;
根据所述通过时间,查询所述拖体航行至此处母船的经度lon、母船的纬度lat、所述航向角、所述拖体深度以及所述拖体距水底高度;
通过坐标变换的方式,将所述母船的经度lon和所述母船的纬度lat,转换为所述拖体UTM坐标。
7.如权利要求6所述的基于合成孔径声纳的水底目标埋深深度计算方法,其特征在于,所述拖体UTM坐标的计算过程具体包括:
对所述母船的经度lon和所述母船的纬度lat进行坐标变换,计算母船UTM坐标;
计算拖体与定位装置沿航行方向的距离,并根据所述距离和所述航向角,计算拖体UTM坐标修正量;
根据所述拖体UTM坐标修正量和所述母船UTM坐标,计算所述拖体UTM坐标。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110896058.9A CN113608225B (zh) | 2021-08-05 | 2021-08-05 | 一种基于合成孔径声纳的水底目标埋深深度计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110896058.9A CN113608225B (zh) | 2021-08-05 | 2021-08-05 | 一种基于合成孔径声纳的水底目标埋深深度计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113608225A true CN113608225A (zh) | 2021-11-05 |
CN113608225B CN113608225B (zh) | 2024-02-13 |
Family
ID=78307044
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110896058.9A Active CN113608225B (zh) | 2021-08-05 | 2021-08-05 | 一种基于合成孔径声纳的水底目标埋深深度计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113608225B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114200458A (zh) * | 2021-12-06 | 2022-03-18 | 广东智能无人系统研究院 | 一种基于自主水下机器人的水下掩埋目标探查方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN206531959U (zh) * | 2017-01-13 | 2017-09-29 | 苏州桑泰海洋仪器研发有限责任公司 | 声纳系统信号采集传输控制装置 |
CN107678034A (zh) * | 2017-11-16 | 2018-02-09 | 中科探海(苏州)海洋科技有限责任公司 | 一种掩埋目标高效三维探测声呐 |
CN110794412A (zh) * | 2019-11-05 | 2020-02-14 | 苏州桑泰海洋仪器研发有限责任公司 | 一种轻便型掩埋管缆路由探测系统 |
CN110907937A (zh) * | 2018-09-18 | 2020-03-24 | 中国科学院声学研究所 | 一种基于“t”型阵的掩埋物合成孔径三维成像方法 |
CN110907938A (zh) * | 2018-09-18 | 2020-03-24 | 中国科学院声学研究所 | 一种近场的快速下视合成孔径三维成像方法 |
CN111190185A (zh) * | 2020-04-08 | 2020-05-22 | 汕头大学 | 一种基于群体智能的水下掩埋雷自主探测方法及系统 |
-
2021
- 2021-08-05 CN CN202110896058.9A patent/CN113608225B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN206531959U (zh) * | 2017-01-13 | 2017-09-29 | 苏州桑泰海洋仪器研发有限责任公司 | 声纳系统信号采集传输控制装置 |
CN107678034A (zh) * | 2017-11-16 | 2018-02-09 | 中科探海(苏州)海洋科技有限责任公司 | 一种掩埋目标高效三维探测声呐 |
CN110907937A (zh) * | 2018-09-18 | 2020-03-24 | 中国科学院声学研究所 | 一种基于“t”型阵的掩埋物合成孔径三维成像方法 |
CN110907938A (zh) * | 2018-09-18 | 2020-03-24 | 中国科学院声学研究所 | 一种近场的快速下视合成孔径三维成像方法 |
CN110794412A (zh) * | 2019-11-05 | 2020-02-14 | 苏州桑泰海洋仪器研发有限责任公司 | 一种轻便型掩埋管缆路由探测系统 |
CN111190185A (zh) * | 2020-04-08 | 2020-05-22 | 汕头大学 | 一种基于群体智能的水下掩埋雷自主探测方法及系统 |
Non-Patent Citations (1)
Title |
---|
陈强;刘国祥;李永树;丁晓利;: "干涉雷达永久散射体自动探测――算法与实验结果", 测绘学报, no. 02 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114200458A (zh) * | 2021-12-06 | 2022-03-18 | 广东智能无人系统研究院 | 一种基于自主水下机器人的水下掩埋目标探查方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113608225B (zh) | 2024-02-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112505750B (zh) | 一种深拖多道地震拖缆姿态确定方法及处理终端 | |
Guan et al. | An effective method for submarine buried pipeline detection via multi-sensor data fusion | |
Singh et al. | Microbathymetric mapping from underwater vehicles in the deep ocean | |
CN104133217B (zh) | 一种水下运动目标与水流的三维速度联合测定方法及装置 | |
JP2003114275A (ja) | バイスタティックレーダシステム | |
CN113640808B (zh) | 浅水区海底电缆埋深探测方法及装置 | |
Xin et al. | A TOA/AOA underwater acoustic positioning system based on the equivalent sound speed | |
RU2012153734A (ru) | Способ съемки рельефа дна акватории и устройство для съемки рельефа дна акватории | |
Godin | The calibration of shallow water multibeam echo-sounding systems | |
JP3515751B2 (ja) | 海底三次元構造の再構成方法 | |
Okino et al. | Measurement of seabed topography by multibeam sonar using CFFT | |
CN113608225B (zh) | 一种基于合成孔径声纳的水底目标埋深深度计算方法 | |
US20180003809A1 (en) | System and method for calibration of echo sounding systems and improved seafloor imaging using such systems | |
CN111880185A (zh) | 一种水下目标勘测处理方法及系统 | |
CN112147578B (zh) | 一种高精度深水发射阵及多元垂直接收阵阵元定位系统与方法 | |
CN112904428B (zh) | 海洋浅地层剖面探测系统及方法 | |
KR20040092508A (ko) | Gps와 gpr을 이용한 하상 지반 조사 시스템 | |
CN115032695B (zh) | 一种基于声学侧反射的海底管缆追踪式探测方法 | |
RU2529207C1 (ru) | Система навигации буксируемого подводного аппарата | |
RU2480790C1 (ru) | Способ определения местоположения измеренных глубин звуковыми сигналами | |
CN114018224A (zh) | 一种海图水深数据检核系统和方法 | |
JPH01216288A (ja) | 海底地形調査方法 | |
CN203358842U (zh) | 海底沉船姿态测量系统 | |
RU2736231C1 (ru) | Способ определения распределения скорости звука | |
RU2771434C1 (ru) | Устройство для съемки подводной поверхности айсберга |
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 |