CN106570770A - 果园水肥一体化追肥量精准估算方法 - Google Patents

果园水肥一体化追肥量精准估算方法 Download PDF

Info

Publication number
CN106570770A
CN106570770A CN201610998312.5A CN201610998312A CN106570770A CN 106570770 A CN106570770 A CN 106570770A CN 201610998312 A CN201610998312 A CN 201610998312A CN 106570770 A CN106570770 A CN 106570770A
Authority
CN
China
Prior art keywords
fertilizer
fruit
amount
fruit tree
bud
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
CN201610998312.5A
Other languages
English (en)
Other versions
CN106570770B (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.)
Shandong Agricultural University
Original Assignee
Shandong Agricultural University
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 Shandong Agricultural University filed Critical Shandong Agricultural University
Priority to CN201610998312.5A priority Critical patent/CN106570770B/zh
Publication of CN106570770A publication Critical patent/CN106570770A/zh
Application granted granted Critical
Publication of CN106570770B publication Critical patent/CN106570770B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/02Agriculture; Fishing; Forestry; Mining
    • AHUMAN NECESSITIES
    • A01AGRICULTURE; FORESTRY; ANIMAL HUSBANDRY; HUNTING; TRAPPING; FISHING
    • A01CPLANTING; SOWING; FERTILISING
    • A01C21/00Methods of fertilising, sowing or planting
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/70Denoising; Smoothing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10024Color image
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20024Filtering details
    • G06T2207/20032Median filtering
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20036Morphological image processing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30181Earth observation
    • G06T2207/30188Vegetation; Agriculture

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Soil Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Business, Economics & Management (AREA)
  • Mining & Mineral Resources (AREA)
  • Tourism & Hospitality (AREA)
  • General Health & Medical Sciences (AREA)
  • Human Resources & Organizations (AREA)
  • Marketing (AREA)
  • Primary Health Care (AREA)
  • Strategic Management (AREA)
  • Economics (AREA)
  • Health & Medical Sciences (AREA)
  • General Business, Economics & Management (AREA)
  • Marine Sciences & Fisheries (AREA)
  • Animal Husbandry (AREA)
  • Agronomy & Crop Science (AREA)
  • Environmental Sciences (AREA)
  • Image Analysis (AREA)

Abstract

果园水肥一体化追肥量精准估算方法,包括果树产量预估模块和追肥量计算模块,所述的果树产量预估模块通过图像采集、图像处理和数据分析,利用面积聚类估测果树花芽数量,结合往年疏花经验和商品果划分单果重标准,预估果树产量;所述的追肥量计算模块根据果树产量与施肥量之间的经验关系,得出果树需肥量,结合土壤含肥量、肥料利用率、肥料养分含量,根据公式:实际追肥量=(需肥量‑土壤含肥量)/肥料利用率×肥料养分含量,计算果树实际追肥量。本发明通过果树产量预估、往年经验运用、土壤含肥量检测以及果树实际追肥量计算结合,做到适时、适量追肥,以便充分发挥肥效,满足苹果安全、丰产、优质的需要。

Description

