CN109798875B - 一种在矿区视频***中实时标注采矿范围线的方法 - Google Patents
一种在矿区视频***中实时标注采矿范围线的方法 Download PDFInfo
- Publication number
- CN109798875B CN109798875B CN201910190638.9A CN201910190638A CN109798875B CN 109798875 B CN109798875 B CN 109798875B CN 201910190638 A CN201910190638 A CN 201910190638A CN 109798875 B CN109798875 B CN 109798875B
- Authority
- CN
- China
- Prior art keywords
- camera
- coordinates
- mining area
- coordinate system
- mining
- 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.)
- Active
Links
Images
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
- Studio Devices (AREA)
Abstract
本发明公开了一种在矿区视频***中实时标注采矿范围线的方法,包括:步骤1、建立矿区摄影测量坐标系:步骤2、将采矿许可证中矿区范围拐点坐标转换到矿区摄影测量坐标系中,根据采矿许可证中矿区范围拐点坐标通过GPS实地测量出拐点对应的高程;步骤3、确定摄像机的内方位元素和畸变系数;步骤4、通过GPS获得摄像机外方位元素的坐标,并将其转换到矿区摄影测量坐标系中;步骤5、计算出矿区范围拐点坐标对应的像平面坐标系上的坐标并连线,形成采矿范围线,并在矿区视频***中实时显示采矿范围线。本发明克服了物理标志桩容易被损坏及移位等缺点,同时大大减小了工作人员识别矿山越界开采的难度。
Description
技术领域
本发明涉及测绘技术领域,尤其涉及一种在矿区视频***中实时标注采矿范围线的方法。
背景技术
矿山安全事故是我国矿山普遍存在严重问题,该问题严重威胁人们生命和生活,引起了国家和社会的普遍关注。矿山业主的越界开采是引发矿山事故的重要原因之一,但是现在一直没有成熟的越界开采的实时监测技术。黄潘等人提出了利用天然地震确定震源位置方法,通过计算出***震源时空参数,来判断是否进行了越界开采。汪云甲等人的专利方法:《多源数据监测矿区形变的时空特性及越界开采识别方法》,利用合成孔径雷达判断矿区沉降与形变,识别越界开采区域。现有的方法或专利技术均要求高成本的数据监测,同时需要监测人员有着过硬的相关专业知识,并不能真正的做到实时监测越界开采发生与否的要求。
在矿山的现有监控设备与数据基础上,如果能够做到在监控画面中标明开采范围,将可以辅助监控人员实时、准确的判断越界开采情况发生与否。因此,提出一种经济实用的在矿区视频***实时标注采矿范围线的方法,是本领域极有价值的问题。
发明内容
本发明要解决的技术问题在于针对现有技术中的缺陷:目前采矿许可证中仅规定了矿区范围线拐点坐标,为了实时的发现及识别矿区存在的越界开采现象,提供一种在矿区视频***中实时标注采矿范围线的方法。
本发明解决其技术问题所采用的技术方案是:
本发明提供一种在矿区视频***中实时标注采矿范围线的方法,在矿区内设置实时监控矿区画面的矿区视频***,包括多台摄像机;该方法包括以下步骤:
步骤1、建立矿区摄影测量坐标系:根据矿区的实际地理位置和地形地貌特征选取坐标系原点,以天顶方向为坐标系z轴,以地理东方向为x轴,地理北方向为y轴建立右手空间直角坐标系;
步骤2、将采矿许可证中矿区范围拐点坐标转换到矿区摄影测量坐标系中,根据采矿许可证中矿区范围拐点坐标通过GPS实地测量出拐点对应的高程;
步骤3、通过直接线性变换解法建立摄像机像点坐标和相应物点在矿区摄影测量坐标系下的直接的线性关系,并由此确定摄像机的内方位元素和畸变系数;
步骤4、通过GPS获得摄像机外方位元素的坐标,并将其转换到矿区摄影测量坐标系中;
步骤5、根据摄像机外方位元素在矿区摄影测量坐标系中的坐标,计算出矿区范围拐点坐标对应的像平面坐标系上的坐标并连线,形成采矿范围线,并在矿区视频***中实时显示采矿范围线。
进一步地,本发明的步骤3中确定摄像机的内方位元素和畸变系数的方法具体为:
步骤3.1、建立控制场并测定控制点坐标:直接线性变换检校方法是借助实验场中坐标己知的标志点进行检校,建立控制场是基础;在控制场中布设一定数量的基本控制点,控制点保证稳定;通过建立高精度控制网的方法精密测定基本控制点的坐标,运用工程测量前方交会方法,精确测定控制点的物方空间坐标;
步骤3.2、拍摄标志点照片并获取控制点像平面坐标:在控制场中不同位置,相机以不同旋转角度拍摄控制点照片,每张照片拍摄6个以上的控制点;将照片传输到计算机中,通过像点量测软件获得标志点像平面坐标;
步骤3.3、DLT求解检校参数值:利用获得的标志点的物方空间坐标和像平面坐标,依据DLT模型计算出模型的11个l系数值,通过l系数值和相机内方位元素的关系计算出相机的内方位元素和畸变系数。
进一步地,本发明的步骤4中通过GPS获得摄像机外方位元素的坐标,并将其转换到矿区摄影测量坐标系中的具体方法为:
通过GPS获得摄像机外方位元素并将其转换到矿区摄影测量坐标系(XO,YO,ZO)中:将相机固定在工作区域;将GPS接收机天线相位中心和摄像机成像平面在水平方向对齐,记录下GPS接收机返回的位置坐标,减去GPS接收机相对摄像机的相对高度即得到摄像机的位置坐标(XC,YC,ZC);摄像机位置坐标在矿区摄影测量坐标系下的坐标为(YC-YO,XC-XO,ZG-ZO),记作(XcYc,Zc);相机的初始外方位角元素通过两台GPS接收机解算,在不考虑相机感光元件旋转的情况下,得到相机外方位角元素为
进一步地,本发明的步骤4中相机外方位角元素随着相机姿态的变化而实时更新。
进一步地,本发明的步骤5中计算像平面坐标系上的坐标的具体方法为:
步骤5.1、确定摄像机旋转矩阵:由步骤4得到的摄像机外方位元素的角元素得到对应的摄像机旋转矩阵,记作R;
步骤5.2、根据步骤5.1所确定的摄像机旋转矩阵确定各个控制点所对应的共线条件方程式;
步骤5.3、寻找位于摄影平面前方的矿区范围拐点坐标:由于位于摄影平面后方的范围点坐标一定不成像,所以要区分位于摄影平面前方的点;首先记待区分点位置矢量为g=(Xg,Yg,Zg),记相机位矢为c=(Xc,Yc,Zc),设摄像机成像面在像空间坐标系下的法向量为i=(0,0,1),则摄像机成像面在摄影测量坐标系下的法向量为R·i;若R·i·(g-c)>0,则位矢g对应的矿区范围拐点在摄影平面前方;
步骤5.4、计算矿区范围拐点坐标对应的像平面坐标系上的坐标并连线:将步骤5.3所得到的位于摄影平面前方的矿区范围拐点坐标代入步骤5.2中所计算的共线方程中,得到待求的矿区范围拐点坐标对应的像点坐标(x,y),将其按照顺序连接并显示。
进一步地,本发明的步骤5.2中的共线条件方程式的形式为:
其中,dx,dy为相机畸变修正项,相机的内方位元素为x0,y0,fx,fy。
本发明产生的有益效果是:本发明的在矿区视频***中实时标注采矿范围线的方法,现有矿区范围线标注需要在范围线拐点处设置标志桩,本发明无需设置任何物理标志桩,通过在视频中实时生成虚拟标线的方式替代物理标志桩,可克服物理标志桩容易被损坏及移位等缺点;另外,本发明的方法还具有以下优点:1.利用摄影测量技术在视频中实时生成虚拟边界线线来替代物里边界线,无需设置任何物理标志,可克服物理标志桩容易被损坏及移位等缺点。2.虚拟边界线降低了购置埋设物理标志物的成本。3.在视频中实时标注边界线,大大减小了工作人员识别矿山越界开采的难度。
附图说明
下面将结合附图及实施例对本发明作进一步说明,附图中:
图1是本发明实施例的总流程图;
图2是本发明实施例的步骤5的流程图;
图3是本发明实施例的实际应用场景图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
如图1所示,本发明实施例以露天采石矿为例对本发明的流程进行一个具体的阐述,该方法包括以下步骤:
步骤1,建立矿区摄影测量坐标系:根据矿区的实际地理位置和地形地貌特征选取坐标系原点,以天顶方向为坐标系z轴,以地理东方向为x轴,地理北方向为y轴建立右手空间直角坐标系。其坐标系原点在大地坐标系中的位置记作(XO,YO,ZO)。
建立摄影测量坐标系应以方便后续计算为原则,最好将坐标系原点设置在已有的地面控制点上。若无较好的地面控制点,可以考虑选择一空旷区域,使用GPS测量坐标并埋设地面控制点。
步骤2,将采矿许可证中矿区范围拐点坐标转换到矿区摄影测量坐标系(XO,YO,ZO)中:注意到实施例中采矿许可证中矿区范围拐点坐标为1980西安大地坐标系。矿区范围拐点坐标为大地坐标系,而大地坐标系是左手系。矿区范围拐点平面坐标记作(XG,YG),则转换后的矿区范围拐点在矿区摄影测量坐标系下的平面坐标为(YG-YO,XG-XO)。
根据采矿许可证中矿区范围拐点坐标通过GPS实地测量出拐点对应的高程,记作ZG。则矿区范围拐点在矿区摄影测量坐标系下的坐标为(YG-YO,XG-XO,ZG-ZO),我们记作(Xg,Yg,Zg)。
在实施例中,步骤1中的坐标系原点坐标是利用GPS实地测量得到,则需要先将矿区摄影测量坐标系原点坐标由WGS-84坐标系转化为1980西安大地坐标系。
实施例过程涉及矿区范围拐点坐标数据如下表:
步骤3,通过直接线性变换解法建立摄像机像点坐标和相应物点在矿区摄影测量坐标系下的直接的线性关系,并由此确定摄像机的内方位元素和畸变系数。
在本实施例中,摄像机的内方位元素包括像主点坐标(x0,y0),和摄像机的x方向主距fx,和y方向主距fy,畸变系数包括坐标轴不垂直误差dβ和比例尺不一致误差ds。
在所述的步骤3中,相机检校的具体操作方法如下。
步骤3.1,建立控制场并测定控制点坐标:直接线性变换检校方法是借助实验场中坐标己知的标志点进行检校,建立控制场是基础。控制场中要求布设一定数量的基本控制点。控制点须保证相当的稳定。通过一些建立高精度控制网的方法精密测定基本控制点的坐标,运用工程测量前方交会方法,精确测定控制点的物方空间坐标。
步骤3.2,拍摄标志点像片并获取控制点像平面坐标:在控制场中不同位置,相机以不同旋转角度拍摄控制点像片,每张相片需至少拍摄6个以上的控制点。由于视频监控领域的相机为固定焦距相机,所以不存在焦距改变的情况。将数字图像传输到计算机中,通过像点量测软件获得标志点像平面坐标。
步骤3.3,DLT求解检校参数值。利用获得的标志点的物方空间坐标和像平面坐标,依据DLT模型计算出模型的11个l系数值。通过l系数值和相机内方位元素的关系计算出相机的内方位元素和畸变系数。
步骤4,通过GPS获得摄像机外方位元素并将其转换到矿区摄影测量坐标系(XO,YO,ZO)中:将相机固定在工作区域。将GPS接收机天线相位中心和摄像机成像平面在水平方向对齐,记录下GPS接收机返回的位置坐标,减去GPS接收机相对摄像机的相对高度即得到摄像机的位置坐标(XC,YC,ZC)。摄像机位置坐标在矿区摄影测量坐标系下的坐标为(YC-YO,XC-XO,ZG-ZO),记作(XcYc,Zc)。相机的初始外方位角元素可通过两台GPS接收机解算。在不考虑相机感光元件旋转的情况下,可以得到相机外方位角元素为相机外方位角元素需要随着相机姿态的变化而实时更新。
在本实施例中,摄像机为旋转半球面相机,可以通过旋转更新姿态。旋转半球面相机在旋转过程中可以实时更新相机的外方位角元素。
步骤5,根据相机内外方位元素,计算出矿区范围拐点坐标对应的像平面坐标系上的坐标并连线:
如图2所示,在所述的步骤5中,具体操作方法如下:
步骤5.1,确定摄像机旋转矩阵:由步骤四得到的摄像机外方位元素的角元素可以得到对应的摄像机旋转矩阵,记作R。其中R矩阵中的各个元素为:
步骤5.2,根据步骤5.1所确定的摄像机旋转矩阵确定各个控制点所对应的共线条件方程式。共线条件方程是求解边界点像平面坐标的基础。
对于本例而言,加入了相机坐标轴不垂直性畸变和比例尺不一致性畸变的共线条件方程的形式为:
其中x0,y0,fx,a1,a2,a3,b1,b2,b3,c1,c2,c3,ds,dβ,(XC,YC,ZC)已知,待求的像点坐标(x,y)和该像点对应的摄影测量坐标(Xg,Yg,Zg)为未知数。
步骤5.3,寻找位于摄影平面前方的矿区范围拐点坐标:由于位于摄影平面后方的范围点坐标一定不成像,所以要区分位于摄影平面前方的点。首先记待区分点位置矢量为g=(Xg,Yg,Zg)。记相机位矢为c=(Xc,Yc,Zc)设摄像机成像面在像空间坐标系下的法向量为i=(0,0,1),则摄像机成像面在摄影测量坐标系下的法向量为R·i。
由矢量运算法则知,若R·i·(g-c)>0,则位矢g对应的矿区范围拐点在摄影平面前方。
对于本实施例来说,在半球相机改变其外方位角元素的时候,都要进行一次步骤5.2和5.3。而对于某些固定姿态的摄像机,则在初始化阶段进行一次步骤5.2和5.3即可。
步骤5.4,计算矿区范围拐点坐标对应的像平面坐标系上的坐标并连线:将步骤5.3所得到的位于摄影平面前方的矿区范围拐点坐标代入步骤5.2中所计算的共线方程中。得到待求的矿区范围拐点坐标对应的像点坐标(x,y),将其按照顺序连接并显示即可。具体显示效果见附图3。
应当理解的是,对本领域普通技术人员来说,可以根据上述说明加以改进或变换,而所有这些改进和变换都应属于本发明所附权利要求的保护范围。
Claims (1)
1.一种在矿区视频***中实时标注采矿范围线的方法,其特征在于,在矿区内设置实时监控矿区画面的矿区视频***,包括多台摄像机;该方法包括以下步骤:
步骤1、建立矿区摄影测量坐标系:根据矿区的实际地理位置和地形地貌特征选取坐标系原点,以天顶方向为坐标系z轴,以地理东方向为x轴,地理北方向为y轴建立右手空间直角坐标系;
步骤2、将采矿许可证中矿区范围拐点坐标转换到矿区摄影测量坐标系中,根据采矿许可证中矿区范围拐点坐标通过GPS实地测量出拐点对应的高程;
步骤3、通过直接线性变换解法建立摄像机像点坐标和相应物点在矿区摄影测量坐标系下的直接的线性关系,并由此确定摄像机的内方位元素和畸变系数;
步骤4、通过GPS获得摄像机外方位元素的坐标,并将其转换到矿区摄影测量坐标系中;
步骤5、根据摄像机外方位元素在矿区摄影测量坐标系中的坐标,计算出矿区范围拐点坐标对应的像平面坐标系上的坐标并连线,形成采矿范围线,并在矿区视频***中实时显示采矿范围线;
步骤5中计算像平面坐标系上的坐标的具体方法为:
步骤5.1、确定摄像机旋转矩阵:由步骤4得到的摄像机外方位元素的角元素得到对应的摄像机旋转矩阵,记作R;
步骤5.2、根据步骤5.1所确定的摄像机旋转矩阵确定各个控制点所对应的共线条件方程式;
步骤5.3、寻找位于摄影平面前方的矿区范围拐点坐标:由于位于摄影平面后方的范围点坐标一定不成像,所以要区分位于摄影平面前方的点;首先记待区分点位置矢量为g=(Xg,Yg,Zg),记相机位矢为c=(Xc,Yc,Zc),设摄像机成像面在像空间坐标系下的法向量为i=(0,0,1),则摄像机成像面在摄影测量坐标系下的法向量为R·i;若R·i·(g-c)>0,则位矢g对应的矿区范围拐点在摄影平面前方;
步骤5.4、计算矿区范围拐点坐标对应的像平面坐标系上的坐标并连线:将步骤5.3所得到的位于摄影平面前方的矿区范围拐点坐标代入步骤5.2中所计算的共线方程中,得到待求的矿区范围拐点坐标对应的像点坐标(x,y),将其按照顺序连接并显示;
步骤5.2中的共线条件方程式的形式为:
其中,dx,dy为相机畸变修正项,相机的内方位元素为x0,y0,fx,fy;
步骤3中确定摄像机的内方位元素和畸变系数的方法具体为:
步骤3.1、建立控制场并测定控制点坐标:直接线性变换检校方法是借助实验场中坐标己知的标志点进行检校,建立控制场是基础;在控制场中布设一定数量的基本控制点,控制点保证稳定;通过建立高精度控制网的方法精密测定基本控制点的坐标,运用工程测量前方交会方法,精确测定控制点的物方空间坐标;
步骤3.2、拍摄标志点照片并获取控制点像平面坐标:在控制场中不同位置,相机以不同旋转角度拍摄控制点照片,每张照片拍摄6个以上的控制点;将照片传输到计算机中,通过像点量测软件获得标志点像平面坐标;
步骤3.3、DLT求解检校参数值:利用获得的标志点的物方空间坐标和像平面坐标,依据DLT模型计算出模型的11个l系数值,通过l系数值和相机内方位元素的关系计算出相机的内方位元素和畸变系数;
步骤4中通过GPS获得摄像机外方位元素的坐标,并将其转换到矿区摄影测量坐标系中的具体方法为:
通过GPS获得摄像机外方位元素并将其转换到矿区摄影测量坐标系(XO,YO,ZO)中:将相机固定在工作区域;将GPS接收机天线相位中心和摄像机成像平面在水平方向对齐,记录下GPS接收机返回的位置坐标,减去GPS 接收机相对摄像机的相对高度即得到摄像机的位置坐标(XC,YC,ZC);摄像机位置坐标在矿区摄影测量坐标系下的坐标为(YC-YO,XC-XO,ZG-ZO),记作(Xc,Yc,Zc);相机的初始外方位角元素通过两台GPS接收机解算,在不考虑相机感光元件旋转的情况下,得到相机外方位角元素为
步骤4中相机外方位角元素随着相机姿态的变化而实时更新。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910190638.9A CN109798875B (zh) | 2019-03-13 | 2019-03-13 | 一种在矿区视频***中实时标注采矿范围线的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910190638.9A CN109798875B (zh) | 2019-03-13 | 2019-03-13 | 一种在矿区视频***中实时标注采矿范围线的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109798875A CN109798875A (zh) | 2019-05-24 |
CN109798875B true CN109798875B (zh) | 2021-08-31 |
Family
ID=66563483
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910190638.9A Active CN109798875B (zh) | 2019-03-13 | 2019-03-13 | 一种在矿区视频***中实时标注采矿范围线的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109798875B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116337012B (zh) * | 2023-03-20 | 2023-10-20 | 重庆地质矿产研究院 | 基于倾斜摄影测量技术的露天矿山界桩位置布设优化方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104618688B (zh) * | 2015-01-19 | 2017-09-29 | 荣科科技股份有限公司 | 一种可视化监控防护方法 |
CN205899722U (zh) * | 2016-08-19 | 2017-01-18 | 中国有色金属长沙勘察设计研究院有限公司 | 一种用于露天矿山开采的电子围栏*** |
-
2019
- 2019-03-13 CN CN201910190638.9A patent/CN109798875B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN109798875A (zh) | 2019-05-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103557841B (zh) | 一种提高多相机合成影像摄影测量精度的方法 | |
US11887273B2 (en) | Post capture imagery processing and deployment systems | |
US8897482B2 (en) | Stereo photogrammetry from a single station using a surveying instrument with an eccentric camera | |
US20160187130A1 (en) | Method for determining a position and orientation offset of a geodetic surveying device and such a surveying device | |
Sanz‐Ablanedo et al. | Reducing systematic dome errors in digital elevation models through better UAV flight design | |
KR101346323B1 (ko) | 지상기준점과 부가매개변수를 이용한 비측량용 디지털 카메라 촬영 시스템의 자체검정 방법 | |
CN109523471B (zh) | 一种地面坐标和广角摄像机画面坐标的转换方法、***以及装置 | |
US20130113897A1 (en) | Process and arrangement for determining the position of a measuring point in geometrical space | |
CN102034238A (zh) | 基于光学成像测头和视觉图结构的多摄像机***标定方法 | |
CN103411587B (zh) | 定位定姿方法及*** | |
CN106408601A (zh) | 一种基于gps的双目融合定位方法及装置 | |
CN113012292B (zh) | 一种基于无人机航拍的ar远程施工监控方法及*** | |
CN110736447B (zh) | 一种集成式图像采集设备竖直方向水平位置检校方法 | |
CN106767895B (zh) | 利用标志点定向求解相机外方位元素的***和方法 | |
CN110986888A (zh) | 一种航空摄影一体化方法 | |
KR101346192B1 (ko) | 항공측량 전용 촬영이미지의 실시간 정사 보정용 항공측량시스템 | |
CN115760999A (zh) | 一种基于gis辅助的单目相机校准和目标地理位置提取方法 | |
CN109798875B (zh) | 一种在矿区视频***中实时标注采矿范围线的方法 | |
CN111127564B (zh) | 一种基于几何定位模型的视频影像校正方法 | |
KR101346206B1 (ko) | 정사영상의 정밀 처리를 위한 gps 기반 항공측량시스템 | |
CN117308915A (zh) | 一种测绘工程中特殊地形的测绘*** | |
KR100952136B1 (ko) | 경사영상의 외부표정요소 보정방법 및 이를 이용한 구조물의 3차원 정보를 측정하는 방법 | |
CN111561867A (zh) | 一种飞机表面形貌数字化测量方法 | |
Zheng et al. | A novel measurement method based on silhouette for chimney quasi-static deformation monitoring | |
WO2022126339A1 (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |