首页 / 专利库 / 视听技术与设备 / 图像分割 / 分水岭分割 / 一种基于OCT图像处理的枇杷早期瘀伤鉴别方法

一种基于OCT图像处理的枇杷早期瘀伤鉴别方法

阅读:926发布:2020-05-14

专利汇可以提供一种基于OCT图像处理的枇杷早期瘀伤鉴别方法专利检索,专利查询,专利分析的服务。并且本 发明 公开了一种基于OCT 图像处理 的枇杷早期瘀伤 鉴别 方法。采集枇杷果实的具有细胞图像细节的SD‑OCT图像,使用双三次插值 算法 ,对图像进行降 采样 缩小图像 分辨率 ,进行高斯模糊降噪处理,提取枇杷果实目标和背景的分界线,取分界线的最高点并作为参考点,根据参考点使得分界线 变形 为沿参考点的一条直线,对图像进行均值滤波,对图像进行二值化处理,将二值图像处理获得各个细胞对应的细胞区域,通过对细胞区域进行分析计算获得瘀伤鉴别的结果。本发明方法实现了枇杷果肉的早期瘀伤的全自动化检测,并完成瘀伤组织的皮下细胞标识和判别,提高了检测效率,为枇杷内在品质在线检测奠定技术 基础 。,下面是一种基于OCT图像处理的枇杷早期瘀伤鉴别方法专利的具体信息内容。

1.一种基于OCT图像处理的枇杷早期瘀伤鉴别方法,其特征在于:包括以下步骤:
1)采集枇杷果实的具有细胞图像细节的SD-OCT图像;
2)使用双三次插值算法,对图像进行降采样,缩小图像分辨率
3)对步骤2)获得的SD-OCT图像进行高斯模糊降噪处理;
4)提取枇杷果实目标和背景的分界线;
5)取分界线的最高点并作为参考点,计算分界线与参考点的纵向坐标差并作为位移量,对分界线上除参考点所在列外的每一列进行上下位移,使得分界线变形为沿参考点的一条直线;
6)取3×3模板,对图像进行均值滤波;
7)设置阈值,对图像进行二值化处理,获得二值图像;
8)将二值图像处理获得各个细胞对应的细胞区域,通过对细胞区域进行分析计算获得瘀伤鉴别的结果;所述步骤8)具体为:
8.1)针对二值图像的每个像素,计算像素的最短距离:若所在像素为零像素,最短距离为所在像素与最近的非零像素之间的距离;若所在像素为非零像素,最短距离为零;
8.2)使用分岭算法,使用步骤8.1)所述的最短距离,对图像根据细胞不同进行分割,得到分割后的各个细胞区域;
8.3)在分割后的各个细胞区域中进行筛选,去除表皮细胞的细胞区域,保留果肉细胞的细胞区域;
8.4)计算每个细胞区域的费雷特直径和等效直径,保留费雷特直径满足费雷特直径下限阈值≤R1≤费雷特直径上限阈值,且最大等效直径小于等效直径阈值的细胞区域;
8.5)由步骤8.4)获得的所有细胞区域计算总区域表面积、平均区域面积、平均费雷特直径、平均等效圆直径和单位面积细胞数:
总区域表面积:定义为所有细胞区域的面积之和;
平均区域面积=总细胞区域表面积/细胞区域数;
平均费雷特直径=所有细胞区域的费雷特直径之和/细胞区域数;
平均等效圆直径=所有细胞区域的等效圆直径之和/细胞区域数;
单位面积细胞数=细胞区域数/OCT图像占有的面积;
8.6)设正常和瘀伤枇杷的标准样本集合,分别计算标准样本的总区域表面积、平均区域面积、平均费雷特直径、平均等效圆直径、单位面积细胞数阈值,通过聚类分析判断获得瘀伤鉴别的结果。
2.根据权利要求1所述的一种基于OCT图像处理的枇杷早期瘀伤鉴别方法,其特征在于:所述的SD-OCT图像的图像清晰度达到肉眼能够清晰分辨枇杷表皮和果肉的细胞。
3.根据权利要求1所述的一种基于OCT图像处理的枇杷早期瘀伤鉴别方法,其特征在于:所述步骤4)具体为:
4.1)设置滤波模板为[-1,1],对OCT图像进行第一次滤波;
4.2)设置滤波模板为[1,-1],对OCT图像进行第二次滤波;
4.3)对滤波后的图像进行归一化处理;
4.4)对图像进行二值化变换:设立限,设置图像大于等于门限的像素为1,图像小于门限的像素为0;
4.5)对图像进行闭运算操作;
4.6)对图像进行开运算操作;
4.7)对二值化后图像中的每一列像素点,自上向下搜索该列出现第一个灰度值为1的像素点并记录作为目标和背景的分界线;
4.根据权利要求1所述的一种基于OCT图像处理的枇杷早期瘀伤鉴别方法,其特征在于:所述步骤7)具体为:设置分割限,设置图像中大于等于0且小于等于分割限的像素为1,图像中大于分割限的像素为0。

