CN112364802A - 一种崩塌滑坡灾害体的变形监测方法 - Google Patents
一种崩塌滑坡灾害体的变形监测方法 Download PDFInfo
- Publication number
- CN112364802A CN112364802A CN202011300788.XA CN202011300788A CN112364802A CN 112364802 A CN112364802 A CN 112364802A CN 202011300788 A CN202011300788 A CN 202011300788A CN 112364802 A CN112364802 A CN 112364802A
- Authority
- CN
- China
- Prior art keywords
- crack
- image
- deformation
- monitoring
- coordinate system
- 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
- 238000012544 monitoring process Methods 0.000 title claims abstract description 81
- 238000000034 method Methods 0.000 title claims abstract description 51
- 230000000007 visual effect Effects 0.000 claims abstract description 36
- 238000004364 calculation method Methods 0.000 claims abstract description 7
- 230000008859 change Effects 0.000 claims description 22
- 238000000605 extraction Methods 0.000 claims description 11
- 238000005311 autocorrelation function Methods 0.000 claims description 9
- 230000003287 optical effect Effects 0.000 claims description 8
- 238000003708 edge detection Methods 0.000 claims description 7
- 230000011218 segmentation Effects 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 5
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 238000004458 analytical method Methods 0.000 claims description 3
- 239000011541 reaction mixture Substances 0.000 claims description 2
- YHXISWVBGDMDLQ-UHFFFAOYSA-N moclobemide Chemical compound C1=CC(Cl)=CC=C1C(=O)NCCN1CCOCC1 YHXISWVBGDMDLQ-UHFFFAOYSA-N 0.000 claims 8
- 238000012545 processing Methods 0.000 abstract description 10
- 238000003384 imaging method Methods 0.000 description 9
- 206010017076 Fracture Diseases 0.000 description 6
- 238000006073 displacement reaction Methods 0.000 description 6
- 208000010392 Bone Fractures Diseases 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 239000011435 rock Substances 0.000 description 4
- 238000000691 measurement method Methods 0.000 description 3
- 230000006378 damage Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000035515 penetration Effects 0.000 description 2
- 239000002689 soil Substances 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 238000003648 Ljung–Box test Methods 0.000 description 1
- 229910001294 Reinforcing steel Inorganic materials 0.000 description 1
- 229910000831 Steel Inorganic materials 0.000 description 1
- 208000013201 Stress fracture Diseases 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 239000011083 cement mortar Substances 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000003706 image smoothing Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000002156 mixing Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 238000011155 quantitative monitoring Methods 0.000 description 1
- 230000010076 replication Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 239000010959 steel Substances 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 239000002352 surface water Substances 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
- 239000002023 wood Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/50—Context or environment of the image
- G06V20/52—Surveillance or monitoring of activities, e.g. for recognising suspicious objects
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
- G06Q50/265—Personal security, identity or safety
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/20—Image preprocessing
- G06V10/26—Segmentation of patterns in the image field; Cutting or merging of image elements to establish the pattern region, e.g. clustering-based techniques; Detection of occlusion
- G06V10/267—Segmentation of patterns in the image field; Cutting or merging of image elements to establish the pattern region, e.g. clustering-based techniques; Detection of occlusion by performing operations on regions, e.g. growing, shrinking or watersheds
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/40—Extraction of image or video features
- G06V10/44—Local feature extraction by analysis of parts of the pattern, e.g. by detecting edges, contours, loops, corners, strokes or intersections; Connectivity analysis, e.g. of connected components
-
- G—PHYSICS
- G08—SIGNALLING
- G08B—SIGNALLING OR CALLING SYSTEMS; ORDER TELEGRAPHS; ALARM SYSTEMS
- G08B21/00—Alarms responsive to a single specified undesired or abnormal condition and not otherwise provided for
- G08B21/02—Alarms for ensuring the safety of persons
- G08B21/10—Alarms for ensuring the safety of persons responsive to calamitous events, e.g. tornados or earthquakes
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Multimedia (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Economics (AREA)
- General Business, Economics & Management (AREA)
- Development Economics (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Entrepreneurship & Innovation (AREA)
- Health & Medical Sciences (AREA)
- Computer Security & Cryptography (AREA)
- Educational Administration (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Game Theory and Decision Science (AREA)
- General Health & Medical Sciences (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- Emergency Management (AREA)
- Pit Excavations, Shoring, Fill Or Stabilisation Of Slopes (AREA)
Abstract
本发明提供了一种崩塌滑坡灾害体的变形监测方法,属于地质灾害监测领域。所述崩塌滑坡灾害体的变形监测方法,构建世界坐标系及视觉传感器坐标系,并采用激光可视化传感器对灾害体进行实时侦测,采集基于坐标系的图像数据,关键裂缝识别模块读取图像数据后,提取目标特征,获得关键裂缝图像数据,计算所述关键裂缝的形状及实际变化量,并基于世界坐标系恢复所述关键裂缝的三维地理坐标;基于时间序列构建变形监测模型,分析关键裂缝的变形,结合三维地理坐标进行灾害预警。本发明在有效时间内对崩塌滑坡灾害体的变形趋势进行判断和及时预警,提高了崩塌滑坡灾害体变形监测的实时性、准确性和稳定性,数据处理计算量小、适应性强,节约了监测资源。
Description
技术领域
本发明属于地质灾害监测领域,具体涉及一种崩塌滑坡灾害体的变形监测方法。
背景技术
危岩地段存在滑坡等灾害隐患,危险发生时会使建筑物、甚至使整个居民点遭到毁坏,使公路和铁路被掩埋。由危岩带来的损失,不只是建筑物毁坏的直接损失,并且常因此而使交通中断,给运输带来重大损失,还会给人民的生命财产及安全带来较大危险。滑坡灾害监测包括地表变形监测、深部位移监测、力学参数监测以及环境影响因素(地表水、地下水、降雨量等)监测和宏观地质现象监测。
在无有效方法监测深部大变形的情况下,地表形变监测是滑坡、崩塌滑坡灾害体监测预警的重要内容和有效手段。对于地表形变监测,采用的方法包括以全站仪、GPS技术为代表的离散单点式测量方法,以InSAR技术为代表的光学影像面式测量方法。
其中,以全站仪、GPS技术为代表的离散单点式测量方法,虽然可同时测出单个测量点的三维(X、Y、Z)位移量,不受视通条件限制,且精度在不断提高,但是,只能对布设了待测点的特定点进行监测,无法监测复杂环境或设备及人员无法到达的区域,监测点数少、工作效率低、监测周期长、成本相对较高,同时经常受到天气和时间的限制,一旦破坏监测点,资料的连续性会受到严重影响。以InSAR技术为代表的光学影像面式测量方法,是一种主动式遥感成像方式全天候对地观测成像,同时对植被具有一定的穿透能力,覆盖范围大、空间分辨率高,但是该方法对于地表覆盖的变化、卫星轨道参数的误差和大气参数的变化非常敏感;干涉像对之间基线的选择也有一定的限制条件;高山地区成像时存在雷达阴影和迭掩现象;数据的时间分辨率有时满足不了地表变化监测的要求;对于地表出现的剧烈变化无法监测出来;虽然覆盖范围大,但是精度较低。
对于宏观地质现象监测,主要是利用简单实用的工具来测量,在崩塌滑坡灾害体裂缝两侧(或上下)插入木筋、钢筋、埋桩法标记,使用钢卷尺进行裂缝测量,还可以在裂缝上黏贴水泥砂浆片来测量裂缝变形情况。但是不能实时获取裂缝变化情况,具有一定的滞后性。
对于深部位移监测,通常采用地声监测法,利用地声发射仪、地音探测仪等,采集岩体变形微破裂或破坏时释放出的应力波强度、频度等信号资料,分析、判断崩滑体变形情况。仪器一般应设置在崩滑体应力集中部位,地表、地下均可,灵敏度较高,可连续监测,但仅适用于岩质崩滑体或斜坡的变形监测,且在崩滑体匀速变形阶段不宜使用。
对于力学参数监测,通常采用地应力监测,利用埋设于钻孔、平硐、竖井内的地应力计监测岩质崩滑体内不同部位的应力变化,分析、判断崩滑体变形情况。也可在地表安设水平应力计,监测地表应力变化情况,分辨拉力区、压力区等。另外,利用差动传递式土压力计、应变式压力计,可监测土质崩滑体地表应力变化情况。
由此可见,传统的崩塌滑坡灾害体变形监测在进行监测工作前需布控变形监测网,且还要在监测体上布置一定数量的监测点,然后才可以用相应的测绘仪器对监测体进行变形监测。布网和布点工作复杂繁琐,需要投入大量的工作时间及人力物力,而且布设好的控制点和监测点很容易受到外界的影响与破坏,无法满足当前监测技术要求。
另外,近年兴起的遥感(RS)法和近景摄影法,利用地球卫星或飞机,周期性的拍摄崩滑体的变形,适用于大范围、区域性崩滑体监测,虽然精度较高,但是费用也相对较高,并不适合24小时实时动态监测,且不能实时测出崩滑体的具体变化,具有一定的滞后性。
上述监测方法均无法满足需要越来越高的安全需求的崩滑体监测。
发明内容
鉴于上述问题,本发明实施例提供了一种崩塌滑坡灾害体的变形监测方法,通过对崩塌滑坡灾害体关键裂缝位移变形的量化监测,实时监测崩塌滑坡灾害体的变形,预测灾害体的运动趋势并进行预警,提高对崩塌滑坡灾害体变形的识别能力及变形趋势判别的准确性。
为了实现上述目的,本发明实施例所采用的技术方案如下:
本发明实施例提供了一种崩塌滑坡灾害体的变形监测方法,包括如下步骤:
步骤S1,构建世界坐标系,并以世界坐标系为基准坐标系,构建视觉传感器坐标系;
步骤S2,采用激光可视化视觉传感器对当前崩塌滑坡灾害体进行实时侦测,并基于视觉传感器坐标系采集灾害体中的图像数据,以供关键裂缝识别模块调用;
步骤S3,关键裂缝识别模块从所述视觉传感器中读取基于视觉传感器坐标系的灾害体图像数据,对基于图像数据的二维数字信号进行目标特征提取,获得关键裂缝的图像数据,计算所述关键裂缝的形状及实际变化量,并基于世界坐标系恢复所述关键裂缝的三维地理坐标;
步骤S4,基于时间序列构建变形监测模型,分析关键裂缝的变形,结合三维地理坐标进行灾害预警。
作为本发明的一个优选实施例,所述步骤S3中所述对基于图像数据的二维数字信号进行目标特征提取,获得所述灾害体图像数据中的关键裂缝的图像数据,计算所述裂缝的形状及实际变化量,具体包括:
步骤S301,将采集到的彩色图像转化为灰度图像;
步骤S302,采用中值滤波器对图像进行滤波处理,把数字图像或数字序列中一点的值用该点的一个领域中各点值的中值代替,让周围的像素值接近真实值,消除噪声点;
步骤S303,从滤波后的图像数据中提取裂缝特征,并根据裂缝特征对图像进行分割,标定出关键裂缝的图像数据;
步骤S304,采用Canny算子进行亚像素级边缘检测,提取出裂缝与背景图像中的交界线;
步骤S305,根据所述裂缝与背景图像中的交界线,获取裂缝形状与变化率。
作为本发明的一个优选实施例,所述步骤S301,进一步为:利用加权平均法将彩色图像进行灰度转化,根据彩色图像中R、G、B三个分量进行不同权值的加权平均,采用式(5)得到灰度图像:
f(i,j)=0.30R(i,j)+0.59G(i,j)+0.11B(i,j) (5)
式(5)中,(i,j)是图像中的像素坐标。
作为本发明的一个优选实施例,所述步骤S303中对图像进行分割,采用Otsu阈值分割算法来有效标定图像中的裂缝特征,通过在图像的灰度取值范围内选择一灰度阈值,并对全部像素点应用该阈值。
作为本发明的一个优选实施例,所述阈值确定过程如下:
用C0表示裂缝图像,C1表示背景,图像像素总数为N,灰度区间[0,L-1],C0和C1类方差分别为:
其中T为阈值,pi为像素点中值,w0,w1,u0,u1的计算分别如下:
pi=ni/N (10)
其中,ni代表在灰度区间内灰度级为i的像素点数;
C0和C1类的类间方差为:
C0和C1类的类间方差为:
通过让T在[0,L-1]范围内依次取值,使式(18)中T取得最大值,解得阈值T。
作为本发明的一个优选实施例,所述步骤S305裂缝形状与变化率计算过程如下:
设Xi为第ti时刻裂缝区域的像素总值,则Xi+1为第ti+1时刻裂缝区域的像素总值,在第i帧的裂缝目标区域像素总值设为Xi,第i+1帧的像素总值为Xi+1,令第i帧的裂缝面积为Li,第i+1帧区域面积为Li+1,在连续视觉图像帧下的裂缝面积变化率:
式(19)中,j表示第j帧,i+j表示第i+j帧,ΔB表示裂缝的面积变化率,限定于在时间范围ti到ti+1这一间断时间内。
作为本发明的一个优选实施例,所述步骤S4中构建变形监测及预警模型,具体包括如下步骤:
步骤S401,构建自回归滑动平均ARIMA模型;
步骤S402,计算自相关函数和偏自相关函数并进行特点分析,优化ARIMA模型;
步骤S403,基于优化后的ARIMA模型,采用最小二乘法计算ARIMA模型参数;
步骤S404,采用Ljung-Box对优化后的ARIMA模型进行检验,获得最终变形监测与预警模型。
作为本发明的一个优选实施例,所构建的ARIMA模型为:
式(20)中,是p个自回归系数,et是满足均值为零,方差为的白噪声序列;{Xt}是平稳、正态、零均值的时间序列,p和q代表自起始开始的步数;对所述模型进行d阶差分,获得非平稳时间序列模型ARIMA(p,d,q)。
作为本发明的一个优选实施例,所述步骤S4进行灾害预警,采用优化并检验后的ARIMA(p,d,q)模型,根据历史获取的裂缝监测数据X1,X2,···,Xt-1,Xt推测未来T期的Xt+T,用Xt(T)代表预测值,最小均方误差预测公式如下:
Xt(T)=E(Xt+T|X1,X2,···,Xt-1,Xt) (26)。
作为本发明的一个优选实施例,所述视觉传感器坐标系是固定在传感器上的直角坐标系,原点定义在传感器光心,光心到图像平面的距离为视觉传感器有效焦距。
本发明实施例技术方案具有如下有益效果:
采用机器视觉与图像处理相结合的方式量化崩塌滑坡灾害体关键裂缝位移变形,对崩塌滑坡灾害体的裂缝变形实时动态监测,利用图像处理和特征提取算法,对崩塌滑坡灾害体裂缝变形监测数据进行快速分析与处理,精确、可靠、有效地计算出崩塌滑坡灾害体关键裂缝的具体变形量,再通过构建适合崩塌滑坡灾害体变形监测预警的时间序列数学模型,确保了在有效的时间内对崩塌滑坡灾害体的运动趋势做出判断,提高了崩塌滑坡灾害体变形监测的实时性、准确性、和时效性,为崩塌滑坡灾害体运动趋势预测提供有效的理论依据和数据支撑,具有重要的工程应用价值。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施方式提供的崩塌滑坡灾害体的变形监测方法流程图;
图2为本发明实施方式变形监测方法中构建坐标体系示意图;
图3为本发明实施方式变形监测方法中视觉传感器布局示例图。
具体实施方式
下面通过参考示范性实施例,并结合附图,对本发明技术问题、技术方案和优点进行详细阐明。以下所述示范性实施例仅用于解释本发明,而不能解释为对本发明的限制。本技术领域技术人员可以理解,除非另外定义,这里使用的所有术语(包括技术术语和科学术语)具有与本发明所属领域中的普通技术人员的一般理解相同的意义。还应该理解的是,诸如通用字典中定义的那些术语应该被理解为具有与现有技术的上下文中的意义一致的意义,并且除非在这里进行定义,否则不会用理想化或过于正式的含义来解释。
本发明实施方式提供了一种崩塌滑坡灾害体的变形监测方法,通过适合裂缝监测的视觉传感器标定及透视模型获取崩塌滑坡灾害体关键裂缝的基础数据,再基于图像处理和特征提取识别计算出关键裂缝变形位移量,利用亚像素级边缘检测算法,结合Otsu阈值分割算法及Canny算子,有效提取崩塌滑坡灾害体重点变形区物理图像像素,进而计算出重点变形区裂缝变化的实际物理量值,对关键裂缝进行量化监测、分析处理,再基于监测预警的时间序列模型,在有效时间内对崩塌滑坡灾害体的变形趋势进行判断和及时预警,提高崩塌滑坡灾害体变形监测的实时性、准确性、稳定性以及时效性,数据处理计算量小、适应性强,节约监测资源。
图1示出了本发明实施方式提供的崩塌滑坡灾害体的变形监测方法流程图。如图1所示,所述崩塌滑坡灾害体的变形监测方法包括如下步骤:
步骤S1,构建世界坐标系,并以世界坐标系为基准坐标系,构建视觉传感器坐标系。
如图2所示,本步骤中,世界坐标系(P(Xw,Yw,Zw))是一个基准坐标系,用来描述视觉传感器和关键裂缝在灾害体中的具体位置。通过建立图像坐标系、成像平面坐标系、视觉传感器坐标系、世界坐标系来构建出裂缝透视模型。基于世界坐标系,构建视觉传感器坐标系,则视觉传感器基于自身坐标系采集的图像信息,可转化到世界坐标系下;视觉传感器坐标系OcXcYcZc是固定在传感器上的直角坐标系,原点定义在传感器光心,光心到图像平面的距离称之为视觉传感器有效焦距。基于视觉传感器坐标系,还可以构建图像坐标系、成像平面坐标系,其中,图像坐标系Ouv以数组形式存在,每个元素(像素)以坐标形式出现,(u,v)分别代表像素在数组中的列和行,这里的坐标即为像素点在视觉传感器坐标系中的坐标值;图像坐标只表示像素的列和行,并没有实际表示像素在图像中的物理位置;在图像坐标系基础上建立以物理量单位(如毫米)表示的成像平面坐标系O1xy(x,y),表示以物理单位度量的成像平面坐标系。
这里将视觉传感器下图像像素的坐标值,及加了单位量后的坐标值,将坐标集的集合,分别称为图像坐标系、成像平面坐标系,所述图像坐标系和成像平面坐标系对应于相应的图像,方便图像数据的调用。
步骤S2,采用激光可视化视觉传感器对当前崩塌滑坡灾害体进行实时侦测,并基于视觉传感器坐标系采集灾害体中的图像数据,以供关键裂缝识别模块调用。
如图3所示,本步骤中,采用激光可视化视觉传感器,具备激光、透雾功能,从而可以实现夜晚、阴雨天、雾天等灾害多发环境下对灾害体的观测。所述基于视觉传感器坐标系和图像坐标系的图像数据,图像的每一点都具有视觉传感器坐标系的坐标值,坐标值可以反映出图像在视觉传感器坐标系中的相对位置。
步骤S3,关键裂缝识别模块从所述视觉传感器中读取基于视觉传感器坐标系的灾害体图像数据,对基于图像数据的二维数字信号进行目标特征提取,获得所述灾害体图像数据中的关键裂缝的图像数据,计算所述裂缝的形状及实际变化量,并基于世界坐标系恢复所述关键裂缝的三维地理坐标。
本步骤中,图像数据是基于视觉传感器坐标系的二维数字信号,通过视觉传感器坐标系与世界坐标系间的对应关系,计算二维图像与三维场景之间的变换关系,标定主点坐标、比例尺度、焦距、旋转矩阵、平移矢量等关键信息。
同时,对所述二维图像,采用快速边缘检测算法进行目标特征提取,抑制地质灾害现场图像带来的高斯白噪声污染;再根据标定的关键信息,结合Canny算子和Otsu阈值分割算法,计算出崩塌滑坡灾害体整体或局部裂缝的具体变形量。上述计算过程,计算量小,易于实现对裂缝的动态监测。
其中,所述对基于图像数据的二维数字信号进行目标特征提取,获得所述灾害体图像数据中的关键裂缝的图像数据,计算所述裂缝的形状及实际变化量,具体包括:
步骤S301,将采集到的彩色图像转化为灰度图像,有效节约图像的计算量。
优选地,本步骤中利用加权平均法将彩色图像进行灰度转化,根据彩色图像中R、G、B三个分量进行不同权值的加权平均,采用式(5)可以得到灰度图像:
f(i,j)=0.30R(i,j)+0.59G(i,j)+0.11B(i,j) (5)
式(5)中,(i,j)是图像中的像素坐标。
步骤S302,采用中值滤波器对图像进行滤波处理,把数字图像或数字序列中一点的值用该点的一个领域中各点值的中值代替,让周围的像素值接近真实值,这样可以有效消除鼓励噪声点;中值滤波器将图像中每一点的像素值由对应滤波区域R内像素值的中值代替。
I′(u,v)←median{I(u+i,v+j)|(i,j)∈R} (6)
2K+1个像素点pi的中值定义为:
步骤S303,提取裂缝特征,并根据裂缝特征对图像进行区域划分,标定出关键裂缝的图像数据。
采用Otsu阈值分割算法来有效标定图像中的裂缝特征,通过在图像的灰度取值范围内选择一灰度阈值,并对全部像素点应用该阈值。选取阈值是为了能够区分裂缝与图像背景,通过灰度方差拉开距离,这样就可以精确的提取裂缝信息。
具体地,用C0表示裂缝图像,C1表示背景,假设图像像素总数为N,灰度区间[0,L-1],那么C0和C1类方差分别为:
其中T为阈值,w0,w1,u0,u1的计算分别如下:
pi=ni/N (10)
其中,ni代表在灰度区间内灰度级为i的像素点数;
C0和C1类的类间方差为:
C0和C1类的类间方差为:
通过让T在[0,L-1]范围内依次取值,使式(18)中T取得最大值,求解Otsu算法最佳阈值T;以所获得的阈值T为裂缝特征进行目标特征提取,从而获得裂缝信息。
步骤S304,采用Canny算子进行亚像素级边缘检测,提取出裂缝与背景图像中的交界线。
具体包括:
步骤S3041,用高斯滤波器进行图像平滑处理;
步骤S3042,用一阶偏导数的有限差分来计算梯度的幅值和方向;
步骤S3043,对梯度复制应用非极大值抑制;
步骤S3044,用双阈值算法检测和连接边缘。
步骤S305,根据所述裂缝与背景图像中的交界线,获取裂缝形状与变化率。
本步骤中,依据当前帧经过边缘检测后的裂缝像素数目和此段确定的时间后的另外一帧裂缝像素数目的差来计算裂缝在一定时间内的变化率。在裂缝视觉图像中,每一帧的图像的裂缝形体去除背景区域后所占区域都有所不同,裂缝的形态变化都有不一样的像素总值,通过像素的变化来反映裂缝的变化情况,再利用相应的算法处理来计算出裂缝的具体变化物理量。
具体过程包括:
设Xi为第ti时刻裂缝区域的像素总值,则Xi+1为第ti+1时刻裂缝区域的像素总值,连续的视觉数据是以前后帧数来计算,因此我们就可以得出在第i帧的裂缝目标区域像素总值可设为Xi,那么第i+1帧的像素总值为Yi+1,令第i帧的裂缝面积为Li,第i+1帧区域面积为Li+1,因此就可以得出在连续视觉图像帧下的裂缝面积变化率:
式(19)中,j表示第j帧,i+j表示第i+j帧,ΔB表示裂缝的面积变化率,限定于在时间范围ti到ti+1这一间断时间内。
步骤S4,基于时间序列构建,分析关键裂缝的变形,并进行灾害预警。
本步骤中,采用如下过程构建变形监测及预警模型,进行灾害预警:
步骤S401,构建自回归滑动平均(ARIMA)模型。
对于平稳、正态、零均值的时间序列{Xt},若{Xt}的取值不仅与前p步的各个取值相关,而且还与前q步的各个干扰项相关,构建ARIMA模型为:
经过图像处理和特征提取后的裂缝视觉数据不具有平稳性,属于非平稳时间序列模型(ARIMA(p,d,q)),再采用差分法对序列进行d阶差分形成稳定的时间序列模型,差分公式如下:
步骤S402,计算自相关函数和偏自相关函数及特点分析,选取p,d和q值,优化ARIMA(p,d,q)模型。
所述自相关函数和偏自相关函数计算公式如下:
不断修正p,d和q值使自相关系数和偏相关系数随时滞k衰减呈正弦波或指数并趋于零,称为“拖尾”,则认为{Xt}序列是ARIMA序列。
步骤S403,采用最小二乘法计算ARIMA模型参数,将(20)式改写为如下:
将(24)式表示为向量形式:X=F(Y,β)+α
其中,Y=[Xp+1,Xp+2,···,XN]T,F=[fp+1,fp+2,···,fN]T,a=[ep+1,ep+2,···,eN]T;Y表示裂缝数据值X1,X2,···,XN;β表示待求参数。
步骤S404,采用Ljung-Box对优化后的优化ARIMA(p,d,q)模型进行检验。
本步骤中检验公式如下:
步骤S405,采用优化并检验后的ARIMA(p,d,q)模型,根据历史获取的裂缝监测数据,即X1,X2,···,Xt-1,Xt来推测未来T期的Xt+T,用Xt(T)代表预测值,最小均方误差预测公式如下:
Xt(T)=E(Xt+T|X1,X2,···,Xt-1,Xt) (26)
其中,E代表的是最小均方差公式。
将Xt(T)的预测值设置不同的预警预报等级,当达到预设的预警阈值时,可以根据相应的预警等级发出相应的预警信号。在本发明的一个具体实施例中,当预测值为0mm≤Xt(T)<10mm,发出蓝色预警;当预测值为10mm≤Xt(T)<25mm,发出黄色预警;当预测值为25mm≤Xt(T)<40mm,发出橙色预警;当预测值为40mm≤Xt(T)<50mm,发出红色预警。
由以上技术方案可以看出,本发明实施方式所提供的崩塌滑坡灾害体的变形监测方法,在获得灾害体图像数据的基础上,利用亚像素级边缘检测算法,结合Otsu阈值分割算法及Canny算子,有效提取崩塌滑坡灾害体重点变形区物理图像像素,进而计算出关键裂缝的实际变化量,对关键裂缝进行量化监测、分析处理;再基于变形监测预警的时间序列模型,并结合坐标系间的转换关系,在有效时间内对崩塌滑坡灾害体的变形趋势进行判断和及时预警,提高了崩塌滑坡灾害体变形监测的实时性、准确性和稳定性,数据处理计算量小、适应性强,节约了监测资源。
以上所述是本发明的优选实施方式,应当指出,本发明并不受限于以上所公开的示范性实施例,说明书的实质仅仅是帮助相关领域技术人员综合理解本发明的具体细节。对于本技术领域的普通技术人员来说,在不脱离本发明所述原理的前提下,在本发明揭露的技术范围做出的若干改进和润饰、可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。
Claims (10)
1.一种崩塌滑坡灾害体的变形监测方法,其特征在于,包括如下步骤:
步骤S1,构建世界坐标系,并以世界坐标系为基准坐标系,构建视觉传感器坐标系;
步骤S2,采用激光可视化视觉传感器对当前崩塌滑坡灾害体进行实时侦测,并基于视觉传感器坐标系采集灾害体中的图像数据,以供关键裂缝识别模块调用;
步骤S3,关键裂缝识别模块从所述视觉传感器中读取基于视觉传感器坐标系的灾害体图像数据,对基于图像数据的二维数字信号进行目标特征提取,获得关键裂缝的图像数据,计算所述关键裂缝的形状及实际变化量,并基于世界坐标系恢复所述关键裂缝的三维地理坐标;
步骤S4,基于时间序列构建变形监测模型,分析关键裂缝的变形,结合三维地理坐标进行灾害预警。
2.根据权利要求1所述的崩塌滑坡灾害体的变形监测方法,其特征在于,所述步骤S3中所述对基于图像数据的二维数字信号进行目标特征提取,获得所述灾害体图像数据中的关键裂缝的图像数据,计算所述裂缝的形状及实际变化量,具体包括:
步骤S301,将采集到的彩色图像转化为灰度图像;
步骤S302,采用中值滤波器对图像进行滤波处理,把数字图像或数字序列中一点的值用该点的一个领域中各点值的中值代替,让周围的像素值接近真实值,消除噪声点;
步骤S303,从滤波后的图像数据中提取裂缝特征,并根据裂缝特征对图像进行分割,标定出关键裂缝的图像数据;
步骤S304,采用Canny算子进行亚像素级边缘检测,提取出裂缝与背景图像中的交界线;
步骤S305,根据所述裂缝与背景图像中的交界线,获取裂缝形状与变化率。
3.根据权利要求2所述的崩塌滑坡灾害体的变形监测方法,其特征在于,所述步骤S301,进一步为:利用加权平均法将彩色图像进行灰度转化,根据彩色图像中R、G、B三个分量进行不同权值的加权平均,采用式(5)得到灰度图像:
f(i,j)=0.30R(i,j)+0.59G(i,j)+0.11B(i,j) (5)
式(5)中,(i,j)是图像中的像素坐标。
4.根据权利要求2所述的崩塌滑坡灾害体的变形监测方法,其特征在于,所述步骤S303中对图像进行分割,采用Otsu阈值分割算法来有效标定图像中的裂缝特征,通过在图像的灰度取值范围内选择灰度阈值,并对全部像素点应用该灰度阈值。
7.根据权利要求1至6任一项所述的崩塌滑坡灾害体的变形监测方法,其特征在于,所述步骤S4中构建变形监测及预警模型,具体包括如下步骤:
步骤S401,构建自回归滑动平均ARIMA模型;
步骤S402,计算自相关函数和偏自相关函数并进行特点分析,优化ARIMA模型;
步骤S403,基于优化后的ARIMA模型,采用最小二乘法计算ARIMA模型参数;
步骤S404,采用Ljung-Box对优化后的ARIMA模型进行检验,获得最终变形监测与预警模型。
9.根据权利要求8所述的崩塌滑坡灾害体的变形监测方法,其特征在于,所述步骤S4进行灾害预警,采用优化并检验后的ARIMA(p,d,q)模型,根据历史获取的裂缝监测数据X1,X2,…,Xt-1,Xt推测未来T期的Xt+T,用Xt(T)代表预测值,最小均方误差预测公式如下:
Xt(T)=E(Xt+T|X1,X2,…,Xt-1,Xt) (26)。
10.根据权利要求1所述的崩塌滑坡灾害体的变形监测方法,其特征在于,所述视觉传感器坐标系是固定在传感器上的直角坐标系,原点定义在传感器光心,光心到图像平面的距离为视觉传感器有效焦距。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011300788.XA CN112364802B (zh) | 2020-11-19 | 2020-11-19 | 一种崩塌滑坡灾害体的变形监测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011300788.XA CN112364802B (zh) | 2020-11-19 | 2020-11-19 | 一种崩塌滑坡灾害体的变形监测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112364802A true CN112364802A (zh) | 2021-02-12 |
CN112364802B CN112364802B (zh) | 2021-08-03 |
Family
ID=74533689
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011300788.XA Active CN112364802B (zh) | 2020-11-19 | 2020-11-19 | 一种崩塌滑坡灾害体的变形监测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112364802B (zh) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113705108A (zh) * | 2021-09-02 | 2021-11-26 | 南方科技大学 | 实时滑坡灾害监测预警方法及系统 |
CN114494979A (zh) * | 2022-03-29 | 2022-05-13 | 杭州定川信息技术有限公司 | 一种视频识别生态流量泄放的方法 |
CN114676907A (zh) * | 2022-01-17 | 2022-06-28 | 中国地质大学(北京) | 一种区域地质灾害预警方法、装置、存储介质及设备 |
CN114782845A (zh) * | 2022-06-20 | 2022-07-22 | 青岛理工大学 | 爆破后建筑物新生裂缝识别方法、装置、设备及介质 |
CN114812491A (zh) * | 2022-01-19 | 2022-07-29 | 国网甘肃省电力公司经济技术研究院 | 基于长时间序列分析的输电线路地表形变预警方法及装置 |
CN115035182A (zh) * | 2022-06-06 | 2022-09-09 | 桂林理工大学 | 一种山区滑坡灾害预警方法及系统 |
CN115130356A (zh) * | 2022-09-01 | 2022-09-30 | 北京云庐科技有限公司 | 一种基于数字孪生技术的崩塌监测系统与方法 |
CN115410342A (zh) * | 2022-08-26 | 2022-11-29 | 安徽省地质矿产勘查局332地质队 | 一种基于裂缝计实时监测的滑坡灾害智能预警方法 |
CN115938095A (zh) * | 2023-02-22 | 2023-04-07 | 湖北通达数科科技有限公司 | 一种基于集成融合模型的山体滑坡监测预警方法及系统 |
CN116757335A (zh) * | 2023-08-17 | 2023-09-15 | 四川省华地建设工程有限责任公司 | 基于北斗卫星的崩塌预测方法及系统 |
CN117493833A (zh) * | 2023-12-29 | 2024-02-02 | 江西飞尚科技有限公司 | 一种滑坡变形阶段识别方法、系统、存储介质及计算机 |
CN117765689A (zh) * | 2024-02-22 | 2024-03-26 | 山东光合云谷大数据有限公司 | 一种地质灾害监测预警系统 |
CN117877212A (zh) * | 2024-03-12 | 2024-04-12 | 四川省自然资源勘察设计集团有限公司 | 一种补充式视觉化灾害监测报警系统 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104732519A (zh) * | 2015-01-20 | 2015-06-24 | 中国科学院半导体研究所 | 鲁棒的全局阈值分割方法 |
CN105336119A (zh) * | 2015-11-27 | 2016-02-17 | 中国地质调查局水文地质环境地质调查中心 | 一种崩塌体裂缝识别报警装置及其识别方法 |
CN106126483A (zh) * | 2016-06-21 | 2016-11-16 | 湖北天明气和网络科技有限公司 | 一种天气预测的方法及装置 |
CN107179322A (zh) * | 2017-06-15 | 2017-09-19 | 长安大学 | 一种基于双目视觉的桥梁底面裂缝检测方法 |
CN108416985A (zh) * | 2018-04-20 | 2018-08-17 | 北京航天控制仪器研究所 | 一种基于图像识别的地质灾害监测预警系统及预警方法 |
CN109376919A (zh) * | 2018-10-12 | 2019-02-22 | 西安科技大学 | 一种煤矿综采工作面瓦斯涌出量的预测方法 |
CN109584240A (zh) * | 2018-12-20 | 2019-04-05 | 成都理工大学 | 滑坡后缘裂缝位移图像识别方法 |
CN110514113A (zh) * | 2019-06-13 | 2019-11-29 | 杭州电子科技大学 | 一种基于单目视觉摄像头的滑坡裂缝监测方法 |
CN111473776A (zh) * | 2020-05-11 | 2020-07-31 | 中晋环境科技有限公司 | 一种基于单像近景摄影测量的滑坡裂缝监测方法 |
-
2020
- 2020-11-19 CN CN202011300788.XA patent/CN112364802B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104732519A (zh) * | 2015-01-20 | 2015-06-24 | 中国科学院半导体研究所 | 鲁棒的全局阈值分割方法 |
CN105336119A (zh) * | 2015-11-27 | 2016-02-17 | 中国地质调查局水文地质环境地质调查中心 | 一种崩塌体裂缝识别报警装置及其识别方法 |
CN106126483A (zh) * | 2016-06-21 | 2016-11-16 | 湖北天明气和网络科技有限公司 | 一种天气预测的方法及装置 |
CN107179322A (zh) * | 2017-06-15 | 2017-09-19 | 长安大学 | 一种基于双目视觉的桥梁底面裂缝检测方法 |
CN108416985A (zh) * | 2018-04-20 | 2018-08-17 | 北京航天控制仪器研究所 | 一种基于图像识别的地质灾害监测预警系统及预警方法 |
CN109376919A (zh) * | 2018-10-12 | 2019-02-22 | 西安科技大学 | 一种煤矿综采工作面瓦斯涌出量的预测方法 |
CN109584240A (zh) * | 2018-12-20 | 2019-04-05 | 成都理工大学 | 滑坡后缘裂缝位移图像识别方法 |
CN110514113A (zh) * | 2019-06-13 | 2019-11-29 | 杭州电子科技大学 | 一种基于单目视觉摄像头的滑坡裂缝监测方法 |
CN111473776A (zh) * | 2020-05-11 | 2020-07-31 | 中晋环境科技有限公司 | 一种基于单像近景摄影测量的滑坡裂缝监测方法 |
Non-Patent Citations (1)
Title |
---|
周迪 等: ""一种改进的OTSU阈值分割算法"", 《中国计量大学学报》 * |
Cited By (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113705108A (zh) * | 2021-09-02 | 2021-11-26 | 南方科技大学 | 实时滑坡灾害监测预警方法及系统 |
CN114676907A (zh) * | 2022-01-17 | 2022-06-28 | 中国地质大学(北京) | 一种区域地质灾害预警方法、装置、存储介质及设备 |
CN114812491A (zh) * | 2022-01-19 | 2022-07-29 | 国网甘肃省电力公司经济技术研究院 | 基于长时间序列分析的输电线路地表形变预警方法及装置 |
CN114812491B (zh) * | 2022-01-19 | 2023-10-20 | 国网甘肃省电力公司经济技术研究院 | 基于长时间序列分析的输电线路地表形变预警方法及装置 |
CN114494979A (zh) * | 2022-03-29 | 2022-05-13 | 杭州定川信息技术有限公司 | 一种视频识别生态流量泄放的方法 |
CN114494979B (zh) * | 2022-03-29 | 2022-07-22 | 杭州定川信息技术有限公司 | 一种视频识别生态流量泄放的方法 |
CN115035182A (zh) * | 2022-06-06 | 2022-09-09 | 桂林理工大学 | 一种山区滑坡灾害预警方法及系统 |
CN115035182B (zh) * | 2022-06-06 | 2024-03-29 | 桂林理工大学 | 一种山区滑坡灾害预警方法及系统 |
CN114782845A (zh) * | 2022-06-20 | 2022-07-22 | 青岛理工大学 | 爆破后建筑物新生裂缝识别方法、装置、设备及介质 |
CN115410342B (zh) * | 2022-08-26 | 2023-08-11 | 安徽省地质矿产勘查局332地质队 | 一种基于裂缝计实时监测的滑坡灾害智能预警方法 |
CN115410342A (zh) * | 2022-08-26 | 2022-11-29 | 安徽省地质矿产勘查局332地质队 | 一种基于裂缝计实时监测的滑坡灾害智能预警方法 |
CN115130356A (zh) * | 2022-09-01 | 2022-09-30 | 北京云庐科技有限公司 | 一种基于数字孪生技术的崩塌监测系统与方法 |
CN115938095A (zh) * | 2023-02-22 | 2023-04-07 | 湖北通达数科科技有限公司 | 一种基于集成融合模型的山体滑坡监测预警方法及系统 |
CN116757335A (zh) * | 2023-08-17 | 2023-09-15 | 四川省华地建设工程有限责任公司 | 基于北斗卫星的崩塌预测方法及系统 |
CN116757335B (zh) * | 2023-08-17 | 2023-10-20 | 四川省华地建设工程有限责任公司 | 基于北斗卫星的崩塌预测方法及系统 |
CN117493833A (zh) * | 2023-12-29 | 2024-02-02 | 江西飞尚科技有限公司 | 一种滑坡变形阶段识别方法、系统、存储介质及计算机 |
CN117493833B (zh) * | 2023-12-29 | 2024-04-09 | 江西飞尚科技有限公司 | 一种滑坡变形阶段识别方法、系统、存储介质及计算机 |
CN117765689A (zh) * | 2024-02-22 | 2024-03-26 | 山东光合云谷大数据有限公司 | 一种地质灾害监测预警系统 |
CN117765689B (zh) * | 2024-02-22 | 2024-05-10 | 山东光合云谷大数据有限公司 | 一种地质灾害监测预警系统 |
CN117877212A (zh) * | 2024-03-12 | 2024-04-12 | 四川省自然资源勘察设计集团有限公司 | 一种补充式视觉化灾害监测报警系统 |
CN117877212B (zh) * | 2024-03-12 | 2024-06-14 | 四川省自然资源勘察设计集团有限公司 | 一种补充式视觉化灾害监测报警系统 |
Also Published As
Publication number | Publication date |
---|---|
CN112364802B (zh) | 2021-08-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112364802B (zh) | 一种崩塌滑坡灾害体的变形监测方法 | |
CN112598881B (zh) | 一种地质灾害监测方法、装置及计算机设备 | |
CN113611082B (zh) | 一种无人机铁路边坡监控与预警系统及方法 | |
Lubowiecka et al. | Historic bridge modelling using laser scanning, ground penetrating radar and finite element methods in the context of structural dynamics | |
Tarpanelli et al. | Toward the estimation of river discharge variations using MODIS data in ungauged basins | |
CN114353876A (zh) | 一种黄土公路边坡健康监测方法 | |
Wu et al. | Rapid safety monitoring and analysis of foundation pit construction using unmanned aerial vehicle images | |
CN108663017A (zh) | 一种监测城市地铁沿线地表沉降的方法 | |
CN113885025A (zh) | 一种滑坡形变监测方法与可视化服务平台 | |
CN114117913A (zh) | 一种坝后露天渗流量监测系统及测值归因分析方法 | |
CN106526125A (zh) | 双变坡侵蚀槽和壤中流测定仪器及其测定方法 | |
CN116429070A (zh) | 采煤沉陷积水区水资源调查方法及系统 | |
CN111473770A (zh) | 一种基于综合属性测量的边坡状况评定方法 | |
CN112270675B (zh) | 一种基于极化雷达遥感影像的城市内涝区域检测方法 | |
CN104914483B (zh) | 一种基于多源观测数据的南极接地线质量评价方法 | |
Abd Elbasit et al. | Potential of low cost close‐range photogrammetry system in soil microtopography quantification | |
Barbarella et al. | Multi-temporal terrestrial laser scanning survey of a landslide | |
Liu et al. | Architecture planning and geo-disasters assessment mapping of landslide by using airborne LiDAR data and UAV images | |
Sui et al. | Application of 3D laser scanner for monitoring of landslide hazards | |
Liu et al. | Development of building height data in Peru from high-resolution SAR imagery | |
CN113625241A (zh) | 差异沉降监控预警方法 | |
CN113744393A (zh) | 一种多层级边坡滑坡变化监测方法 | |
CN113670454A (zh) | 一种自动观测气象卫星热红外通道在轨场地辐射定标方法 | |
Nadi et al. | Analysis image-based automated 3D crack detection for post-disaster bridge assessment in flyover Mall Boemi Kedaton | |
CN117516487B (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 |