技术领域
[0001] 本
发明属于矿用掘进机技术领域，具体涉及一种矿用悬臂式掘进机可视化辅助导引方法和系统。
背景技术
[0002] 悬臂式掘进机被广泛地应用于
煤矿巷道、地
铁隧道以及铁路工程、洞穴的掘进施工中，但是，悬臂式掘进机截割煤岩时，由于司机所处
位置的原因，司机只能看到一侧的截割，另一侧往往需要有一位副司机专
门看着并用灯光指示司机截割，同时，由于截割产尘和喷雾洒
水等原因，导致巷道形成
质量不好，当巷道断面为半圆拱形或三心拱形时，司机更是难以把握。目前巷道的掘进工艺主要有两种，一是靠司机的经验把握，另一种就是根据位移
传感器测量油缸的伸缩量进行计算得到截割头的
位姿。第一种方式不可靠，容易出现超挖欠挖；第二种方式相较于第一种有了很大的提升，但是，由于井下环境非常复杂和振动的影响较大，易导致
接触式传感器接触不良，使得检测结果可靠性大大降低。
发明内容
[0003] 针对
现有技术存在的缺点，本发明的提供了一种矿用悬臂式掘进机可视化辅助导引方法和系统，解决现有的悬臂式掘进机工作过程中受人为和环境因素的影响难以保证巷道掘进效率和质量的问题。
[0004] 为了实现上述目的，本发明采用如下技术方案予以实现：
[0005] 矿用悬臂式掘进机可视化辅助导引方法，该方法中的掘进机悬臂上设置有标靶，标靶的边上设置有多个
光源点，该方法包括以下步骤：
[0006] 步骤1：采集掘进机运行状态信息以及掘进机悬臂上的标靶图像，所述的掘进机运行状态信息包括掘进机
机身的
横滚角、
俯仰角、
偏航角、掘进机
机体在巷道中的位置以及掘进机机身与两侧煤壁的距离、掘进机机身距离前煤壁的距离；
[0007] 步骤2：通过视觉测量的方法确定掘进机截割头相对于机体的位姿；
[0008] 步骤2.1：提取步骤1采集的标靶图像中的光斑，对光斑进行排序、直线拟合，将每两条直线的交点作为特征点；
[0009] 步骤2.2：利用四点透视成像方法确定所有特征点在摄像机
坐标系中的坐标；
[0010] 步骤2.3：利用对偶四元数的误差模型，根据特征点在标靶坐标系中的坐标以及步骤2.2得到的特征点在摄像机坐标系中的坐标，计算截割头的水平摆角θ1和截割头垂直摆角θ2；
[0011] 步骤2.4：根据步骤2.3确定的θ1和θ2，通过式(1)计算截割头相对于掘进机机体的位姿
[0012]
[0013] 其中，θ1为截割头水平摆角，θ2为截割头垂直摆角，d为油缸的伸缩距离，b1为掘进机抬升关节与伸缩关节之间的高度差；a1为掘进机回转台中心与抬升关节之间水平距离，a2为掘进机抬升关节与伸缩关节之间的距离，a3为掘进机伸缩关节与截割头之间的水平距离；
[0014] 步骤3：根据步骤1获得的掘进机运行状态信息，通过式(2)得到掘进机相对于巷道的位姿
[0015]
[0016] 式中，γ为掘进机机身的横滚角，β为掘进机机身的俯仰角，α为掘进机机身的偏航角，D1为机身到侧方煤壁的距离，L为掘进机到前煤壁的距离；D3为掘进机机身宽度，b2为抬升关节与地面之间的距离；
[0017] 步骤4：根据步骤2得到的截割头相对于机体位姿数据和步骤3得到的掘进机机身相对于巷道断面的位姿数据，通过式(3)可得到截割头相对于巷道断面的位姿[0018]
[0019] 步骤5：根据掘进机工作面巷道断面类型、截割头直径、巷道尺寸、走向以及截割工艺确定掘进机截割头的导引轨迹，并将导引轨迹在上位机上进行图形化界面显示；
[0020] 步骤6：将步骤4得到截割头相对于巷道断面的位姿数据与上位机进行数据交互，将截割头在上位机上图形化界面显示，并根据步骤5确定的导引轨迹进行截割头移动过程中的轨迹导引。
[0021] 具体的，所述的步骤2.3中截割头水平摆角θ1和截割头垂直摆角θ2的计算过程为：
[0022] 步骤2.3.1：利用对偶四元数误差模型构建矩阵A，式(4)所示，
[0023]
[0024]
[0025] 式中，Q(PiW)表示第i个特征点在标靶坐标系下坐标的
正交矩阵；W(PiC)表示第i个特征点在摄像机坐标系下坐标的正交矩阵； 表示第i个特征点在标靶坐标系中的坐标， 表示第i个特征点在摄像机坐标系下的坐标，i＝1,2,...,N，N表示特征点的个数；
[0026] 当矩阵A的特征值最大时，模型的误差最小，求取此时矩阵A的
特征向量r，其中， 表示特征向量r里的四个元素；
[0027] 步骤2.3.2：根据步骤2.3.1得到的特征向量r，通过式(5)确定摄像机坐标系与标靶坐标系之间的旋转矩阵R，
[0028]
[0029] 步骤2.3.3：根据步骤2.3.2获得的旋转矩阵R，通过式(6)中相等矩阵对应元素相等原理计算得到θ1和θ2，
[0030]
[0031] 式中，R为步骤2.3.2计算得到的旋转矩阵，R*为θ1＝0°、θ2＝-90°时，根据步骤2.3.2计算的特征点在摄像机坐标系下的坐标 所对应的旋转矩阵R。
[0032] 具体的，所述的步骤1中，利用捷联惯导传感器获取悬臂式掘进机机体的横滚角、俯仰角和偏航角以及掘进机机体在巷道中的位置信息；利用
超声波传感器测量掘进机机身到两侧煤壁的距离和掘进机机身到前煤壁的距离。
[0033] 本发明还公开了矿用悬臂式掘进机可视化辅助导引系统，该系统中所述的掘进机悬臂上设置有标靶，标靶的边上设置有多个光源点，该系统包括：
[0034] 通信模
块，采集掘进机运行状态信息以及掘进机悬臂上的标靶图像，所述的掘进机运行状态信息包括掘进机机身的横滚角、俯仰角、偏航角以及掘进机机身与两侧煤壁的距离、掘进机机身距离前煤壁的距离；
[0035] 视觉测量模块，用于测量截割头相对于掘进机机体的位姿，具体包括：
[0036] 提取通信模块采集的标靶图像中的光斑，对光斑进行排序、直线拟合，将每两条直线的交点作为特征点；
[0037] 利用四点透视成像原理确定所有特征点在摄像机坐标系中的坐标；
[0038] 利用对偶四元数的误差模型，根据特征点在标靶坐标系中的坐标以及特征点在摄像机坐标系中的坐标，计算截割头的水平摆角θ1和截割头垂直摆角θ2；
[0039] 根据θ1和θ2，通过式(1)计算截割头相对于掘进机机体的位姿
[0040]
[0041] 其中，θ1为截割头水平摆角，θ2为截割头垂直摆角，d为油缸的伸缩距离，b1为掘进机抬升关节与伸缩关节之间的高度差；a1为掘进机回转台中心与抬升关节之间水平距离，a2为掘进机抬升关节与伸缩关节之间的距离，a3为掘进机伸缩关节与截割头之间的水平距离；
[0042] 掘进机位姿检测模块，用于掘进机机身相对于巷道的位姿数据，根据通信模块获得的掘进机运行状态信息，利用公式(2)得到掘进机机身相对于巷道的位姿[0043]
[0044] 式中，γ为掘进机机体的横滚角，β为掘进机机体的俯仰角，α为掘进机机体的偏航角，D1为机身到侧方煤壁的距离，L为掘进机到前煤壁的距离；D3为掘进机机身宽度，b2为抬升关节与地面之间的距离；
[0045]
数据处理模块，用于确定截割头相对于巷道断面的位姿数据，根据视觉测量模块得到的截割头相对于机体位姿数据和掘进机位姿检测模块得到的掘进机机身相对于巷道的位姿数据，通过式(3)得到截割头相对于巷道断面的位姿数据
[0046]
[0047] 轨迹显示模块，用于显示掘进机截割头的导引轨迹，根据掘进机工作面巷道断面类型、截割头直径、巷道尺寸、走向以及截割工艺确定掘进机截割头的导引轨迹，并将导引轨迹在上位机上进行图形化界面显示；
[0048] 轨迹导引模块，用于截割头移动过程中的轨迹导引，将数据处理模块得到的截割头相对于巷道断面的位姿数据与上位机进行数据交互，将截割头在上位机上图形化界面显示，并根据轨迹显示模块确定的导引轨迹进行截割头移动过程中的轨迹导引。
[0049] 具体的，所述的视觉测量模块中截割头的水平摆角θ1和截割头垂直摆角θ2的确定过程为：
[0050] 利用对偶四元数误差模型构建矩阵A，如式(4)所示，
[0051]
[0052]
[0053] 式中，Q(PiW)表示第i个特征点在标靶坐标系下坐标的正交矩阵；W(PiC)表示第i个特征点在摄像机坐标系下坐标的正交矩阵； 表示第i个特征点在标靶坐标系中的坐标， 表示第i个特征点在摄像机坐标系下的坐标，i＝1,2,...,N，N表示特征点的个数；
[0054] 当矩阵A的特征值最大时，模型的误差最小，求取此时矩阵A的特征向量r，其中， 表示特征向量r里的四个元素；
[0055] 根据得到的特征向量r，通过式(5)确定摄像机坐标系与标靶坐标系之间的旋转矩阵R，
[0056]
[0057] 根据获得的旋转矩阵R，通过式(6)中相等矩阵对应元素相等原理计算得到θ1和θ2，
[0058]
[0059] 式中，R为旋转矩阵，R*为θ1＝0°、θ2＝-90°时特征点在摄像机坐标系下的坐标所对应的旋转矩阵R。
[0060] 具体的，所述的通信模块中，利用捷联惯导传感器获取悬臂式掘进机机体的横滚角、俯仰角和偏航角以及掘进机机体在巷道中的位置信息；利用
超声波测距传感器和激光测距传感器测量机身到两侧煤壁的距离和掘进机到前煤壁的距离。
[0061] 与现有技术相比，本发明的有益效果是：
[0062] (1)本发明的通过视觉测量的方法计算出的掘进机机身和截割头空间实时位置，并实现图形化的界面实时显示，可实时了解截割头运动的真实移动轨迹，解决了煤矿井下掘进机工作环境粉尘大、噪音大，操作人员难以掌握掘进断面具体状况。掘进机司机可以通过界面显示以及构建的截割路线的提示，完成可视化的截割轨迹引导，即保证截割质量的要求，同时也避免巷道超挖和欠挖得发生，同时大大提高了巷道掘进的工作效率及工作质量。
[0063] (2)本发明采用视觉测量的方法测量截割头相对于掘进机机体的位姿，并结合捷联惯性导航技术进行组合
定位确定截割头在巷道中的位姿，提高了位姿检测
精度。
附图说明
[0064] 图1为四点透视成像方法原理图。
[0065] 图2为掘进机各个坐标系示意图。
[0068] 图5为仿真实验中掘进机及导引轨迹的可视化界面显示图。
[0069] 图6为仿真实验中掘进机可视化截割过程可视化界面显示图。
具体实施方式
[0070] 本发明的标靶一般是指方形板，并在方形板的四个边上均匀设置有16个红外光源，所述的光斑即为光源在图像中形成的点。
[0071] 本发明中的摄像机坐标系是指以摄像机的
光心为原点X轴沿水平方向，Y轴沿竖直方向，Z轴垂直Y、X轴指向标靶方向。标靶坐标系是指以标靶中心为原点，X轴平行于标靶上下边指向右，Y平行于标靶左右边指向下，Z轴垂直于标靶平面。
[0072] 本发明的矿用悬臂式掘进机可视化辅助导引方法，该方法中的掘进机悬臂上设置有标靶，标靶的边上设置有多个光源点，该方法包括以下步骤：
[0073] 步骤1：采集掘进机运行状态信息以及掘进机悬臂上的标靶图像，具体，使用安装在掘进机机身上的摄像机采集安装在掘进机悬臂上的标靶图像，使用捷联惯导传感器获得掘进机在工作过程中掘进机机身的位姿，即横滚角、俯仰角、偏航角，以及掘进机机体在巷道中的位置信息，使用安装在机身两侧
超声波传感器获得掘进机与两侧煤壁的距离、使用安装在掘进机机身前侧的激光测距传感器测量机身距离前煤壁的距离。利用捷联惯导获得的位置信息确定掘进机的地理位置，以确保掘进机在工作时的方向，防止巷道打偏；再利用横滚角、俯仰角、偏航角以及超声波等传感器确定掘进机相对于巷道的位置信息，进行精确定位。
[0074] 步骤2：确定截割头相对于掘进机机体的位姿数据；
[0075] 步骤2.1：对采集的标靶图像进行降噪处理，优选的使用高斯滤波进行降噪，提取降噪后图像上每个光斑的中心点，利用最小二乘法对标靶图像上的每个边上的光斑点进行拟合，将每两条直线的交点作为特征点；
[0076] 步骤2.2：利用四点透视成像方法(即P4P)确定所有特征点在摄像机坐标系中的坐标，具体包括：将步骤2.1得到的特征点的坐标为
像素坐标，将像素坐标转换到图像坐标系下，如图1所示为共面P4P原理图，像素坐标系与图像坐标系转换关系如下：
[0077]
[0078] 其中，(Xp,Yp)表示图像坐标系下的坐标，(u0,v0)为图像坐标系的原点在像素坐标系下的坐标，(u,v)表示像素坐标系下的坐标；
[0079] 在不考虑成像畸变的情况下，图像坐标系与摄像机坐标系之间的关系：
[0080]
[0081] 其中，f为摄像机的有效焦距，(XC,YC,ZC)表示摄像机坐标系下的坐标；根据以上的转换关系便可以将像素坐标系下的特征点的坐标转换得到在摄像机坐标系下的坐标；
[0082] 如图1所示，四个共面的特征点P1～P4为正方形标靶的四个
顶点，其在像面坐标系中形成的点为C1～C4，C1～C4的坐标可以通过处理图像得到像素坐标，在利用像素坐标系与像面坐标系的关系(公式(7))求得像面坐标系下坐标 再利用P4P方法计算P1～P4在摄像机坐标系中的坐标
[0083] 步骤2.3：根据特征点在标靶坐标系中的坐标以及步骤2.2得到的特征点在摄像机坐标系中的坐标，计算截割头的水平摆角θ1和截割头垂直摆角θ2，即截割头的空间
姿态角；本发明利用对偶四元数的误差模型来确定，进一步的，还可使用基于三维视觉的平行透视模型位姿测量方法，优选对偶四元数模型方法，具体利用对偶四元数的误差模型的计算过程包括：
[0084] 步骤2.3.1：利用对偶四元数误差模型确定摄像机坐标系与标靶坐标系之间的旋转矩阵R；
[0085] ①对偶四元数由两部分构成： 其中ε为对偶运算符，且对偶运算符满足ε2＝0，ε≠0；
[0086] 对偶四元数和纯四元数都具有一个相似的表示：
[0087] 对偶向量 表示坐标系绕其旋转和平移的空间曲线， u为单位向量，表明旋转和平移方向；为旋转和平移的对偶角；r和s分别为纯四元数的
实部和对偶部，其数学表达式为：r＝[r0 r1 r2 r3]T，s＝[s0 s1 s2 s3]T；
[0088] 由对偶四元数的性质可知：rTs＝sTr＝0，rTr＝1。
[0089] 两坐标系间的旋转矩阵和平移量R和T可以表示为
[0090]
[0091]
[0092] 其中，Q(r)、W(r)为正交矩阵
[0093]
[0094]
[0095] ②基于误差模型的最优解
[0096] 由坐标系之间的转换可得 其中 为特征点在在摄像机坐标系中的理论坐标， 为特征点在标靶坐标系中的坐标，T为标靶坐标系到摄像机坐标的平移量。
[0097] 在实际计算中，由于提取中心点误差、计算误差等导致特征点在摄像机坐标系中的理论坐标与实际坐标之间存在误差
[0098] 通过式 和 建立误差模型：
[0099]
[0100] 式中，N为特征点个数；
[0101] 由式 以及式(9)和式(10)可得：
[0102]
[0103] 其中，
[0104]
[0105]
[0106]
[0107] 式中，Q(PiW)表示第i个特征点在标靶坐标系下坐标的正交矩阵；W(PiC)表示第i个特征点在摄像机坐标系下坐标的正交矩阵； 表示第i个特征点在标靶坐标系中的坐标， 为第i个特征点在摄像机坐标系下的坐标，i＝1,2,...,N，N表示特征点的个数；
[0108] 利用对偶四元数的约束条件构造拉格朗函数：
[0109] L(r，s，λ1，λ2)＝rTG1r+sTG2r+G3+4NsTs+λ1(rTr-1)+λ2sT[0110] 其中，λ1和λ2为拉格朗日乘数。
[0111] 分别对r、s求偏导并计算可以求出λ2＝0，且
[0112]
[0113] 把式(15)代入式(14)中，得指标函数
[0114] F(r,s)＝(G3-λ)/N (15)
[0115] 设定 λ表示矩阵A的特征值，要使得指标函数F(r,s)最小，即误差最小，则λ最大；求取此时矩阵A的特征向量r，即对偶四元数的实部，[0116] 步骤2.3.2：根据确定的特征向量r，通过式(5)即可摄像机坐标系与标靶坐标系之间的旋转矩阵R，
[0117]
[0118] 步骤2.3.3：根据步骤2.3.2获得的旋转矩阵R，根据相等矩阵对应元素相等原理，通过式(6)中计算得到θ1和θ2，
[0119]
[0120] 式中，R为步骤2.3.2计算得到的旋转矩阵，R*为θ1＝0°、θ2＝-90°时根据步骤2.3.2计算的特征点在摄像机坐标系下的坐标 所对应的旋转矩阵R。一般通过给定θ1和θ2的方法确定R*，本发明通过设定θ1＝0°、θ2＝-90°即掘进机悬臂水平且位于机体中间的位置，确定了该条件下特征点在摄像机坐标系下的坐标 再根据步骤步骤2.3.2计算出当前位置时摄像机坐标系与标靶坐标系之间的旋转矩阵R，设此时的R为*
R。
[0121] 步骤2.4：对悬臂式掘进机进行运动学分析，如图2所示建立悬臂式掘进机坐标系，基坐标系O0X0Y0Z0的原点与截割臂抬升关节的旋
转轴线在同一高度，X0轴在掘进机的
中轴线上，Z0轴在回转台的
旋转轴线上；回转关节坐标系O1X1Y1Z1的原点与基坐标系原点重合；抬升关节坐标系O2X2Y2Z2的原点在截割臂抬升关节的中心，Z2轴在截割臂抬升关节的轴线上；伸缩关节坐标系O3X3Y3Z3的原点在伸缩油缸前关节的中心，Y3轴在关节的轴线上，Z3轴在油缸
活塞杆的轴线上；截割头坐标系O4X4Y4Z4的原点在截割头最小截割半径的中心，Z4轴在截割头锥台的轴线上；摄像机坐标系OCXCYCZC以摄像机的光心为原点；标靶坐标系OWXWYWZW以正方形标靶中心为原点，ZW轴垂直于红外标靶；巷道断面坐标系OXXXYXZX以巷道断面的左下角点为原点，ZX坐标轴垂直地面，YX坐标轴平行与巷道地面水向左。根据建立的坐标系得到截割头到掘进机机体的变化矩阵，即式(1)；再根据步骤2.3确定的θ1和θ2，计算截割头相对于掘进机机体的位姿数据
[0122]
[0123] 其中，θ1为截割头水平摆角，θ2为截割头垂直摆角，d为油缸的伸缩距离，b1为掘进机抬升关节与伸缩关节之间的高度差；a1为掘进机回转台中心与抬升关节之间水平距离，a2为掘进机抬升关节与伸缩关节之间的距离，a3为掘进机伸缩关节与截割头之间的水平距离。
[0124] 步骤3：根据步骤1获得的掘进机运行状态信息，通过式(2)得到掘进机机身相对于巷道的位姿数据
[0125]
[0126] 式中，γ为掘进机机身的横滚角，β为掘进机机身的俯仰角，α为掘进机机身的偏航角，D1为机身到侧方煤壁的距离，L为掘进机到前煤壁的距离；D3为掘进机机身宽度，b2为抬升关节与地面之间的距离；
[0127] 步骤4：融合惯性导航与视觉测量技术进行组合定位，即根据步骤2得到的截割头相对于机体位姿数据和步骤3得到的掘进机机身相对于巷道断面的位姿数据，通过式(3)可得到截割头相对于巷道断面的位姿数据
[0128]
[0129] 步骤5：根据掘进机工作面巷道断面类型、截割头直径、巷道尺寸、走向以及截割工艺(即截割路径组合方式)确定掘进机截割头的导引轨迹，并将导引轨迹在上位机上进行图形化界面显示；
[0130] 步骤6：将步骤4得到截割头相对于巷道断面的位姿数据与上位机进行数据交互，将截割头在上位机上图形化界面显示，并根据步骤5确定的导引轨迹进行截割头移动过程中的轨迹导引。可用于获取超挖欠挖及进尺查询和超挖欠挖报警等功能。
[0131] 本发明还公开了矿用悬臂式掘进机可视化辅助导引系统，该系统包括：
[0132] 通信模块，用于采集掘进机运行状态信息以及掘进机悬臂上的标靶图像，其中，掘进机运行状态信息包括掘进机机身的横滚角、俯仰角、偏航角以及掘进机机身与两侧煤壁的距离、掘进机机身距离前煤壁的距离；
[0133] 视觉测量模块，用于测量截割头相对于掘进机机体的位姿，具体包括：
[0134] 提取通信模块采集的标靶图像中的光斑，对光斑进行排序、直线拟合，将每两条直线的交点作为特征点；
[0135] 利用四点透视成像原理确定所有特征点在摄像机坐标系中的坐标，具体包括：利用像素坐标系与像面坐标系的关系(公式(7))求得像面坐标系下坐标 再利用P4P方法计算P1～P4在摄像机坐标系中的坐标
[0136] 利用对偶四元数的误差模型，根据特征点在标靶坐标系中的坐标以及特征点在摄像机坐标系中的坐标，计算截割头的水平摆角θ1和截割头垂直摆角θ2，具体包括：
[0137] 利用对偶四元数误差模型构建矩阵A，式(4)所示，
[0138]
[0139]
[0140] 式中，Q(PiW)表示第i个特征点在标靶坐标系下坐标的正交矩阵；W(PiC)表示第i个特征点在摄像机坐标系下坐标的正交矩阵； 表示第i个特征点在标靶坐标系中的坐标， 表示第i个特征点在摄像机坐标系下的坐标，i＝1,2,...,N，N表示特征点的个数；
[0141] 当矩阵A的特征值最大时，模型的误差最小，求取此时矩阵A的特征向量r，其中， 表示特征向量r里的四个元素；
[0142] 然后根据得到的特征向量r，通过式(5)确定摄像机坐标系与标靶坐标系之间的旋转矩阵R，
[0143]
[0144] 最后根据获得的旋转矩阵R，通过式(6)中相等矩阵对应元素相等原理计算得到θ1和θ2，
[0145]
[0146] 式中，R为步骤2.3.2计算得到的旋转矩阵，R*为θ1＝0°、θ2＝-90°时，根据步骤2.3.2计算的特征点在摄像机坐标系下的坐标 所对应的旋转矩阵R。
[0147] 根据θ1和θ2，通过式(1)计算截割头相对于掘进机机体的位姿数据[0148]
[0149] 其中，θ1为截割头水平摆角，θ2为截割头垂直摆角，d为油缸的伸缩距离，b1为掘进机抬升关节与伸缩关节之间的高度差；a1为掘进机回转台中心与抬升关节之间水平距离，a2为掘进机抬升关节与伸缩关节之间的距离，a3为掘进机伸缩关节与截割头之间的水平距离；
[0150] 掘进机位姿检测模块，用于掘进机机身相对于巷道的位姿数据，根据通信模块获得的掘进机运行状态信息，利用公式(2)得到掘进机机身相对于巷道的位姿数据[0151]
[0152] 式中，γ为掘进机机体的横滚角，β为掘进机机体的俯仰角，α为掘进机机体的偏航角，D1为机身到左侧煤壁的距离，L为掘进机到前煤壁的距离；D3为掘进机机身宽度，b2为抬升关节与地面之间的距离；
[0153] 数据处理模块，用于确定截割头相对于巷道断面的位姿数据，根据视觉测量模块得到的截割头相对于机体位姿数据和掘进机位姿检测模块得到的掘进机机身相对于巷道的位姿数据，通过式(3)得到截割头相对于巷道断面的位姿数据
[0154]
[0155] 轨迹显示模块，用于显示掘进机截割头的导引轨迹，根据掘进机工作面巷道断面类型、截割头直径、巷道尺寸、走向以及截割工艺确定掘进机截割头的导引轨迹，并将导引轨迹在上位机上进行图形化界面显示；
[0156] 轨迹导引模块，用于截割头移动过程中的轨迹导引，将数据处理模块得到的截割头相对于巷道断面的位姿数据与上位机进行数据交互，将截割头在上位机上图形化界面显示，并根据轨迹显示模块确定的导引轨迹进行截割头移动过程中的轨迹导引。
[0157] 仿真实验：
[0158] 本仿真实验中掘进机悬臂上的标靶为正方形标靶，标靶的四个边上均匀设置有16个红外光源灯，掘进机引导过程包括：
[0159] 步骤1：根据可视化辅助截割的需要，对掘进机运行状态信息进行采集，使用捷联惯导传感器获得掘进机在工作过程中掘进机机身的位姿、使用安装在机身两侧超声波传感器获得掘进机与两侧煤壁的距离、使用安装在掘进机机身前侧的激光测距传感器测量机身距离前煤壁的距离、以及使用防爆摄像机安装在掘进机机身上的方形标靶进行
图像采集；将获得的数据通过串口通信传输到机载计算机中，将各传感器的数据处理与计算并存入
数据库。在本仿真实验中，掘进机截割头在轨迹引导的初始阶段，各传感器实时采集的截割头的横滚角γ为0.120°，俯仰角β为-0.069°，偏航角α为0.172°，掘进机机身距离煤矿右侧的距离为32cm。
[0160] 步骤2：通过视觉测量的方法确定掘进机截割头相对于机体的位姿；具体如图2所示，步骤2.1：利用高斯滤波对提取的图像进行降噪处理，提取降噪后图像中的光斑中心点，对方形标靶上四个边的每个边上的光斑点进行排序、拟合最后确定计算特征点；
[0161] 步骤2.2：利用四点透视成像方法确定四个特征点在摄像机坐标系中的坐标；
[0162] 步骤2.3：利用对偶四元数这一数学工具，建立误差模型，建立截割头的位姿误差模型，利用拉格朗日乘子法求最小值的方法，确定相机坐标系与地面坐标系之间的旋转矩阵R；再通过式(6)所示的旋转矩阵R与θ1和θ2的关系，确定了θ1和θ2；
[0163] 步骤2.4：利用式(1)计算截割头相对于掘进机机体的位姿04T，并存入数据库；
[0164] 步骤3：根据步骤1中获得的掘进机运行状态信息，通过式(2)计算得到掘进机机身相对于巷道的位姿数据
[0165] 步骤4：通过式(3)可得到截割头相对于巷道断面的位姿数据
[0166] 步骤5：根据掘进机工作面巷道断面类型、截割头直径、巷道尺寸、走向以及截割工艺(即截割路径组合方式)确定掘进机截割头的导引轨迹，并将导引轨迹在上位机上进行图形化界面显示。如图5所示为掘进机可视化显示图，图中右上部分的弯折线表示导引轨迹，弯折曲线上的原点表示截割头的；图中左上角的视图表示掘进机在巷道中的视图，其中方框表示巷道。
[0167] 步骤6：将步骤4得到截割头相对于巷道断面的位姿数据与上位机进行数据交互，将截割头在上位机上图形化界面显示，如图5所示的原点为表示截割头；并根据步骤5确定的导引轨迹进行截割头移动过程中的轨迹导引。如图6所示，图上的截割头
覆盖部分即为截割头的行走路径。从图6可以看出，当掘进机行走到此位置时，各传感器实时采集的截割头的横滚角γ为0.008°，俯仰角β为-0.013°，偏航角α为0.132°，掘进机机身距离煤矿右侧的距离为21.1cm。
[0168] 设定掘进机的初始状态为：截割头水平摆角为0°，垂直摆角为-90°。在掘进机实际运行过程中，测量得到掘进机运行过程中某一时刻截割头抬升25.369°，向右转动15.814°；通过本仿真实验获得掘进机在该时刻截割头水平摆角15.562°和截割头垂直摆角为-
115.685°，则截割头水平摆角和截割头垂直摆角分别变化为15.562°及25.685°，与实际测量的得到的截割头的抬升角度和向右转动角度很接近，说明本发明的方法和系统具有很高的精确度。
[0169] 本发明通过视觉测量的方法计算出的截割头空间实时位置，并实现图形化的界面实时显示，掘进机司机可以通过界面显示以及构建的截割路线的提示，完成可视化的截割轨迹引导。即保证截割质量的要求，同时也避免巷道超挖和欠挖得发生。在巷道掘进过程中如果超挖和欠挖发生，可视化界面立即做出提示，掘进机司机可以根据情况做出补救措施。另外还可通过查询页面，可以对掘进距离及历史报警情况进行查看。
[0170] 需要说明的是本发明并不局限于以上具体实方式中，凡在本
申请技术方案
基础上做的等同变换均落入本发明的保护范围。