说明书全文

一种基于OCT图像处理的枇杷早期瘀伤鉴别方法

技术领域

[0001] 本发明属于果内部品质自动化检测领域,涉及OCT图像处理方法,尤其是涉及了一种基于OCT图像处理的枇杷早期瘀伤鉴别方法。

背景技术

[0002] 枇杷是我国的特有水果特产之一,其在采摘、销售、运输、存储过程中内部品质的无损快速检测方法是枇杷产业发展面临的主要技术问题。枇杷在各产业销售过程中,极易收到外损伤,造成后期的腐败变质。枇杷的瘀伤可能出现在采摘、存储、运输、包装等各个环节,在销售早期不易被察觉。瘀伤后的枇杷货架期大大缩短,由于细胞结构的破损,组织逐步褐变,严重影响了消费者的满意度和回购率。
[0003] 在无损检测枇杷内部结构上,一般使用光谱法或者高光谱成像,需要大型设备保证光谱信息的全面采集,耗费比较多的检测时间和成本,且对检测人员有一定技术要求,且高光谱图像很难真正反映其内部情况,光谱特征随着枇杷种类的变化有一定的偏移。谱域光学相干层析成像(SD-OCT)通过测量物质的光学干涉特性展现内部其结构形态和分布,目前SD-OCT图像已经用于人体多个组织的识别、定量测量、定性鉴定,报道表明图像可以清楚展现生物组织的层次结构。目前OCT图像法在农业、养殖领域主要应用有:观察苹果的表皮结构、区别海水有核珍珠与淡水无核珍珠、观察种子的内部细胞结构、观察植物叶片的生长缺陷等。该方法用于枇杷产业,具有较广阔的应用前景。
[0004] 由于在产业应用中,枇杷OCT图像反差小,特征不明显,采用人工判别,基本无法判别早期瘀伤的情况。在枇杷OCT图像应用过程中,也没有成系统地报道图像处理算法,各项研究尚处于起步阶段,现有技术缺少能进行枇杷早期瘀伤鉴别的方法。

发明内容

