技术领域
[0001] 本发明属于植物叶片面积测量领域，具体涉及一种叶根、叶中段和叶尖分别测量的叶面积精确计算方法。本申请是2022年2月11日提交的专利申请2022101302446《采用最小二乘法拟合参数的叶面积分类计算方法》的分案申请。
相关背景技术
[0002] 叶片是植物进行光合作用和蒸腾作用的主要器官，有关植物生长情况的诸多信息都与叶片的形态特征有着密不可分的联系。针对植物叶片参数进行的快速、精准的测量是研究植物生长发育等生命活动的重要手段。
[0003] 植物叶面积是叶片测量中的主要参数之一，叶面积的大小影响着植物光合物质的积累，叶面积的变化也体现着植物生长的状态，是研究植物生理生化、遗传育种、作物栽培技术等方面的重要指标。因此植物叶面积的测量对于调整群体结构、充分利用光热资源、指导作物栽培密度及合理施肥等有重要意义。
[0004] 现有的叶面积测量方法中，传统的手工测量方法如复印称重法、直尺测量法等测量速度慢，且测量精度低。利用计算机技术测量叶面积方法大体可分为图像测量法、光电扫描法、三维点云法。通过摄像机采集图像，对图像数据进行数字化处理的图像测量法是较为便捷的方法，但现有的叶面积图像测量方法精确度不高，易受光照等环境因素的干扰，且自身算法复杂，适应性较差，一般作为叶面积的粗略测量。使用光电器件，如扫描仪对叶片进行扫描成像，获取的图像数据质量较高，特别是图像背景为纯白色，运用较为简单的算法也能获得高精度的测量结果，但此类设备形态类似扫描仪，体积较大，难以用于室外测量叶面积，一般适用于实验室内使用，而且叶片运输回实验室的途中不易保存，叶片失水造成叶片缩小、卷边，导致测量误差。三维点云法运用kinect或激光雷达设备采集植物的三维点云形态信息，该设备能获取丰富的植物表型数据：但该种方式噪声干扰多，对于叶面积高精度测量影响较大，且算法复杂，操作困难，成本较高，难以普及。
[0005] 因此，研究一种新的准确、高效、便捷的叶面积测量计算方法，且叶面积测量计算过程不接触植物叶片，对植物无损。
具体实施方式
[0043] 实施例一
[0044] 叶根、叶中段和叶尖分别测量的叶面积精确计算方法，包括以下步骤：
[0045] 采集获得叶子图像，从中提取出完整的单片叶子的图像；
[0046] 将单片叶子的图像输入叶子识别模型，得到叶子的类别；
[0047] 沿叶片中心线方向将叶片划分为叶根区1、叶中段2和叶尖区3，将叶根区1划分为第一叶根部111、第二叶根部112和第三叶根部113；
[0048] 针对每种类别的叶片，分别计算叶根区、叶中段、叶尖区对应的三角形或四边形的面积，并计算得到叶片面积；
[0049] 针对不同类别的叶片，将叶片面积作为因变量，将叶根区、叶中段、叶尖区对应的多边形的面积分别作为自变量，采用最小二乘法对各个自变量的系数进行拟合，得到不同类别叶子的叶片面积计算公式；
[0050] 根据叶片面积计算公式计算得到待测叶子的面积，具体包括：
[0051] 将待测叶子的图像输入叶子识别模型，得到待测叶子的类别；
[0052] 分别计算叶根区、叶中段、叶尖区对应的三角形或四边形的面积；
[0053] 根据待测叶子的类别选择对应的叶片面积计算公式，计算得到待测叶子的面积。
[0054] 实施例中，叶子的类别包括桃心形、类椭圆形。
[0055] 桃心形叶片面积计算公式如下
[0056] SP＝α1S111+α2S112+α3S113+α4S2+α5S3+d1
[0057] 式中SP表示桃心形叶片的总面积，S111表示叶根区的第一叶根部111对应的三角形的面积，S112表示叶根区的第二叶根部112对应的四边形的面积，S113表示叶根区的第三叶根部113对应的三角形的面积，S2表示叶中段2对应的四边形的面积，S3表示叶尖区3对应的三角形的面积；α1表示第一叶根部三角形面积的拟合系数，α2表示第二叶根部四边形面积的拟合系数，α3表示第三叶根部三角形面积的拟合系数，α4表示叶中段四边形面积的拟合系数，α5表示叶尖区三角形面积的拟合系数，d1表示桃心形叶片面积拟合的截距。
[0058] 对于桃心形叶片，如图1所示，A点、B点为叶根区边缘点，叶尖为C点，线段CQ与线段AB垂直，垂足为Q点，线段CQ即为叶片的高，即LCQ＝h；线段EF与线段CQ垂直，垂足为R点，线段GH与线段CQ垂直，垂足为T点；R点、T点将线段CQ三等分，即
[0059]
[0060] 分别自A点、B点向线段EF作垂线，垂足为U、V点，LEU＝w1，LUV＝w2，LVF＝w3，LGH＝w4；
[0061]
[0062]
[0063]
[0064]
[0065]
[0066] 实施例中，提取出完整的单片叶子的图像后，对单片叶子的图像进行整体旋转变换，使得叶子图像中的线段AB与图像坐标系的Y轴平行，以图像坐标系为基准，获取桃心形叶片中的w1、w2、w3、w4和h的像素长度，并根据像素坐标系与世界坐标系的转换关系，得到实际的w1、w2、w3、w4和h的长度。
[0067] 以测取柑橘叶片的高度h为例，将拍摄叶片图像的相机相对柑橘树在最佳拍摄距离的位置通过三角架固定，并测量相机镜头距离拍摄柑橘叶片的距离即物距u，对柑橘叶片进行拍照得到柑橘叶片图像，保持相机焦距不变，将画有2cm×2cm正方形的空白纸放置在柑橘树的叶片位置并利用相机拍照，得到柑橘叶片位置的正方形图像。利用相机成像公式：其中u表示物距即叶片或正方形离相机镜头的距离，v表示像距即叶片的图
像或正方形的图像离相机镜头的距离，l像素表示相机拍摄的图像中的像素尺寸，l实际表示图像对应的实体的实际尺寸。焦距f可从拍摄的柑橘叶片图片的属性参数中读取，根据相机拍摄的正方形图像的像素尺寸l实际以及正方形的实际尺寸，计算得到像距v，并计算像距v与物距u的比值，即得到此相机的像素尺寸与拍摄的实体实际尺寸的比例 利用现有技术的图像处理技术，得到柑橘叶片图像中叶片高度h的像素大小，再结合比例 即得到柑橘叶片的高度h的实际大小。
[0068] 类椭圆形叶片面积计算公式如下
[0069] SE＝β1S1+β2S2+β3S3+d2
[0070] 式中SE表示类椭圆形叶片的总面积，S1表示叶根区对应的三角形面积，S2表示叶中段对应的四边形的面积，S3表示叶尖区对应的三角形的面积；β1表示叶根区三角形面积的拟合系数，β2表示叶中段四边形面积的拟合系数，β3表示叶尖区三角形面积的拟合系数，d2表示类椭圆形叶片面积拟合的截距。
[0071] 如图2所示，M点、N点为类椭圆叶片两端的端点，线段JK与线段MN垂直，垂足为Y点，线段OP与线段MN垂直，垂足为Z点；Y点、Z点将线段MN三等分；L表示类椭圆叶片的长度；线段JK的长度LJK＝b1，线段OP的长度LOP＝b2；
[0072]
[0073]
[0074]
[0075] 与上述的测取桃心形叶片的尺寸数据相类似，对类椭圆形叶片的图像进行整体旋转变换，使得叶子图像中的线段MN与图像坐标系的X轴平行，以图像坐标系为基准，获取类椭圆形叶片中的b1、b2和L的像素长度，并根据图像坐标系与世界坐标系的转换关系，得到实际的b1、b2和L的长度。
[0076] 实施例的步骤5中，采摘多片叶子作为样本，测量得到叶子样本的真实面积，利用叶子样本的真实面积对叶片面积计算公式中自变量的系数进行拟合。如图3所示，测量得到叶子样本的真实面积，具体过程如下：
[0077] 1)在空白纸上画边长为2cm的正方形，将其作为比对标尺；
[0078] 2)将叶子样本放置在正方形旁边，利用摄像设备获取叶片和正方形的图像，如图4所示；
[0079] 3)将获取的叶片和正方形的图像导入到r2v矢量软件，得到叶片和正方形的矢量图像，如图5所示；
[0080] 4)将叶片和正方形的矢量图像输入到AUTO CAD软件，如图6所示，利用AUTO CAD软件得到叶片的矢量图像的面积SCAD以及正方形的矢量图像的面积SQ；
[0081] 5)计算得到叶子样本的真实面积，
[0082]
[0083] 式中SR表示叶子样本的真实面积。
[0084] 实施例中，叶子识别模型采用现有技术的卷积神经网络CNN，如2021年第10期《中国图象图形学报》刊登的张珂、冯晓晗等人的论文“图像分类的深度卷积神经网络模型综述”公开的卷积神经网络。实施例的叶面积计算公式中叶片各部位面积的系数的拟合采用现有技术的最小二乘法，如2019年第4期《无线互联科技》刊登的莫小琴的论文“基于最小二乘法的线性与非线性拟合”公开的最小二乘法。
[0085] 如图7所示，以桃心形叶片的叶片面积计算公式为例，利用实验室采集的多个桃心形叶片样本，拟合得到叶片面积计算公式中自变量的系数的具体过程包括：
[0086] 1)提取出完整的单片叶子的图像后，对单片叶子的图像进行整体旋转变换，使得叶子图像中的线段AB与图像坐标系的Y轴平行，以图像坐标系为基准，获取桃心形叶片中的w1、w2、w3、w4和h的像素长度；
[0087] 2)并根据图像坐标系与世界坐标系的转换关系，得到各个叶片样本实际的w1、w2、w3、w4和h的长度；
[0088] 3)根据得到的w1、w2、w3、w4和h的实际长度，算出叶片样本的S111、S112、S113、S2、S3，利用EXCEL对多个叶片样本的面积进行数据整合；
[0089] 4)将S1、S2、S3、S111、S112、S113作为变量以及对应的不同叶片样本的数据导入模型，得到变量的序号，并依次计算变量参数的最小二乘估计，判断是否完成所有变量参数的估计计算，若完成所有变量参数的计算即模型中仅剩截距时，结束，否则执行步骤5)；
[0090] 5)随机进行方差分析以及模型拟合检验和各变量T检验，找出数字变量的偏回归方和的最小值和单变量的最大值；计算用于变量检验的p‑value值大小，进行显著性水平检验，根据p‑value值大小判断是否剔除变量。当无变量可剔除时则结束，最终得到回归方程及其系数。
[0091] 实施例二
[0092] 如图3所示，柑橘叶片面积测量方法，包括以下步骤：
[0093] 1)在空白纸上画边长为2cm的正方形，将其作为比对标尺；
[0094] 2)将叶子样本放置在正方形旁边，利用摄像设备获取叶片和正方形的图像，如图4所示；
[0095] 3)将获取的叶片和正方形的图像导入到r2v矢量软件，得到叶片和正方形的矢量图像，如图5所示；
[0096] 4)将叶片和正方形的矢量图像输入到AUTO CAD软件，如图6所示，利用AUTO CAD软件得到叶片的矢量图像的面积SCAD以及正方形的矢量图像的面积SQ；
[0097] 5)计算得到叶子样本的真实面积，
[0098]
[0099] 式中SR表示叶子样本的真实面积。
[0100] 实施例三
[0101] 如图1所示，桃心型叶片的面积测量方法，包括以下步骤：
[0102] 步骤1：沿叶片中心线方向将叶片划分为叶根区1、叶中段2和叶尖区3，将叶根区1划分为第一叶根部111、第二叶根部112和第三叶根部113；
[0103] 步骤2：分别计算叶根区、叶中段、叶尖区对应的三角形或四边形的面积，并计算得到叶片面积；
[0104] 步骤3：将叶片面积作为因变量，将叶根区、叶中段、叶尖区对应的多边形的面积分别作为自变量，采用最小二乘法对各个自变量的系数进行拟合，得到不同类别叶子的叶片面积计算公式；
[0105] 步骤4：根据叶片面积计算公式计算得到待测叶子的面积。
[0106] 桃心形叶片面积计算公式如下
[0107] SP＝α1S111+α2S112+α3S113+α4S2+α5S3+d1
[0108] 式中SP表示桃心形叶片的总面积，S111表示叶根区的第一叶根部对应的三角形的面积，S112表示叶根区的第二叶根部对应的四边形的面积，S113表示叶根区的第三叶根部对应的三角形的面积，S2表示叶中段对应的四边形的面积，S3表示叶尖区对应的三角形的面积；α1表示第一叶根部三角形面积的拟合系数，α2表示第二叶根部四边形面积的拟合系数，α3表示第三叶根部三角形面积的拟合系数，α4表示叶中段四边形面积的拟合系数，α5表示叶尖区三角形面积的拟合系数，d1表示桃心形叶片面积拟合的截距。
[0109] 对于桃心形叶片，如图1所示，A点、B点为叶根区边缘点，叶尖为C点，线段CQ与线段AB垂直，垂足为Q点，线段CQ即为叶片的高，即LCQ＝h；线段EF与线段CQ垂直，垂足为R点，线段GH与线段CQ垂直，垂足为T点；R点、T点将线段CQ三等分，即
[0110]
[0111] 分别自A点、B点向线段EF作垂线，垂足为U、V点，LEU＝w1，LUV＝w2，LVF＝w3，LGH＝w4；
[0112]
[0113]
[0114]
[0115]
[0116]
[0117] 实施例中，提取出完整的单片叶子的图像后，对单片叶子的图像进行整体旋转变换，使得叶子图像中的线段AB与图像坐标系的Y轴平行，以图像坐标系为基准，获取桃心形叶片中的w1、w2、w3、w4和h的像素长度，并根据图像坐标系与世界坐标系的转换关系，得到实际的w1、w2、w3、w4和h的长度。
[0118] 以测取柑橘叶片的高度h为例，将拍摄叶片图像的相机相对柑橘树在最佳拍摄距离的位置通过三角架固定，并测量相机镜头距离拍摄柑橘叶片的距离即物距u，对柑橘叶片进行拍照得到柑橘叶片图像，保持相机焦距不变，将画有2cm×2cm正方形的空白纸放置在柑橘树的叶片位置并利用相机拍照，得到柑橘叶片位置的正方形图像。利用相机成像公式：其中u表示物距即叶片或正方形离相机镜头的距离，v表示像距即叶片的图
像或正方形的图像离相机镜头的距离，l像素表示相机拍摄的图像中的像素尺寸，l实际表示图像对应的实体的实际尺寸。焦距f可从拍摄的柑橘叶片图片的属性参数中读取，根据相机拍摄的正方形图像的像素尺寸l实际以及正方形的实际尺寸，计算得到像距v，并计算像距v与物距u的比值，即得到此相机的像素尺寸与拍摄的实体实际尺寸的比例 利用现有技术的图像处理技术，得到柑橘叶片图像中叶片高度h的像素大小，再结合比例 即得到柑橘叶片的高度h的实际大小。
[0119] 测量得到叶子样本的真实面积的方法与实施例一中的测量得到叶子样本的真实面积的方法相同。
[0120] 实施例四
[0121] 如图2所示，类椭圆型叶片的面积测量方法，包括以下步骤：
[0122] 步骤1：沿叶片中心线方向将叶片划分为叶根区1、叶中段2和叶尖区3，将叶根区1划分为第一叶根部111、第二叶根部112和第三叶根部113；
[0123] 步骤2：分别计算叶根区、叶中段、叶尖区对应的三角形或四边形的面积，并计算得到叶片面积；
[0124] 步骤3：将叶片面积作为因变量，将叶根区、叶中段、叶尖区对应的多边形的面积分别作为自变量，采用最小二乘法对各个自变量的系数进行拟合，得到不同类别叶子的叶片面积计算公式；
[0125] 步骤4：根据叶片面积计算公式计算得到待测叶子的面积。
[0126] 类椭圆形叶片面积计算公式如下
[0127] SE＝β1S1+β2S2+β3S3+d2
[0128] 式中SE表示类椭圆形叶片的总面积，S1表示叶根区对应的三角形面积，S2表示叶中段对应的四边形的面积，S3表示叶尖区对应的三角形的面积；β1表示叶根区三角形面积的拟合系数，β2表示叶中段四边形面积的拟合系数，β3表示叶尖区三角形面积的拟合系数，d2表示类椭圆形叶片面积拟合的截距。
[0129] 如图2所示，M点、N点为类椭圆叶片两端的端点，线段JK与线段MN垂直，垂足为Y点，线段OP与线段MN垂直，垂足为Z点；Y点、Z点将线段MN三等分；L表示类椭圆叶片的长度；线段JK的长度LJK＝b1，线段OP的长度LOP＝b2；
[0130]
[0131]
[0132]
[0133] 与上述的测取桃心形叶片的尺寸数据相类似，对类椭圆形叶片的图像进行整体旋转变换，使得叶子图像中的线段MN与图像坐标系的X轴平行，以图像坐标系为基准，获取类椭圆形叶片中的b1、b2和L的像素长度，并根据图像坐标系与世界坐标系的转换关系，得到实际的b1、b2和L的长度。
