CN109738440A - 一种基于智能手机的岩体结构面产状非接触测量方法 - Google Patents
一种基于智能手机的岩体结构面产状非接触测量方法 Download PDFInfo
- Publication number
- CN109738440A CN109738440A CN201910005679.6A CN201910005679A CN109738440A CN 109738440 A CN109738440 A CN 109738440A CN 201910005679 A CN201910005679 A CN 201910005679A CN 109738440 A CN109738440 A CN 109738440A
- Authority
- CN
- China
- Prior art keywords
- rock mass
- relational expression
- discontinuity
- target
- smart phone
- 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
Landscapes
- Length Measuring Devices By Optical Means (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了一种基于智能手机的岩体结构面产状非接触测量方法,包括:在岩体露头或岩体开挖面上布置一个标靶,使用手机内置电子罗盘测量其产状;利用智能手机对待测量的岩体进行人工移动拍照;将获取的手机照片导入专业摄影测量软件生成三维岩体点云模型;基于标靶产状进行标定,得到相机坐标系与世界坐标系的转换矩阵;采用RANSAC算法提取岩体结构面,获取相机坐标系下结构面单位法向矢量;采用转换矩阵得到世界坐标系下结构面的单位法向矢量,从而计算得到结构面的产状。本发明方法适用于天然岩体露头及边坡、地下隧洞开挖面的岩体结构面测量,其优点是:现场操作简便易行,对复杂场地的适应性强。
Description
技术领域
本发明属于工程地质测量技术领域,具体涉及一种基于智能手机的岩体结构面产状非接触测量方法。
背景技术
岩体工程的稳定性主要受结构面控制,正确测量岩体的结构面产状对于建立岩体工程地质模型、开展岩体工程稳定性评价与指导支护设计具有重要意义。
传统岩体结构面产状的测量主要有三类方法。
第一类方法是接触式人工测量,这种方法是利用地质罗盘近距离直接测量结构面的产状。由于岩体露头范围大,特别是对于高位岩体,测量人员往往很难到达测量位置,另一方面,对于一些危岩体,对人身危险性较大,不宜近距离开展接触测量。此外,人工接触测量耗时长,效率低,受主观影响大。
第二类方法是非接触式摄影测量技术,采用相机按照一定的要求对被测对象连续拍摄图像,同时使用全站仪或GPS测量一定数量的控制点(3个以上),通过控制点坐标来标定摄影测量得到的点云坐标,测量过程中需要使用照相机和全站仪等多个仪器配合,从而导致现场测量操作复杂,效率低下。
第三类方法是使用三维激光扫描技术,利用激光扫描仪获取测量岩体的三维点云数据,由于测量得到的点云数据是相对坐标,该方法同样需要使用全站仪或GPS测量一定的控制点,再实现点云数据从相对坐标到世界坐标的转换。由于激光扫描仪造价比较高,仪器笨重,还需要配合全站仪使用,从而导致现场测量操作复杂,实施困难。
结合上述三类方法存在的不足之处,寻求一种简易、轻便的岩体结构面产状测量方法就显得很有必要。
发明内容
本发明的目的在于提供一种基于智能手机的岩体结构面产状非接触测量方法,能够测量难以近距离接触的结构面产状,克服了传统方法仪器笨重、需要多个仪器配合、高位岩体接触困难的弊端,为正确认识地质构造提供依据。
本发明的原理在于:在边坡岩体露头面或开挖面上放置一个矩形标靶,使得标靶一条边保持水平,采用手机进行非接触移动拍摄,将获取的照片导入摄影测量软件中建立岩体三维点云模型;采用RANSAC算法提取岩体结构面,获取相机坐标系下结构面单位法向矢量;同时利用标靶产状得到相机坐标系与世界坐标系的转换矩阵;基于转换矩阵得到世界坐标系下结构面的单位法向矢量,从而计算得到结构面的产状。
一种边坡岩体结构面产状测量方法,包括如下步骤:
1)在待测量岩体露头或开挖面上布设一个标靶,使得标靶的一条边保持水平,利用手机内置的电子罗盘测量标靶面的产状,记录其倾向和倾角;
2)对待测量岩体进行拍照,并确保两幅照片之间达到一定的重叠率,且拍摄的照片中含有标靶;
3)将拍摄获取的照片导入专业的摄影测量软件,生成基于相机坐标系的三维岩体点云模型;
4)在三维岩体点云模型中找到标靶,提取标靶面三个角点A,B,C的坐标,运用关系式1计算向量AB与向量AC的单位矢量ab和ac,运用关系式2计算标靶面单位法向矢量ad;
关系式1为:
关系式2为:
ad=ab×ac
关系式1和关系式2中,xA、yA、zA为点A的坐标,xB、yB、zB为点B的坐标,xC、yC、zC为点C的坐标;
5)根据罗盘测量得到的标靶面产状,运用关系式3计算世界坐标系三个轴的单位矢量ox、oy、oz,其中ox、oy轴分别和正东、正北方向保持一致,oz轴竖直向上;
关系式3为:
关系式3中,α0为标靶面倾向,β0为标靶面倾角;
6)从相机坐标系下的三维岩体点云模型中提取岩体结构面,得到结构面的单位法向矢量,运用关系式4将其转换为世界坐标系下的单位法向矢量;
关系式4为:
关系式4中,(nx,ny,nz)为结构面在世界坐标系下的单位法向矢量,为结构面在相机坐标系下的单位法向矢量;
7)采用关系式5计算岩体结构面的产状
关系式5为:
关系式5中,α为结构面倾向,β为结构面倾角;
8)重复步骤6)和7),提取岩体露头或岩体开挖坡面出现的所有结构面的产状信息。
进一步的,所述步骤1)中的标靶为矩形。
进一步的,所述步骤2)中两幅照片之间的重叠率大于60%。
进一步的,所述步骤2)中利用手机对待测量岩体进行拍照。
进一步的,所述步骤6)中采用RANSAC算法从相机坐标系下的三维岩体点云模型中提取岩体结构面。
与现有技术相比,本发明方法的优点是:
1.本发明通过布设一个矩形标靶对测量区域的正北方向进行标定,基于坐标转换得到相机坐标与世界坐标的转换矩阵,不需要其它仪器的配合来实现点云坐标的标定。
2.本发明所提供的测量方法在整个野外操作过程中,仅需要一部手机即可完成全部外业测量,操作简单易行,不需要使用笨重的仪器,减轻了负重,提高了效率。
3.本发明无需获取所拍照片的GPS信息,通过标靶面的产状即可推算出坐标转换关系,对于无法接受GPS信号或GPS信号不良的地方,本发明具有较强的适用性。
附图说明
图1为本发明实施例的流程图;
图2为本发明实施例的现场照片图;
图3为本发明实施例的结构面提取图,其中A,B,C为标靶上三个角点;
图4为本发明实施例从标靶面坐标系转换到世界坐标系的示意图。
具体实施方式
下面结合附图及具体实施例,对本发明进行进一步说明。
如图1所示一种基于智能手机的岩体结构面产状非接触测量方法,包括如下步骤:
1)在待测量岩体露头或开挖面上布设一个矩形标靶,使得标靶的一条边保持水平,使用手机内置电子罗盘测量标靶面的产状,记录其倾向和倾角;
2)利用智能手机对待测量岩体进行拍照,确保两幅照片之间有60%以上的重叠率,且拍摄的照片中含有标靶;
3)将获取的多幅手机照片导入专业的摄影测量软件,生成基于相机坐标系的三维岩体点云模型;
4)在三维岩体模型中找到标靶,提取标靶面三个角点A,B,C的坐标,运用关系式1计算向量AB与向量AC的单位矢量ab和ac,运用关系式2计算标靶面单位法向矢量ad。
关系式1为:
关系式2为:
ad=ab×ac
关系式1和关系式2中,xA、yA、zA为点A的坐标,xB、yB、zB为点B的坐标,xC、yC、zC为点C的坐标;
5)根据罗盘测量得到的标靶面产状,运用关系式3计算世界坐标系三个轴的单位矢量ox、oy、oz,其中ox、oy轴分别和正东、正北方向保持一致,oz轴竖直向上。
关系式3为:
关系式3中,α0为标靶面倾向,β0为标靶面倾角;
6)采用RANSAC算法从相机坐标系下的三维岩体点云模型中提取岩体结构面,得到结构面单位法向矢量,运用关系式4将其转换为世界坐标系下的单位法向矢量。
关系式4为:
关系式4中,(nx,ny,nz)为结构面在世界坐标系下的单位法向矢量,为结构面在相机坐标系下的单位法向矢量;
7)采用关系式5计算岩体结构面的产状信息。
关系式5为:
关系式5中,α为结构面倾向,β为结构面倾角;
8)重复步骤6)和7),提取岩体露头或岩体开挖坡面出现的所有结构面的产状信息。
如图2~4所示,下面通过一个具体实施例对本发明方法的技术方案进行详细说明。
(1)在武汉市南望山某一岩石露头面上布设一个矩形标靶(图2),使得标靶的一条边保持水平,用手机内置电子罗盘测量标靶面的产状,测得标靶面的倾向α0为197°,倾角β0为52°;
(2)利用智能手机对待测量岩体进行拍摄,确保两幅照片之间有60%以上的重叠率,且拍摄的照片中含有标靶,共拍摄25照片;
(3)将拍摄的25张手机照片导入专业的摄影软件中,生成基于相机坐标系的三维岩体点云模型,如图3所示;
(4)在三维岩体点云模型中测得标靶面A,B,C三个角点的坐标分别为(-0.326,1.553,18.631),(0.695,1.783,18.075),(-0.166,0.379,18.438),导入关系式1计算出向量AB的单位矢量ab=(0.862,0.194,-0.469),向量AC的单位矢量ac=(0.133,-0.978,-0.160),然后把ab、ac导入关系式2中计算出标靶面的单位法向矢量ab=(-0.490,0.076,-0.869);
(5)如图4所示,将测得标靶面的倾向α0=197°、倾角β0=52°、单位矢量ab、ac、ad导入关系式3中得到世界坐标系三个方向轴的单位矢量:ox=(0.085,-0.988,0.131),oy=(-0.177,0.115,0.977),oz=(-0.980,-0.106,-0.165);
(6)采用RANSAC算法从相机坐标系下的三维岩体点云模型中提取岩体结构面P1的单位法向矢量,运用关系式4得到结构面在世界坐标系下的单位法向矢量为(-0.030,0.688,-0.725)。
(7)采用关系式5得到岩体结构面P1的倾向为177.5°,倾角为43.5°。
(8)重复步骤(6)和(7),提取该岩体露头出现的结构面P2-P24,得到这些结构面的产状信息,见表1。
表1结构面产状
本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代,但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。
Claims (5)
1.一种基于智能手机的岩体结构面产状非接触测量方法,其特征在于,包含以下步骤:
1)在待测量岩体露头或开挖面上布设一个标靶,使得标靶的一条边保持水平,利用手机内置的电子罗盘测量标靶面的产状,记录其倾向和倾角;
2)对待测量岩体进行拍照,并确保两幅照片之间达到一定的重叠率,且拍摄的照片中含有标靶;
3)将拍摄获取的照片导入专业的摄影测量软件,生成基于相机坐标系的三维岩体点云模型;
4)在三维岩体点云模型中找到标靶,提取标靶面三个角点A,B,C的坐标,运用关系式1计算向量AB与向量AC的单位矢量ab和ac,运用关系式2计算标靶面单位法向矢量ad;
关系式1为:
关系式2为:
ad=ab×ac
关系式1和关系式2中,xA、xA、zA为点A的坐标,xB、yB、zB为点B的坐标,xC、yC、zC为点C的坐标;
5)根据罗盘测量得到的标靶面产状,运用关系式3计算世界坐标系三个轴的单位矢量ox、oy、oz,其中ox、oy轴分别和正东、正北方向保持一致,oz轴竖直向上;
关系式3为:
关系式3中,α0为标靶面倾向,β0为标靶面倾角;
6)从相机坐标系下的三维岩体点云模型中提取岩体结构面,得到结构面的单位法向矢量,运用关系式4将其转换为世界坐标系下的单位法向矢量;
关系式4为:
关系式4中,(nx,ny,nz)为结构面在世界坐标系下的单位法向矢量,为结构面在相机坐标系下的单位法向矢量;
7)采用关系式5计算岩体结构面的产状
关系式5为:
关系式5中,α为结构面倾向,β为结构面倾角;
8)重复步骤6)和7),提取岩体露头或岩体开挖坡面出现的所有结构面的产状信息。
2.根据权利要求1所述基于智能手机的岩体结构面产状非接触测量方法,其特征在于:所述步骤1)中的标靶为矩形。
3.根据权利要求1所述基于智能手机的岩体结构面产状非接触测量方法,其特征在于:所述步骤2)中两幅照片之间的重叠率大于60%。
4.根据权利要求1所述基于智能手机的岩体结构面产状非接触测量方法,其特征在于:所述步骤2)中利用手机对待测量岩体进行拍照。
5.根据权利要求1所述基于智能手机的岩体结构面产状非接触测量方法,其特征在于:所述步骤6)中采用RANSAC算法从相机坐标系下的三维岩体点云模型中提取岩体结构面。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910005679.6A CN109738440B (zh) | 2019-01-03 | 2019-01-03 | 一种基于智能手机的岩体结构面产状非接触测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910005679.6A CN109738440B (zh) | 2019-01-03 | 2019-01-03 | 一种基于智能手机的岩体结构面产状非接触测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109738440A true CN109738440A (zh) | 2019-05-10 |
CN109738440B CN109738440B (zh) | 2020-05-19 |
Family
ID=66363256
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910005679.6A Active CN109738440B (zh) | 2019-01-03 | 2019-01-03 | 一种基于智能手机的岩体结构面产状非接触测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109738440B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110426742A (zh) * | 2019-08-09 | 2019-11-08 | 浙江岩创科技有限公司 | 结构面产状的测量方法 |
CN110426020A (zh) * | 2019-08-20 | 2019-11-08 | 中国科学院地理科学与资源研究所 | 一种便携式地形摄影测量***及其测量方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103697854A (zh) * | 2013-12-10 | 2014-04-02 | 广西华锡集团股份有限公司 | 一种非接触式结构面产状测量方法 |
CN103983234A (zh) * | 2014-05-21 | 2014-08-13 | 重庆大学 | 一种基于智能移动设备的岩层产状测量方法 |
CN104183017A (zh) * | 2014-07-29 | 2014-12-03 | 浙江大学 | 基于地面三维激光点云实现地质体产状自动提取的方法 |
CN105180890A (zh) * | 2015-07-28 | 2015-12-23 | 南京工业大学 | 融合激光点云和数字影像的岩体结构面产状测量方法 |
CN106595571A (zh) * | 2016-12-20 | 2017-04-26 | 黄河勘测规划设计有限公司 | 基于手持式激光测距仪的地下洞室地质编录方法 |
CN106595567A (zh) * | 2016-12-20 | 2017-04-26 | 黄河勘测规划设计有限公司 | 地质结构面产状测量方法 |
CN108489403A (zh) * | 2018-06-08 | 2018-09-04 | 绍兴文理学院 | 基于三维激光扫描的露天矿山边坡岩体节理产状快速精细取值方法 |
-
2019
- 2019-01-03 CN CN201910005679.6A patent/CN109738440B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103697854A (zh) * | 2013-12-10 | 2014-04-02 | 广西华锡集团股份有限公司 | 一种非接触式结构面产状测量方法 |
CN103983234A (zh) * | 2014-05-21 | 2014-08-13 | 重庆大学 | 一种基于智能移动设备的岩层产状测量方法 |
CN104183017A (zh) * | 2014-07-29 | 2014-12-03 | 浙江大学 | 基于地面三维激光点云实现地质体产状自动提取的方法 |
CN105180890A (zh) * | 2015-07-28 | 2015-12-23 | 南京工业大学 | 融合激光点云和数字影像的岩体结构面产状测量方法 |
CN106595571A (zh) * | 2016-12-20 | 2017-04-26 | 黄河勘测规划设计有限公司 | 基于手持式激光测距仪的地下洞室地质编录方法 |
CN106595567A (zh) * | 2016-12-20 | 2017-04-26 | 黄河勘测规划设计有限公司 | 地质结构面产状测量方法 |
CN108489403A (zh) * | 2018-06-08 | 2018-09-04 | 绍兴文理学院 | 基于三维激光扫描的露天矿山边坡岩体节理产状快速精细取值方法 |
Non-Patent Citations (2)
Title |
---|
NA CHEN等: "Automatic extraction of blocks from 3D point clouds of fractured rock", 《COMPUTERS AND GEOSCIENCES》 * |
周春霖等: "双目***的岩体结构面产状非接触测量方法", 《岩石力学与工程学报》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110426742A (zh) * | 2019-08-09 | 2019-11-08 | 浙江岩创科技有限公司 | 结构面产状的测量方法 |
CN110426020A (zh) * | 2019-08-20 | 2019-11-08 | 中国科学院地理科学与资源研究所 | 一种便携式地形摄影测量***及其测量方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109738440B (zh) | 2020-05-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104964673B (zh) | 一种可定位定姿的近景摄影测量***和测量方法 | |
Gallay et al. | Large-scale and high-resolution 3-D cave mapping by terrestrial laser scanning: a case study of the Domica Cave, Slovakia | |
CN100450152C (zh) | 地质数码影像编录***及其使用方法 | |
CN101709964B (zh) | 大型洞室仪测成像可视化地质编录方法 | |
CN106323176B (zh) | 一种露天矿边坡的三维位移监测方法 | |
CN107525475B (zh) | 一种目标物高度测量方法 | |
CN107024174A (zh) | 基于三维激光扫描技术的粉状堆料体积测量装置及方法 | |
CN107655459A (zh) | 一种野外岩石结构面粗糙度的量测及计算方法 | |
CN105953777B (zh) | 一种基于深度图的大比例尺倾斜影像测图方法 | |
CN109738440A (zh) | 一种基于智能手机的岩体结构面产状非接触测量方法 | |
CN204963858U (zh) | 一种可定位定姿的近景摄影测量*** | |
CN106846477B (zh) | 一种编录野外地质影像的地质标记解译建模方法 | |
CN109631849A (zh) | 一种基于倾斜摄影的高陡边坡危岩测量方法 | |
CN103808306B (zh) | 基于windows的平板式施工地质可视化快速编录方法 | |
CN101770656B (zh) | 一种基于立体正射影像对的大场景立体模型生成及量测方法 | |
CN103162670B (zh) | 一种大比例尺地形图测绘内业屋檐改正方法 | |
CN106123853B (zh) | 一种隧道掌子面与边墙作地质体界面的测绘方法 | |
CN109357635A (zh) | 一种基于图像的特色建筑立面图测绘方法 | |
CN101298985A (zh) | 一种摄影测量物方控制杆及其使用方法 | |
CN107917692A (zh) | 一种基于无人机的块体识别方法 | |
JPWO2015071940A1 (ja) | 情報処理装置、情報処理方法及びプログラム | |
US11487044B2 (en) | Trench cross-section reference line setting device and trench cross-section analysis method using the same | |
CN106441222A (zh) | 基于数码照相的线路横断面测绘装置及方法 | |
CN105737802B (zh) | 基于体感摄影技术的堆积剖面空间结构信息分析方法 | |
CN114046779B (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 |