基于波段比方法的高光谱反射图像中玉米胚部的分割方法技术

本发明专利技术公开了一种基于波段比方法的高光谱反射图像中玉米胚部的分割方法。胚面朝上，采集玉米种子的高光谱图像；提取胚乳部和胚部的多个像素点及其光谱，归一化处理；并计算其任意两波段的光谱比值作为波段比，选取最优波段比；选取波段子集，计算所有像素在波段子集的方差，进行阈值分割获得种子二值图像，掩模运算获得去背景后的玉米种子高光谱图像；由去背景后的玉米种子高光谱图像计算其最优波段比下的灰度图像，对于灰度图像采用区域增长法分割得到胚部二值图像，再进行图像处理获得种子胚部二值图像，再次掩模运算得到玉米种子胚部的高光谱图像。本发明专利技术方法实现对可见近红外高光谱图像中玉米种子胚部的有效分割，分割效果好，实用性较强。

【技术实现步骤摘要】

本专利技术涉及了一种基于波段比方法的高光谱反射图像中玉米胚部的分割方法，可应用于各种黄色玉米种子胚部质量检测，尤其是对胚部化学特性分析等领域的应用。
技术介绍
种子的胚部是种子中最重要的部分，可以发育成植物的根、茎和叶。如果种子的胚部受损对植物的后续发育影响最大。玉米种子的胚部大，体积占籽粒的1/3左右，重量占籽粒的10％-20％。玉米胚中含有较多亲水基，比胚乳更易吸湿。同时玉米胚部集中了玉米籽粒中84％的脂肪、83％的无机盐、65％的糖和22％的蛋白质，营养丰富，适宜霉菌的生长和繁殖，容易霉变。胚部脂肪含量高，吸湿性强，容易酸败，导致玉米种子的生活力降低。研究玉米种子的胚部特性，对玉米种业有重要意义。高光谱图像融合了传统图像技术和光谱技术的优势，可以同时提取出对象的图像和光谱信息，能够反映对象的颜色、形态和化学成分差异。光谱技术用于谷物化学成分检测的研究较多，但并不能研究单个种子的化学成分空间差异等特性。高光谱图像技术可以提取图像中每个像素点的光谱进行分析，可以有效的分析种子各部位的化学成分指标。
技术实现思路
本专利技术的目的在于，针对单独分析玉米种子胚部特性的需求，提出一种基于波段比方法的高光谱反射图像中玉米胚部的分割方法，克服了在传统高光谱图像分析中需要手工选择种子胚部区域后在再人工提取胚部光谱的问题，实现玉米种子胚部在高光谱图像中的准确、快速分割。本专利技术的技术方案如图1所示，具体步骤如下：1)将玉米种子胚面朝上，通过高光谱反射图像采集系统采集N粒玉米种子在x个可见近红外波段的高光谱图像；x一般具体实施可取400、1000或者更大的数值，使得每段可见近红外波段的波段范围小于2nm。2)如图2所示，从高光谱反射图像中选取胚乳部和胚部的多个像素点(至少200个像素点)，提取各像素点的光谱，对所有波段的高光谱反射图像进行归一化处理；3)对于步骤2)中提取得到的各像素点，计算其任意两波段的光谱比值作
为波段比，采用KW(Kruskal and Wallis)检验法选取最优波段比；4)如图2所示，选取500-580nm波段子集，选取涵盖500-580nm波段的所有可见近红外波段作为波段子集，计算高光谱图像所有像素在波段子集的方差，并构建波段方差图并进行阈值分割，获得种子二值图像；5)图像背景分割：利用阈值分割后得到的二值图像，对玉米种子的各可见近红外波段的高光谱图像进行掩模运算实现背景分割，获得去背景后的玉米种子高光谱图像；6)胚部分割及掩模：如图3所示，由去背景后的玉米种子高光谱图像计算其最优波段比下的灰度图像，对于灰度图像采用区域增长法分割得到胚部二值图像，再对胚部二值图像和步骤4)获得的种子二值图像进行图像处理操作获得种子胚部二值图像，用种子胚部二值图像对步骤5)获得的去背景后的玉米种子高光谱图像进行掩模运算，得到玉米种子胚部的高光谱图像，实现种子高光谱图像中的胚部分割。所述的步骤2)中，归一化计算方法采用下式方法计算： I λ i ′ = I λ i m a x ( I λ 0 , ... , I λ X ) ]]>其中，λi为第i个波段的波长，I′λi表示第i个波段下归一化后的图像像素反射值，Iλi表示第i个波段下采集到的图像像素反射值。所述步骤3)采用KW(Kruskal and Wallis)检验法选取最优波段比具体采用以下方式：3.1)提取的所有像素点分为位于胚乳部和位于胚部的两类，以像素点的任意两波段的波段比为特征，根据波段比对像素点进行编号，形成特征集合为Xj；3.2)对于胚乳部和胚部，分别计算每类像素点编号的平均值，记为3.3)采用以下公式计算第j个像素点的波段比的统计量KWj，选取统计量KWj最大的波段比作为最优波段比： KW j = 12 N ( N + 1 ) Σ k = 1 M N k ( R k ‾ - N + 1 2 ) 2 ]]>其中，k表示分类的序号，Nk表示k类像素点数量，N表示提取获得的胚乳部和位于胚部两类所有像素点的总数，M表示k类像素点的总数；统计量KWj越大，说明该特征下的类间离差越大，该特征的分类能力越强。所述根据波段比对像素点进行编号具体为：对于胚乳部和胚部所有提取得到的像素点，按波段比从小到大排序编号；如果多个像素点的波段比相同，则这几个像素点按随机顺序排序编号；如果所有像素点的波段比全部相同，则该
波段比下的KWj值为0。如果存在多个统计量KWj最大的波段比，则分别计算胚乳部的像素点光谱和胚部的像素点光谱在统计量KWj最大的波段比下的平均值，取二者平均值差异最大值所在的波段比作为最优波段比。所述步骤6)中对灰度图像采用区域增长法进行区域增长分割时，种子点取为0，差分相似度阈值取为1。所述步骤6)中种子胚部二值图像采用以下方式图像处理获得：先将胚部二值图像和种子二值图像相乘，然后依次进行填充空洞、腐蚀、删除白色连通区域和膨胀得到玉米种子胚部二值图。删除白色连通区域是指删除含有小于100个白色像素的连通区域。所述的光谱反射图像采集系统在采集前经过黑白校正，具体采用以下方式进行校正：首先将标准聚四氟乙烯白板(Spectralon,Labsphere Inc.，USA)放于成像区域，调整曝光时间，采集本标准白板的全白反射标定图像；然后盖上镜头盖，关闭光源，获取暗场全黑标定图像；接着采用本文档来自技高网...

