CN111024053A - 一种凸球面镜反射的激光扫描数据获取与解算方法 - Google Patents

一种凸球面镜反射的激光扫描数据获取与解算方法 Download PDF

Info

Publication number
CN111024053A
CN111024053A CN201911152397.5A CN201911152397A CN111024053A CN 111024053 A CN111024053 A CN 111024053A CN 201911152397 A CN201911152397 A CN 201911152397A CN 111024053 A CN111024053 A CN 111024053A
Authority
CN
China
Prior art keywords
point cloud
laser
convex spherical
spherical mirror
point
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
Application number
CN201911152397.5A
Other languages
English (en)
Other versions
CN111024053B (zh
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.)
Wuhan University WHU
Original Assignee
Wuhan University WHU
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 Wuhan University WHU filed Critical Wuhan University WHU
Priority to CN201911152397.5A priority Critical patent/CN111024053B/zh
Publication of CN111024053A publication Critical patent/CN111024053A/zh
Application granted granted Critical
Publication of CN111024053B publication Critical patent/CN111024053B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C15/00Surveying instruments or accessories not provided for in groups G01C1/00 - G01C13/00
    • G01C15/002Active optical surveying means
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/24Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures

Landscapes

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

Abstract

本发明属于测绘技术领域,公开了一种凸球面镜反射的激光扫描数据获取与解算方法,将凸球面镜放置在作业空间的死角处,通过三维激光扫描仪对物体进行扫描,获得激光点云扫描数据,将激光点云扫描数据分离为直视表面区域点云、镜面定位圈区域点云、镜面反射区域点云,根据镜面定位圈区域点云获得旋转矩阵,根据旋转矩阵、镜面反射区域点云,解算凸球面镜反射激光点云中激光像点对应的物体上的物点点云坐标;将直视表面区域点云与凸球面镜反射激光点云中激光像点对应的物体上的物点点云坐标进行配准合成,得到物体的完整点云。本发明解决了现有技术中扫描对象非直视区域出现点云数据缺失和空洞的问题,能够得到物体的完整点云。

Description

