基于计算机图像处理技术的叶面特征参数无损测量研究进展

来源 :安徽农业科学 | 被引量 : 0次 | 上传用户:COMMA87730030
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
  摘要相比传统方法计算叶片特征参数,无损图像处理技术无需采摘叶片,在保持植株完整株型的前提下为下一次数据测量提供了保障并应用先进的计算技术准确计算植株的各种特征值,无损图像处理技术逐渐成为国内外学者研究植物生长形态的重要手段。介绍基于数码技术的叶面积计算方法和基于图像处理的叶片轮廓建模方法,概述数码技术应用于叶面积计算的国内外研究进展,总结归纳叶片轮廓提取方法和拟合方法的优缺点,并指出应用图像处理技术计算叶面积仍需进一步解决的问题。
  关键词无损测量;叶面积计算;图像处理;轮廓建模
  中图分类号S176文献标识码A文章编号0517-6611(2014)23-08029-03
  基金项目广东省教育部产学研结合项目(2011B090400510);广东省教育部产学研结合项目(2012B091100480)。
  作者简介林雪梅(1989- ),女,广东阳江人,硕士研究生,研究方向:生物数学与微分方程应用。*通讯作者, 副教授,硕士,硕士生导师,从事花卉栽培、能源植物资源与能源转化技术以及草坪生态与环境管理研究。
  收稿日期20140703叶片是植物光合作用的主要器官,其发育状况和叶面積大小直接影响作物生长、抗病性及产量,叶片特征研究是植物生长研究的重要内容。获取叶片特征参数,一般采用两种方法:一种是传统的破坏性测量,即随机选取某一植株,有损地测量其叶面积、叶长等叶片特征参数及其干、鲜重量。此测量方法的结果虽准确但较为耗时繁琐,不能大规模地搜集数据,且不能对同一植株进行连续的生长测量并破坏植株的生长平衡。二是采用接触式传感器测量植株叶面特征参数等。此方法虽然方便,但是成本高,传感器的开发难度大,且在接触植物的过程中不可避免地影响植物的生长状态。因此,如何对植物生长连续快速地进行无损检测,成为国内外学者的热点研究内容。为此,笔者概述了国内外对无损测量中叶面积计算和叶片轮廓提取与模拟等方面的最新研究进展,并进一步指出相关研究还需解决的问题。
  1叶面积无损测量技术
  1.1图片采样技术图片的采样方法主要包括扫描采样法、相机采样法。其中,扫描采样法在计算叶面积的像素法中应用较广泛。自由路等通过扫描技术应用数字图像法建立了测算叶面积的方法[2-6]。数字图像法能快速、便捷、精确地计算叶面积,但应用也存在局限性:①不能一次测出多片叶子的单个面积;②受参照物面积的影响比较大,有时需要对计算结果进行校正才能得到准确的实际叶面积;③由于从获取图片到计算叶面积整个过程需要较多的时间,并不适合大量的叶面积测量工作;④叶片需要离体测量、破坏植物和采集图像,如果出现倾斜则失真等。2002年鲍雅静等首次提出利用普通相机采集照片通过像素法计算草片的实际面积,但需要将叶片剪下来[7]。
  如何解决叶片图像采样过程中存在有损采样、图像失真以及计算技术等缺陷,成为国内外学者研究的热点。
  1.2无损测量叶面积的计算方法
  1.2.1像素法。像素法是根据图片像素在总像素中所占的比例得出图像面积的一种计算方法,该方法源于网格法,主要应用于数字图像。FematIDiaz等采用扫描仪获得叶片的彩色图像,与图像处理方法相结合,用像素法计算叶面积[8]。
  针对叶片图像采样过程中图像存在有损测量的缺陷,曹志刚等利用数码相机和手持标准面积板获取叶片图像,通过运行VB程序分析图像计算叶面积,进行无损测算[9]。范克俊等对此方法进行了改良,利用数码相机拍照获取叶片和直尺的图像,利用Photoshop软件对图像进行处理, 统计叶片的像素数和每平方厘米面积的像素数,从而计算叶片的实际面积[10]。
  在无损测量中,采用数码相机获取图像的方法得到广泛的应用。马彦平等以数字图像处理技术与Visual C++6.0相结合进行区域标记与区域像素数统计,从而获得玉米叶片实际面积[11]。Chen等用数码相机获取叶片的数字图像,通过图像处理技术获得叶片的图像像素点总数,根据叶片像素点总数与参考面像素点的总数之比得出叶面积系数,从而获得单片叶片的总面积[12]。此方法对于表面有缺陷例如洞、烤焦等的叶面同样适用。
  随着科技的发展,应用智能手机计算叶片面积得到发展。Han等通过智能手机采集叶片图像,以边长为3的长方形为参考物,在手机上通过OTSU方法分离叶片区域和背景,最后根据叶片区域与参考物长方形的像素比例求出该叶片的叶面积。试验结果表明,该方法与传统方法数码相机和电脑相比具有同等精度[13]。
  1.2.2CCD测量面积法。CCD法(Charge-coupled Device)又称为CCD图像传感器,是一种计算不规则平面物体面积的有效方法。它能把光学影像转化为数字信号,将植入的微小光敏物质作像素,像素越多,提供的画面分辨率越高。很多研究学者把CCD图像传感器与计算机相结合,得出所要计算图片中图像像素,从而得到图像的面积。
  张全法等以CCD为图像传感器与计算机相结合统计所得叶片图像的像元数,通过测量选定标准物体的面积得到图像面积的比例系数,从而计算出叶片的实际面积。虽然CCD法具有测量精度高、测量范围大、使用方便等特点,但测量精度会受到很多因素影响,如几何畸变等[14]。针对图像存在几何畸变的缺陷,左欣等基于改进的Hough变换校提取控制点,采用双线性空间变换进行几何校正,从而改进了张全法等的CCD测量面积法,更加适应野外植物叶片面积的非破坏性测量的需要[15]。
  1.2.3CAD图形分析法。随着二维绘图软件AutoCAD的快速发展和普及使用,许多学者利用AutoCAD程序中的area命令计算叶片图像的面积。田青等在对阔叶类和针叶类植物叶面积的测量中,采用数码相机获取叶片的数字图像,应用AutoCAD软件计算植物叶面积[16]。刘浩等利用扫描仪获取番茄和青椒叶片图像后,通过AutoCAD计算叶片实际面积[17]。AutoCAD叶片面积图像处理技术精度高,不受叶片大小、形状和厚度等因素的影响,不仅可以测定像青椒等规则叶片,还可以对如番茄等不规则叶片和针叶类等叶片进行准确测定。   1.2.4Mapinfo方法。Mapinfo是美国Mapinfo公司推出的桌面地理信息系统软件,提供了屏幕矢量化功能,广泛应用于众多领域。吴玉德等是我国把该软件应用于叶面积测算的第一人,通过扫描仪参照标准计算纸采集树木叶片栅格图像,通过Mapinfo软件的屏幕矢量化功能得到树木叶片轮廓矢量图,进而通过该软件的内置菜单得到树木叶面积[18]。Mapinfo方法具有严谨的科学性,其显著优点是快速、精确,可用于所有叶片,不受叶片大小、形状、厚度等因素的影响。不足之处是必须采用标准计算纸,叶面积的测定精度受标准计算纸本身精度的影响比较大,因此必须保证标准计算纸的精度。
  1.2.5机器视觉技术法。机器视觉是一种应用计算机模拟目标外显或宏观视觉功能的技术,是用图像创建或回复显示目标,用计算机来局部实现人的视觉功能和在一定场合拓宽或辅助人的视觉功能。将图像处理与机器视觉技术结合后可以实现连续、快速的无损检测,可对植物叶面积、叶柄夹角、茎杆直径等生长特征参数进行有效检测。
  He采用双目立体视觉,对甘薯幼苗进行了无损检测[19]。陈鼎才等将此方法进一步改进后的方法能更好地适应野外活体叶面积的测量計算,具有较高测量精度和较低测量难度,还可以对叶面积数据进行管理和分析,为作物生长检测提供有力的数据支持,从而具有广泛的应用前景[20]。基于机器视觉的无损测量,通过采集大量的植物生长信息,统计并比较出不同品种植物在不同生长条件或不同生长时期的生长速率,进而指导环境控制,提供合理的生长所需,合理筛选出适合培育的品种,帮助指导实际生产,减少盲目性。
  2图像处理在叶片轮廓建模中的应用
  叶片形状是辨别不同植物的另一个重要特征,是植物识别的重要工具,因此,植物叶片轮廓曲线形状研究十分重要。Wang等应用分割与形态学运算可自动控制标点的分水岭算法,分割具有复杂的背景图像,然后通过所提取的二值图像中的形状特征点,得出完整叶片特征点[21]。Caballero等对植物叶片轮廓的提取进行了探索[22]。Li等指出,采用snakes技术与细胞神经网络系统相结合,可灵活准确地提取植物叶片的叶脉[23]。王静文等采用Williams和Shah提出的贪心算法,通过对轮廓点周围的局部近邻进行检索选取具有更小Snake 模型能量的新轮廓,改进了Snake模型。该方法提取标本叶片轮廓,排除了原来的阴影区域,保留了原来叶片图像的细节,比传统模型提取的轮廓更加精确[24]。
  用图像技术中的轮廓检测有以下几项要求:①轮廓的定位要高;② 不同尺度的轮廓应有良好的响应,并尽量减少漏检;③应对噪声不敏感;④检测灵敏度应受轮廓方向影响小。目前,主要的边缘检测提取方法有边界跟踪法和边缘检测算法。
  2.1边界跟踪法边界跟踪法的基本原理是先根据图像特征,大步距地寻找边缘起始点,再从起始点开始,小步距递进式地寻找出各个要素点,最后将这些要素点连成平滑的曲线。Michael和Suen于1981借鉴该思想提出一个较为完善的边界跟踪算法,但算法的描述以及终止的判别条件并不严格。邓小同等对Michael和Suen的算法做了一些改进,提供了较为严格的数学描述,并给出算法正确性以及封闭性的证明[25]。此后,相关边界跟踪法的研究更多。贺忠海等提出的方法是在图形较简单的情况下,可完成与传统边缘检测方法相同的工作[26]。周秀芝等针对传统边界跟踪算法易出现漏跟踪或者跟踪断裂的确定,利用交叉点描述边界跟踪的不同方向,形成以交叉点为结点的目标边界的树结构描述,通过对树的遍历获得目标边界的跟踪序列[27]。
  传统的边界跟踪算法主要包括虫跟随法、光栅扫描法以及T算法[28]。目前,用于叶片轮廓检测的典型轮廓跟踪法是方向邻域法,主要有四邻域[29]、八邻域[30]、十六邻域[31]。
  2.2边缘检测算子一般导数用来描述连续函数的变化。因此,可以使用偏导数描述边缘的检测算子,并且可以用图像函数最大增长方向的梯度来描述图像函数的变化。边缘检测首先检测出图像局部特性的不连续,然后将这些不连续的边缘像素连成完整的边界。边缘检测算子的实质是检测出符合边缘特性的边缘像素点。
  目前,主要的边缘检测算子有一阶微分边缘检测算子和二阶微分边缘检测算子。其中,一阶微分边缘检测算子包括差分边缘检测算子、Roberts交叉微分算子、Sobel微分算子、Prewitt微分算子。二阶微分边缘检测算子包括Laplacian算子、Wallis算子、LOG算子、canny微分算子。
  现有的边缘检测算子十分成熟,但这些算子都有图像质量及环境的限制。例如,Roberts算子边缘定位准,但对噪声敏感,可用于边缘明显而且噪声较少的图像。Sobel算子对噪声具有平滑作用,能够提供较为精确的边缘信息,但由于局部平均的影响,也会检测出许多伪边缘,且得到的边缘较粗,定位不高。Laplacian算子基于二阶导数的零交叉,具有各向同性,即与坐标轴方向无关,坐标轴旋转后梯度结果不变,因此边缘检测效果较好,但是其对噪声比较敏感[32-34]。
  2.3轮廓线的拟合植物叶形的建模是研究虚拟植物的重要组成部分,定量描述植物叶片形态结构可更好地确定植物叶片形状的几何模型。同时,拟合叶片的轮廓曲线可更好地刻画植物叶片形状。目前,关于拟合叶片曲线的研究已取得了一定的进展。
  根据不同叶片轮廓表现出来的不同形状特征,研究者用不同的数学方法解决了相应的数学问题。Chi等将叶片的轮廓表示成径向角的函数,采用三次Bezier曲线作为叶形几何建模的基本元素,既提高了计算量,又复杂了叶片轮廓的几何特征[35]。曾茂林等以植物叶片轮廓的建模为例,先构造出刻画叶片轮廓形态的形状函数,再利用小波分析法快速重构不同叶片的轮廓形状函数,成功地对植物叶片轮廓进行几何建模。该方法具有数据存储量小、建模速度快等优点,并且适用于锯齿形状的叶片[36]。   拟合轮廓线的方法中,应用最为广泛的是多项式样条函数。例如,苏礼楷等通过边缘扫描算法提取边界点,然后利用三次样条曲线建立了植物叶片的几何模型,具有数据量少、几何特征简单等优点,但不能很好地适用于锯齿形叶片[37]。鉴于该思想,谢秋菊等利用边缘检测技术提取出大豆叶片形态数据,通过计算机图形学中的三次参数样条曲线拟合大豆叶片的轮廓曲线,得到了更具有真实感的大豆叶片模型,为静态虚拟植物叶片建模研究提供必要的方法,具有一定的实践价值[38]。除此之外,王静文等提取叶片图像的角点作为轮廓特征点,利用三次均匀B样条曲线进行初步插值拟合,并根据拟合误差采用改进的节点插入算法进行反复迭代,达到所要求的拟合精度。该方法能够更快地达到逼近精度,具有良好的模拟效果,有一定的实际应用价值[39]。杜德彭等结合三次均匀B样条曲线的特征,利用功能特征组合得到基本型值点,拟合出叶片的轮廓曲线,并反求控制点,得到了基本叶形特征的轮廓模型,提高了对叶片轮廓边缘的控制和交互调节[40]。
  3结语
  综上所述,通过数码技术离题采集叶片图像,利用图像处理技术计算叶面积,不仅快速、可靠和准确,而且降低成本,可推广到其他大田作物,应用到实际生产生活中,为农业生产生活及广大农业科研人员提供切实可行的科学方法。但与传统的抛物线法以及复印称重法相比较,用图像处理法计算叶片面积的结果与实际的结果有一定的差距,原因在于图像处理方面还没有很好地解决叶片倾向、叶片变形、叶片颤动和自然光变化对图像带来的阈值分割误差等问题,读取图片过程容易丢失数据,计算叶面积不免存在误差。所以,今后的研究任务是:选取可靠、快捷、准确的彩色图像获取方法,探索出基于图像处理的最佳轮廓提取及其拟合方法,并且把叶面积无损测量与叶片轮廓建模有机地结合起来,把抛物线法计算叶面积的原理应用到叶片面积的计算中,应用积分理论计算叶长,实现特征值的计算机可视化,不仅可以提高计算精度,而且可为计算叶片特征值提供新方法、新思路。
  42卷23期林雪梅等基于计算机图像处理技术的叶面特征参数无损测量研究进展参考文献
  [1] 冯冬霞, 施生锦.叶面积测定方法的研究效果初报[J].中国农学通报, 2005, 121(6):150-152,155.
  [2] 自由路, 楊俐苹.基于图像处理的植物叶面积测定方法[J].农业网络信息 ,2004(1):36-38.
  [3] 肖强, 叶文景, 朱珠, 等.利用数码相机和Photoshop软件非破坏性测定叶面积的简便方法[J].生态学杂志, 2005, 24(6):711-714.
  [4] 杨劲峰, 陈清, 韩晓日, 等.数字图像处理技术在蔬菜叶面积测量中的应用[J].农业工程学报, 2002, 18(4):155-158.
  [5] 石剑飞, 殷璀艳, 冷锁虎, 等.采用数码图像处理法测定油菜叶面积的方法探讨[J].中国油料作物学报, 2010, 32(3):379-382.
  [6] 李宝光, 陶秀花, 倪国平, 等.扫描像素法测定植物叶面积的研究[J].江西农业学报,2006, 18(3):78-81.
  [7] 鲍雅静, 李政海, 张颖.羊草叶面积测量方法的比较[J].内蒙古大学学报:自然科学版,2002, 33(1):62-64.
  [8] FEMATDIAZ A,VARGASVAZQUEZ D,HUERTAMANZANILLA E, etal.Scanner image methodology (SIM) to measure dimensions of leaves for agronomical applications[J].African Journal of Biotechnology, 2013, 10(10):1840-1847.
  [9] 曹志刚, 冯仲科, 龙春玲, 等.基于数字摄影与计算机图形技术的植物叶面积测量方法[J].河北科技师范学院学报, 2004, 18(1):39-41.
  [10] 苑克俊, 刘庆忠, 李圣龙, 等.利用数码相机测定果树叶面积的新方法[J].园艺学报, 2006, 33(4):829-832.
  [11] 马彦平, 白由路, 高祥照, 等.基于数字图像的玉米叶面积测量方法研究[J].中国农学通报, 2009, 25(22):329-334.
  [12] CHEN B,FU Z, PAN Y, et al.Single leaf area measurement using digital camera image[M]//Computer and Computing Technologies in Agriculture IV.Springer Berlin Heidelberg, 2011:525-530.
  [13] HAN D,ZHANG F.Leaf area measurement embeded in smart mobile phone[M]//Advances in Intelligent Systems.Springer Berlin Heidelberg, 2012:207-213.
  [14] 张全法, 冯绚, 何金田, 等.植物叶面积测量系统的设计及应用[J].河南农业大学学报, 2001, 35(4):383-386.
  [15] 左欣, 韩斌, 程嘉林.基于数字图像处理的植物叶面积测量方法[J].计算机工程与应用, 2006, 27(2):194-221.
  [16] 田青, 曹致中, 张睿.基于数码相机和AutoCAD软件测定园林植物叶面积的简便方法[J].草原与草坪, 2008(3):25-28.
  [17] 刘浩, 孙景生, 段爱旺, 等.基于AutoCAD软件确定番茄与青椒叶面积的简易方法[J].中国农学通报, 2009, 25(5):287-293.   [18] 吳玉德, 张鹏.基于Mapinfo的树木叶面积测定方法[J].林业调查规划, 2005, 30(6):23-24.
  [19] HE D.Development of a plug-transplant production system and an image analysis system for use in sweet potato[D].Japan:Chiba University, 2002.
  [20] 陈鼎才, 王定成, 查金水.基于机器视觉的现实叶片面积测量方法的研究[J].计算机应用, 2006, 26(5):1226-1228.
  [21] WANG X F,HUANG D S,DU J X,et al.Classification of plant leaf images with complicated background[J].Applied Mathematics and Computation 2008, 205(2):916-926.
  [22] CABALLERO C,ARANDA M C.Plant species identification using leaf image retrieval[C]//Proceedings of the ACM International Conference on Image and Video Retrieval.ACM, 2010:327-334.
  [23] LI Y F,ZHU Q S,CAO Y K,et al.A leaf vein extraction method based on cnakes technique[C]//China Neural Networks Council(CNNC),IEEE Computational Intelligence Society Beijing Chapter.Proceedings of 2005 International Conference on Neural Networks & Brain(ICNN&B05).China Neural Networks Council(CNNC)、IEEE Computational Intelligence Society Beijing Chapter,2005:4.
