一种玉米果穗性状测量装置和测量玉米穗行数、穗行倾角以及穗缘角的方法

技术领域

本发明涉及一种在农业科研领域内使用的测量装置，具体的说涉及一种玉米果穗性状测量装置，以及使用该装置测量玉米穗行数及穗行倾角和穗缘角的方法。 

背景技术

果穗是玉米重要的植物器官，植物器官的大小取决于细胞的分裂与生长，形状主要依赖于细胞分化，颜色则与显色物质的代谢密切相关，纹理则是细胞分裂、生长、分化与代谢相互作用的最终形态体现。 

在玉米生产与科研活动中，玉米穗行数是最常见的测定性状。 

在玉米种质资源评价和品种权保护工作中，籽粒排列形式是一个重要性状。一般有三种类型：直型、螺旋状和不规则。果穗籽粒穗行倾角是确定果穗籽粒排列形式的直接依据之一，果穗的穗行倾角即为果穗中轴线与果穗籽粒行之间的夹角。籽粒行倾角大于等于某一阈值，则可视为螺旋状；小于该阈值，则可视为直型；没有明显的籽粒行，则视为不规则。 

在玉米种质资源评价和品种权保护工作中，果穗形状是又一个重要性状。一般有三种穗形：圆锥形、圆筒形和中间型。穗缘角是确定果穗形状的最重要依据之一，穗缘角即为果穗中轴线与果穗边缘线之间的夹角。穗缘角大于等于某一阈值，则可视为圆锥形；小于该阈值，而且，果穗主体部分的矩形度大于某一阈值，则可视为筒形；其它情形，可视为中间型。 

目前，上述性状主要靠人工目测计数，但如采用图像处理方法测量植物形态，不仅可以替代部分人工测量，获得更高精度的结果，而且还能够提供新的候选性状。虽然有文献报道计算机图像处理方法，但是，由于方法固有的局限性，鲁棒性差，实际应用不广泛。中国专利申请201010210313.1公开了一种基于机器视觉的玉米果穗形状参数自动检测装置，中国专利申请201010226335.7公开了一种基于机器视觉的玉米考种装置，但这些专利申请中并没有公开具体的玉米果穗性状的测定方法。 

发明内容

本发明所要解决的技术问题就是提供一种通过采用计算机图像处理方法来实现对玉米果穗的穗行数及穗行倾角和穗缘角进行测量的装置和方法。 

本发明采用如下技术方案： 

一种玉米果穗性状测量装置，该测量装置包括： 

用来采集玉米性状图像信息的采集单元； 

与上述的采集单元相连接的，用于处理上述图像信息并获得最终测量结果的处理单元。 

进一步的，所述的采集单元包括： 

一面可打开的封闭灯箱，该封闭灯箱的底面设有转盘，顶面设有钩子，并且在顶面和一侧面开有摄像孔，在底面的四角装有光源； 

通过数据线与处理单元连接的相机，该相机的镜头可以伸入摄像孔内； 

通过数据线与处理单元连接的扫描仪，该扫描仪的盖板敞开，并有遮光罩罩在扫描区域上。 

进一步的，所述的处理单元为计算机。 

进一步的，所述的光源为白炽灯光源。 

一种玉米穗行数测量方法，该方法包括如下步骤： 

(1)将玉米从中部掰断，利用权利要求1所述的测量装置测得玉米果穗横切面的图像； 

(2)将上述图像转化为二值图像BW并作Radon变换，取过圆点的各直径长度，长度的中位数为外圆周半径； 

(3)玉米籽粒环厚度T等于籽粒深度D乘以k，k∈[0.5，1.0]，籽粒深度D为用户估计值； 

(4)内圆周半径等于外圆周半径减去籽粒环厚度； 

(5)以横切面的质心为圆心，并参照上述内外圆周半径，从图像中剪切出籽粒环，并将籽粒环转换为灰度图像后再进行归一化处理得g’，求其补，即1-g’记为g； 

(6)把直角坐标系向极坐标系转化并取720个等间隔的旋转角，将籽粒环灰度图像变换为矩形图像，记为p’； 

(7)对p’纵向求和后，计算和值向量的最小值，并以最小值所在竖直轴为界，将左侧整体移至右侧末尾，得p； 

(8)根据预测的玉米果穗穗行数，即集合{8、10、12、14、16、18、20、22}或其子集，分别将横切面的圆心角进行等分，计算等分线上果穗横断面补色图像的均值，然后对每一预设穗行数计算平均值，然后等间隔计算图像p的纵向和，然后横向平移一个步长再重新计算，重新计算的次数以玉米籽粒宽度为限，步长以单个像素为最小值； 

(9)对于每一预设穗行数，取其最大均值为代表值，将最大代表值对应的穗行数作为初估值，如其不是其它预设穗行数的公约数，则将其作为穗行数输出。 

在上述第(5)步中，对于黄色玉米和白色玉米，转换函数为取红、绿和蓝分量的最大值。 

在上述第(7)步中，需对和值向量进行平滑处理。 

一种玉米穗行倾角的测量方法，该方法包括如下步骤： 

(1)利用权利要求1所述的测量装置测得穗长方向上的玉米果穗图像后将其转化为二值图像BW； 

(2)根据上述二值图像BW，逐点计算果穗粗度，记为W，取W的最大值Wmax，并将整数区间[1，Wmax]划分为6等分或者7等分，计数W落入各亚区间的频数，最大频数对应的区间上限值即为玉米果穗粗度的众数； 

(3)穗粗阈值等于k乘以粗度的众数，k∈[0.5，1.0]，分别自左或者自右开始，遍历W，第一个大于或者等于穗粗阈值的位置，即是玉米果穗中段起始或者终止位置，分别记为L和R； 

(4)根据上述起始和终止位置在图像上剪切出玉米果穗中段，记为M(x，y)，x＝1，2，...，C，C为穗段长度，即C＝R-L+1，将M转换为灰度图像后再对其进行归一化，得g’，求其补，即1-g’，记为g； 

(5)求图像g的Radon变换，结果记为G，旋转角的取值范围为[30°，150°]； 

(6)计算G的最大值点，最大值点的旋转角减去90°，即是穗行倾角，其值域为[-60°，60°]。 

一种玉米穗缘角的测量方法，该方法包括如下步骤： 

(1)利用权利要求1所述的测量装置测得穗长方向上的玉米果穗图像后将其转化为二值图像BW； 

(2)根据上述二值图像BW，逐点计算果穗粗度，记为W，取W的最大值Wmax，并将整数区间[1，Wmax]划分为6等分或者7等分，计数W落入各亚区间的频数，最大频数对应的区间上限值即为玉米果穗粗度的众数； 

(3)穗粗阈值等于k乘以粗度的众数，k∈[0.5，1.0]，分别自左或者自右开始，遍历W，第一个大于或者等于穗粗阈值的位置，即是玉米果穗中段起始或者终止位置，分别记为L和R； 

(4)根据上述起始和终止位置在图像上剪切出玉米果穗中段，记为M(x，y)，x＝1，2，...，C，C为穗段长度，即C＝R-L+1； 

(5)以M左下角为原点，沿X轴逐点计算玉米果穗粗度，即在X坐标轴上的投影和，连同X坐标一起，记为数据对(X，Y)，根据数据对(X，Y)建立Y对X的直线回归方程： y＝a+b*x； 

(6)回归系数b就是直线与X轴夹角δ的正切函数值，所以，δ＝tan^-1(b)，而穗缘角为δ值的一半。 

本发明的有益效果是： 

本发明所公开的玉米果穗性状测量装置，通过灯箱与扫描仪组合获取玉米果穗图像，通过计算机对获取的图像进行性状测量，实现了对玉米果穗性状图像信息的计算机图像处理。果穗的性状信息提取更客观、科学，应用于农业生产与科研领域，提高了生产效率。 