一种凸球面镜反射的激光扫描数据获取与解算方法
技术领域
本发明涉及测绘技术领域,尤其涉及一种凸球面镜反射的激光扫描数据获取与解算方法。
背景技术
由于国家对文化遗产保护的重视,越来越多的高新技术投入其中,包括利用激光扫描技术进行三维重建工作。激光扫描技术是高精度三维重建的重要手段,能够获取精度优于2mm的对象表面密集点云,在精密测量、虚拟地理环境、文物保护及虚拟现实等方面应用广泛。扫描对象的完整性是高精度三维重建的基本要求,为了保证对象数据尽可能完整,目前常用方法是:多站扫描和多种扫描仪联合扫描。然而对结构复杂且不可移动的对象如石窟寺、古建筑、雕塑等进行扫描时,往往因为作业空间受限、遮挡严重等问题造成扫描死角,在扫描对象非直视区域出现点云数据缺失和空洞等问题。
发明内容
本申请实施例通过提供一种凸球面镜反射的激光扫描数据获取与解算方法,解决了现有技术中由于作业空间受限、遮挡严重等问题造成扫描死角进而使扫描对象非直视区域出现点云数据缺失和空洞的问题。
本申请实施例提供一种凸球面镜反射的激光扫描数据获取与解算方法,包括以下步骤:
S1、将凸球面镜放置在作业空间的死角处;
S2、将三维激光扫描仪放置在与凸球面镜保持通视的位置,通过三维激光扫描仪对物体进行扫描,获得激光点云扫描数据;
S3、将所述激光点云扫描数据分离为直视表面区域点云、镜面定位圈区域点云、镜面反射区域点云;
S4、根据所述镜面定位圈区域点云,解算凸球面镜的镜面位置与姿态,获得旋转矩阵;
S5、根据所述旋转矩阵、所述镜面反射区域点云,解算凸球面镜反射激光点云中激光像点对应的物体上的物点点云坐标;
S6、将所述直视表面区域点云与解算得到的所述凸球面镜反射激光点云中激光像点对应的物体上的物点点云坐标进行配准合成,得到物体的完整点云。
优选的,所述步骤S1中,将凸球面镜固定在一个可多角度旋转移动的支架上。
优选的,所述步骤S2中,通过三维激光扫描仪对物体进行扫描时,调整凸球面镜的角度进行多次扫描。
优选的,所述步骤S4中,通过所述镜面定位圈区域点云解算出旋转矩阵,通过旋转平移使镜面的最大圆平面与三维激光扫描仪的坐标系平面平行。
优选的,所述步骤S5具体为:
建立凸球面镜反射激光点云中激光像点与原点的空间直线方程;
联立凸球面镜反射激光点云中激光像点与原点的空间直线方程、三维空间下的球面方程,获得空间直线与凸球面镜的交点坐标;
根据三维空间下的球面函数,获得空间直线与凸球面镜的交点在凸球面上的切平面方程;
根据切平面方程获得凸球面镜反射激光点云中激光像点关于切平面的对称点,所述对称点的坐标为凸球面镜反射激光点云中激光像点对应的物体上的物点坐标。
优选的,所述凸球面镜反射激光点云中激光像点与原点的空间直线方程为:
Figure BDA0002283895670000021
其中,(x1,y1,z1)是坐标原点(0,0,0),(x2,y2,z2)是凸球面镜反射激光点云中激光像点,t是参数;
所述三维空间下的球面方程为:
(X-a)2+(Y-b)2+(Z-c)2=R2 (2)
其中,(a,b,c)是凸球面镜的球心坐标,R是凸球面镜的半径;
联立(1)(2),求得关于参数t的函数方程:
F(t)=(x2t-a)2+(y2t-b)2+(z2t-c)2-R2 (3)
将(x2,y2,z2)以此代入(3)中,求得t值,得到空间直线与凸球面镜的交点(x,y,z);
所述三维空间下的球面函数为:
f(X,Y,Z)=(X-a)2+(Y-b)2+(Z-c)2-R2 (4)
所述空间直线与凸球面镜的交点在凸球面上的切平面为:
fx′(X-x)+fy′(Y-y)+fz′(Z-z)=0
fx′X+fy′Y+fz′Z-(fx′x+fy′y+fz′z)=0 (5)
所述凸球面镜反射激光点云中激光像点关于切平面的对称点为:
Figure BDA0002283895670000031
Figure BDA0002283895670000032
Figure BDA0002283895670000041
得到的(x3,y3,z3)为所述凸球面镜反射激光点云中激光像点对应的物体上的物点坐标。
本申请实施例中提供的一个或多个技术方案,至少具有如下技术效果或优点:
在本申请实施例中,将凸球面镜放置在作业空间的死角处,通过三维激光扫描仪对物体进行扫描,获得激光点云扫描数据,将激光点云扫描数据分离为直视表面区域点云、镜面定位圈区域点云、镜面反射区域点云,根据镜面定位圈区域点云获得旋转矩阵,根据旋转矩阵、镜面反射区域点云,解算凸球面镜反射激光点云中激光像点对应的物体上的物点点云坐标;将直视表面区域点云与凸球面镜反射激光点云中激光像点对应的物体上的物点点云坐标进行配准合成,得到物体的完整点云。本发明能够实现复杂情况下扫描物体的可能性,得到物体的完整点云。
附图说明
为了更清楚地说明本实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本发明的一个实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的一种凸球面镜反射的激光扫描数据获取与解算方法的流程图;
图2为本发明实施例提供的一种凸球面镜反射的激光扫描数据获取与解算方法中采用的球面镜的结构示意图;
图3为本发明实施例提供的一种凸球面镜反射的激光扫描数据获取与解算方法的实验示意图。
具体实施方式
为了更好的理解上述技术方案,下面将结合说明书附图以及具体的实施方式对上述技术方案进行详细的说明。
本实施例提供了一种凸球面镜反射的激光扫描数据获取与解算方法,参照图1-3,包括以下步骤:
S1、放置凸球面镜,固定支架,将凸球面镜放在作业空间的死角处。
具体的,凸球面镜固定在一个可以多角度旋转移动的支架上,凸球面镜的构成:标准球冠、漫反射圈、支架,其中,凸球面镜的半径R、球冠高等几何参数通过预先检校得到。
S2、三维激光扫描仪获取点云数据,将三维激光扫描仪摆放在可以与凸球面镜保持通视的地方,然后对物体(例如佛像)进行扫描,获取点云数据。
具体的,将凸球面镜放在物体的死角处,然后将三维激光扫描仪放置在三脚架上,调整合适的角度,保持镜子与扫描仪通视,调整扫描仪仪器参数,对物体进行扫描,调整镜子角度,进行多次扫描,尽量保证获得的数据完整。
S3、激光扫描点云分离,依据镜面范围,将激光扫描点云分离成直视表面区域点云、镜面定位圈区域点云、镜面反射区域点云。
具体的,激光扫描点云分离,依据镜面范围,分离成:直视表面区域点云、镜面定位圈区域点云、镜面反射区域点云三个部分,将扫描的完整的点云数据导入软件,然后将三部分分离,分别进行保存。
S4、解算镜面位置与姿态,根据镜面定位圈区域点云解算出旋转矩阵对镜面进行定位。
具体的,所用的凸球面镜的镜面是球形单面反射镜,镜面是标准球的球冠,球的半径R已知,步骤S4中对凸球面镜定位,通过镜面定位圈区域点云解算球面镜位置的旋转矩阵,通过旋转平移将镜面的最大圆平面与扫描仪坐标系xy平面平行。
S5、解算凸球面镜反射激光点云中激光像点对应的物体上的物点坐标。
具体的,在扫描仪坐标系中建立球面镜反射激光点云中激光像点与原点的空间直线方程,用两点式表示:
Figure BDA0002283895670000061
其中,(x1,y1,z1)是坐标原点(0,0,0),(x2,y2,z2)是凸球面镜反射激光点云中激光像点,t是参数。
三维空间下任意球面公式:
(X-a)2+(Y-b)2+(Z-c)2=R2 (2)
联立(1)(2),求得关于参数t的函数方程:
F(t)=(x2t-a)2+(y2t-b)2+(z2t-c)2-R2 (3)
将(x2,y2,z2)以此代入(3)中,求得t值,那么可以得到空间直线与凸球面镜的交点(x,y,z)。
三维空间下球面函数为:
f(X,Y,Z)=(X-a)2+(Y-b)2+(Z-c)2-R2 (4)
空间直线与凸球面镜的交点在凸球面上的切平面为:
fx′(X-x)+fy′(Y-y)+fz′(Z-z)=0
fx′X+fy′Y+fz′Z-(fx′x+fy′y+fz′z)=0 (5)
凸球面镜反射激光点云中激光像点关于切平面的对称点为:
Figure BDA0002283895670000062
Figure BDA0002283895670000071
Figure BDA0002283895670000072
得到的(x3,y3,z3)是凸球面镜反射激光点云中激光像点对应的物体上的物点坐标。
S6、直视表面区域点云与解算得到的凸球面镜反射激光点云中激光像点对应的物体上的物点点云做配准合成完整物体。
具体的,判断并选择两部分特征点,选取多对,进行配准合成。
本发明实施例提供的一种凸球面镜反射的激光扫描数据获取与解算方法至少包括如下技术效果:
1、本发明利用球面镜解算点云的方法和球面镜定位的方法,实现了凸球面镜反射的激光扫描数据获取与解算,能够准确地解算球面镜反射激光点云中激光像点对应的物体上的物点坐标。
2、本发明可实现凸球面镜反射激光点云中激光像点对应的物体上的物点坐标的解算,得到物体的完整点云,对后续的三维建模有极大的帮助,实现复杂情况下扫描物体的可能性,有助于实现文化保护的完整性。
最后所应说明的是,以上具体实施方式仅用以说明本发明的技术方案而非限制,尽管参照实例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。

