专利名称：水果内部品质快速检测模型的建模方法
技术领域：
本发明涉及一种水果内部品质检测模型的建模方法，特别涉及一种基于计算机断层扫描快速检测水果内部品质的模型的建模方法，属于农产品无损检测技术领域。
背景技术：
水果的品质检测在日常的生产生活中有很广泛的应用。但是传统的化学检测方法，比如测定糖的蒽酮比色法，要对水果进行复杂的前处理，不但步骤繁琐，费时费力，并且有一定的实验技能的要求。因此对于农产品快速无损检测的方法，国内外的很多学者都正在研究光学特性法由于农产品的内部成分和外部特性不同，在不同波长光线照射下会有不同的吸收、反射或散射特性，即农产品的分光反射率或吸收率在某一特定波长内会比 其他部分大，根据这一特性结合光学检测装置实现农产品品质的无损检测，主要包括紫外光谱法、可见光谱法和近红外光谱检测分析法。该方法利用实验室的光谱仪器，而且受限于物料对于光谱的消耗，导致有效的光谱信息提取困难，影响其检测精度。高效液相色谱法由高效液相色谱对待测样品进行分析，虽然此法可以得到非常高的准确度，但是由于仪器对分析样品的要求非常高，同样需要较为复杂繁琐的处理。前人的研究多是采用同一种水果，但是对于不同产地的苹果，由于其生长环境的极大差异，水果主要成分的含量也有很大的区别。CT技术是利用软X射线穿透物体的原理，通过扫描图像可以得出被摄物体内部确切位置的部分物质的特性信息。在医学，无损检测等方面得到了广泛的应用。

发明内容
本发明的目的在于提供一种基于计算机断层扫描的水果内部品质的快速检测模型的建模方法。该方法可以根据水果的产地，水果的贮藏时间以及水果计算机断层扫描的CT值准确预测出水果的内部品质，以实现水果无损检测。本发明的水果内部品质快速检测模型的建模方法，其步骤如下I)试验样品的准备挑选大小均匀，无伤疤，虫蛀的水果作为待测水果，试验期间将待测水果贮藏于不见光透气的纸箱中，保持室温，空气湿度在509^70%之间；2) CT扫描及数据来源的确定将待测水果放入CT机进行纵向扫描，扫描结束后，选取水果最大剖面的CT图像，作为获取CT值的数据来源。3) CT值的获取在最大水果剖面的CT图像上，均匀选取6个直径Icm的圆形区域作为试验区域并编号，其中3个靠近果核,另外3个远离果核,利用eFilm Workstation软件获取各试验区域CT值。4)水果内部品质参数的测量在待测水果最大剖面对应的6个试验区域，分别切取直径1cm、厚度为2cm的圆形水果块，按相应国家标准分别检测各水果块的可溶性固形物含量、酸度、PH值以及含水率；
5)周期试验按5 7天为一个周期，共测试6 7个周期，每个周期重复上述步骤I) 4);6)水果产地参数的确定根据检测数据分别计算同一产地水果的可溶性固形物含量、酸度、PH值以及含水率的平均值，将此平均值分别作为该水果可溶性固形物含量模型、酸度模型、PH值模型以及含水率模型中产地的参数；7)预测模型的建立将步骤3)中得到的各试验区域CT值、步骤4)中得到的水果可溶性固形物含量、酸度、PH值以及含水率品质参数、步骤5)的周期试验总天数以及步骤
6)中得到的水果产地参数，利用SAS数据处理软件中的最小二乘拟合方法进行拟合，分别建立水果可溶性固形物含量、酸度、pH值以及含水率的预测模型Y=Bc^a1Xja2Xja3Xja4X1Wa5X1Xfa6X1Xja7X2Wa8X2Xja9X32,其中 X1 是水果各试验区域的CT值，X2是水果的贮藏总天数，X3是相应模型的水果产地参数，Y是水果相应的可溶性 固形物含量、酸度、PH值以及含水率的预测值，BtTa9是由最小二乘拟合方法得到的方程各项的系数。本发明中，所说的待测水果为苹果、梨或桃。为了全面提取水果的信息，步骤3)中所述的3个靠近果核和3个远离果核的试验区域通常采用沿水果剖面周向交替选取。本发明中，对水果进行CT扫描的参数设置为，扫描参数扫描方式扇形形束360度旋转扫描，并行连续数据采集；扫描速度2. 0秒；连续X射线扫描层厚5mm ；扫描与重建并行进行方向角0-355度，5度递增；最大扫描范围500mm ;扫描层厚5mm ；进床速度37. 5mm /秒，窗宽460 520，窗位 180 210 ；图象质量参数空间分辨率-15.01p/cmi0%MTF ；密度分辨率5.OmmiO. 3%。本发明方法有如下优点，I)操作简便，耗时短，只需将水果经过计算机断层扫描的CT值和贮藏的天数，就可以预测水果内部品质参数；2)把不同的产地作为预测水果内部品质的一个参数，将离散的变量通过相应品质参数的平均值将其变为连续的参数，提高预测模型的精度；3)灵敏度高，重复性好，可以实现水果无损检测。