本发明所公开的玉米穗行数测量方法从果穗横切面图像剪切出籽粒环，这增强了方法的鲁棒性，能够抵抗各种穗轴芯异常和籽粒间隙异常的干扰。 

本发明所公开的玉米穗行倾角的测量方法根据玉米果穗粗度的众数精确确定剪切位置，而不是简单地将穗长三等分。这有效避免了玉米果穗秃顶和基部缺粒的干扰，增强了方法的鲁棒性。旋转角的取值范围为[30°，150°]，这是根据玉米的生物学特性作出的选择。与取值0度至180度相比，可以大大提高速度，而且避免了个别果穗存在籽粒间隙几乎全部位于同一径向直线上时，误报为籽粒行间隙。 

本发明所公开的玉米穗缘角的测量方法根据玉米果穗粗度的众数精确确定剪切位置，而不是简单地将穗长三等分。这有效避免了玉米果穗秃顶和基部缺粒的干扰，增强了方法的鲁棒性。采用稳健直线回归方法，进一步增加方法的鲁棒性。 

附图说明

图1是本发明的整体结构示意图。 

具体实施方式

为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。 

实施例1，请参照图1，图1示意性的示出了本发明所公开的一种玉米果穗性状测量装置，该测量装置包括用来采集玉米性状图像信息的采集单元，与上述的采集单元相连接的，用于处理上述图像信息并获得最终测量结果的处理单元。前述的采集单元包括，一面可打开的封闭灯箱2，该封闭灯箱2相当于一个暗箱，其底面设有转盘3，顶面设有钩子4，并且在顶面和一侧面开有摄像孔5，在底面的四角装有白炽灯光源1；通过数据线与处理单元连接的相机7，该相机7的镜头可以伸入摄像孔I 5或摄像孔II 6内；通过数据线与处理单元连接的 扫描仪10，该扫描仪10的盖板11敞开，并有遮光罩12罩在扫描区域上。在拍摄玉米侧面图像时，白炽灯光源1打开，可将玉米果穗平放在转盘3上，相机7通过顶面的摄像孔I 5拍摄图像，也可以将玉米果穗的一端挂在钩子4上，玉米果穗自然下垂，相机7通过侧面的摄像孔II6拍摄图像，为了能更全面的拍摄玉米果穗侧面的图像，可将转盘3或钩子4的角度进行旋转，每旋转120度角，拍摄一副图像。在用扫描仪10扫描玉米横断面的图像时，将玉米从中间位置掰断，将断面朝下平放在扫描仪10上，扫描时把扫描仪盖板11敞开或者卸掉，同时，盖上遮光罩12，扫描的图像背景为黑色，以消除背景影响。最后将拍摄的图像和扫描的图像传输给处理单元以备测量特征使用，在本实施例中，处理单元采用普通的计算机8，计算机8与显示器9相连接。 

本发明还公开了一种玉米穗行数测量方法，该方法包括如下步骤： 

(1)将玉米从中部掰断，利用权利要求1所述的测量装置测得玉米果穗横切面的图像； 

(2)将上述图像转化为二值图像BW并作Radon变换，取过圆点的各直径长度，长度的中位数为外圆周半径； 

(3)玉米籽粒环厚度T等于籽粒深度D乘以k，k∈[0.5，1.0]，值越大，籽粒环越厚，籽粒深度D为用户估计值； 

(4)内圆周半径等于外圆周半径减去籽粒环厚度； 

(5)以横切面的质心为圆心，并参照上述内外圆周半径，从图像中剪切出籽粒环，并将籽粒环转换为灰度图像后再进行归一化处理得g’，求其补，即1-g’记为g； 

(6)把直角坐标系向极坐标系转化并取720个等间隔的旋转角，将籽粒环灰度图像变换为矩形图像，记为p’； 