Claims (6)

1.一种凸球面镜反射的激光扫描数据获取与解算方法,其特征在于,包括以下步骤:
S1、将凸球面镜放置在作业空间的死角处;
S2、将三维激光扫描仪放置在与凸球面镜保持通视的位置,通过三维激光扫描仪对物体进行扫描,获得激光点云扫描数据;
S3、将所述激光点云扫描数据分离为直视表面区域点云、镜面定位圈区域点云、镜面反射区域点云;
S4、根据所述镜面定位圈区域点云,解算凸球面镜的镜面位置与姿态,获得旋转矩阵;
S5、根据所述旋转矩阵、所述镜面反射区域点云,解算凸球面镜反射激光点云中激光像点对应的物体上的物点点云坐标;
S6、将所述直视表面区域点云与解算得到的所述凸球面镜反射激光点云中激光像点对应的物体上的物点点云坐标进行配准合成,得到物体的完整点云。
2.根据权利要求1所述的凸球面镜反射的激光扫描数据获取与解算方法,其特征在于,所述步骤S1中,将凸球面镜固定在一个可多角度旋转移动的支架上。
3.根据权利要求1所述的凸球面镜反射的激光扫描数据获取与解算方法,其特征在于,所述步骤S2中,通过三维激光扫描仪对物体进行扫描时,调整凸球面镜的角度进行多次扫描。
4.根据权利要求1所述的凸球面镜反射的激光扫描数据获取与解算方法,其特征在于,所述步骤S4中,通过所述镜面定位圈区域点云解算出旋转矩阵,通过旋转平移使镜面的最大圆平面与三维激光扫描仪的坐标系平面平行。
5.根据权利要求1所述的凸球面镜反射的激光扫描数据获取与解算方法,其特征在于,所述步骤S5具体为:
建立凸球面镜反射激光点云中激光像点与原点的空间直线方程;
联立凸球面镜反射激光点云中激光像点与原点的空间直线方程、三维空间下的球面方程,获得空间直线与凸球面镜的交点坐标;
根据三维空间下的球面函数,获得空间直线与凸球面镜的交点在凸球面上的切平面方程;
根据切平面方程获得凸球面镜反射激光点云中激光像点关于切平面的对称点,所述对称点的坐标为凸球面镜反射激光点云中激光像点对应的物体上的物点坐标。
6.根据权利要求5所述的凸球面镜反射的激光扫描数据获取与解算方法,其特征在于,所述凸球面镜反射激光点云中激光像点与原点的空间直线方程为:
Figure FDA0002283895660000021
其中,(x1,y1,z1)是坐标原点(0,0,0),(x2,y2,z2)是凸球面镜反射激光点云中激光像点,t是参数;
所述三维空间下的球面方程为:
(X-a)2+(Y-b)2+(Z-c)2=R2 (2)
其中,(a,b,c)是凸球面镜的球心坐标,R是凸球面镜的半径;
联立(1)(2),求得关于参数t的函数方程:
F(t)=(x2t-a)2+(y2t-b)2+(z2t-c)2-R2 (3)
将(x2,y2,z2)以此代入(3)中,求得t值,得到空间直线与凸球面镜的交点(x,y,z);
所述三维空间下的球面函数为:
f(X,Y,Z)=(X-a)2+(Y-b)2+(Z-c)2-R2 (4)
所述空间直线与凸球面镜的交点在凸球面上的切平面为:
fx′(X-x)+fy′(Y-y)+fz′(Z-z)=0
fx′X+fy′Y+fz′Z-(fx′x+fy′y+fz′z)=0 (5)
所述凸球面镜反射激光点云中激光像点关于切平面的对称点为:
Figure FDA0002283895660000031
Figure FDA0002283895660000032
Figure FDA0002283895660000033
得到的(x3,y3,z3)为所述凸球面镜反射激光点云中激光像点对应的物体上的物点坐标。
CN201911152397.5A 2019-11-22 2019-11-22 一种凸球面镜反射的激光扫描数据获取与解算方法 Active CN111024053B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911152397.5A CN111024053B (zh) 2019-11-22 2019-11-22 一种凸球面镜反射的激光扫描数据获取与解算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911152397.5A CN111024053B (zh) 2019-11-22 2019-11-22 一种凸球面镜反射的激光扫描数据获取与解算方法

