专利名称：一种快速无损检测苹果内部质量的方法
技术领域：
本发明涉及一种检测苹果内部质量的方法，特别是关于一种针对苹果水心病和褐变的快速无损检测苹果内部质量的方法。
背景技术：
我国是世界第一水果生产大国，其中苹果是最主要的品种，但每年出口量却较少，仅占总产量的2.1％。制约我国苹果出口的一个重要原因是国内对苹果的分选检测能力弱，速度慢，试验环境条件差，分选技术水平达不到国际市场的要求。苹果水心病是其在生长期和贮藏期经常发生的一种内部生理失调病，在我国各苹果产区均有发生，水心病果的果心呈水渍状，严重时病变部分接近果皮，使果皮也出现水渍状透明似蜡，后期病果内病组织可发生败坏。苹果褐变是在贮藏过程中常见的一种病害，除严重水心病会导致苹果褐变外，高CO2、低O2、冷冻害及衰老崩溃等都能引起褐变。果实褐变可发生在果肉、果心和果皮，褐变组织变坚实、湿润或松散干燥，果肉变苦，病变严重时整个果实全部褐变。由于产生苹果水心病、褐变的病果一般与正常果在外观上没有区别，因此原有对苹果内在品质的检验方法，只能是通过观察随机样品的切片来进行，但这种方法属于破坏性抽样检测的方法，不但浪费极大，而且对出口产品分级毫无意义。因此必须采用非破坏性的无损伤检测方法对苹果内部品质进行评价分级，才能适应国际市场的要求。
目前无损伤检测苹果水心病的方法包括以下几种(1)浮力密度法据Fidler等报道，根据苹果水心病果比重发生变化的特征，可通过浮力分离的原理来分离水心病果，利用酒精水溶液可分选78％的水心病果；Cavalieri等报告，根据水心病果的密度可分选100％的严重病果。但这种密度检测法只能分离质量合格的苹果与严重水心病果(4级)，而无法对相邻级别的苹果进行分离。(2)透光图像密度法据Throop报告，使用可视计算机检测苹果的透光率，并根据苹果图像估算面积、计算密度来判断苹果水心病，尽管该方法能正确分离出75％的水心病果，但不足以精确区分水心病的严重程度，原因是病果各级之间密度差异很小；而且Upchurch和Throop的进一步研究发现，检测效果与计算机可视系统的摄像机灵敏度有关。(3)X射线扫描法Shahin和Tollner报道，利用X射线扫描苹果样品，经过对图像的处理可以正确分离75％的轻微病果、40％的中等程度病果和60％的严重病果。(4)核磁共振法Clark和Bieleski报告，可使用核磁共振成像跟踪苹果在贮藏过程中水心病的变化，但该技术成本太高，检测速度慢，无法应用到苹果分选线上。Staoru Tsuchikawa等报告，采用TOE-NIR(time-of-flight Near-infrared)光谱技术，通过检测光通过样品的时间来判断水心病。(5)差分法经典的光学水心病检测方法为760nm和810nm的差分法，检测精度还是不够理想。据报道，褐变果有用600nm和740nm检测的。另外，Bramlage探讨了利用光学方法检测完整苹果的可能性；而Paola等使用time-resolved反射光谱对梨褐心进行了无损伤检测，但只能检测到果皮以下2cm范围内，对果心周围的轻微褐变检测比较困难。
发明内容
针对上述问题，本发明的目的是提供一种针对苹果水心病和褐变的快速无损检测苹果内部质量的方法。
为实现上述目的，本发明采取以下技术方案一种快速无损检测苹果内部质量的方法，它包括以下步骤(1)用光谱范围在638～1295nm的可见-近红外光谱仪对待测苹果逐个进行连续光谱检测；(2)将被测苹果的连续光谱检测结果，带入预先建立的苹果内部质量判别函数，得到苹果内部质量的判定。
所述预先建立的苹果内部质量判别函数为苹果水心病等级的判别函数，其采用以下方法建立(1)采集苹果样品集，用可见-近红外光谱仪对样品集中的样品进行波长范围在650nm～950nm的连续光谱检测，并对检测结果进行预处理，得到在连续光谱范围内OD值变化的连续光谱；(2)将每个经连续光谱检测后的样品沿赤道方向切开，对切片截面进行拍照后输入计算机，计算苹果水心病的病变程度，并按病变程度对样品分级；(3)对样品集中各样品连续光谱与切片后的水心病分级结果，采用TQ analyst 6.2分析软件包进行计算分析后，确定出与步骤(2)中苹果水心病分级相对应的判别函数；(4)将未知被测苹果的波长范围在650nm～950nm的连续光谱检测结果进行预处理后，将要判断的观测值带入到判别函数中，哪一个判别函数值最大，就属于哪一类苹果的等级。
所述步骤(2)中计算苹果水心病病变的程度，是指计算病变面积占整个截面面积的百分比，其包括五级A级＝0％，B级≤1％，C级＝1～3％，D级＝3～10％，E级≥10％以上。
所述预先建立的苹果内部质量判别函数为苹果褐变的判别函数，其采用以下方法建立(1)采集苹果样品集，用可见-近红外光谱仪对样品集中的样品进行波长范围在650nm～950nm的连续光谱检测，并对检测结果进行预处理，得到在连续光谱范围内OD值变化的连续光谱；(2)将每个经连续光谱检测后的样品沿赤道方向切开，对切片截面进行拍照后输入计算机，按照是否发生褐变分为正常果和褐变果两级；(3)将样品集中各样品在(715±5)nm、(750±5)nm、(810±5)nm三个特征波长处的OD值与切片后的苹果褐变分级结果，采用SAS 8.0 discrim进行计算分析后，确定出苹果褐变的判别函数；褐变果＝-741.58928-71.5161OD(715±5)+1508OD(750±5)-1640OD(810±5)；正常果＝-830.51874-81.46691OD(715±5)+1600OD(750±5)-1736OD(810±5)；(4)将未知被测苹果在(715±5)nm、(750±5)nm、(810±5)nm三个特征波长处的OD值，带入判别函数进行计算，两方程中计算结果数值较大的方程，被判定为代表未知被测苹果的检测结果。
所述预处理包括光谱数据的平滑、求导、多元散射校正等。
所述三个特征波长为715nm、750nm和810nm。
一种快速无损检测苹果内部质量的方法，它包括以下步骤(1)用滤光片仪，逐个对待测苹果检测波长在(715±5)nm、(750±5)nm、(810±5)nm三个特征波长处的OD值；(2)将被测苹果的(715±5)nm、(750±5)nm、(810±5)nm三个特征波长处的OD值，带入预先建立的苹果内部质量判别函数，得到苹果内部质量的判定。
所述预先建立的苹果内部质量判别函数为苹果褐变的判别函数，其采用以下方法建立(1)采集苹果样品集，用滤光片仪逐个对待测苹果检测波长在(715±5)nm、(750±5)nm、(810±5)nm三个特征波长处的OD值；(2)将每个检测后的样品沿赤道方向切开，对切片截面进行拍照后输入计算机，按照是否发生褐变分为正常果和褐变果两级；(3)对样品集中各样品在(715±5)nm、(750±5)nm、(810±5)nm三个特征波长处的OD值与切片后的苹果褐变分级结果，采用SAS8.0discrim法进行计算分析后，确定出苹果褐变的判别函数褐变果＝-741.58928-71.5161OD(715±5)+1508OD(750±5)-1640OD(810±5)；正常果＝-830.51874-81.46691OD(715±5)+1600OD(750±5)-1736OD(810±5)；(4)将未知被测苹果在(715±5)nm、(750±5)nm、(810±5)nm三个特征波长处的OD值，带入判别函数进行计算，两方程中计算结果数值较大的方程，被判定为代表未知被测苹果的检测结果。
所述三个特征波长为715nm、750nm和810nm。
本发明由于采取以上技术方案，其具有以下优点1、本发明由于通过预先对苹果样品集进行近红外连续光谱采集，并结合对光谱采集后的样品切片进行观察和相关性分析，建立苹果水心病和褐变程度分级的判别函数，因此无论是在线检测时还是进行抽样检测时，都可以逐个判别每个苹果的等级，保证每个苹果的分选质量。2、本发明检测时，只需对苹果进行一次连续光谱的采集，便可以分别通过预先建立的两大类判别函数，进行水心病和褐变两种病症的同时检测。3、本发明根据水心病的危害程度将其分为五个等级，检测时计算机软件计算被测苹果的连续光谱属于哪一级，便可以得到水心病等级结论。4、本发明考虑到苹果褐变的不可接受性，将其仅分为褐变果和好果两个等级，而且根据实验和分析结果，确定了三个特征波长(715±5)nm、(750±5)nm和(810±5)nm处的OD值作为变量，得到褐变果和正常果的判别函数，因此检测时，计算机软件只需采集连续光谱在这三个特征波长处的OD值，便可以通过判别函数的计算得出是褐变果还是正常果的结论。本发明采用方法设备简单，可以同时检测出苹果水心病和苹果褐变，大量减少因抽样检测所造成的果品浪费，其不但采摘时可以对苹果进行现场检测，而且可以经常性地对贮藏的苹果进行监测，减少贮藏损失。
具体实施例方式
实施例1关于苹果水心病的快速无损检测一、首先建立苹果水心病的判别函数(1)样品采集及检测仪器在同区域、同一采摘时期、不同的果树株棵上、同一果树株棵上的不同部位上，随机采摘具有代表性的红富士苹果样品200个，果径范围52.4mm～93.3mm。
采用CCD光栅型可见-近红外光谱仪，其波长范围为638nm～1295nm，计算机内装有TQ analyst 6.2分析软件包(美国NICOLET公司)，用于对光谱数据进行处理。
(2)对样品进行连续光谱采集用CCD光栅型可见-近红外光谱仪对样品集中的各样品进行连续光谱采集，波长范围为638nm～1295nm，每隔0.35nm照射一个点，每张光谱共照射2047个点，光源位于仪器上部，监测器位于底端，两者成180°，积分时间为1000ms，平均照射次数为10次；可以选择多个测点，取平均光谱用作分析。由于在650nm之前和950nm之后的波长峰值不明显，故仅选择波长范围在650nm～950nm之间的光谱。
(3)将样品切片观察将每个检测后的样品，沿赤道方向切开，对切片的截面进行拍摄，并将照片输入计算机，通过Scion Image图像处理软件对照片进行内部病变面积的计算，根据病变面积占整个截面面积的百分比确定出每个样品水心病的病变程度，并将所有样品苹果按其病变程度人为地划分为五个级别，各级别如下A级＝0％，B级≤1％，C级＝1～3％，D级＝3～10％，E级＝≥10％以上。
(4)分析各样品光谱与切片结果的相关性，建立苹果水心病判别函数根据测试过程的精确记录，利用TQ analyst 6.2分析软件包(美国NICOLET公司)对样品集中各样品的近红外线检测结果与相应样品切片后的内部质量分级结果进行相关性分析。首先对每条光谱进行Savitzky-Golay25点3次平滑预处理，为了消除苹果大小和散射造成的影响，光谱透过苹果内部时的路径变化，利用SNV(标准正态变量)进行校正。程序采用分类方法中的discriminant analysis(判别分析)，对所选光谱范围进行分析判断，计算机分别自动给出五个等级水心病病变对应的判别函数，但无具体判别函数显示，其一般表达形式为组(A、B、C、D、E)＝K0+K1×X1+K2×X2+…+Kn×XnK0-常数项，Kn-变量系数，Xn-波长该程序通过计算各级之间的Mahalanobis(马氏)距离，确定与已知材料级别中最相近的未知级别作为计算级。即将新观测值带入到五个判别函数中，哪个判别函数值最大，该观测就属于哪一类。
二、快速无损检测苹果水心病采用CCD光栅型可见-近红外光谱仪，对待测苹果逐个进行波长范围在650nm～950nm的连续光谱检测，然后将被测苹果的连续光谱检测结果带入上述苹果内部质量判别函数，通过计算机软件的计算结果，可以得到被测苹果更接近哪一级。下面是应用本发明方法进行的一组苹果水心病检测判别及验证结果(如表1所示)表1快速无损检测苹果水心病判别结果
本实施例经对受测苹果切片检验表明，利用本发明方法快速无损检测苹果水心病，其正确判别率如表2所示。
表2正确判别率
实施例2关于苹果褐变的快速无损检测一、首先建立苹果褐变的判别函数(1)样品采集及检测仪器在同区域、同一采摘时期、不同的果树株棵上、同一果树株棵上的不同部位上，随机采摘具有代表性的富士苹果样品184个，果径范围7.33cm～9.5cm。将所有样品苹果分四袋装入聚氯乙烯薄膜(厚0.8mm)密封袋中，分别充入CO2气体并使其含量达到30％，O2含量为5％。当聚氯乙烯薄膜袋自然漏气发扁时即充入空气使其保持一定的气压，并在室内自然状态下放置一个月，促使苹果发生褐变后取出。
采用CCD光栅型可见-近红外光谱仪，其波长范围为638nm～1295nm，计算机内装有TQ analyst 6.2分析软件包(美国NICOLET公司)，用于对光谱数据进行处理。
(2)对样品进行连续光谱照射用CCD光栅型可见-近红外光谱仪对样品集中各样品进行连续光谱采集，波长范围为638nm～1295nm；每隔0.35nm照射一个点，每张光谱共照射2047个点，光源位于仪器上部，监测器位于底端，两者成180°，积分时间为1000ms，平均照射次数为10次。每个苹果可以照射三次，每次水平旋转120度，对三次检测的平均光谱作分析。检测时可以仅从果柄向上朝向光源的一个方向做检测，也可以从两个方向做检测，一个是沿果轴方向，一个是垂直果轴方向。
(3)将样品切片观察将每个检测后的样品，沿赤道方向切开，对切片的截面进行拍摄，并将照片输入计算机，通过Scion Image图像处理软件对照片进行苹果内部病变面积的计算，根据病变面积占整个截面面积的百分比确定出每一个苹果褐变的程度，并将褐变程度划分为六个级别，各级别如下0级＝0％，1级＝0～10％，2级＝10～25％，3级＝25～40％，4级＝40～60％，5级＞60％。
(4)分析光谱与切片结果的相关性建立苹果褐变判别函数苹果发生褐变对消费者来说就是不可接受的，所以检测褐变苹果只需要有两个指标，一个是好一个是坏(褐变)，我们把不同程度的褐变归为一类，与好果做判别分析。选择三个单波长715nm、750nm、810nm的OD值，作为判别指标，715nm是褐变变化最明显的波长，750nm处是吸收峰，810nm吸收最少。通过SAS 8.0discrim一般判别程序计算结果，得到苹果褐变果和正常果的判别函数褐变果＝-741.58928-71.5161OD715+1508OD750-1640OD810正常果＝-830.51874-81.46691OD715+1600OD750-1736OD810将未知被测苹果在715nm、750nm、810nm处的OD值，带入判别函数进行计算，两方程中计算结果数值较大的方程，被判定为代表未知被测苹果的检测结果。
二、快速无损检测苹果褐变检测苹果褐变，可以采用两种方法，一种是与苹果水心病同时检测苹果褐变；另一种是单独检测苹果褐变。与水心病同时检测时，与水心病采用同一个CCD光栅型可见-近红外光谱仪，对待测苹果逐个进行波长范围在650nm～950nm的连续光谱检测，然后将检测结果采用同一个计算机软件进行预处理，再将被测苹果的连续光谱带入水心病判别函数；将连续光谱中波长为715nm、750nm、810nm三点的OD值带入苹果褐变的判别函数，根据两方程计算的数值，按其中数值较大的方程所代表的褐变果或正常果，作为被测苹果质量的判定结果。
如果单独检测苹果褐变，只需检测波长为715nm、750nm、810nm三点的OD值，并将这三点的OD值带入苹果褐变的判别函数，再根据上述同样方法判定褐变果或正常果即可。单独检测苹果褐变可以采用其它更为简单的设备(如滤光片型仪器等)，得到波长为715nm、750nm、810nm三点的OD值。
以下是采用本发明方法进行一组苹果褐变检测判断和验证的结果(如表3所示)表3快速无损检测苹果褐变程度的判别结果
由表3可见，184个被测苹果，实际褐变果有126个，快速无损检测正确判别出褐变果123个，正确率97.63％；误判率为正常果的3个，误判率2.38％。样品苹果实际正常果有58个，快速无损检测正确判别出正常果的53个，正确率91.38％；误判率为褐变果的5个，误判率8.62％。这个结果表明，本发明快速无损检测苹果褐变程度的方法，完全可以用于快速无损检测苹果内部的质量。
上述实施例中，三个特定波长715nm、750nm、810nm的选定是可以有所变化的，但一般在±5nm的范围内变化即可。
由于各地区，各不同树木，及每年的气候情况都会造成苹果内部的情况有所不同，因此在使用时本发明时，可以根据测量结果和实际观察情况，按本发明提供的方法对建立的判别模型进行一些调整，这些调整不应排除在本发明的保护范围之外。
权利要求
1.一种快速无损检测苹果内部质量的方法，它包括以下步骤(1)用光谱范围在638～1295nm的可见-近红外光谱仪对待测苹果逐个进行连续光谱检测；(2)将被测苹果的连续光谱检测结果，带入预先建立的苹果内部质量判别函数，得到苹果内部质量的判定。
2.如权利要求1所述的一种快速无损检测苹果内部质量的方法，其特征在于所述预先建立的苹果内部质量判别函数为苹果水心病等级的判别函数，其采用以下方法建立(1)采集苹果样品集，用可见-近红外光谱仪对样品集中的样品进行波长范围在650nm～950nm的连续光谱检测，并对检测结果进行预处理，得到在连续光谱范围内OD值变化的连续光谱；(2)将每个经连续光谱检测后的样品沿赤道方向切开，对切片截面进行拍照后输入计算机，计算苹果水心病的病变程度，并按病变程度对样品分级；(3)对样品集中各样品连续光谱与切片后的水心病分级结果，采用TQ analyst6.2分析软件包进行计算分析后，确定出与步骤(2)中苹果水心病分级相对应的判别函数；(4)将未知被测苹果的波长范围在650nm～950nm的连续光谱检测结果进行预处理后，将要判断的观测值带入到判别函数中，哪一个判别函数值最大，就属于哪一类苹果的等级。
3.如权利要求2所述的一种快速无损检测苹果内部质量的方法，其特征在于所述步骤(2)中计算苹果水心病病变的程度，是指计算病变面积占整个截面面积的百分比，其包括五级A级＝0％，B级≤1％，C级＝1～3％，D级＝3～10％，E级≥10％以上。
4.如权利要求1所述的一种快速无损检测苹果内部质量的方法，其特征在于预先建立的苹果内部质量判别函数为苹果褐变的判别函数，其采用以下方法建立(1)采集苹果样品集，用可见-近红外光谱仪对样品集中的样品进行波长范围在650nm～950nm的连续光谱检测，并对检测结果进行预处理，得到在连续光谱范围内OD值变化的连续光谱；(2)将每个经连续光谱检测后的样品沿赤道方向切开，对切片截面进行拍照后输入计算机，按照是否发生褐变分为正常果和褐变果两级；(3)将样品集中各样品在(715±5)nm、(750±5)nm、(810±5)nm三个特征波长处的OD值与切片后的苹果褐变分级结果，采用SAS 8.0 discrim进行计算分析后，确定出苹果褐变的判别函数；褐变果＝-741.58928-71.5161OD(715±5)+1508OD(750±5)-1640OD(810±5)正常果＝-830.51874-81.46691OD(715±5)+1600OD(750±5)-1736OD(810±5)(4)将未知被测苹果在(715±5)nm、(750±5)nm、(810±5)nm三个特征波长处的OD值，带入判别函数进行计算，两方程中计算结果数值较大的方程，被判定为代表未知被测苹果的检测结果。
5.如权利要求2或4所述的一种快速无损检测苹果内部质量的方法，其特征在于所述预处理包括光谱数据的平滑、求导、多元散射校正等。
6.如权利要求4所述的一种快速无损检测苹果内部质量的方法，其特征在于所述三个特征波长为715nm、750nm和810nm。
7.一种快速无损检测苹果内部质量的方法，它包括以下步骤(1)用滤光片仪，逐个对待测苹果检测波长在(715±5)nm、(750±5)nm、(810±5)nm三个特征波长处的OD值；(2)将被测苹果的(715±5)nm、(750±5)nm、(810±5)nm三个特征波长处的OD值，带入预先建立的苹果内部质量判别函数，得到苹果内部质量的判定。
8.如权利要求7所述的一种快速无损检测苹果内部质量的方法，其特征在于所述预先建立的苹果内部质量判别函数为苹果褐变的判别函数，其采用以下方法建立(1)采集苹果样品集，用滤光片仪逐个对待测苹果检测波长在(715±5)nm、(750±5)nm、(810±5)nm三个特征波长处的OD值；(2)将每个检测后的样品沿赤道方向切开，对切片截面进行拍照后输入计算机，按照是否发生褐变分为正常果和褐变果两级；(3)对样品集中各样品在(715±5)nm、(750±5)nm、(810±5)nm三个特征波长处的OD值与切片后的苹果褐变分级结果，采用SAS 8.0discrim法进行计算分析后，确定出苹果褐变的判别函数褐变果＝-741.58928-71.5161OD(715±5)+1508OD(750±5)-1640OD(810±5)正常果＝-830.51874-81.46691OD(715±5)+1600OD(750±5)-1736OD(810±5)(4)将未知被测苹果在(715±5)nm、(750±5)nm、(810±5)nm三个特征波长处的OD值，带入判别函数进行计算，两方程中计算结果数值较大的方程，被判定为代表未知被测苹果的检测结果。
9.如权利要求7或8所述的一种快速无损检测苹果内部质量的方法，其特征在于所述三个特征波长为715nm、750nm和810nm。
全文摘要
本发明涉及一种快速无损检测苹果内部质量的方法，它包括以下步骤(1)将光谱范围在638～1295nm的可见－近红外光谱仪对待测苹果逐个进行连续光谱检测；(2)将被测苹果的连续光谱检测结果，带入预先建立的苹果内部质量判别函数，得到苹果内部质量的判定。本发明检测时，只需对苹果进行一次连续光谱的采集，便可以分别通过预先建立的两大类判别函数，进行水心病和褐变两种病症的同时检测。本发明可以大量减少因抽样检测所造成的果品浪费，其不但采摘时可以对苹果进行现场检测，而且可以经常性地对贮藏的苹果进行监测，减少贮藏损失。
文档编号G06F19/00GK1789978SQ20041009862
公开日2006年6月21日 申请日期2004年12月14日 优先权日2004年12月14日
发明者韩东海, 刘新鑫 申请人:中国农业大学