基于椭球腔模型的海域物体定位方法
【技术领域】
[0001] 本发明属于近海生态环境检测技术的领域,特别是指基于椭球腔模型的海域物体 定位方法。
【背景技术】
[0002] 当前全球陆地资源日趋紧张,环境不断恶化,世界各国开始积极开发海洋资源。发 展海洋经济成为沿海国家国民经济的重要支柱,也是可持续发展战略的前沿阵地。然而海 洋开发带来巨大经济效益的同时,也带来一系列资源和生态环境问题。我国的近海生态环 境状况也不理想,造成我国近岸海域污染的主要原因是陆源污染物。海洋垃圾影响海洋景 观,威胁航行安全,甚至会对海洋生态系统的健康产生致命的影响,严重阻碍海洋经济的可 持续发展。因此近海生态环境检测技术的研宄和开发对于我国海洋环境监测和海洋环境保 护事业的健康发展极其必要。近岸海域海洋垃圾电子监测可使我国相关部门和组织掌握近 岸海域海洋垃圾持续存量以及近海生态系统状况,以有效控制、监管和治理近岸海域生态 环境,保障我国海洋经济持续稳定发展。
[0003] 目前,通过超声波检测海洋垃圾研宄在海洋生态环境监测中十分流行。一种是通 过利用无线局域网技术传输数据,采用射频技术传输超声波发送和接收端的同步信号,利 用病态数学和良态数学的理论知道参考点的布置的超声波定位技术在多种局域环境中对 移动物体进行定位。移动物体至参考点距离的误差是影响该技术定位精度的主要因素。另 一种是通过根据超声波的渡越时间(超声波发射器发出的超声波经水质的传播到接收器 的时间)测距以及方向角检查法定位的结合来定位水中被测目标的位置与运动轨迹,该技 术定位精度与被测物体的大小是否合适有着直接的关系。第三种是利用自主实现的超声波 六元传感器阵列进行TDOA测距(利用超声波信号和射频信号获取TDOA所需的时间差),并 采用基于测地距离的多维定标算法(Geodesic Distance MDS)进行无线传感器网络节点定 位技术。该技术的定位精度与网络环境和节点移动有明显的关系,很难应对比较复杂的环 境下的物体定位。
[0004] 在超声波阵列监测过程中,各发射-接收节点对所监测到的物体面,即对应的被 测物体表面反射点,各不相同。以上三种技术都无法很好解决反射点漂移问题,因此它们并 不能良好的应用于海洋物体区域定位中。
【发明内容】
[0005] 本发明的目的在于解决现有该领域的上述的问题,提供一种基于椭球腔模型的海 域物体定位方法,可对海域物体进行准确的立体区域定位。
[0006] 基于椭球腔模型的海域物体定位方法,包括以下步骤:
[0007] 超声波传感器任意一节点发射超声波信号,其它节点接收直接发射的超声波以及 被测物体反射回来的超声波;分别包含有被测物体距发射点最近点和最远点且绕发射点与 任意接收点连线旋转形成的两个椭球面,交集形成"椭球腔";改变接收点取得更多的"椭球 腔","椭球腔"之间进行交集得出质心;所得的质心形成质心群,求出质心群的质心;以质心 群的质心为球心,以距离发射点最近的点到球心的距离为半径形成的球形区域,以该球形 区域作为该物体所在位置以及所在空间的模拟;超声波阵列各节点进行遍历发射,按同样 方法依序获取被测物体所在位置以及所在空间的模拟;融合所有的获取的模拟球形区域, 得到最终的物体检测图。
[0008] 作为优选,以节点I点发射为例,首先取KxpypZi,),J(x2,y 2,z2), K(x3, y3, z3),B(x4, y4, z4)为参考点,0(x, y, z)为待测物体距发射点I (Xp yp Zp )最近点, P(x',y',z')为发射点I (Xl,yi,Zl,)最远可达点,J点接收到I点发射的信号时刻和0 点反射回来的信号时刻分别为&和t/,K点接收到I点发射的信号时刻和0点反射回来 的信号时刻分别为tjP 12',B点接收到I点发射的信号时刻和0点反射回来的信号时刻 分别为t#Pt3';则:
[0009] dioJ= v ? (t /
[0010] diok= v ? (t 2r -t2) +dik (1)
[0011] diob= v (t 3r -t3) +dib
[0012] 则以I,J为焦点,0为椭球上任意一点可得椭球TQIW方程为:
[0014] 同理椭球TQiqk和椭球TQ IQB方程分别为:
[0017] 联立方程(2)、(3)、⑷可得0(x,y,z)坐标;
[0018] 同理,对于点 I (X!,y!,z!,),J(x2, y2, z2),K(x3, y3, z3),B(x4, y4, z4)为参考点,最远 可达点P(x',y',z')为椭球上一点,J点接收到I点发射的信号时刻和P点反射回来 的信号时刻分别为1 4和14',K点接收到I点发射的信号时刻和P点反射回来的信号时刻 分别为tjP 15',B点接收到I点发射的信号时刻和P点反射回来的信号时刻分别为、和 V,则可得反射路径长度分别为:
[0019] dipj=v.(t4' -tJ+dij
[0020] dipk= v ? (t 5' _t5)+dik (5)
[0021] dipb=v.(t6' -t 6) +dib
[0022] 基于三点的椭球TQIW、TQIPK和TQ IPB方程分别为:
[0026] 联立方程(6)、(7)、⑶可得P (x W,^ )坐标;
[0027] 形成以发射点与任意接收点为两焦点,绕两焦点连线旋转的且分别包含有0, P的 内外椭球面,所述内外椭球面交集形成"椭球腔"。
[0028] 作为优选,通过改变接收点取得更多的"椭球腔","椭球腔"之间交集,通过质心算 法求出交集的质心。
[0029] 作为优选,所得的质心形成质心群,通过质心算法得到求出质心群的质心。
[0030] 作为优选,通过改变接收点取得8个类似"椭球腔"区域,"椭球腔"之间进行交集 得出质心。
[0031] 进一步,从中任取3个"椭球腔"的相互交集确定一个封闭的曲面多面体,并通过 质心算法得到交集的质心,总共可得U个交集以及质心。
[0032] 进一步,所得的C丨个质心形成质心群,通过质心算法求出质心群的质心。
[0033] 本发明具有以下有益效果:本发明基于"椭球腔"模型的近岸海域物体立体区域 定位方法,通过超声波阵列和相关融合技术解决了超声波在海洋生态环境监测的应用中, 由于被测物体表面不规则、反射点不相同,无法对海域物体进行准确的立体区域定位的问 题,,能更有效地应用到海洋环境监测中。
【附图说明】
[0034] 图1是基于椭球腔模型的海域物体定位方法的流程图。
【具体实施方式】
[0035] 下面结合具体实施例,并结合附图,对本发明的技术方案作进一步的说明:
[0036] 根据附图1基于椭球腔模型的海域物体定位方法的流程图可知:
[0037] 以超声波传感器各节点按3X3阵列
沿海岸竖直向下安置为例,以海 岸水平线为Z轴(如LIA > )正方向,竖直向下为Y轴正方向(如LMN >),垂直于发 射阵列的法向量方向为X轴正方向。
[0038]以节点I点为发射点,首先I (Xi, y:,z),J(x2, y2, z2),K(x3, y3, z3),B(x4, y4, z4)为 参考点,〇(x,y,z)为待测物体距发射点I(Xl,yi, Zl)最近点,P(x',y',z')为发射点 I(Xl,yi,Zl)最远可达点,J点接收到I点发射的信号时刻和0点反射回来的信号时刻分别为 tJP ti',K点接收到I点发射的信号时刻和0点反射回来的信号时刻分别为t
jP 12',B 点接收到I点发射的信号时刻和〇点反射回来的信号时刻分别为tjP 13';则:
[0039] dioJ= v ? (t /
[0040] diok=v.(t2' -t2) +dik (1)
[0041] diob=v.(t3' -t3) +dib
[0042] 则以I,J为焦点,0为椭球上任意一点可得椭球TQIW方程为:
[0044] 同理椭球TQiqk和椭球TQ IQB方程分别为:
[0047] 联立方程(2)、(3)、⑷可得0(x,y,z)坐标。
[0048] 同理,对于点 I (Xi, y。Zi,),J(x2, y2, z2),K(x3, y3, z3),B(x4, y4, z4)为参考点,最远 可达点P(x',y',z')为椭球上一点,J点接收到I点发射的信号时刻和P点反射回来 的信号时刻分别为14和1 4',K点接收到I点发射的信号时刻和P点反射回来的信号时刻 分别为tjP 15',B点接收到I点发射的信号时刻和P点反射回来的信号时刻分别为、和 V,则可得反射路径长度分别为:
[0049] dipj=v.(t/ -t 4) +(1^-
[0050] dipk= v ? (t 5' _t5)+dik (5)
[0051] dipb=v.(t6' -t 6) +dib
[0052] 基于三点的椭球TQIW、TQIPK和TQ IPB方程分别为:
[0056] 联立方程(6)、(7)、(8)可得 P (x',y',z')坐标。
[0057] 任意焦点相同,与0,P形成的两个椭球即形成一个"椭球腔",如上述方程即可得 到3个"椭球腔"分别为:(TQIW,TQIP;)、(TQiqk,TQIPK)和(TQiqb,TQIPB),则可得三个"椭球腔" 的交集,并求出质心,记为QK Yp Z)。
[0058]依次类推,可得到C83 (56)个质心,分别记为Q (X2, Y2, Z2),……Q (X56, Y56, Z56),通过 平均质心算法即 X〇= (XdX2……X56)/56,Y〇= (Yi+Y2……Y56)/56,Z〇= (Zi+Z2……Z56)/56, 可得质心群质心,即物体模拟区域的球心Q(X(1,%,zQ)。
[0059]根据点I (Xl,yi,Zl,)接收到从物体上反射回来的信号的最早时刻h,可确定距离 发射点最近的物体上的点,并可对其定位,如上述点〇(x, y, z),则以Q、0之间距离 <。为半 径,以、2(|,)为球心的一个球即可得该物体所在位置以及所在空间的模拟。
[0060] 超声波阵列其他各节点,包括A、B、C、J、K、L、M、N进行遍历发射,按同样方法依序 获取被测物体所在位置及所在空间的模拟。
[0061] 融合所有按以上方法获取的模拟球形区域,得到最终的物体检测图。
[0062] 上述【具体实施方式】用来解释说明本发明,而不是对本发明进行限制,在本发明的 精神和权利要求的保护范围内,对本发明做出的任何修改和改变,都落入本发明的保护范 围。
【主权项】
1. 基于椭球腔模型的海域物体定位方法,其特征在于,包括以下步骤: SOl超声波传感器任意一节点发射超声波信号,其它节点接收直接发射的超声波以及 被测物体反射回来的超声波; S02分别包含有被测物体距发射点最近点和最远点且绕发射点与任意接收点连线旋转 形成的两个椭球面,交集形成"椭球腔"; S03改变接收点取得更多的"椭球腔","椭球腔"之间进行交集得出质心; S04所得的质心形成质心群,求出质心群的质心; S05以质心群的质心为球心,以距离发射点最近的点到球心的距离为半径形成的球形 区域,以该球形区域作为该物体所在位置以及所在空间的模拟; S06超声波阵列各节点进行遍历发射,按同样方法依序获取被测物体所在位置以及所 在空间的模拟; S07融合所有的获取的模拟球形区域,得到最终的物体检测图。2. 据权利要求1所述的基于椭球腔模型的海域物体定位方法,其特征在于,步 骤S02通过以下方式实现:以节点I点发射为例,首先取I (X1, Y1, Z1, ),J (x2, y2, Z2), K(x3, y3, z3),B(x4, y4, z4)为参考点,0(x, y, z)为待测物体距发射点I (X1, Y1, Z1,)最近点, P(x<,太,^ )为发射点I U1, yi,Zl,)最远可达点,J点接收到I点发射的信号时刻和〇 点反射回来的信号时刻分别为&和t/,K点接收到I点发射的信号时刻和O点反射回来 的信号时刻分别为tjP 12',B点接收到I点发射的信号时刻和O点反射回来的信号时刻 分别为t#Pt^ ;则: di〇J=v(ti, -t^+dij- diok=v.(t2' _t2)+dik (I) di〇b= V (t 3r -t3) +dib 则以I,J为焦点,〇为椭球上任意一点可得椭球TQiw方程为:同理椭球TQiqk和椭球TQ IQB方程分别为:联立方程(2)、(3)、(4)可得0(x,y,z)坐标; 同理,对于点 I (X1, Y1, Z1,),J(x2, y2, z2),K(x3, y3, z3),B(x4, y4, z4)为参考点,最远可达 点P(x<,太,z^ )为椭球上一点,J点接收到I点发射的信号时刻和P点反射回来的信 号时刻分别为14和14',K点接收到I点发射的信号时刻和P点反射回来的信号时刻分别 为tjP t ^,B点接收到I点发射的信号时刻和P点反射回来的信号时刻分别为tjP t ^, 则可得反射路径长度分别为:基于三点的椭球TQIW、TQipk和TQ IPB方程分别为:联立方程(6)、(7)"8)可得?〇^',7',2')坐标; 形成以发射点与任意接收点为两焦点,绕两焦点连线旋转的且分别包含有〇, P的内外 椭球面,所述内外椭球面交集形成"椭球腔"。3. 据权利要求1所述的基于椭球腔模型的海域物体定位方法,其特征在于:步骤S03 是通过质心算法求出交集的质心。4. 据权利要求1所述的基于椭球腔模型的海域物体定位方法,其特征在于:步骤S04 是通过质心算法求出质心群的质心。5. 据权利要求1所述的基于椭球腔模型的海域物体定位方法,其特征在于:通过改变 接收点取得8个"椭球腔","椭球腔"之间进行交集得出质心。6. 据权利要求5所述的基于椭球腔模型的海域物体定位方法,其特征在于:从所述8 个"椭球腔"中任取3个"椭球腔"交集,并通过质心算法得到交集的质心,总共可得C〗个质 心。7. 据权利要求6所述的基于椭球腔模型的海域物体定位方法,其特征在于:所述的G 个质心形成质心群,通过质心算法求出质心群的质心。
【专利摘要】本发明公开了基于椭球腔模型的海域物体定位方法,包括任一节点发射超声波信号,其它节点接收超声波;得出分别位于绕两个焦点连线旋转形成的两个椭球交集形成“椭球腔”;取得更多“椭球腔”交集得出质心;求出质心群的质心;以此为球心,距离发射点最近的点为半径形成的球形区域,以该球形区域作为该物体所在位置以及所在空间的模拟;各节点遍历发射,融合所有的获取的模拟球形区域。本发明解决了超声波在海洋生态环境监测的应用中,由于被测物体表面不规则、反射点不相同,无法对海域物体进行准确的立体区域定位的问题,更有效地应用到海洋环境监测中。
【IPC分类】G01S15/08, G01S15/46, G01S5/28
【公开号】CN104898092
【申请号】CN201510211941
【发明人】徐斌
【申请人】浙江工商大学
【公开日】2015年9月9日
【申请日】2015年4月29日