【技术保护点】
一种基于波段比方法的高光谱反射图像中玉米胚部的分割方法，其步骤包括：1)将玉米种子胚面朝上，通过高光谱反射图像采集系统采集N粒玉米种子在x个可见近红外波段的高光谱图像；2)从高光谱反射图像中选取胚乳部和胚部的多个像素点，提取各像素点的光谱，对所有波段的高光谱反射图像进行归一化处理；3)对于步骤2)中提取得到的各像素点，计算其任意两波段的光谱比值作为波段比，采用KW(Kruskal and Wallis)检验法选取最优波段比；4)选取涵盖500‑580nm波段的所有可见近红外波段作为波段子集，计算高光谱图像所有像素在波段子集的方差，构建波段方差图并进行阈值分割，获得种子二值图像；5)图像背景分割：利用阈值分割后得到的二值图像，对玉米种子的各可见近红外波段的高光谱图像进行掩模运算实现背景分割，获得去背景后的玉米种子高光谱图像；6)胚部分割及掩模：由去背景后的玉米种子高光谱图像计算其最优波段比下的灰度图像，对于灰度图像采用区域增长法分割得到胚部二值图像，再对胚部二值图像和步骤4)获得的种子二值图像进行图像处理操作获得种子胚部二值图像，用种子胚部二值图像对步骤5)获得的去背景后的玉米种子高光谱图像进行掩模运算，得到玉米种子胚部的高光谱图像，实现种子高光谱图像中的胚部分割。...

【技术特征摘要】
1.一种基于波段比方法的高光谱反射图像中玉米胚部的分割方法，其步骤包括：1)将玉米种子胚面朝上，通过高光谱反射图像采集系统采集N粒玉米种子在x个可见近红外波段的高光谱图像；2)从高光谱反射图像中选取胚乳部和胚部的多个像素点，提取各像素点的光谱，对所有波段的高光谱反射图像进行归一化处理；3)对于步骤2)中提取得到的各像素点，计算其任意两波段的光谱比值作为波段比，采用KW(Kruskal and Wallis)检验法选取最优波段比；4)选取涵盖500-580nm波段的所有可见近红外波段作为波段子集，计算高光谱图像所有像素在波段子集的方差，构建波段方差图并进行阈值分割，获得种子二值图像；5)图像背景分割：利用阈值分割后得到的二值图像，对玉米种子的各可见近红外波段的高光谱图像进行掩模运算实现背景分割，获得去背景后的玉米种子高光谱图像；6)胚部分割及掩模：由去背景后的玉米种子高光谱图像计算其最优波段比下的灰度图像，对于灰度图像采用区域增长法分割得到胚部二值图像，再对胚部二值图像和步骤4)获得的种子二值图像进行图像处理操作获得种子胚部二值图像，用种子胚部二值图像对步骤5)获得的去背景后的玉米种子高光谱图像进行掩模运算，得到玉米种子胚部的高光谱图像，实现种子高光谱图像中的胚部分割。2.根据权利要求1所述的一种基于波段比方法的高光谱反射图像中玉米胚部的分割方法，其特征在于：所述的步骤2)中，归一化计算方法采用下式方法计算： I λ i ′ = I λ i m a x ( I λ 0 , ... , I λ X ) ]]>其中，λi为第i个波段的波长，I′λi表示第i个波段下归一化后的图像像素反射值，Iλi表示第i个波段下采集到的图像像素反射值。3.根据权利要求1所述的一种基于波段比方法的高光谱反射图像中玉米胚部的分割方法，其特征在于：所述步骤3)采用KW(Kruskal and Wallis)检验法选取最优波段比具体采用以下方式：3.1)提取的所有像素点分为位于胚乳部和位于胚部的两类，以像素点的任
\t意两波段的波段比为特征，根据波段比对像素点进行编号，形成特征集合为Xj；3.2)对于胚乳部和胚部，分别计算每类像素点编号的平均值，记为3.3)采用以下公式计算第j个像素点的波段比的统计量KWj，选取统计量KWj最大的波段比作为最优波段比： KW j = 12 N ( N ...

【专利技术属性】
技术研发人员：成芳，杨小玲，
申请(专利权)人：浙江大学，
类型：发明
国别省市：浙江;33