CN110110946B - 基于各向异性Delaunay剖分的水质预测预警系统及其实现方法 - Google Patents
基于各向异性Delaunay剖分的水质预测预警系统及其实现方法 Download PDFInfo
- Publication number
- CN110110946B CN110110946B CN201910448164.3A CN201910448164A CN110110946B CN 110110946 B CN110110946 B CN 110110946B CN 201910448164 A CN201910448164 A CN 201910448164A CN 110110946 B CN110110946 B CN 110110946B
- Authority
- CN
- China
- Prior art keywords
- water quality
- graph
- river
- establishing
- tubular
- 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.)
- Active
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 36
- 238000000034 method Methods 0.000 title claims description 26
- 239000011159 matrix material Substances 0.000 claims abstract description 18
- 238000012544 monitoring process Methods 0.000 claims abstract description 18
- 238000012549 training Methods 0.000 claims abstract description 16
- 238000005259 measurement Methods 0.000 claims abstract description 14
- 238000001914 filtration Methods 0.000 claims abstract description 8
- 238000013527 convolutional neural network Methods 0.000 claims abstract description 6
- 238000010586 diagram Methods 0.000 claims description 12
- 238000013528 artificial neural network Methods 0.000 claims description 4
- 238000010801 machine learning Methods 0.000 claims description 3
- 230000011218 segmentation Effects 0.000 description 8
- 238000000605 extraction Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 230000009466 transformation Effects 0.000 description 3
- 230000002708 enhancing effect Effects 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 239000003550 marker Substances 0.000 description 2
- 238000000638 solvent extraction Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 102000002274 Matrix Metalloproteinases Human genes 0.000 description 1
- 108010000684 Matrix Metalloproteinases Proteins 0.000 description 1
- 230000004913 activation Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000011176 pooling Methods 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
-
- 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
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A20/00—Water conservation; Efficient water supply; Efficient water use
- Y02A20/152—Water filtration
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Tourism & Hospitality (AREA)
- Strategic Management (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- General Health & Medical Sciences (AREA)
- Software Systems (AREA)
- Biomedical Technology (AREA)
- General Engineering & Computer Science (AREA)
- Molecular Biology (AREA)
- Mathematical Physics (AREA)
- Evolutionary Computation (AREA)
- Data Mining & Analysis (AREA)
- Computational Linguistics (AREA)
- Development Economics (AREA)
- Biophysics (AREA)
- General Business, Economics & Management (AREA)
- Computing Systems (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (AREA)
- Marketing (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Entrepreneurship & Innovation (AREA)
- Game Theory and Decision Science (AREA)
- Educational Administration (AREA)
- Primary Health Care (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明属于水质预测预警系统技术领域,尤其为基于各向异性Delaunay剖分的水质预测预警系统,包括通过高程图提取河网信息、通过若干个监测站点之间图结构的建立和图卷积网络的建立、训练与预测,采用海森矩阵滤波用于高程图增强河流结构,得到管状测度(Vesselness)与管状向量场,构建变分模型,将高程图原始灰度值、管状测度及其各向异性结合,提出新的能量模型;根据图结构以及管状度量距离,确定标志节点与其邻域,采用水质历史数据以及河网信息训练图卷积神经网络,对水质数据进行预测。
Description
技术领域
本发明属于水质预测预警系统技术领域,具体涉及基于各向异性Delaunay剖分的水质预测预警系统及其实现方法。
背景技术
目前水质检测过程中河道水文信息不明,污染点源难以获取,企业存在偷排漏排的现象,在对水质预警进行检测过程中水质参数不易获取,针对目前的水质预测预警过程中所暴露的问题,有必要开发出水质预警系统。
发明内容
为解决现有技术中存在的上述问题,本发明提供了基于各向异性Delaunay剖分的水质预测预警系统及其实现方法,具有获取河网信息,便于对河网水质及逆行预测预警的特点。
为实现上述目的,本发明提供如下技术方案:基于各向异性Delaunay剖分的水质预测预警系统,包括高程图提取河网信息、若干个监测站点之间图结构的建立和图卷积网络的建立、训练与预测,其中高程图提取河网信息,采用海森矩阵(HessianMatrix)滤波用于高程图增强河流结构,得到管状测度(Vesselness)与管状向量场,构建变分模型,将高程图原始灰度值、管状测度及其各向异性结合,提出新的能量模型,并应用图分割(Graph-cuts)算法,以分割河道区域;
若干个所述监测站点之间图结构的建立,用于将管状测度度量空间定义下的各向异性Delaunay三角剖分,建立监测站点间的图结构;
图卷积网络的建立、训练与预测,根据图结构以及管状度量距离,确定标志节点与其邻域,采用水质历史数据以及河网信息训练图卷积神经网络,用于预测水质数据。
作为本发明的基于各向异性Delaunay剖分的水质预测预警系统及其实现方法优选技术方案,若干个所述监测站之间图结构的建立的拓扑关系需由空间位置与河网信息进一步抽象提取得出。
作为本发明的基于各向异性Delaunay剖分的水质预测预警系统及其实现方法优选技术方案,所述建立监测站点间的图结构当中边的权值由河道宽度,长度等多维几何特征共同决定。
作为本发明的基于各向异性Delaunay剖分的水质预测预警系统的实现方法,包括以下的步骤:
第一步:对河网信息进行建立,并提取河网信息;
第二步:建立监测站点图结构;
第三步:通过卷积神经网络进行机器学习,对图卷积网络进行建立、训练与预测;
第四步:选取部分节点作为卷积滤波的中心点,对标志节点进行选取;
第五步:通过经典图卷积网络方法对节点邻域与卷积网络结构进行确定;
与现有技术相比,本发明的有益效果是:通过对河网信息进行提取,建立监测站点结构之后再通过算法对河网信息进行学习、训练以及预测,根据图结构以及管状度量距离,确定标志节点与其邻域,采用水质历史数据以及河网信息训练图卷积神经网络,对水质数据进行预测。
附图说明
附图用来提供对本发明的进一步理解,并且构成说明书的一部分,与本发明的实施例一起用于解释本发明,并不构成对本发明的限制。在附图中:
图1为本发明的高程图示意图;
图2为本发明中的管状测度结构示意图;
图3为本发明中的河网分割结构示意图;
图4为本发明中的像素灰度结构示意图;
图5为本发明中的像素识别结构示意图;
图6为本发明中的权重结构示意图;
图7为本发明中的站点位置结构示意图;
图8为本发明中的各向异性Delaunay剖分结构示意图;
图10为本发明中的经典Dijkstra算法结构示意图;
图11为本发明中的聚类结构示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
实施例
请参阅图1-11,本发明提供以下技术方案:基于各向异性Delaunay剖分的水质预测预警系统,包括高程图提取河网信息、若干个监测站点之间图结构的建立和图卷积网络的建立、训练与预测,其中高程图提取河网信息,采用海森矩阵(HessianMatrix)滤波用于高程图增强河流结构,得到管状测度(Vesselness)与管状向量场,构建变分模型,将高程图原始灰度值、管状测度及其各向异性结合,提出新的能量模型,并应用图分割(Graph-cuts)算法,以分割河道区域;
若干个所述监测站点之间图结构的建立,用于将管状测度度量空间定义下的各向异性Delaunay三角剖分,建立监测站点间的图结构;
图卷积网络的建立、训练与预测,根据图结构以及管状度量距离,确定标志节点与其邻域,采用水质历史数据以及河网信息训练图卷积神经网络,用于预测水质数据。
作为本发明的基于各向异性Delaunay剖分的水质预测预警系统及其实现方法优选技术方案,若干个所述监测站之间图结构的建立的拓扑关系需由空间位置与河网信息进一步抽象提取得出。
作为本发明的基于各向异性Delaunay剖分的水质预测预警系统及其实现方法优选技术方案,所述建立监测站点间的图结构当中边的权值由河道宽度,长度等多维几何特征共同决定。
在图分割的经典应用当中,设置N-link的边权重,见式(8):
以相邻像素的灰度值差值衡量分割这两像素,引进海森矩阵滤波得到的管状测度以及其各向异性至能量泛函并影响边的权重设置,见式(9):
引入,衡量管状测度之间的差值,测算两个
像素的河道走向的内积,用以衡量两个像素的河道走向是否一致。该各向异性
的引进,有效的过滤分离位置靠近但走向不一致的两根河道。而这类情况在复杂的河网结
构中经常发生。
引入各向异性至图分割算法中,分割中减少该类现象。再执行最大流/最小分割(max-flow/min-cut)算法后,获取河道分割结果,二维分割结果如图3所示;
第二步:建立监测站点图结构;
式(11)中最后一项中为欧氏距离,则距离可通过变换原空间,并在新空间中
求解欧氏距离得到,该变换被称为拉伸变换。求解度量下的各向异性Delaunay三角
剖分。对于点集,其在度量下的Delaunay三角剖分的求解如下:
二维空间中用于生成Delaunay三角剖分的经典方法有逐点插入法与对角线交换
法。本项目拟采用对角线交换法,该算法对给定点集的任一三角剖分,遍历其所有凸四边
形,检查两种可能的对角线,选择较短之一,获得Delaunay三角剖分,见图7、图8以及图9所
示,该算法的单步操作。在度量下对角线的长度定义为距离。
图结构的拓扑获得后,图中边的权值可以直接由距离定义,代表在该
地形上两站点之间的远近关系,也可以通过计算,;点间在管状测度上的最短距离定义,
该最短距离可通过经典Dijkstra算法得到,相较更为精确。建立站点间图结构的
整个过程如图10所示;
第三步:通过卷积神经网络进行机器学习,对图卷积网络进行建立、训练与预测;
第四步:选取部分节点作为卷积滤波的中心点,对标志节点进行选取;
参阅图11,对中所有节点进行有序标记,排序并选取最高的个节
点作为标志节点;在本实施例中,除了代表图的节点,本身还具有几何坐标,大大降低了
这一步骤的难度。标志节点可以通过对全体节点进行组聚类,得到对全体节点的划分,
再进一步对每一组提取一个标志节点。图11即为经典K-means算法对所有节点进行聚类。对
第组节点,选取标志节点的方法如下:见式(14);
第五步:通过经典图卷积网络方法对节点邻域与卷积网络结构进行确定;
经典图卷积网络方法对于给定的感受野的长度,往往会通过迭代寻找邻域(1-
邻域,2-邻域,…)的方式,确定某标志节点的邻域。这样获取的邻域节点的数目往往会大于,还需进一步对图结构进行归一化,使得邻域数目固定。而本项目中节点自带的坐标和距
离信息,可以很好的解决这一问题。对于标志节点,选取度量下离最近的个节点,即
组成了的邻域;见式(15);
利用Delaunay剖分的几何特性(每个节点总是与最近的其他节点相连),更高效的
计算的邻域。确定了个标志节点的长度邻域后,通过假设节点的输入属性为
维向量,边的输入属性为维向量,则可生成维度与维度的张量。节
点的输入为该监测站点水质历史数据,边的输入为度量下的距离。通过设计一维的卷积
层,步长分别为与,作用于这两个张量。其余的池化层与全连接层则与经典的卷积网络
结构类似,最后输出层如果不使用激活函数,即可作为回归学习的工具,用于预测标志站点
的未来水质参数。
最后应说明的是:以上所述仅为本发明的优选实施例而已,并不用于限制本发明,尽管参照前述实施例对本发明进行了详细的说明,对于本领域的技术人员来说,其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (4)
1.基于各向异性Delaunay剖分的水质预测预警系统的实现方法,其特征在于,包括高程图提取河网信息、若干个监测站点之间图结构的建立和图卷积网络的建立、训练与预测,其中;
高程图提取河网信息,采用海森矩阵(HessianMatrix)滤波用于高程图增强河流结构,得到管状测度(Vesselness)与管状向量场,构建变分模型,将高程图原始灰度值、管状测度及其各向异性结合,提出新的能量模型,并应用图分割(Graph-cuts)算法,以分割河道区域;
若干个所述监测站点之间图结构的建立,用于将管状测度度量空间定义下的各向异性Delaunay三角剖分,建立监测站点间的图结构;
图卷积网络的建立、训练与预测,根据图结构以及管状度量距离,确定标志节点与其邻域,采用水质历史数据以及河网信息训练图卷积神经网络,用于预测水质数据。
2.根据权利要求1所述的基于各向异性Delaunay剖分的水质预测预警系统的实现方法,其特征在于:所述建立监测站点间的图结构当中边的权值由河道宽度,长度多维几何特征共同决定。
3.根据权利要求1-2任一项所述的基于各向异性Delaunay剖分的水质预测预警系统的实现方法,包括以下的步骤:
第一步:对河网信息进行建立,并提取河网信息;
第二步:建立监测站点图结构;
第三步:通过卷积神经网络进行机器学习,对图卷积网络进行建立、训练与预测;
第四步:选取部分节点作为卷积滤波的中心点,对标志节点进行选取;
第五步:通过经典图卷积网络方法对节点邻域与卷积网络结构进行确定。
4.根据权利要求1所述的基于各向异性Delaunay剖分的水质预测预警系统的实现方法的应用系统,其特征在于:若干个所述监测站之间图结构的建立的拓扑关系需由空间位置与河网信息进一步抽象提取得出。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910448164.3A CN110110946B (zh) | 2019-05-27 | 2019-05-27 | 基于各向异性Delaunay剖分的水质预测预警系统及其实现方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910448164.3A CN110110946B (zh) | 2019-05-27 | 2019-05-27 | 基于各向异性Delaunay剖分的水质预测预警系统及其实现方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110110946A CN110110946A (zh) | 2019-08-09 |
CN110110946B true CN110110946B (zh) | 2021-09-07 |
Family
ID=67492489
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910448164.3A Active CN110110946B (zh) | 2019-05-27 | 2019-05-27 | 基于各向异性Delaunay剖分的水质预测预警系统及其实现方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110110946B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111458471B (zh) * | 2019-12-19 | 2023-04-07 | 中国科学院合肥物质科学研究院 | 一种基于图神经网络的水域检测预警方法 |
CN113034425B (zh) * | 2019-12-25 | 2024-05-28 | 阿里巴巴集团控股有限公司 | 数据处理方法、设备及存储介质 |
CN112529296B (zh) * | 2020-12-11 | 2023-01-31 | 联通(浙江)产业互联网有限公司 | 水质预测方法、装置及服务器 |
CN112651665B (zh) * | 2021-01-14 | 2022-12-20 | 浙江鸿程计算机系统有限公司 | 一种基于图神经网络的地表水水质指标预测方法和装置 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102156779B (zh) * | 2011-04-13 | 2013-03-20 | 北京石油化工学院 | 地下水流仿真与预测分析方法 |
CN104376380B (zh) * | 2014-11-17 | 2017-07-21 | 北京工业大学 | 一种基于递归自组织神经网络的氨氮浓度预测方法 |
-
2019
- 2019-05-27 CN CN201910448164.3A patent/CN110110946B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN110110946A (zh) | 2019-08-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110110946B (zh) | 基于各向异性Delaunay剖分的水质预测预警系统及其实现方法 | |
CN103310453B (zh) | 一种基于子图像角点特征的快速图像配准方法 | |
CN104616278B (zh) | 三维点云兴趣点检测方法和系统 | |
CN101980250B (zh) | 基于降维局部特征描述子和隐条件随机场的目标识别方法 | |
CN109165540B (zh) | 一种基于先验候选框选择策略的行人搜索方法和装置 | |
CN103456022B (zh) | 一种高分辨率遥感图像特征匹配方法 | |
CN103745459B (zh) | 一种非结构化点云特征点检测方法及其提取方法 | |
CN106529591A (zh) | 一种基于改进的mser图像匹配算法 | |
CN100545865C (zh) | 一种对图像初始分割边界进行优化的自动分割方法 | |
CN108629809B (zh) | 一种精确高效的立体匹配方法 | |
CN103679734A (zh) | 基于svm和pde的有眼台风二维表面风场反演方法 | |
Abdelsamea et al. | A SOM-based Chan–Vese model for unsupervised image segmentation | |
CN104835168B (zh) | 基于全局凸优化变分模型的快速多相图像分割方法 | |
CN104732546A (zh) | 区域相似性和局部空间约束的非刚性sar图像配准方法 | |
CN109741358B (zh) | 基于自适应超图学习的超像素分割方法 | |
CN105976364A (zh) | 基于简化带权无向图的统计平均模型构建方法 | |
CN101630407A (zh) | 基于两视几何和图分割的伪造区域定位方法 | |
CN108846845B (zh) | 基于缩略图与分层模糊聚类的sar图像分割方法 | |
CN105912977A (zh) | 基于点聚类的车道线检测方法 | |
CN110930413A (zh) | 一种基于弱监督多核分类优化合并的图像分割方法 | |
CN111553296B (zh) | 一种基于fpga实现的二值神经网络立体视觉匹配方法 | |
CN115035164A (zh) | 一种运动目标识别方法及装置 | |
CN109191503A (zh) | 基于条件随机场的遥感影像变化检测方法及系统 | |
CN107564013A (zh) | 融合局部信息的场景分割修正方法与系统 | |
CN105809657A (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 |