Publications (2)

Publication Number Publication Date
CN111024053A true CN111024053A (zh) 2020-04-17
CN111024053B CN111024053B (zh) 2021-01-01

Family

ID=70206721

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911152397.5A Active CN111024053B (zh) 2019-11-22 2019-11-22 一种凸球面镜反射的激光扫描数据获取与解算方法

Country Status (1)

Country Link
CN (1) CN111024053B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103234480A (zh) * 2013-04-16 2013-08-07 北京理工大学 一种环形凸非球面的快速面形检测方法
CN107830812A (zh) * 2017-09-14 2018-03-23 同济大学 一种适用于隧道内定位和形变分析的激光反射片实现方法
CN108051796A (zh) * 2017-12-29 2018-05-18 合肥嘉东光学股份有限公司 一种基于tof的小型化同轴式激光雷达系统
US10094916B1 (en) * 2017-06-09 2018-10-09 Waymo Llc LIDAR optics alignment systems and methods
CN108692656A (zh) * 2017-04-07 2018-10-23 中国人民解放军信息工程大学 一种激光扫描数据获取方法及装置

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103234480A (zh) * 2013-04-16 2013-08-07 北京理工大学 一种环形凸非球面的快速面形检测方法
CN108692656A (zh) * 2017-04-07 2018-10-23 中国人民解放军信息工程大学 一种激光扫描数据获取方法及装置
US10094916B1 (en) * 2017-06-09 2018-10-09 Waymo Llc LIDAR optics alignment systems and methods
CN107830812A (zh) * 2017-09-14 2018-03-23 同济大学 一种适用于隧道内定位和形变分析的激光反射片实现方法
CN108051796A (zh) * 2017-12-29 2018-05-18 合肥嘉东光学股份有限公司 一种基于tof的小型化同轴式激光雷达系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
张帆: "基于激光扫描的隧道管片结构变形数据处理及可视化分析方法", 《现代隧道技术》 *

