CN105205855A - 一种水土流失情况测定的方法 - Google Patents
一种水土流失情况测定的方法 Download PDFInfo
- Publication number
- CN105205855A CN105205855A CN201510582002.0A CN201510582002A CN105205855A CN 105205855 A CN105205855 A CN 105205855A CN 201510582002 A CN201510582002 A CN 201510582002A CN 105205855 A CN105205855 A CN 105205855A
- Authority
- CN
- China
- Prior art keywords
- block
- image
- dimensional point
- same place
- dimensional
- 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 Processing (AREA)
Abstract
本发明涉及一种水土流失情况测定的方法。方法包括:图像采集、对采集的图像进行分块处理、对各区块分别进行三维重建、拼接三维重建后的各区块得到该时刻和根据不同时刻的三维点云计算水土流失量。本发明的方法能够在降雨过程中观测。通过拍摄、解算降雨过程中下垫面的数字影像,提取土壤侵蚀过程中下垫面变化的时空信息,更加准确地描述土壤侵蚀过程。同时本发明也能够应用于类似的体积变化的量测。
Description
技术领域
本发明涉及图像三维重建技术和水土流失情况测定方法。
背景技术
土壤侵蚀过程是侵蚀力与下垫面相互作用的过程,下垫面既是土壤侵蚀作用的对象,也是土壤侵蚀作用的结果,在土壤侵蚀过程中的每一时刻,都使得下垫面的土体发生位移,下垫面的形态发生变化。围绕观测下垫面土体的位移量和下垫面形态的变化,发展了多种多样的观测技术和方法:
⑴接触式观测法
通过皮尺、卷尺、游标卡尺等常规的测量工具,对下垫面的形态变化进行简单测量。优点:简单易行;缺点:精度偏低,且只能观测次降雨条件下土壤侵蚀的总量,而不能观测土壤侵蚀过程中在不同时段或不同空间的变化量;
⑵示踪观测法
利用137Cs、7Be等天然放射性核素、稳定性稀土元素等按照指定要求添加到土壤里面,对下垫面的形态进行监测。优点:可观测土壤侵蚀速率、空间分布特征;缺点:观测的时间间隔太长、空间定位精度一般在米级以上;
⑵维激光扫描观测法
通过三维激光扫描仪对下垫面形态进行测量;原理:利用激光测距技术,通过记录被测物体表面大量的密集的点的三维坐标、反射率和纹理等信息,可快速测量被测物表面的三维点数据。优点:空间定位精度高;缺点:仪器不防水,不能在降雨过程中观测,且在沟蚀阶段会在沟道出现漏测,最大缺陷是测量效率太慢,不能监测下垫面形态动态演变过程;
⑷传统摄影测量观测法
用于监测河道侵蚀和重力侵蚀,其观测时间间隔分别为次降雨、天、年甚或更长,适合长时段侵蚀总量的观测,而不能观测短时段内相应分量的变化;
⑸径流泥沙采样法
在降雨过程当中,利用容器对下垫面侵蚀产生的泥沙按照不同时间段分别采样,对下垫面的土壤侵蚀量监测。优点:简单易行;缺点:泥沙的空间来源无法区分;
那么,如何在降雨条件下,从时空两个维度上观测下垫面形态演化以及表征土壤侵蚀各物理量。亦即,如何随着土壤侵蚀的发生发展,能够在设定的时间结点上,观测并量化表达下垫面的空间形态,进而根据下垫面空间形态计算土壤侵蚀量、输沙率、侵蚀速率等,以上的观测技术和方法都未能很好地解决这个问题。
发明内容
针对现有技术的缺陷或不足,本发明提供了一种基于立体视觉的图像三维重建方法。
本发明提供的基于立体视觉的图像三维重建方法包括:
步骤一,采集图像
采集目标区域某一时间窗口的图像;
步骤二,区域分块
按采集时间先后对所有图像进行序列标号;每连续多帧图像为一任务区块;相邻区块与区块之间有重叠图像,重叠部分的图像有相同序列标号;
步骤三,三维重建
对步骤二得到的各区块分别采用立体视觉技术进行三维重建,得到各区块的三维点云;
对各区块三维点云的各三维点进行编号,得到各三维点的编号特征;任意三维点A的编号特征为[(A1,a1),(A2,a2),…,(An,an),…,(AN-1,aN-1),(AN,aN)],其中:n=1,2,3,…,N;,N为三维点A对应的原始图像的张数,N取正整数;An为三维点A对应的所有原始图像中某一图像n的ID号;an为三维点A在ID号为An的原始图像上的SIFT特征点的ID号,(An,an)为编号特征的编号元组;
步骤四,拼接
(1)确定邻接区块
①寻找有重叠图像的相邻块与块之间的同名点,相邻块与块之间的两个三维点的编号特征交集包括至少一个编号元组时,两个三维点为同名点;
②两个块的同名点的个数占其中一个块中总三维点数的10%以上时,这两个区块为邻接区块;
③同名点个数占区块中三维点个数比例最高的区块为母站;
(2)邻接区块拼接
以母站为基准对其他区块进行拼接:
(2.1)确定稳定的同名点
①构建同名边
对于两个邻接区块,分别在各区块中利用该区块中的同名点组合构建边网,并计算每条边的长度;其中一个区块中有n个同名点,则该区块中组合构建的边网有条边;
②确定稳定的同名点
邻接区块之间的各组同名边中两条边长度的比值构成比例数组,比例数组包括个数值;采用基于拉依达准则(3σ准则)的奇异数据滤波法对个数值进行检验,得到可靠数值点,,可靠数值点对应的同名点作为稳定的同名点;
③基于稳定的同名点做两个三维区块之间的旋转、平移、缩放参数的计算,从而完成两个区块之间的拼接;
完成所有邻接区块之间的拼接后得到该时刻下测量区域的点云集合,去除点云集合中观测区域之外的点云,得到该时刻下测量区域的三维点云。
本发明还提供了一种水土流失情况测定方法,方法包括:
步骤A,构建参考面:采集目标区域下雨前的图像,执行权利要求步骤二至步骤四,利用所得三维点云中的多个三维点拟合得到的平面为参考面;
步骤B,获取不同时刻目标区域下雨时的三维点云:每一时刻的图像采集和处理执行权利要求1所述步骤一至步骤四;
步骤C,统一坐标系:将下雨时各时刻三维点云的坐标统一到同一坐标系下;
步骤D,计算参考面各时刻参考面以下的DEM体积量:将某时刻的三维点云内插得到对应的DEM数据,计算该时刻参考面以下的DEM体积量;
步骤E,相邻时刻参考面以下的DEM体积量差值为相邻时刻时间段内的土壤侵蚀量。
步骤A和步骤B的图像采集包括:采集目标区域某一时间窗口的土壤图像,所采集图像的重叠度至少为5度;单位平方米内,在10-12秒内采集图像的数量为120-240帧,且单位平方米至少采集120帧。
所述目标区域设有多个标靶。
所述时间窗口时长:15-180s。
步骤二按采集时间先后对所有图像进行序列标号,每连续50-230帧图像为一任务区块;相邻区块与区块之间有图像重叠,重叠图像的帧数占两个相邻区块中每个区块内图像帧数的10%-20%。
本发明的技术效果为:
1、本发明的方法能够在降雨过程中观测。通过拍摄、解算降雨过程中下垫面的数字影像,提取土壤侵蚀过程中下垫面变化的时空信息,更加准确地描述土壤侵蚀过程。同时本发明也能够应用于类似的体积变化的量测。
2、本发明的方法观测精度高,方法可观测一定时间段内自然界土壤侵蚀过程,对土壤侵蚀观测的空间和时间分布率更高。
3、本发明的方法对计算机硬件单元适用范围宽,可以根据实际的硬件环境,计算过程中的多个区块可通过联网进行并行计算或通过单机分批次计算。
4、本发明在确定邻接区块时采用SIFT特征算子的ID号构造编号组集合,通过区块之间的集合交集计算同名点,较之传统特征匹配算法计算速度至少提高了4倍。
附图说明
图1A为实施例中降雨前下垫面的示意图;图1B为实施例中降雨后下垫面的示意图;
图2为实施例中步骤一所采集图像的示意图;
图3为实施例中步骤三各区块的点云示意图;
图4为步骤四中三维重建后的示意图;
图5为步骤五中的参考面及计算示意图;
图6A为激光扫描仪观测得到的数字点云图;图6B为本发明的方法得到的数字点云图。
具体实施方式
以下是发明人提供的具体实施例,以对本发明的技术效果作进一步解释说明。
实施例:
该实施例于2014年7月9日在黄土高原土壤侵蚀与旱地农业国家重点实验室人工模拟降雨厅进行。
实验中人工模拟降雨的强度为90毫米每小时,所构造的下垫面如图1A、图1B所示,其中图1A为未发生侵蚀的(降雨前拍摄)下垫面图片,图1B为发生侵蚀的(降雨后拍摄)下垫面图片,钢槽底部为斜面,钢槽内部放置有待测定土壤,土壤类型为黄绵土,土壤装填容重为1.3g/cm3,下垫面装置的土槽尺寸为1米宽,5米长的立方体容器,下垫面底面与水平面构成15度倾斜角度。
分时段进行摄影并进行数字影像的解算,计算机视觉的下垫面土壤侵蚀观测法:从降雨开始至降雨结束对实验下垫面进行全覆盖的数字影像采集,采集时间分别在降雨开始后50min、70min、90min、100min,第100min后每隔10min采集一次,采集器距离下垫面的平均高度为80±5cm,数字影像的采集密度为150-170帧/m2;
收集水沙全样并换算流失量:在降雨开始后按照计算机视觉的下垫面土壤侵蚀观测的相同时间间隔收集实验过程中的水沙全样,计算样品中所含的土壤重量,并按土槽装填时1.3g/cm3的容重换算相应的体积量。
图像采集和图像处理过程如下:
一、获取不同时刻下测量区域的三维点云
步骤一,采集图像
区域采集面积:5平方米,本次采集得到1013张照片,采集过程历时62秒,如图2所示;
步骤二,对步骤一所采集的所有图像进行区域任务分块
按采集时间先后对步骤一所有图像进行序列标号,将步骤一中得到的1013张数字影像,以每块180张、相邻区块重叠65张的任务分割模式分割为9个任务子区,其中第9子区具有93张照片;
步骤三,三维重建
①对步骤二得到的各区块分别采用立体视觉技术进行三维重建,得到各区块的数字点云。
②特征点编号
对步骤三所得到的各个区块中三维点进行特征编号:任意三维点A的编号特征为[(A1,a1),(A2,a2),…,(An,an),…,(AN-1,aN-1),(AN,aN)],其中:n=1,2,3,…,N;,N为三维点A对应的原始图像的张数,N取正整数;An为三维点A对应的所有原始图像中某一图像n的ID号;an为三维点A在ID号为An的原始图像上的SIFT特征点的ID号,(An,an)为编号特征的编号元组;比如点云中的P点,该点在4张原始图像成像,其四张图像的ID号码是100,101,104,106,分别在其上的SIFT特征点的ID号码是:52,78,99,48.则编号特征表.[(100,52),(101,78),(104,99),(106,48)];如图3所示。
步骤四,拼接
(1)确定邻接区块
①有同名图像的区块为邻接区块;寻找有同名图像块与块之间的同名点,三维点与三维点的编号组集合的交集包括至少两个元素时为同名点;
②采用sift描述子匹配技术确定没有同名序列区块之间的同名特征点,两个区块的同名点的个数占其中一个区块中总三维点数的10%以上时,这两个区块为邻接区块;
定义同名点个数占区块中点云个数比例高的区块为母站;
(2)邻接区块拼接
以母站为基准对其他区块进行拼接:
(2.1)确定稳定的同名点
①构建同名边
对于两个邻接区块,分别在各区块中利用该区块中的同名点组合构建边网,并计算每条边的长度;例如:设其中一个区块中有n个同名点,则该区块中组合构建的边网有条边;
②确定稳定的同名点
邻接区块之间的各组同名边长度的比值构成比例数组,包括个数值;
采用基于拉依达准则(3σ准则)的奇异数据滤波法对个数值进行检验,得到可靠数值点,可靠数值点对应的同名点作为稳定的同名点:
对比例数组中的数值按大小排序,接着去掉排序后数组中的前、后各占总数值个数3%-5%的数值,得到中间数组,然后求取中间数组的平均值和方差,以3倍方差作为置信区间,中间数组中的数值落在置信区间以内的为可靠数值点,可靠数值点对应的同名点作为稳定的同名点;
基于稳定的同名点做两个三维区块之间的旋转、平移、缩放参数的计算,从而完成两个区块之间的拼接;
完成所有邻接区块之间的拼接后得到该时刻下测量区域的点云集合,去除点云集合中观测区域之外的点云,得到该时刻下测量区域的三维点云,如图4所示;
以两个标靶之间的实际尺寸为参考,确定三维点云的实际尺寸;
重复步骤一至步骤四,得到其他时刻下测量区域的三维点云;
二,计算土壤侵蚀量
(1)定义参考面:以未发生侵蚀的(降雨前)三维点云为基础拟合参考平面,具体为:选取150到200个三维点,这些三维点经过拟合得到一张平面,该平面即定义为整个测量过程中的参考面。
(2)统一降雨时各时刻三维点云的坐标系:某一时刻的三维点云作为一个基准,统一按照均匀分布在下垫面的周边的4-6个标靶的空间位置进行空间相似变换,由此,三维点云的坐标系统纳入到标靶坐标系统中,标靶在整个降雨过程中空间位置是稳定不变的,因此所有的点云数据都具有了统一的坐标系。
(3)参考面以下体积的计算:将某一时刻三维点云内插得到对应的DEM数据。人机交互在DEM数据上指定方向参考点,如图5中十字交叉点,计算参考面以下的DEM体积量。
(4)相邻时刻体积量的较差即为该时段的侵蚀量。
该实施例参考面相关数据为:图5中格网为参考平面;十字为方向参考点位置;参考面左顶点:(1.319944,160.862244);参考面右底点:(59.419609,-15.33840);DEM步距:0.1毫米;方向参考点坐标:x=20.219984,y=79.760979,z=-21.923674。
表1本发明的方法与径流泥沙采集法观测结果对比表
与传统的径流泥沙采集法平行观测结果表明全程平均准确度为1.05%。
发明人将本发明的检测方法与三维激光扫描仪检测观测方法进行了对比:
采用实施例最后一次的观测结果与三维激光扫描仪观测结果进行对比,在实验观测降雨结束后,用三维激光扫描仪和本方法观测并计算下垫面的三维数字点云,结果如表2和图6所示。
表2本发明的方法与激光扫描观测结果对比
与激光扫描仪平行观测表明:本发明的方法克服了激光扫描仪不能够在降雨过程中观测下垫面演化的缺陷;二是下垫面形态的数字点云的密度提高了27.8%,即本发明的方法测量每个点云所代表的物点尺寸要更小,对下垫面三维形态的表达要更准确;三是本发明的方法克服了激光扫描仪观测时沟道底部激光线不能投射到位而造成的漏测现象(图6A中黑色部分区域是漏测区域),实现了全覆盖的数字影像采集。
Claims (6)
1.一种基于立体视觉的图像三维重建方法,其特征在于,方法包括:
步骤一,采集图像
采集目标区域某一时间窗口的图像;
步骤二,区域分块
按采集时间先后对所有图像进行序列标号;每连续多帧图像为一任务区块;相邻区块与区块之间有重叠图像,重叠部分的图像有相同序列标号;
步骤三,三维重建
对步骤二得到的各区块分别采用立体视觉技术进行三维重建,得到各区块的三维点云;
对各区块三维点云的各三维点进行编号,得到各三维点的编号特征;任意三维点A的编号特征为[(A1,a1),(A2,a2),…,(An,an),…,(AN-1,aN-1),(AN,aN)],其中:n=1,2,3,…,N;,N为三维点A对应的原始图像的张数,N取正整数;An为三维点A对应的所有原始图像中某一图像n的ID号;an为三维点A在ID号为An的原始图像上的SIFT特征点的ID号,(An,an)为编号特征的编号元组;
步骤四,拼接
(1)确定邻接区块
①寻找有重叠图像的相邻区块与区块之间的同名点:相邻区块与区块之间的某两个三维点的编号特征交集包括至少一个编号元组时,两个三维点为同名点;
②两个区块的同名点的个数占其中一个区块中总三维点数的10%以上时,这两个区块为邻接区块;
③同名点个数占区块中三维点个数比例最高的区块为母站;
(2)邻接区块拼接
以母站为基准对其他区块进行拼接:
(2.1)确定稳定的同名点
①构建同名边
对于两个邻接区块,分别在各区块中利用该区块中的同名点组合构建边网,并计算每条边的长度;其中一个区块中有n个同名点,则该区块中组合构建的边网有条边,两个邻接区块之间具有组同名边,同名边对应同名点;
②确定稳定的同名点
邻接区块之间的各组同名边中两条边长度的比值构成比例数组,比例数组包括个数值;采用基于拉依达准则(3σ准则)的奇异数据滤波法对个数值进行检验,得到可靠数值点,可靠数值点对应的同名点作为稳定的同名点;
③基于稳定的同名点做两个邻接区块之间的旋转、平移、缩放参数的计算,从而完成两个区块之间的拼接;
完成所有邻接区块之间的拼接后得到该时刻下测量区域的点云集合,去除点云集合中观测区域之外的点云,得到该时刻下测量区域的三维点云。
2.一种水土流失情况测定方法,其特征在于,方法包括:
步骤A,构建参考面:采集目标区域降雨前的图像,执行权利要求1所述步骤二至步骤四,利用所得三维点云中的多个三维点拟合得到的平面为参考面;
步骤B,获取不同时刻目标区域下雨时的三维点云:每一时刻的图像采集和处理执行权利要求1所述步骤一至步骤四;
步骤C,统一坐标系:将降雨时各时刻三维点云的坐标统一到同一坐标系下;
步骤D,计算各时刻参考面以下的DEM体积量:将某时刻的三维点云内插得到对应的DEM数据,计算该时刻参考面以下的DEM体积量;
步骤E,相邻时刻参考面以下的DEM体积量差值为相邻时刻时间段内的土壤侵蚀量。
3.权利要求2所述的土壤流失情况测定方法,其特征在于,步骤A和步骤B的图像采集包括:采集目标区域某一时间窗口的图像,所采集图像的重叠度至少为5度;单位平方米内,在10-12秒内采集图像的数量为120-240帧,且单位平方米至少采集120帧。
4.权利要求2所述的土壤流失情况测定方法,其特征在于,所述目标区域设有多个标靶。
5.如权利要求3所述的土壤流失情况测定方法,其特征在于,所述时间窗口时长:15-180s。
6.如权利要求3所述的土壤流失情况测定方法,其特征在于,步骤二按采集时间先后对所有图像进行序列标号,每连续50-230帧图像为一任务区块;相邻区块与区块之间有图像重叠,重叠图像的帧数占两个相邻区块中每个区块内图像帧数的30%-45%。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510582002.0A CN105205855B (zh) | 2015-09-14 | 2015-09-14 | 一种水土流失情况测定的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510582002.0A CN105205855B (zh) | 2015-09-14 | 2015-09-14 | 一种水土流失情况测定的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105205855A true CN105205855A (zh) | 2015-12-30 |
CN105205855B CN105205855B (zh) | 2017-11-21 |
Family
ID=54953513
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510582002.0A Expired - Fee Related CN105205855B (zh) | 2015-09-14 | 2015-09-14 | 一种水土流失情况测定的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105205855B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108180897A (zh) * | 2018-01-06 | 2018-06-19 | 中国科学院、水利部成都山地灾害与环境研究所 | 坡耕地土壤水力侵蚀速率快速估算方法 |
CN108332719A (zh) * | 2018-01-06 | 2018-07-27 | 中国科学院、水利部成都山地灾害与环境研究所 | 坡耕地土壤耕作侵蚀速率快速估算方法 |
CN109141403A (zh) * | 2018-08-01 | 2019-01-04 | 上海航天控制技术研究所 | 一种星敏感器小窗口访问的图像处理系统及其方法 |
CN110363857A (zh) * | 2019-06-14 | 2019-10-22 | 西北农林科技大学 | 降雨过程中土壤侵蚀地表形态dem构建、观测方法及装置 |
CN110390687A (zh) * | 2019-07-29 | 2019-10-29 | 四川大学 | 一种基于三维激光扫描的河道冲淤测量方法 |
CN110763194A (zh) * | 2019-11-01 | 2020-02-07 | 中国矿业大学 | 一种地面三维激光扫描无标靶监测矿区地表沉陷的方法 |
CN112509025A (zh) * | 2020-12-03 | 2021-03-16 | 山东省科学院海洋仪器仪表研究所 | 一种基于三维欧氏距离的岩石空间结构距离图的计算方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101793543A (zh) * | 2010-04-02 | 2010-08-04 | 郑粉莉 | 利用三维激光扫描技术动态监测坡面尺度侵蚀发育过程的研究方法 |
US20140077969A1 (en) * | 2012-09-19 | 2014-03-20 | The Boeing Company | Forest Sensor Deployment and Monitoring System |
CN104181548A (zh) * | 2014-09-01 | 2014-12-03 | 中国电建集团成都勘测设计研究院有限公司 | 基于三维激光扫描技术对水土保持状况的监测方法 |
-
2015
- 2015-09-14 CN CN201510582002.0A patent/CN105205855B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101793543A (zh) * | 2010-04-02 | 2010-08-04 | 郑粉莉 | 利用三维激光扫描技术动态监测坡面尺度侵蚀发育过程的研究方法 |
US20140077969A1 (en) * | 2012-09-19 | 2014-03-20 | The Boeing Company | Forest Sensor Deployment and Monitoring System |
CN104181548A (zh) * | 2014-09-01 | 2014-12-03 | 中国电建集团成都勘测设计研究院有限公司 | 基于三维激光扫描技术对水土保持状况的监测方法 |
Non-Patent Citations (3)
Title |
---|
J POESEN ET AL: "Gully erosion and environmental change:importance and research needs", 《CATENA》 * |
张姣 等: "利用三维激光扫描技术动态监测沟蚀发育过程的方法研究", 《水土保持通报》 * |
李瑾杨: "基于点云数据的冲沟溯源侵蚀过程动态可视化研究与实现", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108180897A (zh) * | 2018-01-06 | 2018-06-19 | 中国科学院、水利部成都山地灾害与环境研究所 | 坡耕地土壤水力侵蚀速率快速估算方法 |
CN108332719A (zh) * | 2018-01-06 | 2018-07-27 | 中国科学院、水利部成都山地灾害与环境研究所 | 坡耕地土壤耕作侵蚀速率快速估算方法 |
CN109141403A (zh) * | 2018-08-01 | 2019-01-04 | 上海航天控制技术研究所 | 一种星敏感器小窗口访问的图像处理系统及其方法 |
CN110363857A (zh) * | 2019-06-14 | 2019-10-22 | 西北农林科技大学 | 降雨过程中土壤侵蚀地表形态dem构建、观测方法及装置 |
CN110363857B (zh) * | 2019-06-14 | 2022-09-20 | 西北农林科技大学 | 降雨过程中土壤侵蚀地表形态dem构建、观测方法及装置 |
CN110390687A (zh) * | 2019-07-29 | 2019-10-29 | 四川大学 | 一种基于三维激光扫描的河道冲淤测量方法 |
CN110390687B (zh) * | 2019-07-29 | 2021-06-01 | 四川大学 | 一种基于三维激光扫描的河道冲淤测量方法 |
CN110763194A (zh) * | 2019-11-01 | 2020-02-07 | 中国矿业大学 | 一种地面三维激光扫描无标靶监测矿区地表沉陷的方法 |
CN112509025A (zh) * | 2020-12-03 | 2021-03-16 | 山东省科学院海洋仪器仪表研究所 | 一种基于三维欧氏距离的岩石空间结构距离图的计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105205855B (zh) | 2017-11-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105205855A (zh) | 一种水土流失情况测定的方法 | |
CN111623722B (zh) | 一种基于多传感器的边坡变形立体监测系统和方法 | |
CN109870108B (zh) | 输电线路覆冰检测方法及装置 | |
CN112434617B (zh) | 一种基于多源遥感数据的森林生物量变化监测方法及系统 | |
CN105004611A (zh) | 一种研究混凝土断裂过程区材料力学行为方法 | |
CN106597416A (zh) | 一种地面GPS辅助的LiDAR数据高程差的误差修正方法 | |
CN114998536A (zh) | 基于新型基础测绘的模型生成方法、装置及存储介质 | |
Xue et al. | Novel SfM-DLT method for metro tunnel 3D reconstruction and visualization | |
CN103886555A (zh) | 一种基于海量三维激光扫描点云数据的处理方法 | |
CN104062653B (zh) | 基于三维微地貌重建的雷达遥感多尺度面粗糙度测量方法 | |
CN113869629A (zh) | 一种基于激光点云的输电线路安全风险分析判断评估方法 | |
CN109345617A (zh) | 一种基于长条带多站点云的链式高精度拼接与平差方法 | |
Corradetti et al. | Quantitative analysis of folds by means of orthorectified photogrammetric 3D models: a case study from Mt. Catria, Northern Apennines, Italy | |
CN116011291A (zh) | 一种隧道洞口多源信息融合的边坡稳定性数字化评估平台 | |
CN112825190A (zh) | 一种精度评估方法、系统、电子设备和存储介质 | |
CN111982822A (zh) | 一种长时间序列高精度植被指数改进算法 | |
CN114140691B (zh) | 基于无人机可见光图像预测作物倒伏角度和形态的方法 | |
Motayyeb et al. | Fusion of UAV-based infrared and visible images for thermal leakage map generation of building facades | |
CN102865814A (zh) | 植物群体三维重建误差测量方法 | |
CN114004950A (zh) | 一种基于BIM与LiDAR技术的路面病害智能识别与管理方法 | |
CN112504156A (zh) | 一种基于前景网格的结构表面应变测量系统及测量方法 | |
Marino et al. | Measuring rock slope damage on rubble mound breakwater through digital photogrammetry | |
Lascelles et al. | Automated digital photogrammetry: a valuable tool for small‐scale geomorphological research for the non‐photogrammetrist? | |
Shen et al. | A review of terrestrial laser scanning (TLS)-based technologies for deformation monitoring in engineering | |
CN110363857B (zh) | 降雨过程中土壤侵蚀地表形态dem构建、观测方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20171121 Termination date: 20190914 |