CN115879695A - 一种基于K-means算法的灌溉喷头布局方法 - Google Patents
一种基于K-means算法的灌溉喷头布局方法 Download PDFInfo
- Publication number
- CN115879695A CN115879695A CN202211418096.4A CN202211418096A CN115879695A CN 115879695 A CN115879695 A CN 115879695A CN 202211418096 A CN202211418096 A CN 202211418096A CN 115879695 A CN115879695 A CN 115879695A
- Authority
- CN
- China
- Prior art keywords
- uniformity
- irrigation
- center
- land
- grid points
- 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
- 230000002262 irrigation Effects 0.000 title claims abstract description 55
- 238000003973 irrigation Methods 0.000 title claims abstract description 55
- 239000007921 spray Substances 0.000 title claims abstract description 51
- 238000004422 calculation algorithm Methods 0.000 title claims abstract description 24
- 238000000034 method Methods 0.000 title claims abstract description 11
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 23
- 238000005507 spraying Methods 0.000 claims abstract description 4
- 239000002699 waste material Substances 0.000 abstract description 5
- 238000009434 installation Methods 0.000 abstract description 3
- 238000001556 precipitation Methods 0.000 description 12
- 238000010586 diagram Methods 0.000 description 4
- 230000001788 irregular Effects 0.000 description 4
- 239000003621 irrigation water Substances 0.000 description 4
- 238000012544 monitoring process Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000010801 machine learning Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 239000002689 soil Substances 0.000 description 2
- 206010063385 Intellectualisation Diseases 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
Images
Classifications
-
- 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
- Y02A40/00—Adaptation technologies in agriculture, forestry, livestock or agroalimentary production
- Y02A40/10—Adaptation technologies in agriculture, forestry, livestock or agroalimentary production in agriculture
- Y02A40/22—Improving land use; Improving water use or availability; Controlling erosion
Landscapes
- Spray Control Apparatus (AREA)
- Nozzles (AREA)
Abstract
本发明公开了一种基于K‑means算法的灌溉喷头布局方法,包括步骤:获取灌溉土地的顶点坐标将灌溉土地构成多边形;获取喷头喷洒水量的3D分布数据;获取将多边形包围的最小矩形,并产生等间距的网格点;使用多边形内外点判别算法,选出多边形内的所有网格点;将多边形内的所有网格点作为灌溉土地的数据集;随机选取数据集中的多个样本作为初始中心,并产生多个空簇的集合;使用K‑means算法,计算多个喷头的安装位置;计算多个喷头下的克里斯琴森均匀度、分布均匀度和统计均匀度,再计算综合均匀度;选取综合均匀度最优的喷头个数及其安装位置,作为灌溉土地的喷头布局。本发明得到的喷头合理布局,充分发挥了喷灌设备的性能,避免了水资源浪费。
Description
技术领域
本发明涉及机器学习和农业灌溉技术领域,特别是一种基于K-means算法的灌溉喷头布局方法。
背景技术
农业灌溉用水占据了农业用水的绝大部分,农业灌溉用水的高效利用尤为重要。喷灌技术在当前农业灌溉中被大量运用,是农业节水灌溉的一种重要技术措施,具有节水、增产、增收的优点。喷灌方式有全移动式和半固定式、固定式,其中固定式具有操作方便,效率高,占地少,便于实现灌溉智能化等优点。喷灌水量分布均匀性是衡量喷灌质量和喷头水力性能的重要指标,是喷灌系统规划设计中的重要参数。
然而,现有的灌溉喷头布局主要依靠人工凭经验完成,当土地面积很大或者土地形状不规则的时候很难做到较优的布局,并且不能对喷头布局的喷灌效果进行有效评估,这样会造成人力资源的浪费、水资源浪费及喷灌设备不能充分发挥其性能的浪费。
发明内容
本发明的目的是提供一种基于K-means算法的灌溉喷头布局方法。
实现本发明目的的技术方案如下:
一种基于K-means算法的灌溉喷头布局方法,包括:
步骤1:获取灌溉土地的顶点坐标p1、p2、p3、…、pn,将灌溉土地构成多边形P=p1p2p3…pn;获取喷头喷洒水量的3D分布数据;
步骤2:获取将P包围的最小矩形R,并在R中产生等间距的网格点;遍历R内的所有网格点,使用多边形内外点判别算法,选出P内的所有网格点;将P内的所有网格点作为灌溉土地的数据集X={x1,x2,…,xi,…,xn};其中,xi为网格点,对应一棵植株;n为P内网格点的数量;
步骤3:随机选取X中的k个样本作为初始中心Z={z1,z2,…,zj,…,zk},并产生k个空簇c的集合C={c1,c2,…,cj,…,ck};其中,zj表示第j个簇cj的中心,即第j个喷头的安装位置;cj表示能被第j个喷头灌溉的植株所构成的簇;
步骤4:使用K-means算法,计算k个喷头的安装位置;
4.1根据当前中心Z,遍历X中所有的xi,将xi分配到与其欧式距离最近的簇中心zj所对应的簇cj中;
4.2对簇cj内的所有网格点分别计算质心得到簇cj的新中心z′j;遍历集合C中所有的簇,得到新中心Z’;
4.3如果新中心Z’相对于当前中心Z没有发生变化,则令当前中心Z为聚类中心Z={z1,z2,…,zj,…,zk},得到k个喷头的安装位置;否则,令Z=Z’,并将集合C内所有的簇重置为空簇,返回4.1;
步骤5:根据聚类中心Z={z1,z2,…,zj,…,zk}和喷头喷洒水量的3D分布数据,分别计算得到k个喷头下的克里斯琴森均匀度CUk、分布均匀度DUk和统计均匀度SUk,再计算得到综合均匀度Uk=(CUk+DUk+SUk)/3;
步骤6:令喷头个数k=1,2,…,m,分别按照步骤3至步骤5的方法,计算得到综合均匀度U1,U2,…,Um;
步骤7:选取U1,U2,…,Um中最大值所对应的喷头个数k,与其聚类中心Z={z1,z2,…,zj,…,zk},作为灌溉土地的喷头布局。
本发明的有益效果在于,通过对任意形状和任意规模大小的灌溉土地进行智能划分,基于K-means算法和综合均匀度指标对灌溉土地的喷头合理布局,充分发挥了喷灌设备的性能,避免了水资源浪费。
附图说明
图1是本发明的框架结构图。
图2a至图2d是灌溉土地在不同喷头个数下的喷头布局图。其中,图2a中喷头个数为3,图2b中喷头个数为4,图2c中喷头个数为5,图2d中喷头个数为8。
图3a是实施例中灌溉喷头水量的截面分布图。
图3b是实施例中灌溉喷头水量的3D分布图。
图4a是在规则土地情况下水量3D分布图。
图4b是在规则土地情况下灌溉性能指数图。
图5a是在不规则土地情况下水量分布热力图。
图5b是在不规则土地情况下灌溉性能指数图。
具体实施方式
一种基于K-means算法的灌溉喷头布局方法,包括:
步骤1:输入灌溉土地的坐标信息(即多边形P=p1p2p3…pkpk+1…pn的顶点坐标),灌溉喷头的喷洒水量的3D分布数据。
步骤2:获取能够将多边形P包围的最小矩形R,并在矩形R中产生等间距的网格点。遍历该矩形R内的所有网格点,使用多边形内外点判别算法判断这些点是否是多边形P内的点,将多边形P内的所有网格点作为所要灌溉区域的模拟土地数据集X={x1,x2,…,xi,…,xn}。其中数据集X中每一个网格点xi就代表一棵植株,n代表模拟土地数据集中网格点的个数(即植株的数量)。
步骤3:随机选取数据集X={x1,x2,…,xi,…,xn}中的k个样本作为初始聚类中心Z={z1,z2,…,zj,…,zk}(即喷头安装的位置),并产生k个空簇集合C={c1,c2,…,cj,…,ck},其中zj表示第j个簇cj的中心点(即第j个喷头的安装位置),cj表示被第j个喷头所灌溉的植株构成的簇。
步骤4:遍历数据集X中所有的点xi,将数据对象xi分配到与之距离最近的簇中心zj所对应的簇cj中,其距离计算选用欧式距离。
步骤5:遍历簇C中所有的簇cj,并对每个簇内的所有元素分别计算质心得到新的中心点z′j,并将该点作为该簇新的中心点,即zj=z′j。
步骤6:重复(4)(5)步,当聚类中心Z不再发生变化时K-means算法结束。此时可以得到当喷头个数为k时喷头的最佳安装位置为Z={z1,z2,…,zj,…,zk}。
步骤7:针对步骤6所得到的喷头最佳安装位置可以得到喷灌区域X的喷头布局,根据该布局得到评价灌溉性能的克里斯琴森均匀度CUk、分布均匀度DUk和统计均匀度SUk,并计算得到喷灌性能的综合均匀度Uk=(CUk+DUk+SUk)/3。
步骤8:分别选取不同的喷头个数k=1,2,…,m,重复进行步骤(3)至步骤(7),可以得到对应的综合均匀度U={U1,U2,…,Uk,…,Um}。
步骤9:输出到综合均匀度U中值最大所对应的喷头个数k和对应的聚类中心Z={z1,z2,…,zj,…,zk}。
其中步骤3到步骤6采用的是K-means算法。
K-means算法是机器学习中经典的聚类算法,通过不断迭代聚类,将一个数据集划分为若干个簇(组或类),其中划分得到的簇具有簇内元素相似度高,簇间相似度低的特点。
下面结合附图和具体实施例对本发明进一步说明。
一、基于K-means算法的灌溉喷头布局:
本发明的整体框架结构如图1所示。首先是获取土地的边界信息,比如顶点坐标,通过内外点判别算法利用计算机仿真技术产生均匀分布的网格点得到模拟土地S,将K-means算法作用于模拟土地S上进行土地划分,根据设置不同簇的个数K得到聚类的中心点(簇)及喷头对应的喷灌区域。然后根据划分的结果计算出克里斯琴均匀度、分布均匀度和统计均匀度,最后选出最佳均匀度所对应的喷头个数、各喷头的位置以及最佳喷灌均匀度的值。
二、模拟土地的灌溉喷头布局:
选取一个规则的正方形土地作为对象,利用多边形内外点判别算法得到这个正方形内均匀分布的点,即模拟土地。利用K-means算法对原模拟土地进行聚类划分,当k分别取3,4,5,8时得到对应的土地划分2a,2b,2c,2d。图中的黑色点表示喷头安装的位置,以喷头为中心不同符号标记的区域表示该喷头所灌溉的区域。
三、喷头水量分布
本实施例采用型号为Nelson R33的喷头,在0.35MPa工作压力下获得实测的降水深度,通过三次样条插值得到相对平滑的降水深度截面图,详见图3a。其中,x轴表示离喷头的距离,y轴表示降水的深度。对应的喷头水量3D分布详见图3b。
四、灌溉性能指数
喷灌水量分布均匀性是衡量喷灌质量和喷头水力性能的重要指标,是喷灌系统规划设计中的重要参数。目前比较主流的均匀度主要包括克里斯琴森均匀度、分布均匀度和统计均匀度,本实施例采用这三个指标来检验喷灌布局的性能,以下是对以上三个均匀度的详细介绍。
(1)克里斯琴森均匀度
克里斯琴森均匀度是常用且经典的喷灌水量均匀度评价指标,其计算公式如下:
(2)分布均匀度
灌溉喷头喷出的水应对作物尽可能全覆盖,若某些区域不能得到一定降水量则会严重影响作物的生长,克里斯琴森均匀度未能很好地反应这种情况,但分布均匀度能够较好地体现这一点,其计算公式如下:
(3)统计均匀度
PEREIRA等人提出使用统计均匀度SU作为喷灌均匀度的评价指标,其公式为:
图4a是在规则土地情况下喷头个数为4的喷灌水量的3D分布,图4b是对应的喷灌性能指数,包括克里斯琴森均匀度、分布均匀度和统计均匀度。可以发现,随着喷头个数的增加喷灌的性能指数先是有明显的提升,当喷头个数达到4的时候喷灌效果最好;当喷头个数超过4个的时候,随着喷头个数的增加其喷灌性反而有降低的趋势。因此,在实际生产应用当中安装喷头的个数并不是越多越好。
图5a是在不规则土地情况下喷头个数为20的喷头水量的热力图,颜色越深则表示水量越大。图5b是对应的喷头灌溉性能指数,可以发现当喷头个数设置为20到25这个区间的时候喷灌效果比较理想。通过喷灌性能指数图可以得到喷头个数为不同值时对应的灌溉性能指数,利用喷灌性能指数可以在实际工作生产中提供最佳喷头分布的方案。
Claims (1)
1.一种基于K-means算法的灌溉喷头布局方法,其特征在于,包括:
步骤1:获取灌溉土地的顶点坐标p1、p2、p3、…、pn,将灌溉土地构成多边形P=p1p2p3…pn;获取喷头喷洒水量的3D分布数据;
步骤2:获取将P包围的最小矩形R,并在R中产生等间距的网格点;遍历R内的所有网格点,使用多边形内外点判别算法,选出P内的所有网格点;将P内的所有网格点作为灌溉土地的数据集X={x1,x2,…,xi,…,xn};其中,xi为网格点,对应一棵植株;n为P内网格点的数量;
步骤3:随机选取X中的k个样本作为初始中心Z={z1,z2,…,zj,…,zk},并产生k个空簇c的集合C={c1,c2,…,cj,…,ck};其中,zj表示第j个簇cj的中心,即第j个喷头的安装位置;cj表示能被第j个喷头灌溉的植株所构成的簇;
步骤4:使用K-means算法,计算k个喷头的安装位置;
4.1根据当前中心Z,遍历X中所有的xi,将xi分配到与其欧式距离最近的簇中心zj所对应的簇cj中;
4.2对簇cj内的所有网格点分别计算质心得到簇cj的新中心z′j;遍历集合C中所有的簇,得到新中心Z’;
4.3如果新中心Z’相对于当前中心Z没有发生变化,则令当前中心Z为聚类中心Z={z1,z2,…,zj,…,zk},得到k个喷头的安装位置;否则,令Z=Z’,并将集合C内所有的簇重置为空簇,返回4.1;
步骤5:根据聚类中心Z={z1,z2,…,zj,…,zk}和喷头喷洒水量的3D分布数据,分别计算得到k个喷头下的克里斯琴森均匀度CUk、分布均匀度DUk和统计均匀度SUk,再计算得到综合均匀度Uk=(CUk+DUk+SUk)/3;
步骤6:令喷头个数k=1,2,…,m,分别按照步骤3至步骤5的方法,计算得到综合均匀度U1,U2,…,Um;
步骤7:选取U1,U2,…,Um中最大值所对应的喷头个数k,与其聚类中心Z={z1,z2,…,zj,…,zk},作为灌溉土地的喷头布局。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211418096.4A CN115879695B (zh) | 2022-11-14 | 2022-11-14 | 一种基于K-means算法的灌溉喷头布局方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211418096.4A CN115879695B (zh) | 2022-11-14 | 2022-11-14 | 一种基于K-means算法的灌溉喷头布局方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115879695A true CN115879695A (zh) | 2023-03-31 |
CN115879695B CN115879695B (zh) | 2024-02-09 |
Family
ID=85759814
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211418096.4A Active CN115879695B (zh) | 2022-11-14 | 2022-11-14 | 一种基于K-means算法的灌溉喷头布局方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115879695B (zh) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060144957A1 (en) * | 2004-12-02 | 2006-07-06 | Jacobsen Stephen C | Intelligent sprinkler irrigation system |
CN1908615A (zh) * | 2006-08-11 | 2007-02-07 | 西北农林科技大学 | 一种基于连续系统的喷灌均匀度测量方法 |
CN103506237A (zh) * | 2013-10-11 | 2014-01-15 | 江苏大学 | 一种正方形喷洒域喷头的布置方法 |
CN105160096A (zh) * | 2015-08-31 | 2015-12-16 | 江苏大学 | 一种变转速喷头三维动态水量分布模型建立方法 |
CN106577180A (zh) * | 2016-12-13 | 2017-04-26 | 中国科学院地理科学与资源研究所 | 一种快速确定微喷灌喷灌强度及微喷头选择布局的方法 |
CN107907312A (zh) * | 2017-07-31 | 2018-04-13 | 内蒙古自治区水利科学研究院 | 一种快速监测滴灌灌水系统灌水均匀度的方法和冲洗方法 |
CN110472317A (zh) * | 2019-07-31 | 2019-11-19 | 中国石油大学(华东) | 一种用于校园绿地的洒水喷头空间布局优化的方法 |
CN114610828A (zh) * | 2022-03-15 | 2022-06-10 | 河北工程大学 | 一种基于K-means算法的微地形智能分区管理方法 |
-
2022
- 2022-11-14 CN CN202211418096.4A patent/CN115879695B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060144957A1 (en) * | 2004-12-02 | 2006-07-06 | Jacobsen Stephen C | Intelligent sprinkler irrigation system |
CN1908615A (zh) * | 2006-08-11 | 2007-02-07 | 西北农林科技大学 | 一种基于连续系统的喷灌均匀度测量方法 |
CN103506237A (zh) * | 2013-10-11 | 2014-01-15 | 江苏大学 | 一种正方形喷洒域喷头的布置方法 |
CN105160096A (zh) * | 2015-08-31 | 2015-12-16 | 江苏大学 | 一种变转速喷头三维动态水量分布模型建立方法 |
CN106577180A (zh) * | 2016-12-13 | 2017-04-26 | 中国科学院地理科学与资源研究所 | 一种快速确定微喷灌喷灌强度及微喷头选择布局的方法 |
CN107907312A (zh) * | 2017-07-31 | 2018-04-13 | 内蒙古自治区水利科学研究院 | 一种快速监测滴灌灌水系统灌水均匀度的方法和冲洗方法 |
CN110472317A (zh) * | 2019-07-31 | 2019-11-19 | 中国石油大学(华东) | 一种用于校园绿地的洒水喷头空间布局优化的方法 |
CN114610828A (zh) * | 2022-03-15 | 2022-06-10 | 河北工程大学 | 一种基于K-means算法的微地形智能分区管理方法 |
Non-Patent Citations (2)
Title |
---|
DO PRODO等: "Water distribution from medium-size sprinkler in solid set sprinkler systems", 《REVISTA BRASILEIRA DE ENGENHARIA AGRICOLA E ABIENTAL》, vol. 20, no. 3, 31 December 2016 (2016-12-31), pages 195 - 201 * |
劳东青等: "喷头水量分布仿真及组合优化软件系统研究", 《节水灌溉》, no. 1, 5 January 2010 (2010-01-05), pages 42 - 46 * |
Also Published As
Publication number | Publication date |
---|---|
CN115879695B (zh) | 2024-02-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104992072B (zh) | 基于空间网格剖分的作业地块自动识别与面积统计方法 | |
CN105760611B (zh) | 一种低压均匀喷洒喷头空间流道的优化设计方法 | |
CN102789579A (zh) | 基于计算机视觉技术的温室作物水肥胁迫状态识别方法 | |
CN102608620B (zh) | 一种依据反射强度与地形形态的激光扫描点云植被过滤方法 | |
Sanchez et al. | The effects of pressure, nozzle diameter and meteorological conditions on the performance of agricultural impact sprinklers | |
CN107274384B (zh) | 一种小麦出苗均匀度测定方法 | |
Chang et al. | Enforced self‐organizing map neural networks for river flood forecasting | |
CN105160096B (zh) | 一种变转速喷头三维动态水量分布模型建立方法 | |
CN109993420B (zh) | 基于泰森多边形网格的电网负荷可视化方法 | |
Dwomoh et al. | Field performance characteristics of fluidic sprinkler | |
CN105184048B (zh) | 一种基于多源数据的逐年土地利用空间数据的制备方法 | |
CN107845099A (zh) | 一种面向轮胎点云的数据分割方法 | |
Ge et al. | Comparisons of spray characteristics between vertical impact and turbine drive sprinklers—A case study of the 50PYC and HY50 big gun-type sprinklers | |
CN108959746A (zh) | 区域农业用水量核定与分解方法 | |
CN115879695A (zh) | 一种基于K-means算法的灌溉喷头布局方法 | |
Stambouli et al. | Performance of new agricultural impact sprinkler fitted with plastic nozzles | |
Liu et al. | Estimation model of canopy stratification porosity based on morphological characteristics: A case study of cotton | |
CN103439299B (zh) | 一种作物群体光能空间分布量化方法 | |
Khadra et al. | Optimization of on demand pressurized irrigation networks and on-farm constraints | |
CN107369206A (zh) | 一种玉米群体三维模型构建方法及系统 | |
CN105468819B (zh) | 一种利用响应面法优化垂直摇臂式喷头工作参数的方法 | |
CN108510010A (zh) | 一种基于预筛选的密度峰值聚类方法及系统 | |
CN108470015A (zh) | 一种平移式喷灌机组合均匀度的计算方法 | |
Dabbous | A study of sprinkler uniformity evaluation methods | |
CN101826098A (zh) | 一种基于ab直方图的空间查询选择率估计方法 |
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 |