(7)对p’纵向求和后，计算和值向量的最小值，并以最小值所在竖直轴为界，将左侧整体移至右侧末尾，得p； 

(8)根据预测的玉米果穗穗行数，即集合{8、10、12、14、16、18、20、22}或其子集，分别将横切面的圆心角进行等分，计算等分线上果穗横断面补色图像的均值，然后对每一预设穗行数计算平均值，然后等间隔计算图像p的纵向和，然后横向平移一个步长再重新计算，重新计算的次数以玉米籽粒宽度为限，步长以单个像素为最小值； 

(9)对于每一预设穗行数，取其最大均值为代表值，将最大代表值对应的穗行数作为初估值，如其不是其它预设穗行数的公约数，则将其作为穗行数输出。 

在上述第(5)步中，对于黄色玉米和白色玉米，转换函数为取红、绿和蓝分量的最大值。对于其它颜色的玉米，寻找适当的变换函数，使籽粒与籽粒间隙的对比最大。 

在上述第(7)步中，需对和值向量进行平滑处理。 

本发明还公开了一种玉米穗行倾角的测量方法，该方法包括如下步骤： 

(1)利用权利要求1所述的测量装置测得穗长方向上的玉米果穗图像后将其转化为二值图像BW； 

(2)根据上述二值图像BW，逐点计算果穗粗度，记为W，取W的最大值Wmax，并将整数区间[1，Wmax]划分为6等分或者7等分，计数W落入各亚区间的频数，最大频数对应的区间上限值即为玉米果穗粗度的众数； 

(3)穗粗阈值等于k乘以粗度的众数，k∈[0.5，1.0]，值越大，起止距离越短，分别自左或者自右开始，遍历W，第一个大于或者等于穗粗阈值的位置，即是玉米果穗中段起始或者终止位置，分别记为L和R； 

(4)根据上述起始和终止位置在图像上剪切出玉米果穗中段，记为M(x，y)，x＝1，2，...，C，C为穗段长度，即C＝R-L+1，将M转换为灰度图像后再对其进行归一化，得g’，求其补，即1-g’，记为g； 

(5)求图像g的Radon变换，结果记为G，旋转角的取值范围为[30°，150°]，实际应用时，可以根据实际情况，进一步缩小旋转角的取值范围； 

(6)计算G的最大值点，最大值点的旋转角减去90°，即是穗行倾角，其值域为[-60°，60°]。 

在上述第(4)步中，对于黄色玉米和白色玉米，转换函数为取红、绿和蓝分量的最大值。对于其它颜色的玉米，寻找适当的变换函数，使籽粒与籽粒间隙的对比最大。 

一种玉米穗缘角的测量方法，该方法包括如下步骤： 

(1)利用权利要求1所述的测量装置测得穗长方向上的玉米果穗图像后将其转化为二值图像BW； 

(2)根据上述二值图像BW，逐点计算果穗粗度，记为W，取W的最大值Wmax，并将整数区间[1，Wmax]划分为6等分或者7等分，计数W落入各亚区间的频数，最大频数对应的区间上限值即为玉米果穗粗度的众数； 

(3)穗粗阈值等于k乘以粗度的众数，k∈[0.5，1.0]，值越大，起止距离越短，分别自左或者自右开始，遍历W，第一个大于或者等于穗粗阈值的位置，即是玉米果穗中段起始或者终止位置，分别记为L和R； 

(4)根据上述起始和终止位置在图像上剪切出玉米果穗中段，记为M(x，y)，x＝1，2，...，C， C为穗段长度，即C＝R-L+1； 

(5)以M左下角为原点，沿X轴逐点计算玉米果穗粗度，即在X坐标轴上的投影和，连同X坐标一起，记为数据对(X，Y)，根据数据对(X，Y)建立Y对X的直线回归方程：y＝a+b*x； 

(6)回归系数b就是直线与X轴夹角δ的正切函数值，所以，δ＝tan^-1(b)，而穗缘角为δ值的一半。