CN106595599A - 用于尾矿库干滩长度与滩顶高程的自动巡航算法 - Google Patents

用于尾矿库干滩长度与滩顶高程的自动巡航算法 Download PDF

Info

Publication number
CN106595599A
CN106595599A CN201611193413.1A CN201611193413A CN106595599A CN 106595599 A CN106595599 A CN 106595599A CN 201611193413 A CN201611193413 A CN 201611193413A CN 106595599 A CN106595599 A CN 106595599A
Authority
CN
China
Prior art keywords
beach
dry
cruise
point
dlt
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.)
Pending
Application number
CN201611193413.1A
Other languages
English (en)
Inventor
王建国
张�诚
李磊
高美娟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xian University of Technology
Xian Technological University
Original Assignee
Xian Technological University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Xian Technological University filed Critical Xian Technological University
Priority to CN201611193413.1A priority Critical patent/CN106595599A/zh
Publication of CN106595599A publication Critical patent/CN106595599A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C11/00Photogrammetry or videogrammetry, e.g. stereogrammetry; Photographic surveying
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C3/00Measuring distances in line of sight; Optical rangefinders
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C5/00Measuring height; Measuring distances transverse to line of sight; Levelling between separated points; Surveyors' levels

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Multimedia (AREA)
  • Electromagnetism (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

用于尾矿库干滩长度与滩顶高程的自动巡航算法,包括以下步骤:步骤1、将干滩测量仪置于尾矿库附近的高点,测量并记录干滩测量仪的高程;步骤2、预设干滩长度巡航轨迹和滩顶高程巡航轨迹,根据干滩长度巡航轨迹寻找干滩长度测量点,根据滩顶高程巡航轨迹寻找滩顶高程测量点;步骤3、记录干滩测量仪与干滩长度测量点的距离、水平巡航角、垂直巡航角;在坝体上选择两个参考点,记录干滩测量仪分别与两个坝体参考点的距离、水平巡航角、垂直巡航角,计算干滩长度;步骤4、记录干滩测量仪与滩顶高程测量点的距离、水平巡航角、垂直巡航角,计算滩顶高程。解决了现有技术中存在的人工测量受环境制约准确率低的问题。

Description

用于尾矿库干滩长度与滩顶高程的自动巡航算法
技术领域
本发明属于矿冶安全监控方法技术领域,涉及一种用于尾矿库干滩长度与滩顶高程的自动巡航算法。
背景技术
尾矿干排,是近年来国内逐渐兴起的一项尾矿处置工艺,是指选矿流程输出的尾矿浆经多级浓缩后,再经脱水振动筛等高效脱水设备进行脱水处理,形成含水小、易沉淀固化、运输方便易于堆存的矿渣。尾矿干排能够提高尾矿利用率,减少环境污染,现在国家已经明文要求,必须对矿浆进行尾矿干排处理。
尾矿库是指尾矿坝拦截谷口或围地构成的,用以堆存尾矿干排过程中形成的矿渣的场所。在尾矿干排过程中,筑坝拦截矿浆中的矿渣和水,矿渣由于水的冲击在筑坝拦截面沉积,沉积体露出水面的部分称为干滩。干滩长度是指水滩分界点到坝顶方向线的最小垂直距离,滩顶高程是指干滩中的最高点的高程,往往会随着放矿、降雨等因素发生变化。
尾矿库是一个具有高势能的人造泥石流危险源,存在溃坝危险,一旦失事,容易造成重特大事故,干滩长度和滩顶高程是检测尾矿安全的重要指标,传统的测量方法是通过人工定期或不定期采集数据,实现对尾矿库的安全监测,但是,人工采集不仅会存在一些误差和错误,而且受到环境、地形、天气、人员素质等因素的制约,会造成一定的安全隐患。
发明内容
本发明的目的是提供一种用于尾矿库干滩长度与滩顶高程的自动巡航算法,解决了现有技术中存在的人工测量受环境制约准确率低的问题。
本发明所采用的技术方案是,用于尾矿库干滩长度与滩顶高程的自动巡航算法,包括以下步骤:
步骤1、将干滩测量仪置于尾矿库附近的高点,测量并记录干滩测量仪的高程;
步骤2、预设干滩长度巡航轨迹和滩顶高程巡航轨迹,根据干滩长度巡航轨迹寻找干滩长度测量点,根据滩顶高程巡航轨迹寻找滩顶高程测量点;
步骤3、记录干滩测量仪到干滩长度测量点的距离,以及干滩测量仪与干滩长度测量点的水平巡航角、垂直巡航角;在坝体上选择两个参考点,记录干滩测量仪分别到两个坝体参考点的距离,以及干滩测量仪分别与两个坝体参考点的水平巡航角、垂直巡航角,计算干滩长度;
步骤4、记录干滩测量仪到滩顶高程测量点的距离,以及干滩测量仪与滩顶高程测量点的水平巡航角、垂直巡航角,计算滩顶高程。
干滩测量仪包括图像处理模块,图像处理模块分别连接有激光测距仪和云台,云台上设置有摄像机。
步骤2中,预设多条干滩长度巡航轨迹,干滩长度巡航轨迹以水面处为起点,以干滩处为终点,并经过水滩分界线;预设一条滩顶高程巡航轨迹,滩顶高程巡航轨迹经过干滩最高点。
步骤2中,寻找干滩长度测量点和滩顶高程测量点的过程中,巡航操作具体如下:
每条巡航轨迹的起点为S,终点为E,在巡航轨迹起点S处,水平巡航 角为αS,垂直巡航角为βS;在巡航轨迹终点E处,水平巡航角为αE,垂直巡航角为βE;巡航过程中,水平巡航角为α,垂直巡航角为β:
对水平巡航角α进行变换:
计算dltα=(αE-αS) (1)
(1)当dltα>0且|dltα|>180°,Kα=360°-|dltα|;在αS的基础上每次减少0.1°(若αS≤0°时,αS=αS+360°);
(2)当dltα>0且|dltα|<180°时,Kα=|dltα|;如果dltα<0,在αS的基础上每次增加0.1°;
(3)当dltα<0且|dltα|>180°,Kα=360°-|dltα|;在αS的基础上每次增加0.1°(若αS≥360°时,αS=αS-360°);
(4)当dltα<0且|dltα|<180°时,Kα=|dltα|;如果dltα<0,在αS的基础上每次减少0.1°;
对垂直巡航角β进行变换:
计算dltβ=(βE-βS) (2)
(1)当dltβ>0且|dltβ|>180°,Kβ=360°-|dltβ|,K=Kβ/Kα;在βS的基础上每次减少K*0.1°(若βS≤0°时,βS=βS+360°);
(2)当dltβ>0且|dltβ|<180°时,Kβ=|dltβ|,K=Kβ/Kα;如果dltβ<0,在βS的基础上每次增加K*0.1°;
(3)当dltβ<0且|dltβ|>180°,Kβ=360°-|dltβ|,K=Kβ/Kα;在βS的基础上每次增加K*0.1°(若βS≥360°时,βS=βS-360°);
(4)当dltβ<0且|dltβ|<180°时,Kβ=|dltβ|,K=Kβ/Kα;如果dltβ<0,在βS的基础上每次减少K*0.1°。
寻找干滩长度测量点的过程中,当激光测量仪照射到水面时,无距离返 回,当激光测量仪照射到干滩时,有距离返回,水滩分界点,即为干滩长度测量点;
沿一条干滩长度巡航轨迹寻找水滩分界点的过程中,转动云台,每次对水平巡航角α与垂直巡航角β进行变换,用激光测距仪测量距离,若连续三次转云台后,都测量出了距离,即认为找到了水滩分界点W,此时停止对该条干滩长度巡航轨迹的巡航,记录滩测量仪与干滩长度测量点W的距离SW、水平巡航角αW、垂直巡航角βW
当巡航角度超过了终点时,即超出了干滩长度巡航轨迹,认为该条干滩长度轨迹未能测出干滩长度。
计算每条干滩长度巡航轨迹测得的干滩长度的步骤为:
(1)计算水滩分界点在以干滩测量仪W点为原点的坐标系下的坐标为(XW,YW,ZW):
Xw=SW cosβW sinαW (3)
Yw=SW cosβW cosαW (4)
Zw=SW sinβW (5)
其中,SW、αW、βW干分别为滩测量仪与干滩长度测量点W的距离、滩测量仪与干滩长度测量点W的水平巡航角、滩测量仪与干滩长度测量点W的垂直巡航角;
计算坝体参考点P1与坝面参考点P2在W坐标系下的坐标分别为(X1w,Y1w,Z1w)与(X2w,Y2w,Z2w):
X1w=S1cosβ1sinα1 (6)
Y1w=S1cosβ1cosα1 (7)
Z1w=S1sinβ1 (8)
X2w=S2cosβ2sinα2 (9)
Y2w=S2cosβ2cosα2 (10)
Z2w=S2sinβ2 (11)
其中,(X1,Y1,Z1)为坝体参考点P1的坐标,(X2,Y2,Z2)为坝体参考点P2的坐标,S1、α1、β1分别为干滩测量仪与坝体参考点P1的距离、干滩测量仪与坝体参考点P1的水平巡航角、干滩测量仪与坝体参考点P1的垂直巡航角,S2、α2、β2分别为干滩测量仪与坝体参考点P2的距离、干滩测量仪与坝体参考点P2的水平巡航角、干滩测量仪与坝体参考点P2的垂直巡航角;
(2)将以水滩分界点为原点的W坐标系平移至坝体参考点P1为原点的M坐标系,得到水滩分界点在M坐标系下的坐标为(Xm,Ym,Zm):
即:
Xm=Xw-X1w (13)
Ym=Yw-Y1w (14)
Zm=Zw-Z1w (15)
得到坝体参考点P2在M坐标系下的坐标(X2m,Y2m,Z2m):
即:
X2m=X2w-X1w (17)
Y2m=Y2w-Y1w (18)
Z2m=Z2w-Z1w (19)
(3)将M坐标系旋转,使Y坐标轴与两个坝体参考点P1、P2的连线重合,计算旋转的角度θ:
θ=tan-1(Y2m/X2m) (20)
(4)计算干滩长度:
其中,Xc、Yc、Zc分别是以P1、P2所在直线为Y轴的坐标系的X、Y、Z轴,干滩长度即为干滩长度测量点(水滩分界点)距离Yc的最小距离,即,
Xc=Xm cosθ+Ym sinθ (22)
对每条干滩长度巡航轨迹算出的干滩长度Xc取均值,得出干滩长度d。
在寻找滩顶高程的测量点的过程中,通过转动云台,增加或减少水平巡航角α、垂直巡航角β,进行测量;
在测量的距离连续出现之后,出现了测量值的突变(例如某一次的测量值与上一次测量值的差大于10),即认为滩顶高程测量点为突变前的点G;
根据干滩测量仪的高程、干滩测量仪与测点的垂直巡航角βG,运用勾股定理计算出滩顶高程测量点与干滩测量仪的相对高程,滩顶高程测量点与干滩测量仪的相对高程加上干滩测量仪的高程,即为滩顶高程。
测量干滩测量仪的高程的方法如下:选择一个高程已知点,通过摄像机找到该点,根据激光测距仪测出该点与干滩测量仪的距离,以及该点与干滩 测量仪的垂直夹角,根据勾股定理计算出干滩测量仪相对于高程已知点的高程;干滩测量仪相对于测点的高程加上高程已知点的高程,即为干滩测量仪的高程。
本发明的有益效果是:
(1)将干滩测量仪置于尾矿库附近的高点(如附近的山体上),避免了干滩因放矿、降雨等因素变化引起干滩测量仪置位置的变化,实现了全天候采集;
(2)先预设干滩长度巡航轨迹和滩顶高程巡航轨迹,干滩测量仪通过巡航轨迹寻找干滩长度测量点和滩顶高程测量点,克服了人工采集时,受天气、环境的制约,以及测量人员自身因素产生的错误和误差,测量人员足不出户即可对干滩长度和滩顶高程进行测量;
(3)计算方法简便,准确率高,根据干滩长度测量点和滩顶高程测量点以及参考点的相关参数,通过基本几何定理及坐标变就能够计算出干滩长度和滩顶高程;
(4)操作简便,成本低,人工采集需要较多的人力物力,并具有一定的危险性,干滩自动巡航算法实现简单,没有昂贵的硬件要求,在软件实现时,还可以通过设定采集周期实现周期采集或定时采集,在周期采集的同时,也可以进行一次单次采集。
附图说明
图1是尾矿库干滩长度的自动巡航算法原理图;
图2是尾矿库干滩长度的计算原理图;
图3是干滩自动巡航软件流程图。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
用于尾矿库干滩长度与滩顶高程的自动巡航算法,包括以下步骤:
步骤1、将干滩测量仪置于尾矿库附近的高点,测量并记录干滩测量仪的高程;
步骤2、预设干滩长度巡航轨迹和滩顶高程巡航轨迹,根据干滩长度巡航轨迹寻找干滩长度测量点,根据滩顶高程巡航轨迹寻找滩顶高程测量点;
步骤3、记录干滩测量仪到干滩长度测量点的距离,以及干滩测量仪与干滩长度测量点的水平巡航角、垂直巡航角;在坝体上选择两个参考点,记录干滩测量仪分别到两个坝体参考点的距离,以及干滩测量仪分别与两个坝体参考点的水平巡航角、垂直巡航角,计算干滩长度;
步骤4、记录干滩测量仪到滩顶高程测量点的距离,以及干滩测量仪与滩顶高程测量点的水平巡航角、垂直巡航角,计算滩顶高程。
干滩测量仪包括图像处理模块,图像处理模块分别连接有激光测距仪和云台,云台上设置有摄像机。
如图1所示,步骤2中,以预设三条干滩长度巡航轨迹为例:S1,S2,S3分别表示三条干滩长度巡航轨迹的起始点,E1,E2,E3分别表示三条干滩长度巡航轨迹的终止点,干滩长度巡航轨迹以水面处为起点,以干滩处为终点,并经过水滩分界线;预设一条滩顶高程巡航轨迹,滩顶高程巡航轨迹经过干滩最高点。
这是因为水滩分界点的分布并不是一条平行于尾矿坝坝体的直线,因此需要设定多条干滩长度巡航轨迹,采集多个干滩长度,然后取平均值;干滩长度测量点位于坝面与水面的结合处,干滩长度测量点,就是从测量值无距离到测量值有距离的分界处。寻找滩顶高程测量时,因为只会产生一个滩顶 高程,只需要找一条轨迹即可,不需要计算平均值。
如图3所示,步骤2中,寻找干滩长度测量点和滩顶高程测量点的过程中,巡航操作具体如下:
以S1E1巡航轨迹为例,在巡航轨迹起点S1处,水平巡航角为αS1,垂直巡航角为βS1;在巡航轨迹终点E1处,水平巡航角为αE1,垂直巡航角为βE1;巡航过程中,水平巡航角为α’,垂直巡航角为β’:
对水平巡航角α’进行变换:
计算dltα’=(αE1-αS1) (1)
(1)当dltα’>0且|dltα’|>180°,Kα1=360°-|dltα’|;在αS1的基础上每次减少0.1°(若αS1≤0°时,αS1=αS1+360°);
(2)当dltα’>0且|dltα’|<180°时,Kα1=|dltα’|;如果dltα’<0,在αS1的基础上每次增加0.1°;
(3)当dltα’<0且|dltα’|>180°,Kα1=360°-|dltα’|;在αS1的基础上每次增加0.1°(若αS1≥360°时,αS1=αS1-360°);
(4)当dltα’<0且|dltα’|<180°时,Kα1=|dltα’|;如果dltα’<0,在αS1的基础上每次减少0.1°;
对垂直巡航角β’进行变换:
计算dltβ’=(βE1-βS1) (2)
(1)当dltβ’>0且|dltβ’|>180°,Kβ1=360°-|dltβ’|,K=Kβ1/Kα1;在βS1的基础上每次减少K*0.1°(若βS1≤0°时,βS1=βS1+360°);
(2)当dltβ’>0且|dltβ’|<180°时,Kβ1=|dltβ’|,K=Kβ1/Kα1;如果dltβ’<0,在βS1的基础上每次增加K*0.1°;
(3)当dltβ’<0且|dltβ’|>180°,Kβ1=360°-|dltβ’|,K=Kβ1/Kα1;在βS1的基础 上每次增加K*0.1°(若βS1≥360°时,βS1=βS1-360°);
(4)当dltβ’<0且|dltβ’|<180°时,Kβ1=|dltβ’|,K=Kβ1/Kα1;如果dltβ’<0,在βS1的基础上每次减少K*0.1°。
寻找一个干滩长度测量点的过程中,当激光测量仪照射到水面时,无距离返回,当激光测量仪照射到干滩时,有距离返回,水滩分界点W’,即为干滩长度的测量点;
沿一条干滩长度巡航轨迹S1E1寻找水滩分界点W’的过程中,转动云台,每次对水平巡航角α’与垂直巡航角β’进行变换,用激光测距仪测量距离,即认为找到了水滩分界点W’,此时停止对干滩长度巡航轨迹S1E1的巡航,记录滩测量仪到干滩长度测量点W’的距离SW’、水平巡航角αW’、垂直巡航角βW’
当巡航角度超过了终点(以dltα’>0且|dltα’|<180°为例,在αS1的基础上每次增加0.1°,当增加后的αS1超过了αE1)时,即超出了干滩长度巡航轨迹,认为该条干滩长度轨迹未能测出干滩长度。
如图2所示,计算每条干滩长度巡航轨迹测得的干滩长度的步骤为:
以水滩干滩长度测量点W’为例,计算干滩长度巡航轨迹S1E1测得的干滩长度:
(1)计算分界点在以干滩测量仪W’点为原点的坐标系下的坐标为:
Xw'=SW'cosβW'sinαW' (3)
Yw'=SW'cosβW'cosαW' (4)
Zw'=SW'sinβW' (5)
其中,SW’、αW’、βW’分别为干滩测量仪与干滩长度测量点W’的距离、干滩测量仪与干滩长度测量点W’的水平巡航角、干滩测量仪与干滩长度测量 点W’的垂直巡航角;
计算坝体参考点P1与坝面参考点P2在W’坐标系下的坐标分别为(X1w’,Y1w’,Z1w’)与(X2w’,Y2w’,Z2w’):
X1w'=S1cosβ1sinα1 (6)
Y1w'=S1cosβ1cosα1 (7)
Z1w'=S1sinβ1 (8)
X2w'=S2cosβ2sinα2 (9)
Y2w'=S2cosβ2cosα2 (10)
Z2w'=S2sinβ2 (11)
其中,(X1,Y1,Z1)为坝体参考点P1的坐标,(X2,Y2,Z2)为坝体参考点P2的坐标,S1、α1、β1分别为干滩测量仪与坝体参考点P1的距离、干滩测量仪与坝体参考点P1的水平巡航角、干滩测量仪与坝体参考点P1的垂直巡航角,S2、α2、β2分别为干滩测量仪与坝体参考点P2的距离、干滩测量仪与坝体参考点P2的水平巡航角、干滩测量仪与坝体参考点P2的垂直巡航角;
(2)将以水滩分界点为原点的W’坐标系平移至坝体参考点1为原点的M坐标系,得到水滩分界点在M坐标系下的坐标为(Xm,Ym,Zm):
即:
Xm=Xw'-X1w' (13)
Ym=Yw'-Y1w' (14)
Zm=Zw'-Z1w' (15)
得到坝体参考点P2在M坐标系下的坐标(X2m,Y2m,Z2m):
即:
X2m=X2w'-X1w' (17)
Y2m=Y2w'-Y1w' (18)
Z2m=Z2w'-Z1w' (19)
(3)将M坐标系旋转,使Y坐标轴与两个坝体参考点P1、P2的连线重合,计算旋转的θ:
θ=tan-1(Y2m/X2m) (20)
(4)计算干滩长度:
其中,Xc’、Yc’、Zc’分别是以P1、P2所在直线为Y轴的坐标系的X、Y、Z轴,干滩长度即为干滩长度测量点(水滩分界点)距离Yc’的最小距离,即,
Xc'=Xm cosθ+Ym sinθ (22)
对三条干滩长度巡航轨迹算出的干滩长度Xc’取均值,得出干滩长度d。
在寻找滩顶高程的测量点的过程中,通过转动云台,增加或减少水平巡航角α、垂直巡航角β进行测量;
在测量的距离连续出现之后,出现了测量值的突变(例如某一次的测量值与上一次测量值的差大于10),即认为滩顶高程测量点为突变前的点G;
根据干滩测量仪的高程、干滩测量仪与测点的垂直巡航角βG,运用勾股定理计算出滩顶高程测量点与干滩测量仪的相对高程,滩顶高程测量点与干滩测量仪的相对高程加上干滩测量仪的高程,即为滩顶高程。
测量干滩测量仪的高程的方法如下:选择一个高程已知点,通过摄像机找到该点,根据激光测距仪测出该点与干滩测量仪的距离,以及该点与干滩测量仪的垂直夹角,根据勾股定理计算出干滩测量仪相对于高程已知点的高程;干滩测量仪相对于测点的高程加上高程已知点的高程,即为干滩测量仪的高程。
通过上述方式,本发明用于尾矿库干滩长度与滩顶高程的自动巡航算法,避免了干滩因放矿、降雨等因素变化引起干滩测量仪置位置的变化,实现了全天候采集;克服了人工采集时,受天气、环境的制约,以及测量人员自身因素产生的错误和误差,测量人员足不出户即可对干滩长度和滩顶高程进行测量;计算方法简便,准确率高,通过基本几何定理及坐标变就能够计算出干滩长度和滩顶高程;操作简便,成本低,还可以通过设定采集周期实现周期采集或定时采集,在周期采集的同时,也可以进行一次单次采集。

Claims (8)

1.用于尾矿库干滩长度与滩顶高程的自动巡航算法,其特征在于,包括以下步骤:
步骤1、将干滩测量仪置于尾矿库附近的高点,测量并记录干滩测量仪的高程;
步骤2、预设干滩长度巡航轨迹和滩顶高程巡航轨迹,根据干滩长度巡航轨迹寻找干滩长度测量点,根据滩顶高程巡航轨迹寻找滩顶高程测量点;
步骤3、记录干滩测量仪到干滩长度测量点的距离,以及干滩测量仪与干滩长度测量点的水平巡航角、垂直巡航角;在坝体上选择两个参考点,记录干滩测量仪分别到两个坝体参考点的距离,以及干滩测量仪分别与两个坝体参考点的水平巡航角、垂直巡航角,计算干滩长度;
步骤4、记录干滩测量仪到滩顶高程测量点的距离,以及干滩测量仪与滩顶高程测量点的水平巡航角、垂直巡航角,计算滩顶高程。
2.根据权利要求1所述的用于尾矿库干滩长度与滩顶高程的自动巡航算法,其特征在于,干滩测量仪包括图像处理模块,图像处理模块分别连接有激光测距仪和云台,云台上设置有摄像机。
3.根据权利要求2所述的用于尾矿库干滩长度与滩顶高程的自动巡航算法,其特征在于,步骤2中,预设多条干滩长度巡航轨迹,干滩长度巡航轨迹以水面处为起点,以干滩处为终点,并经过水滩分界线;预设一条滩顶高程巡航轨迹,滩顶高程巡航轨迹经过干滩最高点。
4.根据权利要求3所述的用于尾矿库干滩长度与滩顶高程的自动巡航算法,其特征在于,步骤2中,寻找干滩长度测量点和滩顶高程测量点的过程中,巡航操作具体如下:
每条巡航轨迹的起点为S,终点为E,在巡航轨迹起点S处,水平巡航角为αS,垂直巡航角为βS;在巡航轨迹终点E处,水平巡航角为αE,垂直巡航角为βE;巡航过程中,水平巡航角为α,垂直巡航角为β:
对水平巡航角α进行变换:
计算dltα=(αE-αS) (1)
(1)当dltα>0且|dltα|>180°,Kα=360°-|dltα|;在αS的基础上每次减少0.1°(若αS≤0°时,αS=αS+360°);
(2)当dltα>0且|dltα|<180°时,Kα=|dltα|;如果dltα<0,在αS的基础上每次增加0.1°;
(3)当dltα<0且|dltα|>180°,Kα=360°-|dltα|;在αS的基础上每次增加0.1°(若αS≥360°时,αS=αS-360°);
(4)当dltα<0且|dltα|<180°时,Kα=|dltα|;如果dltα<0,在αS的基础上每次减少0.1°;
对垂直巡航角β进行变换:
计算dltβ=(βE-βS) (2)
(1)当dltβ>0且|dltβ|>180°,Kβ=360°-|dltβ|,K=Kβ/Kα;在βS的基础上每次减少K*0.1°(若βS≤0°时,βS=βS+360°);
(2)当dltβ>0且|dltβ|<180°时,Kβ=|dltβ|,K=Kβ/Kα;如果dltβ<0,在βS的基础上每次增加K*0.1°;
(3)当dltβ<0且|dltβ|>180°,Kβ=360°-|dltβ|,K=Kβ/Kα;在βS的基础上每次增加K*0.1°(若βS≥360°时,βS=βS-360°);
(4)当dltβ<0且|dltβ|<180°时,Kβ=|dltβ|,K=Kβ/Kα;如果dltβ<0,在βS的基础上每次减少K*0.1°。
5.根据权利要求4所述的用于尾矿库干滩长度与滩顶高程的自动巡航算法,其特征在于,寻找干滩长度测量点的过程中,当激光测量仪照射到水面时,无距离返回,当激光测量仪照射到干滩时,有距离返回,水滩分界点W,即为干滩长度测量点;
沿一条干滩长度巡航轨迹寻找水滩分界点的过程中,转动云台,每次对水平巡航角α与垂直巡航角β进行变换,用激光测距仪测量距离,若连续三次转云台后,都测量出了距离,即认为找到了水滩分界点W,此时停止对该条干滩长度巡航轨迹的巡航,记录干滩测量仪与干滩长度测量点W的距离SW、水平巡航角αW、垂直巡航角βW
当巡航角度超过了终点时,即超出了干滩长度巡航轨迹,认为该条干滩长度轨迹未能测出干滩长度。
6.根据权利要求5所述的用于尾矿库干滩长度与滩顶高程的自动巡航算法,其特征在于,计算每条干滩长度巡航轨迹得到的干滩长度的步骤为:
(1)计算水滩分界点在以干滩测量仪W点为原点的坐标系下的坐标为(XW,YW,ZW):
Xw=SWcosβWsinαW (3)
Yw=SWcosβWcosαW (4)
Zw=SWsinβW (5)
其中,SW、αW、βW分别为干滩测量仪与干滩长度测量点W的距离、干滩测量仪与干滩长度测量点W的水平巡航角、干滩测量仪与干滩长度测量点W的垂直巡航角;
计算坝体参考点P1与坝面参考点P2在W坐标系下的坐标分别为(X1w,Y1w,Z1w)与(X2w,Y2w,Z2w):
X1w=S1cosβ1sinα1 (6)
Y1w=S1cosβ1cosα1 (7)
Z1w=S1sinβ1 (8)
X2w=S2cosβ2sinα2 (9)
Y2w=S2cosβ2cosα2 (10)
Z2w=S2sinβ2 (11)
其中,(X1,Y1,Z1)为坝体参考点P1的坐标,(X2,Y2,Z2)为坝体参考点P2的坐标,S1、α1、β1分别为干滩测量仪与坝体参考点P1的距离、干滩测量仪与坝体参考点P1的水平巡航角、干滩测量仪与坝体参考点P1的垂直巡航角,S2、α2、β2分别为干滩测量仪与坝体参考点P2的距离、干滩测量仪与坝体参考点P2的水平巡航角、干滩测量仪与坝体参考点P2的垂直巡航角;
(2)将以水滩分界点为原点的W坐标系平移至坝体参考点P1为原点的M坐标系,得到水滩分界点在M坐标系下的坐标为(Xm,Ym,Zm):
即:
Xm=Xw-X1w (13)
Ym=Yw-Y1w (14)
Zm=Zw-Z1w (15)
得到坝体参考点P2在M坐标系下的坐标(X2m,Y2m,Z2m):
即:
X2m=X2w-X1w (17)
Y2m=Y2w-Y1w (18)
Z2m=Z2w-Z1w (19)
(3)将M坐标系旋转,使Y坐标轴与两个坝体参考点P1、P2的连线重合,计算旋转的角度θ:
(4)计算干滩长度:
其中,Xc、Yc、Zc分别是以P1、P2所在直线为Y轴的坐标系的X、Y、Z轴,干滩长度即为干滩长度测量点(水滩分界点)距离Yc的最小距离,即,
Xc=Xmcosθ+Ymsinθ (22)
对每条干滩长度巡航轨迹算出的干滩长度Xc取均值,得出干滩长度d。
7.根据权利要求5所述的用于尾矿库干滩长度与滩顶高程的自动巡航算法,其特征在于,在寻找滩顶高程的测量点的过程中,通过转动云台,增加或减少水平巡航角α、垂直巡航角β,进行测量;
在测量到的距离连续出现之后,出现了测量值的突变(例如某一次的测 量值与上一次测量值的差大于10),即认为滩顶高程测量点为突变前的点G;
根据干滩测量仪的高程、干滩测量仪与测点的垂直巡航角βG,运用勾股定理计算出滩顶高程测量点与干滩测量仪的相对高程,滩顶高程测量点与干滩测量仪的相对高程加上干滩测量仪的高程,即为滩顶高程。
8.根据权利要求7所述的用于尾矿库干滩长度与滩顶高程的自动巡航算法,其特征在于,测量干滩测量仪的高程的方法如下:选择一个高程已知点,通过摄像机找到该点,根据激光测距仪测出该点与干滩测量仪的距离,以及该点与干滩测量仪的垂直夹角,根据勾股定理计算出干滩测量仪相对于高程已知点的高程;干滩测量仪相对于测点的高程加上高程已知点的高程,即为干滩测量仪的高程。
CN201611193413.1A 2016-12-21 2016-12-21 用于尾矿库干滩长度与滩顶高程的自动巡航算法 Pending CN106595599A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611193413.1A CN106595599A (zh) 2016-12-21 2016-12-21 用于尾矿库干滩长度与滩顶高程的自动巡航算法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611193413.1A CN106595599A (zh) 2016-12-21 2016-12-21 用于尾矿库干滩长度与滩顶高程的自动巡航算法

Publications (1)

Publication Number Publication Date
CN106595599A true CN106595599A (zh) 2017-04-26

Family

ID=58602521

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611193413.1A Pending CN106595599A (zh) 2016-12-21 2016-12-21 用于尾矿库干滩长度与滩顶高程的自动巡航算法

Country Status (1)

Country Link
CN (1) CN106595599A (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103837196A (zh) * 2014-03-17 2014-06-04 湖南五舟检测科技有限公司 一种尾矿库干滩参数监测系统及方法
CN103852041A (zh) * 2014-03-05 2014-06-11 北京矿冶研究总院 一种尾矿库干滩长度在线监测的方法
CN104776828A (zh) * 2015-03-19 2015-07-15 中钢集团马鞍山矿山研究院有限公司 尾矿库闭库后滩面坡度倾向观测系统及使用方法
CN104831678A (zh) * 2015-01-30 2015-08-12 王健 一种用于闭库后尾矿库复垦的生态修复方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103852041A (zh) * 2014-03-05 2014-06-11 北京矿冶研究总院 一种尾矿库干滩长度在线监测的方法
CN103837196A (zh) * 2014-03-17 2014-06-04 湖南五舟检测科技有限公司 一种尾矿库干滩参数监测系统及方法
CN104831678A (zh) * 2015-01-30 2015-08-12 王健 一种用于闭库后尾矿库复垦的生态修复方法
CN104776828A (zh) * 2015-03-19 2015-07-15 中钢集团马鞍山矿山研究院有限公司 尾矿库闭库后滩面坡度倾向观测系统及使用方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
张诚: "基于Silverlight的尾矿库安全监测系统", 《工程科技Ⅰ辑》 *

Similar Documents

Publication Publication Date Title
CN111830528A (zh) 基于激光点云的杆塔特征点自动识别及倾斜参数自动测量方法
CN110567680B (zh) 一种基于角度比对的轨道扣件松动检测方法
CN110490918B (zh) 基于三维激光扫描技术的既有铁路道岔岔心提取方法
CN104260751B (zh) 一种多传感器融合的高铁轨道中心线检测系统和检测方法
CN106296694B (zh) 杆塔倾斜智能图像辨识测量方法
CN107833203A (zh) 一种基于图像处理的水平面识别及水位实时测量方法
CN109029283B (zh) 一种基于高度比对的轨道扣件螺栓浮起检测方法
CN112950532B (zh) 一种列车受电弓状态检测方法
CN107685748B (zh) 基于激光位移传感器的列车车轮尺寸参数在线检测方法
CN106157219B (zh) 基于车载扫描系统的道路中心线提取方法和装置
CN108921165B (zh) 基于水尺图像的水位识别方法
Geng et al. UAV-LiDAR-based measuring framework for height and stagger of high-speed railway contact wire
CN105241384B (zh) 一种基于计算机视觉的钢轨焊缝错边量测量方法
CN103852041B (zh) 一种尾矿库干滩长度在线监测的方法
CN103993547A (zh) 带有相对基准点的线激光路面车辙检测系统及方法
CN105571508A (zh) 接触网受电弓的形变检测方法及系统
CN111717244A (zh) 一种列车自动驾驶感知方法和系统
CN108797241B (zh) 一种基于高度比对的轨道扣件螺帽松动检测方法
CN107657144A (zh) 一种基于船舶ais和计程仪数据的近岸流场反演方法
CN110490342A (zh) 一种基于Faster R-CNN的接触网静态几何参数检测方法
CN103552570B (zh) 一种铁路轨道平顺性检测的车载近景摄影测量方法
CN101718527A (zh) 一种基于图像识别计算尾矿库干滩长度的方法
CN204059164U (zh) 一种带有相对基准点的线激光路面车辙检测系统
CN104006804A (zh) 基于观测基准失稳补偿的接触网偏移量检测方法
CN106595599A (zh) 用于尾矿库干滩长度与滩顶高程的自动巡航算法

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
RJ01 Rejection of invention patent application after publication

Application publication date: 20170426

RJ01 Rejection of invention patent application after publication