基于LoG算法的水果热成像检测

韩亚辉, 王琢, 刘佳鑫

韩亚辉, 王琢, 刘佳鑫. 基于LoG算法的水果热成像检测[J]. 红外技术, 2021, 43(7): 709-715.
引用本文: 韩亚辉, 王琢, 刘佳鑫. 基于LoG算法的水果热成像检测[J]. 红外技术, 2021, 43(7): 709-715.
HAN Yahui, WANG Zhuo, LIU Jiaxin. Fruit Thermal Imaging Detection Based on Laplacian of Gaussian Algorithm[J]. Infrared Technology , 2021, 43(7): 709-715.
Citation: HAN Yahui, WANG Zhuo, LIU Jiaxin. Fruit Thermal Imaging Detection Based on Laplacian of Gaussian Algorithm[J]. Infrared Technology , 2021, 43(7): 709-715.

基于LoG算法的水果热成像检测

基金项目: 

中央高校基本科研业务费专项资金资助 2572019CP21

详细信息
    作者简介:

    韩亚辉(1996-),硕士研究生。研究方向:图像处理。E-mail: 2962592004@qq.com

    通讯作者:

    王琢(1979-),硕士,副教授。研究方向:林业工程自动化、检测与测试技术。E-mail: wangzhuo@nefu.edu.cn

  • 中图分类号: TN219