果园水肥一体化追肥量精准估算方法
技术领域
本发明涉及果园管理技术领域,具体涉及一种果园水肥一体化追肥量精准估算方法。
背景技术
施肥是维持土壤肥力,满足果树生长发育、开花结果所需营养元素的重要措施。果树多年生长在同一地点,每年生长、结果都需要从土壤中吸收并消耗大量的无机营养元素。为了满足苹果安全、丰产、优质的需要,果树施肥应根据其需肥规律、树体营养状况、肥料性质、土壤肥力和气候条件而定,尽量做到适时、适量追肥,以便充分发挥肥效。如果盲目追肥,做不到适时适量,不但肥效差,而且会对树体生长和结果造成不良影响(如生长过旺、花芽生成困难、坐果率降低、果实着色不良、果品质量差等)。
果树花前追肥的作用甚为重要,追肥量的多少直接影响果树的开花、坐果和果实发育。传统追肥量的确定是根据经验依据往年果树的产量预测当年的产量,从而确定当年的追肥量,这种确定追肥量的方法缺乏理论依据,无法实现准确科学追肥,影响果树产量。
发明内容
针对上述问题,本发明提供一种果园水肥一体化追肥量精准估算方法,将计算机视觉技术引入果园追肥技术,实现科学精准追肥。
为解决上述技术问题,本发明采取的技术方案为:提供了一种果园水肥一体化追肥量精准估算方法,包括果树产量预估模块和追肥量计算模块,所述的果树产量预估模块通过图像采集、图像处理和数据分析,利用面积聚类估测果树花芽数量,结合往年疏花经验和商品果划分单果重标准,预估果树产量;所述的追肥量计算模块根据果树产量与施肥量之间的经验关系,得出果树需肥量,结合土壤含肥量、肥料利用率、肥料养分含量,根据公式:实际追肥量=(需肥量-土壤含肥量)/肥料利用率×肥料养分含量,计算果树实际追肥量。
具体包括以下步骤:
S1、采集果树图像信息,对果树地上部分做全景采集,得到24位RGB颜色空间的真彩色原始果树图像;
S2、将果树图像由RGB颜色空间转换到Lab颜色空间;
S3、将Lab颜色空间分离成L通道、a通道和b通道;
S4、分别对L、a、b通道求三维像素直方图,选定花芽与背景之间差值明显的通道进行自适应阈值分割,从图像中分离出果树花芽;
S5、通过形态学滤波进行初步降噪,先腐蚀后膨胀处理;
S6、通过中值滤波进一步降噪,即将每个像素的灰度值用其领域的中值代替,去除椒盐噪声;
S7、创建一个动态内存寄存器,对步骤S6处理后的图像进行轮廓检测,找到图像中的所有轮廓,并存放到动态内存寄存器中;
S8、计算所有轮廓中每一个轮廓的面积,根据轮廓面积划分为符合小面积阈值范围的轮廓和符合大面积阈值范围的轮廓;
S9、分别将所有轮廓中符合小面积阈值范围的轮廓绘制到与其外接矩形相同大小的空白图像中,得到单个花芽目标,并统计花芽数量;
S10、分别将所有轮廓中符合大面积阈值范围的轮廓返回至步骤S4,重新进行处理,得到单个花芽目标,并统计花芽数量;
S11、将步骤S9与步骤S10中得到的花芽数量相加得到花芽总数量;
S12、根据步骤S11得到的花芽总数量,结合往年疏花经验得出每个花芽所留花朵数量,从而预估年产果实数量,再根据商品果划分单果重标准规定的单果重量,预估果树产量;
S13、根据果树产量与施肥量之间的经验关系,得出果树需肥量,结合土壤含肥量、肥料利用率、肥料养分含量,根据公式:实际追肥量=(需肥量-土壤含肥量)/肥料利用率×肥料养分含量,计算果树实际追肥量,其中化肥利用率和化肥养分含量均与化肥种类有关。
所述的步骤S2中,通过三维查表插值法将果树图像由RGB颜色空间转换到Lab颜色空间。
所述的步骤S4中,根据三维像素直方图可知,a通道直方图中花芽与土壤、天空等背景之间的差值明显,故选定a通道进行自适应阈值分割,从图像中分离出果树花芽。
所述的步骤S7中,通过遍历的方法对步骤S6处理后的图像进行轮廓检测,找到图像中的所有轮廓。
所述的步骤S8中,符合小面积阈值范围的轮廓为10<轮廓面积<100,符合大面积阈值范围的轮廓为100<轮廓面积<3000;
所述的步骤S12中,商品果划分单果重标准规定的单果重量是指商品果以单果直径划分为75mm、80mm、85mm,其对应的单果重分别为150g、200g、250g。
所述的步骤S13中,通过土壤肥力检测仪检验土壤含肥量,得到土壤中氮、磷、钾肥力的含量信息。
本发明通过果树产量预估、往年经验运用、土壤含肥量检测以及果树实际追肥量计算结合,做到适时、适量追肥,以便充分发挥肥效,满足苹果安全、丰产、优质的需要。
附图说明
图1为本发明的流程图。
具体实施方式
果园水肥一体化追肥量精准估算方法,具体包括以下步骤:
S1、利用数码相机采集果树图像信息,对果树地上部分做全景采集,得到24位RGB颜色空间的真彩色原始果树图像;
S2、通过三维查表插值法将果树图像由RGB颜色空间转换到Lab颜色空间。
S3、将Lab颜色空间分离成L通道、a通道和b通道;
S4、分别对L、a、b通道求三维像素直方图,根据三维像素直方图可知,a通道直方图中花芽与土壤、天空等背景之间的差值明显,故选定a通道进行自适应阈值分割,从图像中分离出果树花芽;
S5、通过形态学滤波进行初步降噪,先腐蚀后膨胀处理;
S6、通过中值滤波进一步降噪,即将每个像素的灰度值用其领域的中值代替,去除椒盐噪声;
S7、创建一个动态内存寄存器,通过遍历的方法对步骤S6处理后的图像进行轮廓检测,找到图像中的所有轮廓,并存放到动态内存寄存器中。
S8、计算所有轮廓中每一个轮廓的面积,根据轮廓面积划分为符合小面积阈值范围的轮廓和符合大面积阈值范围的轮廓,其中符合小面积阈值范围的轮廓为10<轮廓面积<100,符合大面积阈值范围的轮廓为100<轮廓面积<3000;
S9、分别将所有轮廓中符合小面积阈值范围的轮廓绘制到与其外接矩形相同大小的空白图像中,得到单个花芽目标,并统计花芽数量;
S10、分别将所有轮廓中符合大面积阈值范围的轮廓返回至步骤S4,重新进行处理,得到单个花芽目标,并统计花芽数量;
S11、将步骤S9与步骤S10中得到的花芽数量相加得到花芽总数量;
S12、根据步骤S11得到的花芽总数量,结合往年疏花经验得出每个花芽所留花朵数量,即每个花芽留2~3朵花,从而预估年产果实数量,再根据商品果划分单果重标准规定的单果重量,预估果树产量,其中商品果划分单果重标准规定的单果重量是指商品果以单果直径划分为75mm、80mm、85mm,其对应的单果重分别为150g、200g、250g;
S13、根据果树产量与施肥量之间的经验关系,得出果树需肥量,结合土壤含肥量、肥料利用率、肥料养分含量,根据公式:实际追肥量=(需肥量-土壤含肥量)/肥料利用率×肥料养分含量,计算果树实际追肥量,其中通过土壤肥力检测仪检验土壤含肥量,得到土壤中氮、磷、钾肥力的含量信息,化肥利用率和化肥养分含量均与化肥种类有关。

