CN116499364A - 盘煤仪三维激光点云调畸变方法及系统 - Google Patents
盘煤仪三维激光点云调畸变方法及系统 Download PDFInfo
- Publication number
- CN116499364A CN116499364A CN202310786765.1A CN202310786765A CN116499364A CN 116499364 A CN116499364 A CN 116499364A CN 202310786765 A CN202310786765 A CN 202310786765A CN 116499364 A CN116499364 A CN 116499364A
- Authority
- CN
- China
- Prior art keywords
- data
- point cloud
- distortion
- laser point
- coal
- 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 47
- 239000003245 coal Substances 0.000 claims abstract description 99
- 238000012937 correction Methods 0.000 claims abstract description 60
- 239000011159 matrix material Substances 0.000 claims description 46
- 238000012545 processing Methods 0.000 claims description 26
- 230000011218 segmentation Effects 0.000 claims description 9
- 238000006243 chemical reaction Methods 0.000 claims description 5
- 238000010276 construction Methods 0.000 abstract description 5
- 238000005457 optimization Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000012935 Averaging Methods 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 238000009434 installation Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 238000013519 translation Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
-
- 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
- G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
- G01S17/86—Combinations of lidar systems with systems other than lidar, radar or sonar, e.g. with direction finders
-
- 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
- G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
- G01S17/88—Lidar systems specially adapted for specific applications
-
- 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/48—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00
- G01S7/497—Means for monitoring or calibrating
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/80—Analysis of captured images to determine intrinsic or extrinsic camera parameters, i.e. camera calibration
- G06T7/85—Stereo camera calibration
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Electromagnetism (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Theoretical Computer Science (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明涉及盘煤仪技术领域,尤其涉及一种盘煤仪三维激光点云调畸变方法及系统,所述方法包括以下步骤:采集煤堆的3D激光点云数据,并通过煤堆的3D激光点云数据获取畸变位置数据;获取参数标定好的摄像头采集激光标定板的真实位置数据;通过畸变位置数据和真实位置数据对比获取畸变校正函数;通过畸变校正函数对煤堆的3D激光点云数据进行校正,得到畸变校正后的煤堆的3D激光点云数据。本发明还公开了一种盘煤仪三维激光点云调畸变系统。本发明通过利用参数标定好的摄像头来获取真实位置数据,进而对畸变位置数据进行校正,避免了激光扫描进行三维构建中出现点云畸变影响盘煤结果的技术问题。
Description
技术领域
本发明涉及盘煤仪技术领域,尤其涉及一种盘煤仪三维激光点云调畸变方法及系统。
背景技术
目前,针对煤场的煤堆进行盘煤工作,其手段主要包括无人机盘煤或激光盘煤,无人机盘煤大多应用计算机视觉方式进行,目前由于受限于带宽及无人机本身算力所限,通常准确性不能得到很好保证。
而对于激光盘煤技术手段来说,现有技术是将激光扫描仪安装在棚顶上,通过多个激光扫描仪点云数据拼接配准完成整个煤场盘点。但是在激光扫描进行三维构建中会出现点云畸变,从而影响最后的盘煤结果。
上述内容仅用于辅助理解本发明的技术方案,并不代表承认上述内容是现有技术。
发明内容
本发明的主要目的在于提供一种盘煤仪三维激光点云调畸变方法及系统,旨在解决现有技术中激光扫描进行三维构建中出现点云畸变影响盘煤结果的技术问题。
为实现上述目的,本发明提供一种盘煤仪三维激光点云调畸变方法,所述方法包括以下步骤:
采集煤堆的3D激光点云数据,并通过所述煤堆的3D激光点云数据获取畸变位置数据;
获取参数标定好的摄像头采集激光标定板的真实位置数据;
通过所述畸变位置数据和所述真实位置数据对比获取畸变校正函数;
通过所述畸变校正函数对所述煤堆的3D激光点云数据进行校正,得到畸变校正后的煤堆的3D激光点云数据。
可选地,所述通过所述畸变位置数据和所述真实位置数据对比获取畸变校正函数的步骤,包括:
通过旋转矩阵将所述真实位置数据转换成激光雷达数据;
将所述激光雷达数据和所述畸变位置数据进行对比形成偏差矩阵;
通过对所述偏差矩阵进行处理获取畸变校正函数。
可选地,所述通过对所述偏差矩阵进行处理获取畸变校正函数的步骤,包括:
基于垂直点为中心进行半径递增对偏差矩阵中各个点阵偏差数据分别求解系数k1、k2和k3;
通过所述系数k1、k2和k3确定泰勒展开式,所述泰勒展开式为畸变校正函数。
可选地,所述获取参数标定好的摄像头采集激光标定板的真实位置数据的步骤之前,还包括:
通过相机标定板对摄像头进行标定,获取摄像头的内部参数和外部参数,并通过所述内部参数和所述外部参数对所述摄像头进行参数标定。
可选地,所述通过所述畸变校正函数对所述煤堆的3D激光点云数据进行校正,得到畸变校正后的煤堆的3D激光点云数据的步骤之后,还包括:
对所述畸变校正后的煤堆的3D激光点云数据进行配准处理,得到配准之后的煤堆的3D激光点云数据;
对所述配准之后的煤堆的3D激光点云数据进行分割处理,得到煤堆的体积参数。
此外,为实现上述目的,本发明还提供一种盘煤仪三维激光点云调畸变系统,所述系统包括:
数据采集模块,用于采集煤堆的3D激光点云数据,并通过所述煤堆的3D激光点云数据获取畸变位置数据;
参数标定模块,用于获取参数标定好的摄像头采集激光标定板的真实位置数据;
数据比较模块,用于通过所述畸变位置数据和所述真实位置数据对比获取畸变校正函数;
数据校正模块,用于通过所述畸变校正函数对所述煤堆的3D激光点云数据进行校正,得到畸变校正后的煤堆的3D激光点云数据。
可选地,所述数据比较模块,包括:
数据转换模块,用于通过旋转矩阵将所述真实位置数据转换成激光雷达数据;
矩阵形成模块,用于将所述激光雷达数据和所述畸变位置数据进行对比形成偏差矩阵;
函数形成模块,用于通过对所述偏差矩阵进行处理获取畸变校正函数。
可选地,所述函数形成模块,具体用于基于垂直点为中心进行半径递增对偏差矩阵中各个点阵偏差数据分别求解系数k1、k2和k3,通过所述系数k1、k2和k3确定泰勒展开式,所述泰勒展开式为畸变校正函数。
可选地,所述系统还包括:
参数获取模块,用于通过相机标定板对摄像头进行标定,获取摄像头的内部参数和外部参数,并通过所述内部参数和所述外部参数对所述摄像头进行参数标定。
可选地,所述系统还包括:
数据配准模块,用于对所述畸变校正后的煤堆的3D激光点云数据进行配准处理,得到配准之后的煤堆的3D激光点云数据;
数据分割模块,用于对所述配准之后的煤堆的3D激光点云数据进行分割处理,得到煤堆的体积参数。
本发明通过采集煤堆的3D激光点云数据,并通过所述煤堆的3D激光点云数据获取畸变位置数据,获取参数标定好的摄像头采集激光标定板的真实位置数据,通过所述畸变位置数据和所述真实位置数据对比获取畸变校正函数,通过所述畸变校正函数对所述煤堆的3D激光点云数据进行校正,得到畸变校正后的煤堆的3D激光点云数据,从而通过畸变校正后的煤堆的3D激光点云数据获取准确的盘煤结果。
附图说明
图1为本发明盘煤仪三维激光点云调畸变方法第一实施例的流程示意图;
图2为本发明盘煤仪三维激光点云调畸变方法第二实施例的流程示意图;
图3为本发明盘煤仪三维激光点云调畸变方法第三实施例的流程示意图;
图4为本发明盘煤仪三维激光点云调畸变系统第一实施例的功能模块示意图;
图5为本发明盘煤仪三维激光点云调畸变系统第二实施例的功能模块示意图;
图6为本发明盘煤仪三维激光点云调畸变系统第三实施例的功能模块示意图。
本发明目的的实现、功能特点及优点将结合实施例,参照附图做进一步说明。
具体实施方式
应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
参照图1,本发明盘煤仪三维激光点云调畸变方法第一实施例提供一种盘煤仪三维激光点云调畸变方法,所述方法包括:
S10:采集煤堆的3D激光点云数据,并通过所述煤堆的3D激光点云数据获取畸变位置数据;
需要说明的是,用于采集煤堆的煤堆的3D激光点云数据的设备可以是包含一套2d云台,一台测距200m的高灵敏的盘煤激光雷达,并通过云台和激光雷达逐次扫描,形成3d点云信息。
可理解的是,畸变位置数据可以通过激光点采集激光标定板中镂空和不镂空处的特征点获取煤堆的3D激光点云数据,再通过测距和激光与云台的角度换算求出xyz坐标。
S20:获取参数标定好的摄像头采集激光标定板的真实位置数据;
需要说明的是,摄像头参数标定可以是通过相机标定板进行标定,从而获得相机内部参数和外部参数。采集激光标定板的真实位置数据可以是通过采集激光标定板获取相机坐标数据,并通过手动变换获得激光雷达坐标数据坐标,进而得到近似真实世界xyz坐标。
S30:通过所述畸变位置数据和所述真实位置数据对比获取畸变校正函数;
在具体实现中,相机的xyz点云矩阵通过旋转矩阵坐标变换到基于激光雷达的原始坐标上,然后对特征点位置数据和相机点云的近似真实数据对比,形成偏差矩阵,通过对偏差矩阵的处理,基于垂直点为中心进行半径递增对各个点阵偏差数据分别求解k1,k2,k3的系数,进而获得泰勒展开式。后期盘煤仪每次扫描完成后通过泰勒展开式求解真实的xyz坐标。
S40:通过所述畸变校正函数对所述煤堆的3D激光点云数据进行校正,得到畸变校正后的煤堆的3D激光点云数据。
应理解的是,由于相机出厂装配无法百分百标准,相机的透镜面安装和感光cmos阵列面摆放不一定平行和中心点对齐,造成成像既有径向畸变又有切向畸变,可以通过采用泰勒级数展开式对煤堆的3D激光点云数据进行校正,得到畸变校正后的煤堆的3D激光点云数据。
本实施例通过利用参数标定好的摄像头来获取真实位置数据,进而对畸变位置数据进行校正,避免了激光扫描进行三维构建中出现点云畸变影响盘煤结果的技术问题。
参照图2,图2为本发明盘煤仪三维激光点云调畸变方法第二实施例的流程示意图,基于上述图1所示的实施例,提出本发明盘煤仪三维激光点云调畸变方法的第二实施例。
本实施例中,步骤S30包括:
S301:通过旋转矩阵将所述真实位置数据转换成激光雷达数据;
需要说明的是,旋转矩阵通过相机标定的外部参数获得。其中将真实位置数据转换成激光雷达数据函数为:
式中,,/>,/>为激光雷达数据坐标,/>,/>,/>为真实位置数据坐标,/>为旋转矩阵,T为平移向量,t为时间。
S302:将所述激光雷达数据和所述畸变位置数据进行对比形成偏差矩阵;
S303:基于垂直点为中心进行半径递增对偏差矩阵中各个点阵偏差数据分别求解系数k1、k2和k3,通过所述系数k1、k2和k3确定泰勒展开式,所述泰勒展开式为畸变校正函数。
需要说明的是,将所述激光雷达数据和所述畸变位置数据进行比较,可以使用几何关系或图像处理技术来计算其在水平和垂直方向上的偏差值,将计算得到的偏差值组织成一个偏差矩阵。偏差矩阵的每个元素代表相应特征点的偏差值。
应理解的是,基于垂直点为中心进行半径递增可以选择以垂直点为中心,从内到外逐渐增加半径的方式将特征点按照距离图像中心的半径进行分组。对于每个半径分组,可以使用最小二乘法或其他优化算法来拟合或优化偏差数据与理想模型之间的差异,拟合或优化过程会得到每个半径分组对应的系数k1、k2和k3。这些系数可以用于描述相应半径范围内的畸变情况。对于在每个半径分组中获得的系数,可以进行插值或平均操作,以获得连续范围内的系数值。得到系数后,可以使用它们来构建畸变校正函数,例如使用泰勒展开式。这样,校正函数可以在图像采集时应用于每个像素点,以校正畸变。其中,泰勒展开式为:
式中,为畸变校正后的x坐标,/>为畸变校正后的y坐标,x为畸变校正前的x坐标,y为畸变校正前y坐标,r为失真半径,/>,/>,/>为偏差系数。
参照图3,图3为本发明盘煤仪三维激光点云调畸变方法第三实施例的流程示意图,基于上述图2所示的实施例,提出本发明盘煤仪三维激光点云调畸变方法的第三实施例。
本实施例中,步骤S20之前,所述方法还包括:
S200:通过相机标定板对摄像头进行标定,获取摄像头的内部参数和外部参数,并通过所述内部参数和所述外部参数对所述摄像头进行参数标定。
可理解的是,相机的内部参数和外部参数用于进行图像三维场景重构,相机涉及四个坐标系:世界坐标系、相机坐标系、图像物理坐标系、图像像素坐标系。通过相机的内部参数和外部参数将世界坐标系最终转化为图像像素坐标系。
本实施例中,步骤S40之后,所述方法还包括:
S50:对所述畸变校正后的煤堆的3D激光点云数据进行配准处理,得到配准之后的煤堆的3D激光点云数据;
需要说明的是,对煤堆的3D激光点云数据进行配准处理可以采用ICP(IterativeClosest Point)算法进行配准。
在具体实现中,盘煤仪扫描会有重叠部分,而重叠部分会堆出特定高度的煤堆谷条,通过程序在点云中查找谷条,形成特征线,然后通过ICP算法对两片点云特征线进行线-线配准,配准达到一定可信度,将两片点云进行拼接。
需要说明的是,上述“谷条”指的是在盘煤仪扫描交汇的重叠区域中,堆积出的特定高度的煤堆形状。这种堆积的煤堆形状通常被称为“谷条”,因为它们呈现出类似于狭长的沟谷或条带状的外观。这些谷条的高度可以根据具体的需求和操作设置,以便在点云数据处理过程中能够作为特征进行检测和配准。通过在重叠区域堆积谷条,可以在后续的点云数据处理中使用这些特征来进行点云配准、对齐或拼接。这样可以提高点云数据的连续性和准确性,从而获得更完整和准确的煤堆表面模型或数据。
S60:对所述配准之后的煤堆的3D激光点云数据进行分割处理,得到煤堆的体积参数。
在具体实现中,通过对点云所有z轴归零,可以得到垂影面,再对垂影面两两点云间的xy距离,即可求出三角形面积,再对三点的z轴平均值相乘,得到近似的体积,所有的三角体体积积分最终得到煤堆的总体积。
参照图4,本发明盘煤仪三维激光点云调畸变系统第一实施例提供一种盘煤仪三维激光点云调畸变系统,所述系统包括:
数据采集模块10:采集煤堆的3D激光点云数据,并通过所述煤堆的3D激光点云数据获取畸变位置数据;
需要说明的是,用于采集煤堆的煤堆的3D激光点云数据的设备可以是包含一套2d云台,一台测距200m的高灵敏的盘煤激光雷达,并通过云台和激光雷达逐次扫描,形成3d点云信息。
可理解的是,畸变位置数据可以通过激光点采集激光标定板中镂空和不镂空处的特征点获取煤堆的3D激光点云数据,再通过测距和激光与云台的角度换算求出xyz坐标。
参数标定模块20:获取参数标定好的摄像头采集激光标定板的真实位置数据;
需要说明的是,摄像头参数标定可以是通过相机标定板进行标定,从而获得相机内部参数和外部参数。采集激光标定板的真实位置数据可以是通过采集激光标定板获取相机坐标数据,并通过手动变换获得激光雷达坐标数据坐标,进而得到近似真实世界xyz坐标。
数据比较模块30:通过所述畸变位置数据和所述真实位置数据对比获取畸变校正函数;
在具体实现中,相机的xyz点云矩阵通过旋转矩阵坐标变换到基于激光雷达的原始坐标上,然后对特征点位置数据和相机点云的近似真实数据对比,形成偏差矩阵,通过对偏差矩阵的处理,基于垂直点为中心进行半径递增对各个点阵偏差数据分别求解k1,k2,k3的系数,进而获得泰勒展开式。后期盘煤仪每次扫描完成后通过泰勒展开式求解真实的xyz坐标。
数据校正模块40:通过所述畸变校正函数对所述煤堆的3D激光点云数据进行校正,得到畸变校正后的煤堆的3D激光点云数据。
应理解的是,由于相机出厂装配无法百分百标准,相机的透镜面安装和感光cmos阵列面摆放不一定平行和中心点对齐,造成成像既有径向畸变又有切向畸变,可以通过采用泰勒级数展开式对煤堆的3D激光点云数据进行校正,得到畸变校正后的煤堆的3D激光点云数据。
本实施例通过利用参数标定好的摄像头来获取真实位置数据,进而对畸变位置数据进行校正,避免了激光扫描进行三维构建中出现点云畸变影响盘煤结果的技术问题。
参照图5,图5为本发明盘煤仪三维激光点云调畸变系统第二实施例的流程示意图,基于上述图1所示的实施例,提出本发明盘煤仪三维激光点云调畸变系统的第二实施例。
本实施例中,数据比较模块30包括:
数据转换模块301:通过旋转矩阵将所述真实位置数据转换成激光雷达数据;
需要说明的是,旋转矩阵通过相机标定的外部参数获得。其中将真实位置数据转换成激光雷达数据函数为:
式中,,/>,/>为激光雷达数据坐标,/>,/>,/>为真实位置数据坐标,/>为旋转矩阵,T为平移向量,t为时间。
矩阵形成模块302:将所述激光雷达数据和所述畸变位置数据进行对比形成偏差矩阵;
函数形成模块303:基于垂直点为中心进行半径递增对偏差矩阵中各个点阵偏差数据分别求解系数k1、k2和k3,通过所述系数k1、k2和k3确定泰勒展开式,所述泰勒展开式为畸变校正函数。
需要说明的是,将所述激光雷达数据和所述畸变位置数据进行比较,可以使用几何关系或图像处理技术来计算其在水平和垂直方向上的偏差值,将计算得到的偏差值组织成一个偏差矩阵。偏差矩阵的每个元素代表相应特征点的偏差值。
应理解的是,基于垂直点为中心进行半径递增可以选择以垂直点为中心,从内到外逐渐增加半径的方式将特征点按照距离图像中心的半径进行分组。对于每个半径分组,可以使用最小二乘法或其他优化算法来拟合或优化偏差数据与理想模型之间的差异,拟合或优化过程会得到每个半径分组对应的系数k1、k2和k3。这些系数可以用于描述相应半径范围内的畸变情况。对于在每个半径分组中获得的系数,可以进行插值或平均操作,以获得连续范围内的系数值。得到系数后,可以使用它们来构建畸变校正函数,例如使用泰勒展开式。这样,校正函数可以在图像采集时应用于每个像素点,以校正畸变。其中,泰勒展开式为:
式中,为畸变校正后的x坐标,/>为畸变校正后的y坐标,x为畸变校正前的x坐标,y为畸变校正前y坐标,r为失真半径,/>,/>,/>为偏差系数。
参照图6,图6为本发明盘煤仪三维激光点云调畸变系统第三实施例的流程示意图,基于上述图5所示的实施例,提出本发明盘煤仪三维激光点云调畸变系统的第三实施例。
本实施例中,参数标定模块20之前,所述系统还包括:
参数获取模块200:通过相机标定板对摄像头进行标定,获取摄像头的内部参数和外部参数,并通过所述内部参数和所述外部参数对所述摄像头进行参数标定。
可理解的是,相机的内部参数和外部参数用于进行图像三维场景重构,相机涉及四个坐标系:世界坐标系、相机坐标系、图像物理坐标系、图像像素坐标系。通过相机的内部参数和外部参数将世界坐标系最终转化为图像像素坐标系。
本实施例中,数据校正模块40之后,所述系统还包括:
数据配准模块50:对所述畸变校正后的煤堆的3D激光点云数据进行配准处理,得到配准之后的煤堆的3D激光点云数据;
需要说明的是,对煤堆的3D激光点云数据进行配准处理可以采用ICP(IterativeClosest Point)算法进行配准。
在具体实现中,盘煤仪扫描会有重叠部分,而重叠部分会堆出特定高度的煤堆谷条,通过程序在点云中查找谷条,形成特征线,然后通过ICP算法对两片点云特征线进行线-线配准,配准达到一定可信度,将两片点云进行拼接。
需要说明的是,上述“谷条”指的是在盘煤仪扫描交汇的重叠区域中,堆积出的特定高度的煤堆形状。这种堆积的煤堆形状通常被称为“谷条”,因为它们呈现出类似于狭长的沟谷或条带状的外观。这些谷条的高度可以根据具体的需求和操作设置,以便在点云数据处理过程中能够作为特征进行检测和配准。通过在重叠区域堆积谷条,可以在后续的点云数据处理中使用这些特征来进行点云配准、对齐或拼接。这样可以提高点云数据的连续性和准确性,从而获得更完整和准确的煤堆表面模型或数据。
数据分割模块60:对所述配准之后的煤堆的3D激光点云数据进行分割处理,得到煤堆的体积参数。
在具体实现中,通过对点云所有z轴归零,可以得到垂影面,再对垂影面两两点云间的xy距离,即可求出三角形面积,再对三点的z轴平均值相乘,得到近似的体积,所有的三角体体积积分最终得到煤堆的总体积。
需要说明的是,在本文中,术语“包括”、“包含”或者其任何其他变体 意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或 者系统不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还 包括为这种过程、方法、物品或者系统所固有的要素。在没有更多限制的情 况下,由语句“包括一个……”限定的要素,并不排除在包括该要素的过程、 方法、物品或者系统中还存在另外的相同要素。
上述本发明实施例序号仅仅为了描述,不代表实施例的优劣。
通过以上的实施方式的描述,本领域的技术人员可以清楚地了解到上述 实施例方法可借助软件加必需的通用硬件平台的方式来实现,当然也可以通 过硬件,但很多情况下前者是更佳的实施方式。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体 现出来,该计算机软件产品存储在一个存储介质(如ROM/RAM、磁碟、光 盘)中,包括若干指令用以使得一台终端设备(可以是手机,计算机,服务器,空调器,或者网络设备等)执行本发明各个实施例所述的方法。
以上仅为本发明的可选实施例,并非因此限制本发明的专利范围,凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换,或直接或间接运用在其他相关的技术领域,均同理包括在本发明的专利保护范围内。
Claims (10)
1.一种盘煤仪三维激光点云调畸变方法,其特征在于,所述方法包括以下步骤:
采集煤堆的3D激光点云数据,并通过所述煤堆的3D激光点云数据获取畸变位置数据;
获取参数标定好的摄像头采集激光标定板的真实位置数据;
通过所述畸变位置数据和所述真实位置数据对比获取畸变校正函数;
通过所述畸变校正函数对所述煤堆的3D激光点云数据进行校正,得到畸变校正后的煤堆的3D激光点云数据。
2.如权利要求1所述的方法,其特征在于,所述通过所述畸变位置数据和所述真实位置数据对比获取畸变校正函数的步骤,包括:
通过旋转矩阵将所述真实位置数据转换成激光雷达数据;
将所述激光雷达数据和所述畸变位置数据进行对比形成偏差矩阵;
通过对所述偏差矩阵进行处理获取畸变校正函数。
3.如权利要求2所述的方法,其特征在于,所述通过对所述偏差矩阵进行处理获取畸变校正函数的步骤,包括:
基于垂直点为中心进行半径递增对所述偏差矩阵中各个点阵偏差数据分别求解系数k1、k2和k3;
通过所述系数k1、k2和k3确定泰勒展开式,所述泰勒展开式为畸变校正函数。
4.如权利要求1所述的方法,其特征在于,所述获取参数标定好的摄像头采集激光标定板的真实位置数据的步骤之前,还包括:
通过相机标定板对摄像头进行标定,获取摄像头的内部参数和外部参数,并通过所述内部参数和所述外部参数对所述摄像头进行参数标定。
5.如权利要求1~4中任一项所述的方法,其特征在于,所述通过所述畸变校正函数对所述煤堆的3D激光点云数据进行校正,得到畸变校正后的煤堆的3D激光点云数据的步骤之后,还包括:
对所述畸变校正后的煤堆的3D激光点云数据进行配准处理,得到配准之后的煤堆的3D激光点云数据;
对所述配准之后的煤堆的3D激光点云数据进行分割处理,得到煤堆的体积参数。
6.一种盘煤仪三维激光点云调畸变系统,其特征在于,所述系统包括:
数据采集模块,用于采集煤堆的3D激光点云数据,并通过所述煤堆的3D激光点云数据获取畸变位置数据;
参数标定模块,用于获取参数标定好的摄像头采集激光标定板的真实位置数据;
数据比较模块,用于通过所述畸变位置数据和所述真实位置数据对比获取畸变校正函数;
数据校正模块,用于通过所述畸变校正函数对所述煤堆的3D激光点云数据进行校正,得到畸变校正后的煤堆的3D激光点云数据。
7.如权利要求6所述的系统,其特征在于,所述数据比较模块,包括:
数据转换模块,用于通过旋转矩阵将所述真实位置数据转换成激光雷达数据;
矩阵形成模块,用于将所述激光雷达数据和所述畸变位置数据进行对比形成偏差矩阵;
函数形成模块,用于通过对所述偏差矩阵进行处理获取畸变校正函数。
8.如权利要求7所述的系统,其特征在于,所述函数形成模块,具体用于基于垂直点为中心进行半径递增对偏差矩阵中各个点阵偏差数据分别求解系数k1、k2和k3,通过所述系数k1、k2和k3确定泰勒展开式,所述泰勒展开式为畸变校正函数。
9.如权利要求6所述的系统,其特征在于,所述系统还包括:
参数获取模块,用于通过相机标定板对摄像头进行标定,获取摄像头的内部参数和外部参数,并通过所述内部参数和所述外部参数对所述摄像头进行参数标定。
10.如权利要求6~9中任一项所述的系统,其特征在于,所述系统还包括:
数据配准模块,用于对所述畸变校正后的煤堆的3D激光点云数据进行配准处理,得到配准之后的煤堆的3D激光点云数据;
数据分割模块,用于对所述配准之后的煤堆的3D激光点云数据进行分割处理,得到煤堆的体积参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310786765.1A CN116499364B (zh) | 2023-06-30 | 2023-06-30 | 盘煤仪三维激光点云调畸变方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310786765.1A CN116499364B (zh) | 2023-06-30 | 2023-06-30 | 盘煤仪三维激光点云调畸变方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116499364A true CN116499364A (zh) | 2023-07-28 |
CN116499364B CN116499364B (zh) | 2023-09-12 |
Family
ID=87318743
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310786765.1A Active CN116499364B (zh) | 2023-06-30 | 2023-06-30 | 盘煤仪三维激光点云调畸变方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116499364B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117011309A (zh) * | 2023-09-28 | 2023-11-07 | 济宁港航梁山港有限公司 | 基于人工智能及深度数据的自动盘煤系统 |
Citations (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3645253B1 (ja) * | 2004-02-27 | 2005-05-11 | 中電技術コンサルタント株式会社 | 体積算定システム |
CN1687700A (zh) * | 2005-04-15 | 2005-10-26 | 东南大学 | 立体视觉盘煤方法 |
CN107133989A (zh) * | 2017-06-12 | 2017-09-05 | 中国科学院长春光学精密机械与物理研究所 | 一种三维扫描系统参数标定方法 |
CN108020826A (zh) * | 2017-10-26 | 2018-05-11 | 厦门大学 | 多线激光雷达与多路相机混合标定方法 |
KR20190086951A (ko) * | 2018-01-15 | 2019-07-24 | 주식회사 스트리스 | 레이저 관측 장비를 이용한 mms의 캘리브레이션을 위한 장치 및 방법 |
KR102038097B1 (ko) * | 2018-10-24 | 2019-10-29 | 한국중부발전(주) | 라이다를 갖는 드론을 이용한 저탄량 측정 시스템 |
CN111123280A (zh) * | 2019-12-31 | 2020-05-08 | 武汉万集信息技术有限公司 | 激光雷达的定位方法、装置、系统、电子设备及存储介质 |
CN112161568A (zh) * | 2020-09-28 | 2021-01-01 | 山东大学 | 基于多目普通摄像机的大型煤堆体积测量方法 |
CN114565661A (zh) * | 2022-01-20 | 2022-05-31 | 华能汕头海门发电有限责任公司 | 基于图像采集的盘煤系统 |
CN114581619A (zh) * | 2022-02-25 | 2022-06-03 | 浙江工业大学 | 一种基于三维定位和二维建图的煤仓建模方法 |
CN115018902A (zh) * | 2022-06-10 | 2022-09-06 | 北京瓦特曼智能科技有限公司 | 用于确定媒堆库存的方法及处理器 |
CN115063458A (zh) * | 2022-07-27 | 2022-09-16 | 武汉工程大学 | 一种基于三维激光点云的料堆体积计算方法 |
CN115205373A (zh) * | 2022-06-09 | 2022-10-18 | 劢微机器人(深圳)有限公司 | 基于雷达和相机融合的卸货方法、检测装置及存储介质 |
CN115507770A (zh) * | 2022-08-30 | 2022-12-23 | 华能国际电力江苏能源开发有限公司 | 一种煤场盘煤测控方法 |
CN115908539A (zh) * | 2022-11-15 | 2023-04-04 | 南京航空航天大学 | 一种目标体积自动测量方法和装置、存储介质 |
WO2023050188A1 (zh) * | 2021-09-29 | 2023-04-06 | 华能聊城热电有限公司 | 一种基于三维建模技术的数字煤场管理方法及系统 |
-
2023
- 2023-06-30 CN CN202310786765.1A patent/CN116499364B/zh active Active
Patent Citations (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3645253B1 (ja) * | 2004-02-27 | 2005-05-11 | 中電技術コンサルタント株式会社 | 体積算定システム |
CN1687700A (zh) * | 2005-04-15 | 2005-10-26 | 东南大学 | 立体视觉盘煤方法 |
CN107133989A (zh) * | 2017-06-12 | 2017-09-05 | 中国科学院长春光学精密机械与物理研究所 | 一种三维扫描系统参数标定方法 |
CN108020826A (zh) * | 2017-10-26 | 2018-05-11 | 厦门大学 | 多线激光雷达与多路相机混合标定方法 |
KR20190086951A (ko) * | 2018-01-15 | 2019-07-24 | 주식회사 스트리스 | 레이저 관측 장비를 이용한 mms의 캘리브레이션을 위한 장치 및 방법 |
KR102038097B1 (ko) * | 2018-10-24 | 2019-10-29 | 한국중부발전(주) | 라이다를 갖는 드론을 이용한 저탄량 측정 시스템 |
CN111123280A (zh) * | 2019-12-31 | 2020-05-08 | 武汉万集信息技术有限公司 | 激光雷达的定位方法、装置、系统、电子设备及存储介质 |
CN112161568A (zh) * | 2020-09-28 | 2021-01-01 | 山东大学 | 基于多目普通摄像机的大型煤堆体积测量方法 |
WO2023050188A1 (zh) * | 2021-09-29 | 2023-04-06 | 华能聊城热电有限公司 | 一种基于三维建模技术的数字煤场管理方法及系统 |
CN114565661A (zh) * | 2022-01-20 | 2022-05-31 | 华能汕头海门发电有限责任公司 | 基于图像采集的盘煤系统 |
CN114581619A (zh) * | 2022-02-25 | 2022-06-03 | 浙江工业大学 | 一种基于三维定位和二维建图的煤仓建模方法 |
CN115205373A (zh) * | 2022-06-09 | 2022-10-18 | 劢微机器人(深圳)有限公司 | 基于雷达和相机融合的卸货方法、检测装置及存储介质 |
CN115018902A (zh) * | 2022-06-10 | 2022-09-06 | 北京瓦特曼智能科技有限公司 | 用于确定媒堆库存的方法及处理器 |
CN115063458A (zh) * | 2022-07-27 | 2022-09-16 | 武汉工程大学 | 一种基于三维激光点云的料堆体积计算方法 |
CN115507770A (zh) * | 2022-08-30 | 2022-12-23 | 华能国际电力江苏能源开发有限公司 | 一种煤场盘煤测控方法 |
CN115908539A (zh) * | 2022-11-15 | 2023-04-04 | 南京航空航天大学 | 一种目标体积自动测量方法和装置、存储介质 |
Non-Patent Citations (2)
Title |
---|
孔德明;张娜;黄紫双;陈晓玉;沈阅;: "基于激光雷达探测技术的列车车厢载货体积测量方法研究", 燕山大学学报, no. 02 * |
张庆军;何平安;刘震亚;宋家军;: "基于图像校正的测量移动堆料体积的标定", 激光杂志, no. 05 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117011309A (zh) * | 2023-09-28 | 2023-11-07 | 济宁港航梁山港有限公司 | 基于人工智能及深度数据的自动盘煤系统 |
CN117011309B (zh) * | 2023-09-28 | 2023-12-26 | 济宁港航梁山港有限公司 | 基于人工智能及深度数据的自动盘煤系统 |
Also Published As
Publication number | Publication date |
---|---|
CN116499364B (zh) | 2023-09-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107316325B (zh) | 一种基于图像配准的机载激光点云与影像配准融合方法 | |
CN112907676B (zh) | 传感器的标定方法、装置、系统、车辆、设备及存储介质 | |
CN116499364B (zh) | 盘煤仪三维激光点云调畸变方法及系统 | |
CN104835159A (zh) | 用于连续变焦距光学成像系统的数字图像校正方法 | |
CN109727278B (zh) | 一种机载LiDAR点云数据与航空影像的自动配准方法 | |
EP1589482A2 (en) | Three-dimensional image measuring apparatus and method | |
WO2021098448A1 (zh) | 传感器标定方法及装置、存储介质、标定系统和程序产品 | |
JP3850541B2 (ja) | 高度計測装置 | |
CN108629756B (zh) | 一种Kinectv2深度图像无效点修复方法 | |
CN111383279A (zh) | 外参标定方法、装置及电子设备 | |
US8577139B2 (en) | Method of orthoimage color correction using multiple aerial images | |
CN111260539B (zh) | 鱼眼图目标识别方法及其系统 | |
CN111486864B (zh) | 基于立体正八边结构的多源传感器联合标定方法 | |
CN115661262A (zh) | 内外参标定方法、装置及电子设备 | |
CN112927306A (zh) | 拍摄装置的标定方法、装置及终端设备 | |
CN113205563A (zh) | 一种自动驾驶传感器联合标定靶和标定方法 | |
CN111913169A (zh) | 激光雷达内参、点云数据的修正方法、设备及存储介质 | |
CN108917722B (zh) | 植被覆盖度计算方法及装置 | |
CN113534110B (zh) | 一种多激光雷达系统静态标定方法 | |
CN113624142B (zh) | 一种飞机装配紧固件凹凸量检测方法 | |
CN110631555A (zh) | 基于二阶多项式无控制点区域网平差的历史影像正射纠方法 | |
CN107941241B (zh) | 一种用于航空摄影测量质量评价的分辨率板及其使用方法 | |
CN111696141A (zh) | 一种三维全景扫描采集方法、设备及存储设备 | |
CN112419427A (zh) | 用于提高飞行时间相机精度的方法 | |
Dlesk et al. | Structure from motion processing of analogue images captured by Rollei metric camera, digitized with various scanning resolution |
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 | ||
PE01 | Entry into force of the registration of the contract for pledge of patent right | ||
PE01 | Entry into force of the registration of the contract for pledge of patent right |
Denomination of invention: Method and System for 3D Laser Point Cloud Distortion Adjustment of Coal Pan Instrument Effective date of registration: 20231226 Granted publication date: 20230912 Pledgee: Levin commercial bank Limited by Share Ltd. Ji'nan branch Pledgor: JINAN ZUOWEI TECHNOLOGY CO.,LTD. Registration number: Y2023980074468 |