Fruit Thermal Imaging Detection Based on Laplacian of Gaussian Algorithm

  • 摘要: 传统的水果分级与损伤检测大多采用感官评定的方法,随着计算机视觉技术的发展,计算机视觉检测分级技术发展迅速。研究中针对解决水果损伤部位检测的问题提出了一种利用图像处理技术对水果热成像损伤部位进行检测的技术方案。本方案采用Laplacian of Gaussian(LoG)算法对损伤部位进行检测,使用高斯卷积模板抑制噪声,通过设置不同的卷积核尺寸以及σ值获得不同的卷积滤波结果,加强了图像中损伤部位的色彩程度,进而更好地利用边缘检测技术获取损伤部位的边缘信息。采用具有局部损伤的苹果作为研究对象,选取有参考和无参考等5种评价方法,分析卷积过程对于损伤部位边缘检测的影响。结果表明,在水果热成像中LoG算法可以有效地检测水果的损伤部位,卷积核尺寸对于水果损伤部位边缘检测结果的影响远大于σ值,通过增大卷积核尺寸可以有效地加深损坏部分的边缘信息,研究为水果损伤区域检测提供了一种可行的解决方案。
    Abstract: Traditional fruit grading and damage detection mostly use sensory evaluation methods. With the development of computer vision technology, automatic computer vision detection and grading technology developed rapidly. To solve the problem of fruit damage detection, we propose a technical scheme for fruit thermal imaging damage detection using image processing technology. In this scheme, the Laplacian of Gaussian (LoG) algorithm was used to detect the damaged parts; a Gaussian convolution template is used to suppress noise. Different convolution filter results were obtained by varying the convolution kernel sizes and σ values to enhance the color degree of the damaged part in the image. Then, the edge detection technology was used to obtain the edge information of the damaged part. In the experiment, apples with local damage were selected as the research object, and five evaluation methods, including references and non-references, were selected to analyze the influence of the convolution process on the edge detection of damaged parts. The experimental results show that the LoG algorithm can effectively detect the damaged parts of fruits during thermal imaging, and the influence of the convolution kernel size on the edge detection results is far greater than the value of σ. By increasing the size of the convolution kernel, the edge information of the damaged parts can be effectively deepened. This study provides a feasible solution for fruit damage area detection.
  • 近年来,随着红外制冷探测器朝尺寸小、重量轻、功耗低、成本低即低SWaP-C(size, weight and power, cost)方向快速发展,集成小型制冷机的该类中波红外探测器国外已大量应用于武器热瞄镜、便携式手持热像仪、小型无人机、无人车、遥控狙击手和遥控武器站、导弹导引头等空间受限的红外系统[1]。针对该类小型探测器设计一款仅由4片透镜组成的尺寸小、重量轻、成本低的中波红外连续变焦光学系统进而生产结构紧凑、低功耗和低成本的红外变焦热像仪将在手持观瞄具、边防监视系统、小型无人系统等平台得到广泛应用。

    目前采用四片式的红外变焦光学多为双视场变焦系统。文献[2]为四片式非制冷双视场,文献[3]为四片制冷型长波双视场系统,文献[4]为4片制冷型中波双视场系统,目前仅采用4片透镜实现连续变焦的制冷型红外光学系统未见报道。

    从变焦理论分析,一般常用的机械补偿变焦系统由典型的前固定组、变倍组、补偿组、后固定组4组透镜组成。变倍组一般是负透镜,而补偿组可以是正透镜组也可以是负透镜组,前者为正组补偿系统,后者称为负组补偿系统。从光学系统像差校正难易程度、减少透镜数量、降低光学透镜成本考虑,本文采用无后固定组正组补偿变焦系统,即变焦部分由会聚目标光线的前固定组正透镜、变倍负透镜、补偿正透镜构成。为压缩前固定组物镜口径并满足100%冷屏效率,系统采用二次成像方案,利用单片正光焦度透镜将一次像再次中继成像到探测器焦平面。为压缩轴向尺寸,采用二片平面反射镜将光路U型折转,最终实现仅由4片透镜构成的轻小型中波红外连续变焦光学系统。

    机械补偿连续变焦光学系统参数求解就是确定变焦系统在满足像面稳定和焦距在一定范围内连续变化的条件下系统中各组元的焦距、间隔、位移量等参数。通过建立数学模型能方便地计算和分析变焦过程、确定变焦系统高斯光学参数[5]。二组元正组补偿连续变焦系统运动方式如图 1所示。

    图  1  正组补偿光学系统变焦模型
    Figure  1.  Principle diagrams of continuous zoom optical system with mechanical compensation

    变焦系统由于只有运动组才产生像面位移,只需抽出变倍组、补偿组加以分析。

    因变倍组f2′的移动,引起整个运动组分的像面移动为m32(1-m22)dq,因补偿组f3′的移动,引起整个运动组分的像面移动为(1-m32)dΔ。为达到像面稳定,两个运动组像面移动量的代数和必为零。

    $$ m_3^2\left(1-m_2^2\right) \mathrm{d} q+\left(1-m_3^2\right) \mathrm{d} \varDelta=0$$ (1)

    而变倍组f2′、补偿组f3′微分移动量dq、dΔ与其倍率变化dm2、dm3之间的关系为:

    $$ \mathrm{d} q = \frac{{{f_2}'}}{{m_2^2}}\mathrm{d}{m_2} $$ (2)
    $$ \mathrm{d} \varDelta=f_3{ }^{\prime} \mathrm{d} m_3 $$ (3)

    将(2)、(3)代入(1),经整理得到二组元连续变焦微分方程如下:

    $$ \frac{{1 - m_2^2}}{{m_2^2}}{f_2}'\mathrm{d}{m_2} + \frac{{1 - m_3^2}}{{m_3^2}}{f_3}'\mathrm{d}{m_3} = 0 $$ (4)

    式(4)是多变量全微分型微分方程,设U(m2, m3)为原函数,则有dU(m2, m3)=0。

    其通解为:

    $$ U({m_2}, {m_3}) = {f_2}'(\frac{1}{{{m_2}}} + {m_2}) + {f_3}'(\frac{1}{{{m_3}}} + {m_3}) = C $$ (5)

    式中:C为常量;设变倍组f2′、补偿组f3′初始状态都处于系统长焦位置,则:

    $$ m_{2}=m_{2l}\text{;}m_{3}=m_{3l} $$ (6)

    同样有:

    $$ {f_2}'(\frac{1}{{{m_{2l}}}} + {m_{2l}}) + {f_3}'(\frac{1}{{{m_{3l}}}} + {m_{3l}}) = C $$ (7)

    消去常量C,得到方程的特解:

    $$ {f_2}'(\frac{1}{{{m_2}}} - \frac{1}{{{m_{2l}}}} + {m_2} - {m_{2l}}) + {f_3}'(\frac{1}{{{m_3}}} - \frac{1}{{{m_{3l}}}} + {m_3} - {m_{3l}}) = 0 $$ (8)

    将(8)式整理得到补偿组f3′的倍率m3构成的二次方程:

    $$ m_3^2 - b{m_3} + 1 = 0 $$ (9)

    其中

    $$ b = - \frac{{{f_2}'}}{{{f_3}'}}(\frac{1}{{{m_2}}} - \frac{1}{{{m_{2l}}}} + {m_2} - {m_{2l}}) + (\frac{1}{{{m_{3l}}}} + {m_{3l}}) $$ (10)

    解得m3的两根为:

    $$ {m_{31}} = \frac{{b + \sqrt {{b^2} - 4} }}{2} $$ (11)
    $$ {m_{32}} = \frac{{b - \sqrt {{b^2} - 4} }}{2} $$ (12)

    系统参数求解过程如下:

    1)将(2)式积分并整理得到变倍组f2′的倍率m2

    $$ {m_2} = \frac{1}{{\frac{1}{{{m_{2l}}}} + \frac{{{q_2}}}{{{f_2}'}}}} $$ (13)

    2)根据求得的m2按照公式(10)求出系数b,再由(11)、(12)式解得补偿组f3′满足运动方程的两个解m31m32

    3)根据补偿组的两个解求出满足运动方程补偿像面位移的移动量Δ1Δ2

    $$ \varDelta_{1}=f_{3}′(m_{31}-m_{3l})$$ (14)
    $$ \varDelta_{2}=f_{3}′(m_{32}-m_{3l}) $$ (15)

    4)求出系统的总变倍比

    $$ {\varGamma _1} = \frac{{{m_{2l}}{m_{3l}}}}{{{m_2}{m_{31}}}} $$ (16)
    $$ {\varGamma _2} = \frac{{{m_{2l}}{m_{3l}}}}{{{m_2}{m_{32}}}} $$ (17)

    变倍组f2′每移动q2对应着Γ1Γ2,变倍组和补偿组一起同步运动直到预定的总变倍比为止,到达变倍比的要求的最终状态为系统短焦位置。根据上述连续变焦微分模型,利用(10)~(17)式,可解得各组分光焦度分配及间隔位置关系。

    针对复杂的变焦光学系统,建立模型是至关重要的环节,尽管根据(10)~(17)式可推导求解变焦系统高斯参数,但是当光学系统初始参数选择不合适, 会使得系统光焦度分配不合理导致系统像差校正难度大、间隔不合适导致各组件运动中相互碰撞等情况发生。鉴于连续变焦系统的复杂性,系统建模难度大,因此根据连续变焦理论模型编制连续变焦参数计算程序,辅助建立理想光学模型。建模后,设计工作的重点将放在选型以及评价函数的设置和动态修改上,使得设计系统快速收敛[6]。连续变焦光学系统设计流程如图 2所示。

    图  2  光学设计流程
    Figure  2.  Optical design flow chart

    首先,根据设计指标要求需要建立理想光学模型,确定每个组元的参数;再合理选型选材、设定评价函数,进入优化和全局优化;最后根据设计结果评价成像质量。其中函数优化和像质评价环节反复多次迭代,直至达到设计技术指标要求。

    中波制冷型连续变焦光学系统采用昆明物理研究所研制生产的长度方向仅119 mm的小型化制冷型中波红外640×512焦平面探测器组件,该探测器具体参数如表 1所示。连续变焦光学系统主要设计指标见表 2

    表  1  探测器参数
    Table  1.  Parameters of detector
    Detector HgCdTe
    Array size 640×512
    Pixel dimension/μm 15
    NETD/mK ≤22
    Spectral response/μm 3.7−4.8 μm
    Weight ≤380 g
    下载: 导出CSV 
    | 显示表格
    表  2  光学系统设计指标
    Table  2.  Parameters of optical system design
    Working waveband/μm 3.7 to 4.8
    Zoom 10:01
    Field of view 20°×16° to 2.0°×1.6°
    F# 4
    Focal length/mm 27.0~275
    Working temperature/℃ -40 to 60
    下载: 导出CSV 
    | 显示表格

    首先,按照系统指标要求,根据连续变焦理论模型编制连续变焦正组补偿参数计算程序求解变焦光学系统各组元光焦度分配及位置间隔关系。在系统初始参数取值上主要考虑以下几点:

    1)系统采用二次成像,既可以压缩前固定组直径又可以满足100%冷屏效率。中继组初始倍率取为-1×,前端变焦核心按照系统实际变焦范围进行取值,无需缩放;

    2)系统无后固定组,在理论求解中将后固定组倍率取为1×即将后固定组定为无光焦度的虚拟面;

    3)为压缩系统变焦过程中变倍组、补偿组位移量,变倍组长焦初始倍率取较大的倍率值,以符合最速变焦理论;

    4)考虑光路U型折转,补偿组和中继组之间需较大的空间安置两片平面反射镜,则补偿组取较大的焦距值;

    5)系统处于短焦位置时前固定组与变倍组应留出足够的间隔,使两组透镜不至于相碰,初始值设为0.55,补偿组与无光焦度后固定组虚拟面距离初始设为0.55。

    利用计算程序,通过反复调整系统初始参数,观察组元间隔、光焦度分配是否合适,最终确定系统初始值为:

    $$ \begin{gathered} f_{2}′=-1、f_{3}′=1.62、m_{2l}=-1.45、\\m_{3l}=-1.34、d_{12d}=0.55、d_{34d}=0.55 \end{gathered} $$

    表 3为连续变焦光学系统参数简易计算程序按照上述初始值计算得到的5个视场位置的变焦间隔参数分配结果。

    表  3  变焦系统初始间隔参数
    Table  3.  Initial spacing parameters of optical system
    Focal length/mm 275 215 150 78.6 25.9
    f1/f2 spacing/mm 83.49 80.06 74.56 59.44 18.90
    f2/f3 spacing/mm 13.02 24.01 38.57 68.69 125.94
    f3/f4 spacing/mm 67.23 59.67 50.61 35.61 18.90
    下载: 导出CSV 
    | 显示表格

    其次,将上述程序计算的变焦光学系统各组元光焦度及间隔位置参数输入ZEMAX光学设计程序,得到系统近轴光学结构,如图 3所示。通过近轴光学结构分析,系统在各位置的焦距值与计算值吻合、总长一致、间隔布局合理、变焦曲线连续,验证程序计算结果正确,可进入下一步选型工作。

    图  3  近轴光学系统
    Figure  3.  Configuration of paraxial optical system

    再次,考虑各组元材料与调焦镜的选取,由于系统需满足-40℃~+60℃工作环境下成像清晰要求,按照光学系统无热化设计理论,系统的光焦度分配、材料选取、元件间隔都要满足光焦度、消色差、消热差3个方程[7]

    $$ \sum\limits_{i = 1}^j {{h_i}{\varphi _i}} = \varphi $$ (18)
    $$ {\left( {\frac{1}{{{h_1}\varphi }}} \right)^2}\sum\limits_{i = 1}^j {h_i^2{\varphi _i}{\theta _i}} = 0 $$ (19)
    $$ {\left( {\frac{1}{{{h_1}\varphi }}} \right)^2}\sum\limits_{i = 1}^j {h_i^2{\varphi _i}{\chi _i}} = \sum\limits_{i = 1}^j {{\alpha _i}{L_i}} $$ (20)

    式中:hiϕiθiχi分别为系统各透镜组近轴光线高度、光焦度、色差系数及热差系数;h1为第一透镜近轴光线高;ϕ为系统总光焦度;αi为各透镜间隔镜筒材料的线膨胀系数;Li为各间隔镜筒长度。

    该系统采用分步设计技术实现系统连续变焦无热化。首先选取满足上述(18)(19)公式的光焦度和材料分配,实现系统常温状态连续变焦成像清晰及高低温情况下离焦量的线性变化,系统光学材料主要选用硅单晶、锗单晶及硒化锌。其次利用主动补偿技术使系统满足(20)式要求,由于系统透镜数量少,变倍组、补偿组采用凸轮轨道变焦,中继组的轴向移动会产生变倍效果即系统视场焦距发生变化,因此采用前固定组轴向移动来进行主动调焦消热。通过上述无热化设计方法,光学系统在-40℃~+60℃温度范围内保持其性能基本不变。

    最后,设置多重结构,合理设计像差评价函数,不断调整优化。为提升连续变焦系统各视场成像清晰度要求,通过设置多个高次非球面和二元衍射面,以提供更多的自由度,有利于球差、色差、像散等各类像差的校正。

    轻小型中波红外连续变焦光学系统最终设计结果如图 4所示,从上到下依次为长焦275 mm、中焦100 mm、短焦27 mm的系统图。

    图  4  连续变焦光学系统
    Figure  4.  Configuration of continuous zoom optical system

    系统前固定组采用硅单晶材料用于会聚目标景物光线、压缩变倍组透镜尺寸;变倍组采用锗单晶材料,利用锗单晶高折射率、高色散的特性实现大倍率的变焦;补偿组采用硒化锌材料主要利用其较低的温度折射率系数使其在高低温工作环境只产生较小的离焦量以便实时补偿;中继组选择具有低的温度折射率系数、低色散的硅单晶材料平衡前端变焦核的残留色差。

    整个光学系统由4片透镜、两个平面反射镜组成,其中最大透镜为第一透镜其加工直径为71 mm,平面反射镜U型折转后光学系统轴向尺寸长度172 mm,横向尺寸宽度108 mm,光学零件总重量为64 g。该系统光学透镜数量少、重量轻,光路紧凑体积小,冷屏效率100%,适配小型化制冷探测器符合连续变焦光学系统轻小型设计理念。

    由于二元衍射面在消热差、消色差方面的优异特性,系统采用了两个二元衍射面用于减少透镜数量、简化系统设计、提高系统成像质量。

    在补偿组透镜硒化锌材料上引入的二元衍射面参数为Norm Radius=15 mm,A1=-33.416,A2=4.580。经计算得到二元衍射面环带深度随透镜径向的变化如图 5所示。硒化锌二元面环带数为4,最大环带深度2.918 μm,最小环带间隔宽度为2.06 mm。该面型环带间隔宽、加工环带数量少,易于单点金刚石车削加工。

    图  5  硒化锌基底二元面环带与半径的关系
    Figure  5.  Relationship between the ring depth and radius of the ZnSe binary optical element

    在中继组硅单晶材料引入的二元衍射面参数为Norm Radius=16 mm,A1=-60.402,A2=-1.254。其环带深度随透镜径向的变化如图 6所示。硅透镜二元面环带数为9,最大环带深度1.72 μm,最小环带间隔宽度为0.86 mm。该透镜由于材料硬、环带多相对加工难度大,目前昆明物理所光学中心采用单点金刚石车削加工工艺,能制造出满足指标要求的硅基底二元光学元件。

    图  6  硅基底二元面环带与半径的关系
    Figure  6.  Relationship between the ring depth and radius of the silicon binary optical element

    本文应用动态光学理论[8],根据像移补偿公式计算补偿组运动曲线。由于变焦组和补偿组均为沿光轴的一维移动,稳像方程为:

    $$ {\beta _{2m}}{\beta _2}(1 - {\beta _{1m}}{\beta _1}){q_1} + (1 - {\beta _{2m}}{\beta _2}){q_2} = 0 $$ (21)

    式中:β1为变倍组初始位置的垂轴放大率;β1m为变倍组运动后的垂轴放大率;β2为补偿组初始位置的垂轴放大率;β2m补偿组运动后的垂轴放大率;q1为变倍组沿光轴位移量;q2为补偿组沿光轴位移量。

    式中:

    $$ {\beta _{1m}} = \frac{{{\beta _1}{f_1}'}}{{{f_1}' - {\beta _1}{q_1}}} $$ (22)
    $$ {\beta _{2m}} = \frac{{{\beta _2}{f_2}'}}{{{f_2}' + (1 - {\beta _1}{\beta _{1m}}){\beta _2}{q_2} - {\beta _2}{q_2}}} $$ (23)

    由(21)式~(23)式可得出q1q2的运动关系,即:

    $$ Aq_{2}^{2}+Bq_{2}+C=0$$ (24)

    式中:

    $$ A=(f_{1}′-β_{1}q_{1})β_{2}\text{;} $$
    $$ \begin{align} B=&β_{1}β_{2}q_{1}^{2}+[f_{2}′(1-β_{2}^{2})β_{1}-f_{1}′(1-β_{1}^{2})β_{2}]q_{1}- \\&f_{1}′f_{2}′(1-β_{1}^{2}) \end{align}$$
    $$ C=β_{2}^{2}f_{2}′[β_{1}q_{1}-f_{1}′(1-β_{1}^{2})]q_{1} $$

    得到:

    $$ {q_2} = \frac{{ - B \pm \sqrt {{B^2} - 4AC} }}{{2A}} $$ (25)

    根据上述求解公式计算该变焦系统凸轮曲线如图 7所示。变倍组最大行程为28.6 mm、补偿组最大行程35 mm;补偿组曲线变化平滑,有利于凸轮轨道加工。

    图  7  连续变焦凸轮曲线
    Figure  7.  Cam curves of continuous zoom optical system

    系统3个焦距状态下的光学传递函数(MTF)如图 8所示。在3个焦距状态下的系统MTF满足使用要求,光学系统成像质量清晰。

    图  8  光学系统传函曲线
    Figure  8.  MTF curves of continuous zoom optical system

    系统3个焦距状态下的点列图如图 9所示。各视场RMS(root mean square)均小于一个像素,最大弥散斑RMS半径为12.6 μm,小于像元尺寸;最大弥散斑几何半径为24.5 μm,与系统艾里斑半径20.9 μm相当。系统成像质量良好,满足使用要求。

    图  9  光学系统点列图
    Figure  9.  Spot diagrams of continuous zoom optical system

    系统畸变情况如图 10所示,在长焦端小视场位置时,最大畸变量为1.2%,在短焦端大视场位置时的最大畸变量为2.4%,该变焦系统畸变对连续成像无明显影响。

    图  10  光学系统畸变
    Figure  10.  Optical distortion diagrams of continuous zoom optical system

    在高低温工作环境中,系统采用轴向移动前固定组进行主动调焦消热。系统长焦端受环境温度变化,成像质量影响较大,本文主要分析长焦275 mm在高低温下经补偿后的系统成像质量。图 11为系统在高低温下长焦经补偿后的系统调制传递函数。图 12为系统在高低温下长焦经补偿后的系统点列图。从高低温传函图及点列图中看出系统在-40℃~60℃范围内成像质量良好,满足使用要求。

    图  11  系统长焦时高低温下光学传递函数
    Figure  11.  MTF curves of continuous zoom optical system at high and low temperatures
    图  12  系统长焦时高低温下点列图
    Figure  12.  Spot diagrams of continuous zoom optical system at high and low temperatures

    基于小型化制冷中波640×512、像元间距15 μm的焦平面探测器,设计一款具有SWaP-C特征的正组补偿连续变焦光学系统。系统由4片透镜两片平面反射镜组成,F#为4、视场变化范围为20°×16°~2.0°×1.6°,变倍比为10×、U型折叠后系统包络尺寸为172 mm×108 mm、最大物镜口径71 mm、光学零件总重量64 g、零件加工工艺成熟,变焦凸轮曲线平滑,在-40℃~60℃范围内保持较好的成像质量。该轻小型中波红外连续变焦光学系统在导航、搜索、跟踪、警戒、侦察等领域具有广阔的市场前景。

  • 图  1   水果可见光图像

    Figure  1.   Visible light image of fruit

    图  2   损伤部位图像卷积结果

    Figure  2.   Convolution results of damaged image

    图  3   损伤部位边缘检测结果

    Figure  3.   Edge detection results of damaged parts

    图  4   MSE数据图

    Figure  4.   MSE data map

    图  5   PSNR数据图

    Figure  5.   PSNR data map

    图  6   SSIM数据图

    Figure  6.   SSIM data map

    图  7   Laplacian数据图

    Figure  7.   Laplacian data map

    图  8   方差数据图

    Figure  8.   Variance data map

    表  1   滤波结果量化值

    Table  1   Quantized value of filtering results

    MSE σ=5 σ=10 σ=20 σ=30 σ=40 σ=50
    K=7 932.33 1033.55 1159.55 1157.57 1159.55 1159.21
    K=9 1850.98 2007.50 2525.38 2559.60 2561.29 2562.35
    K=11 3331.90 3340.24 4442.35 4816.12 4823.55 4826.10
    K=13 5057.46 5913.20 5876.27 5954.22 5985.96 6002.32
    PSNR K=7 18.4351 17.9875 17.4879 17.4953 17.4879 17.4892
    K=9 15.4568 15.1043 14.1075 14.0491 14.0462 14.0444
    K=11 12.9093 12.8930 11.6547 11.3038 11.2971 11.2948
    K=13 11.0915 10.4126 10.4398 10.3826 10.3595 10.3476
    SSIM K=7 0.961849 0.955389 0.948082 0.948056 0.948080 0.948086
    K=9 0.922108 0.926302 0.913651 0.911872 0.911722 0.911521
    K=11 0.834795 0.842830 0.793065 0.773563 0.773045 0.772895
    K=13 0.732744 0.712764 0.753693 0.755303 0.756641 0.756533
    Laplacian K=7 10.3903 11.1540 10.96140 10.3558 10.2464 10.2421
    K=9 9.63460 9.54343 10.63036 9.35341 9.34649 9.13412
    K=11 6.36219 6.24669 7.62676 7.15089 7.03636 7.00650
    K=13 6.66614 6.64652 7.12676 6.76413 6.57969 6.29320
    Variance K=7 1730.12 1736.66 1722.77 1666.19 1679.57 1679.31
    K=9 1426.23 1456.53 1367.94 1373.06 1366.66 1364.33
    K=11 1319.41 1397.31 1252.20 1274.32 1272.49 1272.55
    K=13 691.399 510.973 563.156 614.730 633.922 630.566
    下载: 导出CSV
  • [1] 李广, 吴限鑫, 王建忠, 等. 水果功能性营养成分及其检测技术研究进展[J]. 农产品质量与安全, 2019(5): 75-80. DOI: 10.3969/j.issn.1674-8255.2019.05.015

    LI Guang, WU Xianxin, WANG Jianzhong, et al. Research progress on functional nutrients in fruit and related testing technology: a reviewv[J]. Quality and Safety of Agro-products, 2019(5): 75-80. DOI: 10.3969/j.issn.1674-8255.2019.05.015

    [2] 车远侠. 水果科普知识在学前教育中的作用——评《水果的秘密》[J]. 中国果树, 2019(6): 28. https://www.cnki.com.cn/Article/CJFDTOTAL-ZGGS201906051.htm

    CHE Yuaxia. The role of fruit science knowledge in preschool education—Comment on "The Secret of Fruits"[J]. China Fruits, 2019(6): 28. https://www.cnki.com.cn/Article/CJFDTOTAL-ZGGS201906051.htm

    [3] 李承龙, 鲁明丽. 基于边缘检测技术的果品分级方法研究[J]. 常熟理工学院学报, 2018, 32(2): 78-82. DOI: 10.3969/j.issn.1008-2794.2018.02.017

    LI Chenglong, LU Mingli. Fruits Grading Algorithm Based on Detection Technology[J]. Journal of Changshu Institute of Technology, 2018, 32(2): 78-82. DOI: 10.3969/j.issn.1008-2794.2018.02.017

    [4] 李娟. 数字处理技术在果品分级中的应用——评《水果品质智能分级技术》[J]. 中国果树, 2019(6): 23. https://www.cnki.com.cn/Article/CJFDTOTAL-ZGGS201906046.htm

    LI Juan. Application of Digital Processing Technology in Fruit Grading——Comment on "Fruit Quality Intelligent Grading Technology"[J]. China Fruits, . 2019(6): 23. https://www.cnki.com.cn/Article/CJFDTOTAL-ZGGS201906046.htm

    [5] 张妍. 果蔬农产品质量安全现状分析[J]. 农业科技与信息, 2019(22): 34-35, 39. DOI: 10.3969/j.issn.1003-6997.2019.22.015

    ZHANG Yan. Analysis on the Status Quo of Quality and Safety of Fruits, Vegetables and Agricultural Products[J]. Agricultural Science-Technology and Information, 2019(22): 34-35, 39. DOI: 10.3969/j.issn.1003-6997.2019.22.015

    [6] 肖付才, 靳雯雯. 数字图像边缘检测算法及在农产品检测中的运用[J]. 热带农业工程, 2019, 43(1): 133-135. https://www.cnki.com.cn/Article/CJFDTOTAL-RDZW201901046.htm

    XIAO Fucai, Jin Wenwen. Algorithm of Digital Image Detection and Its Application in Agricultural Product Processing[J]. Tropical Agricultural Engineering, 2019, 43(1): 133-135. https://www.cnki.com.cn/Article/CJFDTOTAL-RDZW201901046.htm

    [7]

    QIU S S, WANG J. Application of Sensory Evaluation, HS-SPME GC‐MS, E‐Nose, and E‐Tongue for Quality Detection in Citrus Fruits[J]. Journal of Food Science, 2015, 80(10): s2296-s2304. DOI: 10.1111/1750-3841.13012

    [8] 马康凌, 杨绍林, 刀静梅. 中国农产品检测技术现状与展望[J]. 云南农业科技, 2019(2): 60-63. DOI: 10.3969/j.issn.1000-0488.2019.02.024

    MA Kangling YANG Shaolin, DAO Jingmei. Current Status and Prospects of Agricultural Products Testing Technology in China[J]. Yunnan Agricultural Science and Technology, 2019(2): 60-63. DOI: 10.3969/j.issn.1000-0488.2019.02.024

    [9] 徐赛, 孙潇鹏, 张倩倩. 大型厚皮水果的无损检测技术研究[J]. 农产品质量与安全, 2019(5): 30-35, 48. DOI: 10.3969/j.issn.1674-8255.2019.05.007

    XU Sai, SUN Xiaopeng, ZHANG Qianqian. Research progress on nondestructive testing technology applied to large thick-skinned fruit[J]. Quality and Safety of Agro-products, 2019(5): 30-35, 48. DOI: 10.3969/j.issn.1674-8255.2019.05.007

    [10] 门洪, 陈鹏, 邹丽娜, 等. 基于主动热红外技术的苹果损伤检测[J]. 中国农机化学报, 2013, 34(6): 220-224, 229. DOI: 10.3969/j.issn.2095-5553.2013.06.053

    MEN Hong, CHEN Peng, ZOU Lina, et al. Apples'damage detection based on active thermal infrared technique[J]. Journal of Chinese Agricultural Mechanization, 2013, 34(6): 220-224, 229 DOI: 10.3969/j.issn.2095-5553.2013.06.053

    [11] 周建民, 周其显. 基于主动热成像技术的苹果早期机械损伤检测[J]. 农机化研究, 2010, 32(8): 162-165. DOI: 10.3969/j.issn.1003-188X.2010.08.042

    ZHOU Jianmin, ZHOU Qixian. Detection of Early Mechanical Danage inApple Based on Active Thermal Imaging[J]. Journal of Agricultural Mechanization Research, 2010, 32(8): 162-165. DOI: 10.3969/j.issn.1003-188X.2010.08.042

    [12] 马超. 基于多尺度多方向的图像边缘检测算法研究及其应用[D]. 开封: 河南大学, 2019.

    MA Chao. Research and Application of Image Edge Detection Algorithm Based on Multi-scale and Multi-direction[D]. Kaifeng: Henan University, 2019.

    [13] 钟艺晶. 一种绝缘子串红外图像的特征提取方法研究[J]. 机电信息, 2020(21): 86-87. DOI: 10.3969/j.issn.1671-0797.2020.21.042

    ZHONG Yijing. Research on a Feature Extraction Method of Infrar-ed Image of Insulator String[J]. Mechanical and Electrical Information, 2020(21): 86-87. DOI: 10.3969/j.issn.1671-0797.2020.21.042

    [14] 黄玉蕾. 基于形态学滤波结合LOG算法的边缘检测[J]. 计算机测量与控制, 2019, 27(7): 257-260, 284. https://www.cnki.com.cn/Article/CJFDTOTAL-JZCK201907055.htm

    HUANG Yulei. Edge Detection Based on Morphological Filtering Combined with LOG Algorithm[J]. Computer Measurement & Control, 2019, 27(7): 257-260, 284. https://www.cnki.com.cn/Article/CJFDTOTAL-JZCK201907055.htm

    [15] 张阳, 刘缠牢, 卢伟家, 等. 基于LoG算子的双滤波边缘检测算法[J]. 电子测量技术, 2019, 42(4): 95-98. https://www.cnki.com.cn/Article/CJFDTOTAL-DZCL201904017.htm

    ZHANG Yang, LIU Chanlao, LU Weijia, et al. Improved double fi-ltering LoG operator edge detection[J]. Electronic Measurement Technology, 2019, 42(4): 95-98. https://www.cnki.com.cn/Article/CJFDTOTAL-DZCL201904017.htm

    [16] 贺萌. 基于自适应形态学的边缘检测及应用[D]. 长沙: 中南大学, 2013.

    HE Meng. Edge Detection Based on an Adaptive Morphology Algorithm and Application[D]. Changsha: Central South University, 2013.

    [17]

    SLANINA M, RICNY V. Estimating PSNR in High Definition H. 264/AVC Video Sequences Using Artificial Neural Networks[J]. Radio Engineering, 2008, 17(3): 103-108. http://www.oalib.com/paper/2490847

图(8)  /  表(1)
计量
  • 文章访问数:  315
  • HTML全文浏览量:  57
  • PDF下载量:  46
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-08-05
  • 修回日期:  2020-12-03
  • 刊出日期:  2021-06-30

目录

/

返回文章
返回
x 关闭 永久关闭

尊敬的专家、作者、读者:

端午节期间因系统维护,《红外技术》网站(hwjs.nvir.cn)将于2024年6月7日20:00-6月10日关闭。关闭期间,您将暂时无法访问《红外技术》网站和登录投审稿系统,给您带来不便敬请谅解!

预计6月11日正常恢复《红外技术》网站及投审稿系统的服务。您如有任何问题,可发送邮件至编辑部邮箱(irtek@china.com)与我们联系。

感谢您对本刊的支持!

《红外技术》编辑部

2024年6月6日