Claims (8)

1.果园水肥一体化追肥量精准估算方法,其特征在于:包括果树产量预估模块和追肥量计算模块,所述的果树产量预估模块通过图像采集、图像处理和数据分析,利用面积聚类估测果树花芽数量,结合往年疏花经验和商品果划分单果重标准,预估果树产量;所述的追肥量计算模块根据果树产量与施肥量之间的经验关系,得出果树需肥量,结合土壤含肥量、肥料利用率、肥料养分含量,根据公式:实际追肥量=(需肥量-土壤含肥量)/肥料利用率×肥料养分含量,计算果树实际追肥量。
2.根据权利要求1所述的果园水肥一体化追肥量精准估算方法,其特征在于具体包括以下步骤:
S1、采集果树图像信息;
S2、将果树图像由RGB颜色空间转换到Lab颜色空间;
S3、将Lab颜色空间分离成L通道、a通道和b通道;
S4、分别对L、a、b通道求三维像素直方图,选定花芽与背景之间差值明显的通道进行自适应阈值分割,从图像中分离出果树花芽;
S5、通过形态学滤波进行初步降噪,先腐蚀后膨胀处理;
S6、通过中值滤波进一步降噪,去除椒盐噪声;
S7、创建一个动态内存寄存器,对步骤S6处理后的图像进行轮廓检测,找到图像中的所有轮廓,并存放到动态内存寄存器中;
S8、计算所有轮廓中每一个轮廓的面积,根据轮廓面积划分为符合小面积阈值范围的轮廓和符合大面积阈值范围的轮廓;
S9、分别将所有轮廓中符合小面积阈值范围的轮廓绘制到与其外接矩形相同大小的空白图像中,得到单个花芽目标,并统计花芽数量;
S10、分别将所有轮廓中符合大面积阈值范围的轮廓返回至步骤S4,重新进行处理,得到单个花芽目标,并统计花芽数量;
S11、将步骤S9与步骤S10中得到的花芽数量相加得到花芽总数量;
S12、根据步骤S11得到的花芽总数量,结合往年疏花经验得出每个花芽所留花朵数量,从而预估年产果实数量,再根据商品果划分单果重标准规定的单果重量,预估果树产量;
S13、根据果树产量与施肥量之间的经验关系,得出果树需肥量,结合土壤含肥量、肥料利用率、肥料养分含量,根据公式:实际追肥量=(需肥量-土壤含肥量)/肥料利用率×肥料养分含量,计算果树实际追肥量。
3.根据权利要求2所述的果园水肥一体化追肥量精准估算方法,其特征在于:所述的步骤S2中,通过三维查表插值法将果树图像由RGB颜色空间转换到Lab颜色空间。
4.根据权利要求2所述的果园水肥一体化追肥量精准估算方法,其特征在于:所述的步骤S4中,选定a通道进行自适应阈值分割,从图像中分离出果树花芽。
5.根据权利要求2所述的果园水肥一体化追肥量精准估算方法,其特征在于:所述的步骤S7中,通过遍历的方法对步骤S6处理后的图像进行轮廓检测,找到图像中的所有轮廓。
6.根据权利要求2所述的果园水肥一体化追肥量精准估算方法,其特征在于:所述的步骤S8中,符合小面积阈值范围的轮廓为10<轮廓面积<100,符合大面积阈值范围的轮廓为100<轮廓面积<3000。
7.根据权利要求2所述的果园水肥一体化追肥量精准估算方法,其特征在于:所述的步骤S12中,商品果划分单果重标准规定的单果重量是指商品果以单果直径划分为75mm、80mm、85mm,其对应的单果重分别为150g、200g、250g。
8.根据权利要求2所述的果园水肥一体化追肥量精准估算方法,其特征在于:所述的步骤S13中,通过土壤肥力检测仪检验土壤含肥量,得到土壤中氮、磷、钾肥力的含量信息。
CN201610998312.5A 2016-11-11 2016-11-11 果园水肥一体化追肥量精准估算方法 Active CN106570770B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610998312.5A CN106570770B (zh) 2016-11-11 2016-11-11 果园水肥一体化追肥量精准估算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610998312.5A CN106570770B (zh) 2016-11-11 2016-11-11 果园水肥一体化追肥量精准估算方法