其他文献
摘要[目的]探讨安烟1号高产潜力,为其大面积推广提供科学依据。[方法]介绍安烟1号亲本来源与选育过程、主要特征特性(主要植物学性状和农艺性状、主要经济性状、原烟外观质量、烟叶物理特性、烟叶化学成分、原烟感官质量和抗病性)、栽培技术要点及适宜种植区域。[结果]安烟1号是以优质烤烟品种K394的同型细胞质雄性不育系MSK394为母本、以自育优质抗病烤烟新品系3033为父本配制的杂种一代。多年多点试验结
期刊
摘要[目的]探究酒精对设施番茄土壤根结线虫的防治效果。[方法]通过用不同浓度酒精喷施结合地膜覆盖处理种植10年的设施土壤,研究了酒精对设施番茄土壤根结线虫数量、微生物生物量碳的变化特征、番茄的生长发育和土壤养分含量的影响。[结果]酒精处理土壤在短期内对土壤根结线虫有明显的抑制作用,对土壤的养分含量影响不明显,对番茄的生长发育未造成不利影响。[结论]为设施栽培土壤根结线虫的防治提供了理论依据。  关
期刊
摘要 基于JCR 2012年数据,选择由我国主办且在生物学领域内影响因子位于前7名的科技期刊进行分析,分别统计了这7种期刊的主办单位、创刊年、出版商、2008~2012年5年内影响因子的变化情况、出版周期、载文量、篇均面数、栏目设置、参考文献数等以及重要办刊策略的调整,为我国创办高水平的国际期刊提供参考。  关键词 生物学;影响因子;科技期刊  中图分类号S-058文献标识码A文章编号0517-6
期刊
摘要对10头0~12月龄半散养雄性林麝的体重和体尺进行测定,分析了其生长发育规律。结果表明,林麝体重累积生长曲线上升由急剧到平缓。0~4月龄和10~12月龄绝对生长速度呈现明显的上升状态,绝对增重在1月龄达到最大,此后逐渐减慢;0~3月龄生长强度最大,3~6月龄呈直线下降趋势,8月龄开始变得平缓:各项体尺指标在6月龄前增长明显,6~12月龄增长趋缓,6~12月龄肩高变化不大。该研究揭示了半散养林麝
期刊
摘要通过查阅文献资料、解剖图对比、解剖镜观察等方法,对菜豆象、四纹豆象、鹰嘴豆象、绿豆象4种豆象进行了区别鉴定,从外形特征和雄性生殖器等4个方面绘制了识别表,并对该类豆象的防治措施进行了总结。  关键词豆象;识别;防治;检疫措施  中图分类号S433.5文献标识码A文章編号0517-6611(2014)23-07812-02  作者简介朱磊(1979-),男,河南信阳人,农艺师,从事植保植检工作。
期刊
摘要采用超声波技术测定活牛眼肌面积、背膘厚度、眼肌高度和肌间脂肪含量,建立了产肉性能的预测模型,数值的可靠性与准确性高,经方差分析达到极显著差异。该研究为肉牛育种与肉质评定提供了科学依据。  关键词新疆褐牛;眼肌面积;背膘厚度  中图分类号S823文献标识码A文章编号0517-6611(2014)23-07764-05  基金项目自治区科技重大专项课题(20123011610)。  作者简介张杨(
期刊
摘要[目的]了解重要蛀干害虫——云斑天牛的生活习性。[方法]从生物学特性、补充营养、发生规律、综合防治措施4个方面综述了云斑天牛近年来的研究进展。[结果]云斑天牛主要有成虫、卵、幼虫、蛹4种形态;每2~3年发生1代,主要以成虫或幼虫在蛀道中越冬;成蟲食叶和嫩枝皮使新枝枯死,幼虫蛀食枝干皮层和木质部,孵幼虫把皮层蛀成三角形蛀道,木屑和粪便从蛀孔排出,树皮外胀纵裂,是识别云斑天牛为害的重要特征;防治技
期刊
摘要对施秉县主产烟乡镇烟草病虫害进行了调查,并提出了相应的防治对策。结果表明,在施秉地区烟叶生产中常见的病虫害,如烟蚜、烟青虫、TMV、CMV、PVY、青枯病、黑胫病等都有发生。  关键词烟草;病虫害;调查  中图分类号S572文献标识码A文章编号0517-6611(2014)23-07822-03  作者简介陈杰(1986- ),男,湖南衡阳人,硕士,从事烟叶科技推广应用工作。  收稿日期201
期刊
摘要运用 SWOT 方法对云南木薯(Manihot esculenta Crantz)产业的内部条件和外部环境作了详细分析,客观地指出了云南木薯产业发展具有的优势、劣势、机遇与威胁,并在此基础上提出云南木薯产业可持续发展的对策。  关键词木薯产业;现状; SWOT;对策;云南省  中图分类号S632文献标识码A文章编号0517-6611(2014)23-07802-02  作者简介韩本勇(1974
期刊
摘要[目的] 为探索建立邓州市小麦施肥指标体系,指导邓州小麦科学施肥。[方法] 在小麦上开展3414试验、氮肥用量试验、氮、磷、钾丰缺指标试验的基础上,分别建立小麦氮、磷、钾丰缺函数;以土壤有机质含量和土壤有机质丰缺函数为依据,确定耕地肥力水平;以氮肥用量试验,确定不同有机质水平下的氮肥用量;以磷、钾丰缺函数,确定土壤有效磷、速效钾临界值;以3414试验处理④、⑤、⑥、⑦小区产量与对应磷肥用量建立
期刊