[0005] 针对于背景技术中存在的问题,本发明的目的在于提供了一种基于OCT图像处理的枇杷早期瘀伤鉴别方法,能够自动识别OCT图像中枇杷的瘀伤缺陷,并对组织细胞的形态学参数作出了评价,提高了检测效率,配合成像等外观检测方法,为枇杷在线检测奠定技术基础
[0006] 本发明采用的技术方案是包括以下步骤:
[0007] 1)采集枇杷果实的具有细胞细节的SD-OCT图像,所述的SD-OCT图像的图像清晰度达到肉眼能够清晰分辨枇杷表皮和果肉的细胞;
[0008] 2)使用双三次插值算法,对图像进行降采样,缩小图像分辨率
[0009] 3)对步骤2)获得的SD-OCT图像进行高斯模糊降噪处理;
[0010] 4)提取枇杷果实目标和背景的分界线;
[0011] 5)取分界线的最高点并作为参考点,计算分界线与参考点的纵向坐标差并作为位移量,对分界线上除参考点所在列外的每一列进行上下位移,使得分界线变形为沿参考点的一条直线;对于移出图像区域的像素,删除,对于移入图像的新区域,直接补零;
[0012] 6)取3×3模板,对图像进行均值滤波;
[0013] 7)设置阈值,对图像进行二值化处理,获得二值图像,二值图像中的像素为零像素或者非零像素;
[0014] 8)将二值图像处理获得各个细胞对应的细胞区域,通过对细胞区域进行分析计算获得瘀伤鉴别的结果。
[0015] 所述步骤8)具体为:
[0016] 8.1)针对二值图像的每个像素,计算像素的最短距离:若所在像素为零像素,最短距离为所在像素与最近的非零像素之间的距离;若所在像素为非零像素,最短距离为零;
[0017] 8.2)使用分水岭算法,使用步骤8.1)所述的最短距离,对图像根据细胞不同进行分割,得到分割后的各个细胞区域;
[0018] 8.3)在分割后的各个细胞区域中进行筛选,去除表皮细胞的细胞区域,保留果肉细胞的细胞区域;具体是以参考点为基准所在行向下间隔一段距离开始选取余下图像区域进行保留。
[0019] 8.4)计算每个细胞区域的费雷特直径和等效直径,保留费雷特直径满足费雷特直径下限阈值≤R1≤费雷特直径上限阈值,且最大等效直径小于等效直径阈值的细胞区域;
[0020] 8.5)由步骤8.4)获得的所有细胞区域计算总区域表面积、平均区域面积、平均费雷特直径、平均等效圆直径和单位面积细胞数:总区域表面积:定义为所有细胞区域的面积之和;平均区域面积=总细胞区域表面积/细胞区域数;平均费雷特直径=所有细胞区域的费雷特直径之和/细胞区域数;平均等效圆直径=所有细胞区域的等效圆直径之和/细胞区域数;单位面积细胞数=细胞区域数/OCT图像占有的面积;
[0021] 8.6)设正常和瘀伤枇杷的标准样本集合,分别计算标准样本的总区域表面积、平均区域面积、平均费雷特直径、平均等效圆直径、单位面积细胞数阈值,通过聚类分析判断获得瘀伤鉴别的结果。
[0022] 所述步骤4)具体为:
[0023] 4.1)设置滤波模板为[-1,1],对OCT图像进行第一次滤波;
[0024] 4.2)设置滤波模板为[1,-1],对OCT图像进行第二次滤波;
[0025] 4.3)对滤波后的图像进行归一化处理;
[0026] 4.4)对图像进行二值化变换:设立限,设置图像大于等于门限的像素为1,图像小于门限的像素为0;
[0027] 4.5)对图像进行闭运算操作;
[0028] 4.6)对图像进行开运算操作;
[0029] 4.7)对二值化后图像中的每一列像素点,自上向下搜索该列出现第一个灰度值为1的像素点并记录作为目标和背景的分界线。
[0030] 所述步骤7)具体为:设置分割限,设置图像中大于等于0且小于等于分割限的像素为1,图像中大于分割限的像素为0。
[0031] 本发明采用了多次模板匹配滤波方式进行图像处理,其目的就是从OCT图像灰度的变化中提取枇杷大细胞的图像,并在后期计算细胞的多种截面特征,从特征参数判断瘀伤的存在。之前的专利,图像前期处的主要目的是为了提取计算局部组织的光学特性参数,进而用光学参数来判断瘀伤,其检测的技术构思相比本申请以前的现有方法有明显差异。
[0032] 本发明具有的有益效果是:
[0033] 本发明使用OCT图像检测枇杷果实的内部瘀伤缺陷,具有无损、快速、低成本的优点,大大提高了瘀伤判别的效率和准确性。
[0034] 本发明方法采用了细胞结构参数作为评价手段,对不同形状、不同大小、不同厚度、产地、生长环境的瘀伤组织具有普适性,并能自动标记缺陷的细胞位置,较其他方法具有更好定位精度
[0035] 本发明采用细胞颗粒参数作为评价手段,结合拉平变换后的图像,检测效果具有一定的鲁棒性。附图说明
[0036] 图1是本发明方法的流程图
[0037] 图2是典型枇杷样本的OCT图像,其中(a)为正常无缺陷样本,(b)样本果肉组织存在瘀伤缺陷。
[0038] 图3是本发明实施例图象处理过程的原始图像。
[0039] 图4是本发明实施例拉平后获得直线分界线的图像。
[0040] 图5是本发明实施例分割后所有细胞区域的图像。
[0041] 图6是本发明实施例筛选后保留细胞区域的图像。

具体实施方式