Publications (2)

Publication Number Publication Date
CN106570770A true CN106570770A (zh) 2017-04-19
CN106570770B CN106570770B (zh) 2020-12-18

Family

ID=58541742

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610998312.5A Active CN106570770B (zh) 2016-11-11 2016-11-11 果园水肥一体化追肥量精准估算方法

Country Status (1)

Country Link
CN (1) CN106570770B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108077042A (zh) * 2017-12-04 2018-05-29 北京农业智能装备技术研究中心 一种冬小麦全生育期灌溉预警决策方法
CN108154267A (zh) * 2017-12-25 2018-06-12 深圳春沐源控股有限公司 一种施肥机的肥料参数调整方法及装置
CN108684278A (zh) * 2017-04-06 2018-10-23 点豆(山东)网络技术有限公司 智能配肥方法、装置和系统
CN112021016A (zh) * 2020-09-07 2020-12-04 江苏农林职业技术学院 草莓疏花疏果系统
CN115349338A (zh) * 2022-08-17 2022-11-18 江苏省农业机械试验鉴定站 一种基于作物长势的变量追肥控制系统与方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101015262A (zh) * 2007-03-09 2007-08-15 陕西中农瑞丰化肥科技有限责任公司 苹果精准施肥方法
CN103336966A (zh) * 2013-07-15 2013-10-02 山东奥泰机械有限公司 一种应用于农业智能机械的杂草图像辨识方法
WO2013168482A1 (ja) * 2012-05-09 2013-11-14 株式会社クボタ 施肥量設定方法、施肥量設定装置、施肥量設定プログラム
WO2013181069A1 (en) * 2012-06-01 2013-12-05 Rowbot Systems Llc Robotic platform method for performing multiple functions in agricultural systems
CN103646404A (zh) * 2013-12-29 2014-03-19 四川农业大学 一种基于颜色的快速茶树花分割和计数方法
CN103793686A (zh) * 2014-01-16 2014-05-14 中国农业大学 一种果树产量早期预测的方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101015262A (zh) * 2007-03-09 2007-08-15 陕西中农瑞丰化肥科技有限责任公司 苹果精准施肥方法
WO2013168482A1 (ja) * 2012-05-09 2013-11-14 株式会社クボタ 施肥量設定方法、施肥量設定装置、施肥量設定プログラム
WO2013181069A1 (en) * 2012-06-01 2013-12-05 Rowbot Systems Llc Robotic platform method for performing multiple functions in agricultural systems
CN103336966A (zh) * 2013-07-15 2013-10-02 山东奥泰机械有限公司 一种应用于农业智能机械的杂草图像辨识方法
CN103646404A (zh) * 2013-12-29 2014-03-19 四川农业大学 一种基于颜色的快速茶树花分割和计数方法
CN103793686A (zh) * 2014-01-16 2014-05-14 中国农业大学 一种果树产量早期预测的方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
杨成明: "苹果平衡施肥技术", 《果农之友》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108684278A (zh) * 2017-04-06 2018-10-23 点豆(山东)网络技术有限公司 智能配肥方法、装置和系统
CN108684278B (zh) * 2017-04-06 2020-07-17 点豆(山东)网络技术有限公司 智能配肥方法、装置和系统
CN108077042A (zh) * 2017-12-04 2018-05-29 北京农业智能装备技术研究中心 一种冬小麦全生育期灌溉预警决策方法
CN108154267A (zh) * 2017-12-25 2018-06-12 深圳春沐源控股有限公司 一种施肥机的肥料参数调整方法及装置
CN112021016A (zh) * 2020-09-07 2020-12-04 江苏农林职业技术学院 草莓疏花疏果系统
CN115349338A (zh) * 2022-08-17 2022-11-18 江苏省农业机械试验鉴定站 一种基于作物长势的变量追肥控制系统与方法