图I是本发明的水果内部品质快速检测模型的建模方法流程图；图2是最大水果剖面的试验区域选取示意图。
具体实施例方式以下结合附图并通过实施例对本发明作进一步的说明
以不同产地的苹果酸度、可溶性固形物含量的检测模型的建模为例，建模过程参照图I所示流程图。I)试验样品的准备本试验中苹果分别为山东红富士，陕西红富士和新疆阿克苏糖心富士苹果，挑选大小均匀，无伤疤，虫蛀的苹果作为待测水果，试验期间将待测苹果贮藏于不见光透气的纸箱中，保持室温，空气湿度在509^70%之间；2) CT扫描及数据来源的确定将待测苹果放入CT机，进行纵向扫描，扫描结束后，选取苹果最大剖面的CT图像，作为获取CT值的数据来源。3) CT值的获取如图2所示，在最大苹果剖面的CT图像上，均匀选取6个直径Icm的圆形区域作为试验区域并编号，其中3个靠近果核，另外3个远离果核，靠近果核的试验区域和远离果核的试验区域交替选取,编号1，3，5的试验区域远离果核A,编号2，4，6的试验区域靠近果核A,利用eFilm Workstation软件获取各试验区域CT值。4)苹果内部品质参数的测量在待测苹果最大剖面对应的6个试验区域，分别切取直径1cm、厚度为2cm的圆形苹果块，参照国家标准《GB/T 12293-1990水果、蔬菜制品可滴定酸度的测定》检测苹果的酸度，参照国家标准《GB/T12295-1990水果、蔬菜制品可溶性固形物含量的测定折射仪法》检测苹果的可溶性固形物含量。5)周期试验按7天为一个周期，共测试7个周期，每个周期重复上述步骤I)
4)；6)苹果产地参数的确定计算山东红富士苹果可溶性固形物含量平均值为12. 41，酸度平均值为I. 27，将此平均值分别作为山东红富士苹果可溶性固形物含量模型、酸度模型中产地的参数；计算陕西红富士苹果可溶性固形物含量平均值为12. 78，酸度平均值为I. 92，将此平均值分别作为陕西红富士苹果可溶性固形物含量模型、酸度模型中产地的参数；计算新疆阿克苏糖心红富士苹果可溶性固形物含量平均值为12. 62，酸度平均值为0. 68，将此平均值分别作为新疆阿克苏糖心红富士苹果可溶性固形物含量模型、酸度模型中产地的参数；7)预测模型的建立利用SAS数据处理软件中的最小二乘方法进行拟合，将步骤
3)中得到的CT值、步骤4)中得到的苹果品质参数，步骤5)中获得的苹果贮藏天数以及步骤6)中得到的代表苹果产地的数值，建立山东红富士，陕西红富士和新疆阿克苏糖心富士苹果酸度的统一模型，如公式(I)所示Y=O. 253192+0. 008691Xl+0. 003464x2+0. 671846x3+0. 000030405x^+0. 000053048xlX2+0. 000412x22-0. 005817Xix3-0. 020717x2x3+0. 013408x32(I)式中Xl为苹果试验区域的CT值，X2为苹果贮藏的天数，X3为苹果产地参数，Y为苹果的酸度。建立山东红富士，陕西红富士和新疆阿克苏糖心富士苹果可溶性固形物含量含量的统一模型，如公式(2)所不Y=2609. 843001-0. 083779Xl+0. 198376x2_415. 83825xs+0. 000017804x^+0. 000191X1XjO. 000597x22+0. 005451x03-0. 019874x2x3+16. 624252x32(2)式中Xl为苹果试验区域的CT值，X2为苹果贮藏的天数，X3为苹果产地参数，Y为苹果的可溶性固形物含量。
为了便于验证所建模型的准确性，本例步骤7)采取随机选取6个试验区域中的5个试验区域的试验数据进行建模，用剩下的I个试验区域的实测数据进行模型的验证。验证结果见表I表I误差率验证表
权利要求
1.水果内部品质快速检测模型的建模方法，其步骤如下 1)试验样品的准备挑选大小均匀，无伤疤，虫蛀的水果作为待测水果，试验期间将待测水果贮藏于不见光透气的纸箱中，保持室温，空气湿度在509^70%之间； 2)CT扫描及数据来源的确定将待测水果放入CT机进行纵向扫描，扫描结束后，选取水果最大剖面的CT图像，作为获取CT值的数据来源； 3)CT值的获取在最大水果剖面的CT图像上，均匀选取6个直径Icm的圆形区域作为试验区域并编号，其中3个靠近果核，另外3个远离果核,利用eFilm Workstation软件获取各试验区域CT值； 4)水果内部品质参数的测量在待测水果最大剖面对应的6个试验区域，分别切取直径1cm、厚度为2cm的圆形水果块，按相应国家标准分别检测各水果块的可溶性固形物含量、酸度、PH值以及含水率； 5)周期试验按5 7天为一个周期，共测试6 7个周期，每个周期重复上述步骤I) 4)； 6)水果产地参数的确定根据检测数据分别计算同一产地水果的可溶性固形物含量、酸度、PH值以及含水率的平均值，将此平均值分别作为该水果可溶性固形物含量模型、酸度模型、PH值模型以及含水率模型中产地的参数； 7)预测模型的建立将步骤3)中得到的各试验区域CT值、步骤4)中得到的水果可溶性固形物含量、酸度、pH值以及含水率品质参数、步骤5)的周期试验总天数以及步骤6)中得到的水果产地参数，利用SAS数据处理软件中的最小二乘拟合方法进行拟合，分别建立水果可溶性固形物含量、酸度、pH值以及含水率的预测模型A=Bc^a1Xja2Xfa3Xfa4X1Wa5X1X^a6X1Xfa7X2Wa8X2Xfa9X32,其中X1是水果各试验区域的CT值，X2是水果的忙藏总天数，X3是相应模型的水果产地参数，Y是水果相应的可溶性固形物含量、酸度、pH值以及含水率的预测值，B0^a9是由最小二乘拟合方法得到的方程各项的系数。
2.根据权利要求I所述的水果内部品质快速检测模型的建模方法，其特征在于，所说的待测水果为苹果、梨或桃。
3.根据权利要求I所述的水果内部品质快速检测模型的建模方法，其特征在于，对水果进行CT扫描的参数设置为， 扫描参数 扫描方式扇形形束360度旋转扫描，并行连续数据采集； 扫描速度2.0秒； 连续X射线扫描层厚5mm; 扫描与重建并行进行方向角0-355度，5度递增； 最大扫描范围500mm ;扫描层厚5mm ； 进床速度37. 5mm /秒， 窗宽460 520，窗位180 210 ； 图象质量参数 空间高分辨率-8. 51p/cmi50%MTF 密度分辨率5. OmmiO. 3%。
4.根据权利要求I所述的水果内部品质快速检测模型的建模方法，其特征在于，步骤3)中所述的3个靠近果核和3个远离果核的试验区域沿水果剖面周向交替选取 。
全文摘要
本发明公开了一种水果内部品质快速检测模型的建模方法，通过测量水果贮藏期内水果的CT值，进而根据不同水果的产地，水果的贮藏时间和CT值，建立预测模型Y=a0+a1x1+a2x2+a3x3+a4x12+a5x1x2+a6x1x3+a7x22+a8x2x3+a9x32，可用于预测水果内部品质参数，其中x1是水果区域的平均CT值，x2是水果的贮藏天数，x3是相应模型的水果产地参数，Y是相应的品质参数的预测值。采用本发明只需将水果经过计算机断层扫描的CT值和贮藏的天数，就可以预测水果内部品质参数；预测精度高；重复性好，可以实现水果无损检测。
文档编号G01N23/04GK102706908SQ201210179569
公开日2012年10月3日 申请日期2012年5月30日 优先权日2012年5月30日
发明者孙腾, 张京平 申请人:浙江大学