Also Published As

Publication number Publication date
CN111024053B (zh) 2021-01-01

Similar Documents

Publication Publication Date Title
El-Hakim et al. Comparative evaluation of the performance of passive and active 3D vision systems
CN107121062A (zh) 一种机器人三维扫描装置及方法
Sanz‐Ablanedo et al. Reducing systematic dome errors in digital elevation models through better UAV flight design
CN108805976B (zh) 一种三维扫描系统及方法
CN110793464B (zh) 大视场条纹投影视觉三维测量系统和方法
CN111811395A (zh) 基于单目视觉的平面位姿动态测量方法
CN104237868A (zh) 一种多功能实用型激光雷达扫描标靶
CN110360930A (zh) 一种激光位移法线传感器及其测量方法
CN110954017B (zh) 任意曲面镜反射的激光扫描数据获取与解算方法
Cai et al. Calibration method for the rotating axis in panoramic 3D shape measurement based on a turntable
CN114037751A (zh) 结构光系统的标定参数优化方法、存储介质、计算机设备
Amiri Parian et al. Integrated laser scanner and intensity image calibration and accuracy assessment
CN111024053B (zh) 一种凸球面镜反射的激光扫描数据获取与解算方法
CN2650064Y (zh) 一种激光扫描获取近距离物体表面三维数据的测量装置
CN115046498B (zh) 一种单目旋转结构光三维测量系统标定方法
CN105004283A (zh) 一种适用于牙颌模型三维测量的定位系统及其标定方法
CN113432558B (zh) 一种基于激光的不规则物体表面积的测量装置及方法
TW565736B (en) Method for determining the optical parameters of a camera
CN112070885B (zh) 一种基于摄影测量光度立体模型的三维重建方法和系统
MacDonald et al. Determining the coordinates of lamps in an illumination dome
WO2022078438A1 (zh) 一种室内3d信息采集设备
EP1108245B1 (en) Characterization of objects using laser supplemented photogrammetry
CN210154537U (zh) 一种基于数字摄影的金属结构变形测量装置
KR20220107893A (ko) 실제대상물과 가상대상물 사이의 상호작용을 지원하는 장치 및 이를 운영하는 방법
MacDonald et al. 3D reconstruction in an illumination dome

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