CN114882366A - 一种三维场景灾变监控预警方法 - Google Patents

一种三维场景灾变监控预警方法 Download PDF

Info

Publication number
CN114882366A
CN114882366A CN202210580210.7A CN202210580210A CN114882366A CN 114882366 A CN114882366 A CN 114882366A CN 202210580210 A CN202210580210 A CN 202210580210A CN 114882366 A CN114882366 A CN 114882366A
Authority
CN
China
Prior art keywords
point cloud
cloud data
area
monitoring
dem
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
CN202210580210.7A
Other languages
English (en)
Other versions
CN114882366B (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.)
Guangzhou Urban Planning Survey and Design Institute
Original Assignee
Guangzhou Urban Planning Survey and Design Institute
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 Guangzhou Urban Planning Survey and Design Institute filed Critical Guangzhou Urban Planning Survey and Design Institute
Priority to CN202210580210.7A priority Critical patent/CN114882366B/zh
Publication of CN114882366A publication Critical patent/CN114882366A/zh
Application granted granted Critical
Publication of CN114882366B publication Critical patent/CN114882366B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V20/00Scenes; Scene-specific elements
    • G06V20/10Terrestrial scenes
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology
    • G06N3/045Combinations of networks
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/08Learning methods
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/05Geographic models
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/70Arrangements for image or video recognition or understanding using pattern recognition or machine learning
    • G06V10/82Arrangements for image or video recognition or understanding using pattern recognition or machine learning using neural networks
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10Information 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)
  • Theoretical Computer Science (AREA)
  • Software Systems (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Health & Medical Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Computing Systems (AREA)
  • Biophysics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Data Mining & Analysis (AREA)
  • General Engineering & Computer Science (AREA)
  • Computational Linguistics (AREA)
  • Mathematical Physics (AREA)
  • Biomedical Technology (AREA)
  • Molecular Biology (AREA)
  • Geometry (AREA)
  • Multimedia (AREA)
  • Medical Informatics (AREA)
  • Databases & Information Systems (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Remote Sensing (AREA)
  • Computer Graphics (AREA)
  • Image Processing (AREA)

Abstract

本发明公开了一种三维场景灾变监控预警方法,包括:基于改进的种子区域生长算法,对原始图像进行分割处理,得到分割后的原始图像;根据分割后的原始图像对预先构建的泥石流图像识别模型进行训练,得到训练好的泥石流图像识别模型;每隔预设时间段对监测区域进行一次激光扫描,得到当前时刻的点云数据,并基于GNSS导航系统和惯性导航系统,对点云数据进行定位,得到标准激光点云数据;基于标准激光点云数据得到初始DEM;对初始DEM进行滤波得到最终DEM;利用训练好的泥石流图像识别模型和最终DEM对所述监测区域进行地质灾害监测。采用本发明实施例能够实现三维场景灾变的监控预警,并在提高精度的同时降低成本。

Description

一种三维场景灾变监控预警方法
技术领域
本发明涉及图像处理技术领域,尤其涉及一种三维场景灾变监控预警方法。
背景技术
泥石流作为一种大规模地质灾害,极具破坏性,一旦发生,往往会带来巨大的经济损失和人员伤亡。因此,快速准确识别泥石流的发生,对于快速采取有效措施、精准制定应对策略有着重要意义。
现有技术中,一般通过人工法来进行泥石流灾害的监控预警。人工法是一种宏观的测量方法,依靠技术人员根据多次的现场观测情况进行比对实现,技术人员必须具有丰富的经验,人工检测方法单一,数据误差大,人力成本高。
发明内容
本发明提供一种三维场景灾变监控预警方法,通过结合图像识别和激光点云进行泥石流灾害的监控预警,能够在提高精度的同时降低成本。
为实现上述目的,本发明实施例提供了一种三维场景灾变监控预警方法,包括:
获取监测区域的原始图像;
基于改进的种子区域生长算法,对所述原始图像进行分割处理,得到分割后的原始图像;
根据所述分割后的原始图像对预先构建的泥石流图像识别模型进行训练,得到训练好的泥石流图像识别模型;
每隔预设时间段对所述监测区域进行一次激光扫描,得到当前时刻的点云数据,并基于GNSS导航系统和惯性导航系统,对所述当前时刻的点云数据进行定位,得到当前时刻的标准激光点云数据;
基于所述当前时刻的标准激光点云数据得到当前时刻的初始DEM;
对所述当前时刻的初始DEM进行滤波得到当前时刻的最终DEM;
利用所述训练好的泥石流图像识别模型和所述当前时刻的最终DEM对所述监测区域进行地质灾害的监控预警。
作为是上述方案的改进,所述基于改进的种子区域生长算法,对所述原始图像进行分割处理,得到分割后的原始图像,包括:
将所述原始图像进行区域划分,并从中随机选择一个区域作为种子区域;
对所述种子区域周围的每一邻接区域,判断所述邻接区域的平均颜色分量值和所述种子区域的平均颜色分量值的比较结果是否小于预设阈值,若是,将所述邻接区域与所述种子区域合并,得到合并后的区域;
选取一个平均颜色分量值和所述种子区域的平均颜色分量值相同的区域,作为所述种子区域,重复步骤对所述种子区域周围的每一邻接区域,判断所述邻接区域的平均颜色分量值和所述种子区域的平均颜色分量值的比较结果是否小于预设阈值,若是,将所述邻接区域与所述种子区域合并,得到合并后的区域,直到遍历完所述原始图像上所有的区域,得到分割后的原始图像。
作为是上述方案的改进,所述泥石流图像识别模型的构建过程至少包括:卷积操作、池化操作、ReLU变化、Dropout、归一化。
其中,根据下式进行卷积操作:
Figure BDA0003663480980000021
式中,
Figure BDA0003663480980000022
表示第l层卷积层的第j个特征图,f()表示激活函数,Mj表示上一层的输出集合,
Figure BDA0003663480980000023
表示卷积核,
Figure BDA0003663480980000024
表示偏置参数,i表示次数,*表示卷积操作;
根据下式进行池化操作:
Figure BDA0003663480980000025
式中,down()表示采样函数,
Figure BDA0003663480980000026
表示每个输出特征图的多倍偏置参数,
Figure BDA0003663480980000027
表示偏置参数;
根据下式进行ReLU变化:
yijk=max{0,xijk}
式中,xijk表示输入值,yijk表示输出值;
在Dropout中,将50%的隐藏层的输出设置为0;
根据下式进行归一化:
Figure BDA0003663480980000031
Figure BDA0003663480980000032
式中,p表示相邻信道输入特征群,α表示系列参数,k表示第k个核,xijk表示第k个核在位置(i,j)运用ReLU变化的非线性化输出,G(k)表示通道数,beta表示预设参数,G(k′)表示通道数。
作为是上述方案的改进,所述泥石流图像识别模型的训练过程至少包括反向传播。
其中,所述反向传播包括卷积层梯度计算和池化层梯度计算:
根据下式进行卷积层梯度计算:
Figure BDA0003663480980000033
式中,
Figure BDA0003663480980000034
表示节点j在l层的第一个导函数误差,β表示输入系数,°表示每个部分相乘,
Figure BDA0003663480980000035
表示对于节点j的偏导值的激活函数,up()表示采样函数,u表示输入节点,δ表示导函数误差;
Figure BDA0003663480980000036
式中,
Figure BDA0003663480980000037
表示
Figure BDA0003663480980000038
Figure BDA0003663480980000039
的转置,conv2表示一个卷积函数,rot180表示卷积核被转置;
根据下式进行池化层梯度计算:
Figure BDA00036634809800000310
式中,
Figure BDA00036634809800000311
作为是上述方案的改进,所述基于所述当前时刻的标准激光点云数据得到当前时刻的初始DEM,包括:
将所述当前时刻的标准激光点云数据进行自适应分块处理,得到每一分块的点云数据;
基于主成分分析法,对在x方向和y方向的所述每一分块的点云数据进行局部平面拟合处理,并利用反距离加权法和样条加权法对经过局部平面拟合处理过的点云数据进行插值处理,得到初始拟合曲面,对在z方向的所述每一分块的点云数据和所述初始拟合曲面进行滤波处理,得到最终拟合曲面;
对所述最终拟合曲面进行滤波处理,得到当前时刻的初始DEM。
作为是上述方案的改进,所述对所述当前时刻的初始DEM进行滤波得到当前时刻的最终DEM,包括:
采用小波函数来对所述当前时刻的初始DEM进行滤波,得到当前时刻的最终DEM。
与现有技术相比,本发明实施例提供的一种三维场景灾变监控预警方法,通过利用改进的种子区域生长算法识别监测区域发生泥石流的区域、激光点云数据构建表征监测区域的三维场景的数字高程模型(数字高程模型,Digital Elevation Model,DEM),能够实现三维场景灾变监控预警。由此可见,本发明实施例通过结合基于训练好的泥石流图像识别模型的图像识别和激光点云对监测区域进行地质灾害的监控预警,能够在提高精度的同时降低成本。
附图说明
图1是本发明实施例提供的一种三维场景灾变监控预警方法的流程图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
参见图1,图1是本发明实施例提供的一种三维场景灾变监控预警方法的流程图,所述三维场景灾变监控预警方法,包括:
S1、获取监测区域的原始图像;
S2、基于改进的种子区域生长算法,对所述原始图像进行分割处理,得到分割后的原始图像;
S3、根据所述分割后的原始图像对预先构建的泥石流图像识别模型进行训练,得到训练好的泥石流图像识别模型;
S4、每隔预设时间段对所述监测区域进行一次激光扫描,得到当前时刻的点云数据,并基于GNSS导航系统和惯性导航系统,对所述当前时刻的点云数据进行定位,得到当前时刻的标准激光点云数据;
S5、基于所述当前时刻的标准激光点云数据得到当前时刻的初始DEM;
S6、对所述当前时刻的初始DEM进行滤波得到当前时刻的最终DEM;
S7、利用所述训练好的泥石流图像识别模型和所述当前时刻的最终DEM对所述监测区域进行地质灾害的监控预警。
具体地,在步骤S2中,所述基于改进的种子区域生长算法,对所述原始图像进行分割处理,得到分割后的原始图像,包括:
S21、将所述原始图像进行区域划分,并从中随机选择一个区域作为种子区域;
S22、对所述种子区域周围的每一邻接区域,判断所述邻接区域的平均颜色分量值和所述种子区域的平均颜色分量值的比较结果是否小于预设阈值,若是,将所述邻接区域与所述种子区域合并,得到合并后的区域;
可以理解的是,考虑种子区域S0周围的八个邻接区域S,如果S0和S满足生长标准,即满足平均颜色值相同或相近,则将它们合并到同一个区域。
S23、选取一个平均颜色分量值和所述种子区域的平均颜色分量值相同的区域,作为所述种子区域,重复步骤S21对所述种子区域周围的每一邻接区域,判断所述邻接区域的平均颜色分量值和所述种子区域的平均颜色分量值的比较结果是否小于预设阈值,若是,将所述邻接区域与所述种子区域合并,得到合并后的区域,直到遍历完所述原始图像上所有的区域,得到分割后的原始图像。
本发明实施例并不采用将区域中分割出一个种子点作为生长起点进行区域性增长的传统种子区域生长算法,而是选择一个区域作为生长起点,与邻接区域进行合并的方法,能够更改准确的提取区域和位置信息,同时使用平均颜色作为参考值,能够克服图像不均匀和不稳定的特征。
具体地,在步骤S3中,所述泥石流图像识别模型的构建过程至少包括:卷积操作、池化操作、ReLU变化、Dropout、归一化。
具体地,根据下式进行卷积操作:
Figure BDA0003663480980000061
其中,
Figure BDA0003663480980000062
表示第l层卷积层的第j个特征图,f()表示激活函数,Mj表示上一层的输出集合,
Figure BDA0003663480980000063
表示卷积核,
Figure BDA0003663480980000064
表示偏置参数,i表示次数,*表示卷积操作;
根据下式进行池化操作:
Figure BDA0003663480980000065
其中,down()表示采样函数,
Figure BDA0003663480980000066
表示每个输出特征图的多倍偏置参数,
Figure BDA0003663480980000067
表示偏置参数。
通常最大化池化是基于卷积特征的提取,通过卷积隐藏节点,不断降低特征维度,从而降低分类器的负担,通常采用以下方法:
yijk=max{x″ij:i≤i′<i+p,j≤j′≤j+p}
其中,p表示窗口大小,yijk表示(i,j)在采样层的最大输出值,x″ij表示输入值,(i,j)表示节点,i′,j′分别表示在池化过程中的不同节点。
根据下式进行ReLU变化:
yijk=max{0,xijk}
其中,xijk表示ReLU变化的输入值,yijk表示ReLU变化的输出值。
本发明采用ReLU作为激活函数,主要原因是因为sigmoid函数计算误差反向传播梯度需要大量计算,而ReLU更易实现;sigmoid函数在深度网络的反向推导中,容易出现误差梯度消失的问题;ReLU生成的是稀疏网络,减少了参数之间的相互依赖关系从而缓解了网络深度训练的过拟合问题。
Dropout具体包括:
对于50%隐藏层的输出设置为0,并且不对前面和后续有任何影响,采用这种方法得到的参数能够适应不同情况下的网络结构,提高了系统的泛化能力。
根据下式进行归一化:
Figure BDA0003663480980000071
Figure BDA0003663480980000072
其中,p表示相邻信道输入特征群,α表示系列参数,k表示第k个核,xijk表示第k个核在位置(i,j)运用ReLU变化的非线性化输出,G(k)表示通道数,beta为预设参数,G(k′)表示通道数。
在局部响应归一化层进行归一化处理即执行“近抑制”操作,用以对局部输入区域规范化,在相邻信道范围内,归一化算子的形式如上式所示。
具体地,在步骤S3中,所述泥石流图像识别模型的训练过程至少包括反向传播。
具体地,所述反向传播包括卷积层梯度计算和池化层梯度计算。
卷积层梯度计算:
在假设l是通过卷积层l+1进行采样,对于所有的输入节点u,δ代表导函数误差,具体公式如下:
Figure BDA0003663480980000073
其中,
Figure BDA0003663480980000074
表示节点j在l层的第一个导函数误差,β为输入系数,°表示每个部分相乘,
Figure BDA0003663480980000075
表示对于节点j的偏导值的激活函数,up()代表采样函数。
卷积梯度和连接点的梯度相关,计算公式如下:
Figure BDA0003663480980000081
其中,
Figure BDA0003663480980000082
表示
Figure BDA0003663480980000083
Figure BDA0003663480980000084
的转置,conv2表示一个卷积函数,rot表示卷积核被转置。
池化层梯度计算:
和计算卷积层梯度类似,池化层的梯度计算为:
Figure BDA0003663480980000085
其中,
Figure BDA0003663480980000086
具体地,在步骤S4中,每次激光扫描的时候要保证每次飞行的标准一致,作业流程、航带路线保持不变,从而确保每一时刻的点云数据之间坐标基准的统一性。
具体地,在步骤S5中,所述基于所述当前时刻的标准激光点云数据得到当前时刻的初始DEM,包括:
S51、将所述当前时刻的标准激光点云数据进行自适应分块处理,得到每一分块的点云数据;
具体地,在Cyclone6.0中按照Intensity字段来识别不同监测区域的激光点云,使用点云自适应分块算法将标准激光点云数据划分为若干点数近似的子块;
S52、基于主成分分析法,对在x方向和y方向的所述每一分块的点云数据进行局部平面拟合处理,并利用反距离加权法和样条加权法对经过局部平面拟合处理过的点云数据进行插值处理,得到初始拟合曲面,对在z方向的所述每一分块的点云数据和所述初始拟合曲面进行滤波处理,得到最终拟合曲面;
可以理解的是,这一过程是为了对点云数据进行去噪,提高点云数据的精细度。首先是利用主成分分析法来找到两个主成分,在两个主成分(x,y)方向上进行局部平面拟合处理,能够去除较大的噪声,然后利用反距离加权法和样条加权法对初步去噪后的点云数据(即经过局部平面拟合处理过的点云数据)进行插值处理,得到初始拟合曲面,最后基于第三个主成分z方向和初始拟合曲面进行滤波去噪,得到最终拟合曲面。
S53、对所述最终拟合曲面进行滤波处理,得到当前时刻的初始DEM。
具体地,通过滤波算法对所述最终拟合曲面进行滤波,得到当前时刻的初始DEM。
具体地,在步骤S6中,所述对所述当前时刻的初始DEM进行滤波得到当前时刻的最终DEM,包括:
采用小波函数来对所述当前时刻的初始DEM进行滤波,得到当前时刻的最终DEM。
具体地,在步骤S7中,将监测区域的实时图像输入至训练好的卷积神经网络,能够实时识别监测区域是否发生泥石流及发生泥石流的实时位置;将当前时刻的最终DEM和上一时刻的最终DEM进行比较,能够实时识别监测区域是否发生泥石流及发生泥石流的实时位置。
本发明实施例所提供的一种三维场景灾变监控预警方法,通过利用改进的种子区域生长算法识别监测区域发生泥石流的区域、激光点云数据构建表征监测区域的三维场景的数字高程模型,能够实现三维场景灾变监控预警。由此可见,本发明实施例通过结合基于训练好的泥石流图像识别模型的图像识别和激光点云对监测区域进行地质灾害的监控预警,能够在提高精度的同时降低成本。
以上所述是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也视为本发明的保护范围。

Claims (6)

1.一种三维场景灾变监控预警方法,其特征在于,包括:
获取监测区域的原始图像;
基于改进的种子区域生长算法,对所述原始图像进行分割处理,得到分割后的原始图像;
根据所述分割后的原始图像对预先构建的泥石流图像识别模型进行训练,得到训练好的泥石流图像识别模型;
每隔预设时间段对所述监测区域进行一次激光扫描,得到当前时刻的点云数据,并基于GNSS导航系统和惯性导航系统,对所述当前时刻的点云数据进行定位,得到当前时刻的标准激光点云数据;
基于所述当前时刻的标准激光点云数据得到当前时刻的初始DEM;
对所述当前时刻的初始DEM进行滤波得到当前时刻的最终DEM;
利用所述训练好的泥石流图像识别模型和所述当前时刻的最终DEM对所述监测区域进行地质灾害的监控预警。
2.如权利要求1所述的三维场景灾变监控预警方法,其特征在于,所述基于改进的种子区域生长算法,对所述原始图像进行分割处理,得到分割后的原始图像,包括:
将所述原始图像进行区域划分,并从中随机选择一个区域作为种子区域;
对所述种子区域周围的每一邻接区域,判断所述邻接区域的平均颜色分量值和所述种子区域的平均颜色分量值的比较结果是否小于预设阈值,若是,将所述邻接区域与所述种子区域合并,得到合并后的区域;
选取一个平均颜色分量值和所述种子区域的平均颜色分量值相同的区域,作为所述种子区域,重复步骤对所述种子区域周围的每一邻接区域,判断所述邻接区域的平均颜色分量值和所述种子区域的平均颜色分量值的比较结果是否小于预设阈值,若是,将所述邻接区域与所述种子区域合并,得到合并后的区域,直到遍历完所述原始图像上所有的区域,得到分割后的原始图像。
3.如权利要求1所述的三维场景灾变监控预警方法,其特征在于,所述泥石流图像识别模型的构建过程至少包括:卷积操作、池化操作、ReLU变化、Dropout、归一化;
其中,根据下式进行卷积操作:
Figure FDA0003663480970000021
式中,
Figure FDA0003663480970000022
表示第l层卷积层的第j个特征图,f()表示激活函数,Mj表示上一层的输出集合,
Figure FDA0003663480970000023
表示卷积核,
Figure FDA0003663480970000024
表示偏置参数,i表示次数,*表示卷积操作;
根据下式进行池化操作:
Figure FDA0003663480970000025
式中,down()表示采样函数,
Figure FDA0003663480970000026
表示每个输出特征图的多倍偏置参数,
Figure FDA0003663480970000027
表示偏置参数;
根据下式进行ReLU变化:
yijk=max{0,xijk}
式中,xijk表示输入值,yijk表示输出值;
在Dropout中,将50%的隐藏层的输出设置为0;
根据下式进行归一化:
Figure FDA0003663480970000028
Figure FDA0003663480970000029
式中,p表示相邻信道输入特征群,α表示系列参数,k表示第k个核,xijk表示第k个核在位置(i,j)运用ReLU变化的非线性化输出,G(k)表示通道数,beta表示预设参数,G(k′)表示通道数。
4.如权利要求3所述的三维场景灾变监控预警方法,其特征在于,所述泥石流图像识别模型的训练过程至少包括反向传播。
其中,所述反向传播包括卷积层梯度计算和池化层梯度计算:
根据下式进行卷积层梯度计算:
Figure FDA0003663480970000031
式中,
Figure FDA0003663480970000032
表示节点j在l层的第一个导函数误差,β表示输入系数,°表示每个部分相乘,
Figure FDA0003663480970000033
表示对于节点j的偏导值的激活函数,up()表示采样函数,u表示输入节点,δ表示导函数误差;
Figure FDA0003663480970000034
式中,
Figure FDA0003663480970000035
表示
Figure FDA0003663480970000036
Figure FDA0003663480970000037
的转置,conv2表示一个卷积函数,rot180表示卷积核被转置;
根据下式进行池化层梯度计算:
Figure FDA0003663480970000038
式中,
Figure FDA0003663480970000039
5.如权利要求1所述的三维场景灾变监控预警方法,其特征在于,所述基于所述当前时刻的标准激光点云数据得到当前时刻的初始DEM,包括:
将所述当前时刻的标准激光点云数据进行自适应分块处理,得到每一分块的点云数据;
基于主成分分析法,对在x方向和y方向的所述每一分块的点云数据进行局部平面拟合处理,并利用反距离加权法和样条加权法对经过局部平面拟合处理过的点云数据进行插值处理,得到初始拟合曲面,对在z方向的所述每一分块的点云数据和所述初始拟合曲面进行滤波处理,得到最终拟合曲面;
对所述最终拟合曲面进行滤波处理,得到当前时刻的初始DEM。
6.如权利要求1所述的三维场景灾变监控预警方法,其特征在于,所述对所述当前时刻的初始DEM进行滤波得到当前时刻的最终DEM,包括:
采用小波函数来对所述当前时刻的初始DEM进行滤波,得到当前时刻的最终DEM。
CN202210580210.7A 2022-05-26 2022-05-26 一种三维场景灾变监控预警方法 Active CN114882366B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210580210.7A CN114882366B (zh) 2022-05-26 2022-05-26 一种三维场景灾变监控预警方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210580210.7A CN114882366B (zh) 2022-05-26 2022-05-26 一种三维场景灾变监控预警方法

Publications (2)

Publication Number Publication Date
CN114882366A true CN114882366A (zh) 2022-08-09
CN114882366B CN114882366B (zh) 2023-03-21

Family

ID=82678617

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210580210.7A Active CN114882366B (zh) 2022-05-26 2022-05-26 一种三维场景灾变监控预警方法

Country Status (1)

Country Link
CN (1) CN114882366B (zh)

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103529455A (zh) * 2013-10-21 2014-01-22 中铁第四勘察设计院集团有限公司 一种基于机载激光雷达三维的危岩落石调查方法
CN105242279A (zh) * 2015-07-30 2016-01-13 王植 一种基于激光雷达技术的滑坡体变化检测方法
CN107655457A (zh) * 2016-12-23 2018-02-02 航天星图科技(北京)有限公司 一种基于遥感卫星图像的泥石流地质灾害识别方法
CN111340012A (zh) * 2020-05-19 2020-06-26 北京数字绿土科技有限公司 一种地质灾害解译方法、装置、终端设备
US20210048294A1 (en) * 2019-08-15 2021-02-18 China Institute Of Water Resources And Hydropower Research System and method for monitoring deformation of dam slope
CN113723403A (zh) * 2021-08-26 2021-11-30 广东电网有限责任公司广州供电局 一种滑坡监测方法、装置、电子设备及存储介质
CN114049565A (zh) * 2021-11-08 2022-02-15 中国公路工程咨询集团有限公司 基于遥感影像和dem数据的地质灾害判识方法和装置
CN114120129A (zh) * 2021-11-30 2022-03-01 哈尔滨工业大学 基于无人机图像和深度学习的滑坡滑移面的三维识别方法
WO2022061945A1 (zh) * 2020-09-27 2022-03-31 广东电网有限责任公司 一种电力线路安全距离检测方法
CN114267155A (zh) * 2021-11-05 2022-04-01 国能大渡河革什扎水电开发有限公司 一种基于视频识别技术的地质灾害监测预警系统

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103529455A (zh) * 2013-10-21 2014-01-22 中铁第四勘察设计院集团有限公司 一种基于机载激光雷达三维的危岩落石调查方法
CN105242279A (zh) * 2015-07-30 2016-01-13 王植 一种基于激光雷达技术的滑坡体变化检测方法
CN107655457A (zh) * 2016-12-23 2018-02-02 航天星图科技(北京)有限公司 一种基于遥感卫星图像的泥石流地质灾害识别方法
US20210048294A1 (en) * 2019-08-15 2021-02-18 China Institute Of Water Resources And Hydropower Research System and method for monitoring deformation of dam slope
CN111340012A (zh) * 2020-05-19 2020-06-26 北京数字绿土科技有限公司 一种地质灾害解译方法、装置、终端设备
WO2022061945A1 (zh) * 2020-09-27 2022-03-31 广东电网有限责任公司 一种电力线路安全距离检测方法
CN113723403A (zh) * 2021-08-26 2021-11-30 广东电网有限责任公司广州供电局 一种滑坡监测方法、装置、电子设备及存储介质
CN114267155A (zh) * 2021-11-05 2022-04-01 国能大渡河革什扎水电开发有限公司 一种基于视频识别技术的地质灾害监测预警系统
CN114049565A (zh) * 2021-11-08 2022-02-15 中国公路工程咨询集团有限公司 基于遥感影像和dem数据的地质灾害判识方法和装置
CN114120129A (zh) * 2021-11-30 2022-03-01 哈尔滨工业大学 基于无人机图像和深度学习的滑坡滑移面的三维识别方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
孙显辰: "基于深度迁移学习的遥感影像中泥石流沟谷分类", 《中国优秀硕士学位论文全文数据库 基础科学辑》 *
孙鹏 等: "DEM高程插值方法分析", 《城市地理》 *
徐天芝 等: "基于形态学梯度的分水岭彩色图像分割", 《计算机工程与应用》 *

Also Published As

Publication number Publication date
CN114882366B (zh) 2023-03-21

Similar Documents

Publication Publication Date Title
CN108230302B (zh) 一种核电厂冷源海域入侵海洋生物探测及处置方法
WO2020191980A1 (zh) 一种无线传感网络数据漂移盲校准方法
Gomez-Rodriguez et al. Smoke monitoring and measurement using image processing: application to forest fires
CN104680508A (zh) 卷积神经网络和基于卷积神经网络的目标物体检测方法
CN110533665B (zh) 一种抑制扇贝效应和子带拼接效应的sar图像处理方法
CN114332008B (zh) 一种基于多层次特征重构的无监督缺陷检测与定位方法
CN114926797A (zh) 基于边缘约束与特征适应的Transformer双分支道路提取方法及装置
KR101869266B1 (ko) 극한 심층학습 기반 차선 검출 시스템 및 그 방법
CN111488827A (zh) 一种基于多尺度特征信息的人群计数方法及系统
CN114494821A (zh) 基于特征多尺度感知与自适应聚合的遥感影像云检测方法
CN111126185B (zh) 一种针对道路卡口场景的深度学习车辆目标识别方法
CN118247280B (zh) 一种基于多源遥感数据的河流断面形态迁移模拟方法
CN116485885A (zh) 基于深度学习去除视觉slam前端动态特征点的方法
CN112380985A (zh) 变电站入侵异物实时检测方法
Megalingam et al. Computer vision-based measuring method to estimate the diameter of the coconut tree trunk
CN114882366B (zh) 一种三维场景灾变监控预警方法
CN117116096A (zh) 基于多通道交通图像与深度cnn的机场延误预测方法及系统
CN117495935A (zh) 一种基于交叉特征提取窗口与动态特征融合窗口的深度补全方法
CN109948571A (zh) 一种光学遥感图像船舶检测方法
CN113888488A (zh) 一种基于深度残差收缩网络的钢轨缺陷检测方法及系统
CN112036349A (zh) 一种基于遥感影像和人工智能的融雪速率分析方法
Ullah et al. Real-time video dehazing for industrial image processing
AU2021102735A4 (en) Cotton aphid stress monitoring method and system based on high-throughput plant phenotyping platform and deep learning
Mamatov et al. Methods for increasing the contrast of drone agricultural images
Koley et al. A wavelet-based low frequency prior for single-image dehazing

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