Also Published As

Publication number Publication date
CN106570770B (zh) 2020-12-18

Similar Documents

Publication Publication Date Title
CN106570770A (zh) 果园水肥一体化追肥量精准估算方法
CN102968782B (zh) 一种彩色图像中显著对象的自动抠取方法
CN105701829B (zh) 一种套袋绿色果实图像分割方法
CN110796038B (zh) 结合快速区域生长超像素分割的高光谱遥感图像分类方法
CN111798467A (zh) 一种图像分割方法、装置、设备及存储介质
CN102208099A (zh) 一种抗光照变化的作物彩色图像分割方法
CN103778628A (zh) 田间开放环境下的叶片病害程度获取方法及系统
CN105794382A (zh) 一种基于机器视觉技术的果园精准施肥装置及施肥方法
CN111798470B (zh) 一种应用于智能农业的农作物图像实体分割方法及系统
CN108682039A (zh) 一种双目立体视觉测量方法
CN107545550B (zh) 细胞图像偏色校正方法
CN113255434B (zh) 一种融合果实特征与深度卷积神经网络的苹果识别方法
CN104715251A (zh) 一种基于直方图线性拟合的显著目标检测方法
CN102354388B (zh) 一种自适应计算图像底层特征重要性权重的方法
CN106846343A (zh) 一种基于聚类超像素分割的病理图像特征提取方法
CN106666767B (zh) 一种基于视觉技术的高效葵花籽脱壳方法
CN114881127A (zh) 一种基于高分辨率遥感卫星影像的农作物精细分类方法
CN111401121A (zh) 一种基于超像素特征提取实现柑橘分割方法
US8971619B2 (en) Method and a device for extracting color features
CN106952283B (zh) 图像分割方法及装置
Song et al. The study of adaptive multi threshold segmentation method for apple fruit based on the fractal characteristics
CN104899875A (zh) 一种基于积分匹配的快速图像协同显著区域监测方法
CN106558050A (zh) 一种基于自适应三阈值的显著物体分割方法
CN105160677B (zh) 一种适用于空间非合作目标的多特征图像分割定位方法
CN113807129A (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