[0042] 以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
[0043] 本发明的实施例及其实施过程如下:
[0044] 1)使用Thorlabs公司生产的TELSTO 1300V2型SD-OCT成像仪采集枇杷果实的SD-OCT图像样本40个,其中20个含有不同程度的瘀伤缺陷,20个为正常样本;图2为其中2个典型枇杷样本的OCT图像,其中(a)为正常无缺陷样本,(b)样本果肉组织存在瘀伤缺陷。图中可见正常无缺陷组织结构密度较高且紧凑,而存在瘀伤的组织出现了密度较小的稀疏组织。仅使用肉眼,无法判断其瘀伤情况。
[0045] 2)输入如图4所示的原始图像,使用双三次插值算法,对图像进行降采样,将原始图像分辨率(1625*1024)缩小至其1/3。降低采样的目的是便于进行图像的快速处理。
[0046] 3)接着对SD-OCT图像进行高斯模糊降噪处理,去除了散弹噪声。
[0047] 4)提取枇杷果实目标和背景的分界线;
[0048] 4.1)设置模板为[-1,1],对OCT图像进行第一次滤波;
[0049] 4.2)设置模板为[1,-1],对OCT图像进行第二次滤波;
[0050] 二次滤波操作,提取了拉平后图像的垂直突变部分。
[0051] 4.3)对滤波后的图像进行最大最小值归一化处理;
[0052] 4.4)对图像进行二值化变换:设立门限为1.2,设置图像大于等于门限的像素为1,图像小于门限的像素为0;
[0053] 4.5)取3*3方型区域,对图像进行闭运算操作;
[0054] 4.6)取3*3方型区域,对图像进行开运算操作。
[0055] 对于区域进行闭开运算,区分了背景和目标。
[0056] 4.7)对二值化后图像中的每一列像素点,自上向下搜索该列出现第一个灰度值为1的像素点并记录作为目标和背景的分界线。
[0057] 5)取分界线的最高点,并以此作为参考点;计算分界线与参考点的纵向坐标差,并以此为位移量,对分界线上除参考点所在列外的每一列进行上下位移,使得分界线变成经过参考点的一条直线;对于移出图像区域的像素,删除,对于移入图像的新区域,直接补零;拉平变换后如图4所示。
[0058] 6)取3*3模板,对图像进行均值滤波,加大了图像的平滑性。
[0059] 7)设置分割限为80,对图像进行二值化处理,设置图像中大于等于0且小于等于分割限的像素为1,图像中大于分割限的像素为0。
[0060] 8)针对二值图像的每个像素,计算像素的最短距离。对于本是是0的像素,其最短距离定义为与其距离最近的1像素的距离。例如通过边相邻的零像素与最近的非零像素之间的距离为1,通过相邻的零像素与最近的非零像素之间的距离为√2。对于本身是1的像素,其最短距离为0。
[0061] 9)使用分水岭算法,对图像根据细胞不同进行分割,得到分割后的各个细胞区域,如图5所示;
[0062] 10)在分割后的各个细胞区域中进行筛选,去除表皮细胞的细胞区域,保留果肉细胞的细胞区域;本实施例中,以参考点所在的直线分界线为基准,向下截取距离0.07mm到1mm范围内的区域图像,保留的图像如图6所示。
[0063] 11)计算步骤10)所述区域图像中每个区域的费雷特直径和等效直径,选取费雷特直径在30μm到100μm的区域且最大等效直径小于150μm的区域。
[0064] 12)由获得的所有细胞区域计算总区域表面积、平均区域面积、平均费雷特直径、平均等效圆直径、单位面积细胞数。
[0065] 13)设正常和瘀伤枇杷的标准样本集合,分别计算标准样本的总区域表面积、平均区域面积、平均费雷特直径、平均等效圆直径、单位面积细胞数阈值,通过聚类分析判断判断瘀伤的结果。具体实施采用KNN聚类分析方法进行分析判断。
[0066] 在本实施例中,随机取正常和瘀伤的两类样本各10个,计20个作为标准样本,对余下的20个样本进行了分类识别,表1给出了各类参数的统计值,实验结果表明,对于20个样本中的瘀伤样本及正常样本识别率达到了100%。
[0067] 表1大细胞的组织结构参数(95%置信区间)
[0068] 参数 单位 正常组织 瘀伤组织总区域表面积 Mm2 2.08±1.20 1.45±0.07
平均区域面积 Mm2 0.0042±0.0002 0.0043±0.0001
平均费雷特直径 μm 58.63±1.05 59.15±0.85
平均等效圆直径 μm 40.42±0.73 40.62±0.61
单位面积细胞数 - 491.70±25.54 340.15±12.34
[0069] 本发明方法能够用于实现枇杷果肉的早期瘀伤的全自动化检测,完成瘀伤组织的皮下细胞标识和判别,实施通过对不同产地、品种枇杷的瘀伤组织进行检测获得了较强的检测可靠性,提高了检测效率。
[0070] 在本发明实施例中,本领域普通技术人员还可以理解,实现上述实施例方法中的全部或部分步骤是可以通过程序来指令相关的硬件来完成,所述的程序可以在存储于一计算机可读取存储介质中,所述的存储介质,包括ROM/RAM、磁盘、光盘等。
[0071] 以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