本发明属于光谱检测领域，涉及一种桃果实内部硬度空间分布成像的方法。
背景技术：
桃[Prumus persica(L.)Batsch]是属于蔷薇科(Rosaceae)，李属(Prunus L.)桃亚属(Amygdalus L.)的核果类植物。桃原产中国，是我国主要的水果种类之一。软化后的桃果实适口性和风味发生变化，达到理想食用状态。但另一方面软化后桃果实更容易遭受物理伤害和病原浸染，从而缩短贮藏寿命，使桃果产品在市场上的流通和供货期受到一定限制。因此需要深入研究桃果实成熟软化机制，为改进和完善现有桃果实采后物流方式提供理论依据和技术支持。
果实硬度是衡量果实成熟度和贮藏品质的重要指标之一。桃果实的硬度测量通常使用质构仪法，方法是在果实赤道面上选择一点，去皮，选择探头直径穿刺果实一定深度来测定硬度。质构仪法的局限性在于仅能获得被测的小块果实的整体硬度，属于单点测量，无法得知果实内部硬度的空间分布状况，且批量测量时对测定点的要求较高。
目前，高光谱成像技术对桃及其他水果的硬度检测已有报道，一般通过采用高光谱成像技术获取水果表面高光谱图像，从而实现水果硬度的快速无损检测。上述的表面无损检测均停留在果实表面，无法得到果实内部果实硬度状况。无法在果实生长发育和贮藏保鲜阶段对果实内部的硬度变化机制的进一步深入研究、乃至像素级的细化研究提供支持。
技术实现要素：
本发明的目的在于针对现有桃果实硬度检测方法仅在表面检测，无法获得果胶在桃果实内部的空间分布情况的问题，提出一种桃果实内部硬度空间分布成像的方法，通过获取果实内部的高光谱图像，并进行数据处理后得到桃果实内部硬度空间分布分布图，为桃果实品质变化的深入研究提供了数据基础。
本发明解决其技术问题所采用的方案是：一种桃果实内部硬度空间分布成像的方法，其特征在于，包括以下步骤：
步骤1：建立基于桃果实光谱检测硬度含量的定量线性回归方程：
步骤1.1：采集n个桃果实样本分别记为M1、M2、M3、…、Mn；
步骤1.2：对于每一个桃果实样本Mi，1＜i＜n，去除果皮后，环绕果柄轴线将桃果实样本切开形成若干块扇形切块，并将各扇形切块顺着果柄轴线方向、垂直果柄轴线方向分别进行2-4等分切割，切块时候围绕果核进行，并且等分切割时不计算果核尺寸，最终形成m个果实切块，分别记为Mi，j，1＜j＜m；
步骤1.3：采集果实切块Mi，j的每个侧面的高光谱图像；
步骤1.4：采用国家标准方法测定每个果实切块Mi，j的硬度，并分别作为果实切块Mi，j的硬度参考值yi，j；依据国家标准《NY/T 2009-2011水果硬度的测定》；
步骤1.5：选取381nm、384nm、388nm、392nm、394nm、400nm、425nm、456nm、478nm、505nm、537nm、564nm、653nm、679nm、688nm、743nm、811nm、917nm波长为桃果实硬度检测的特征波长，并获取果实切块Mi，j的各个侧面的高光谱图像所有像素点在各个硬度特征波长的光谱平均值，分别记为
步骤1.6：采用多元线性回归将步骤1.5中果实切块Mi，j在各个硬度特征波长的光谱平均值与步骤1.4中的果实切块Mi，j的硬度参考值yi，j进行对应关联拟合，采用以下线性回归方程一进行关联拟合：
建立用于桃果实光谱检测硬度预测的定量线性回归方程二，为：
步骤2：获得待测桃内部果实的硬度分布图：
步骤2.1：选取桃果实待测样本N；
步骤2.2：对桃果实待测样本N进行切块，切块方法按照步骤1.2描述方法进行，共获得p个果实切块，分别记为Nz，1＜z＜p；记录下切块Nz在桃果实待测样本N中的空间坐标；
步骤2.3：采集果实切块Nz的各个侧面的高光谱图像，并分别获得图像中每个像素点在硬度特征波长381nm、384nm、388nm、392nm、394nm、400nm、425nm、456nm、478nm、505nm、537nm、564nm、653nm、679nm、688nm、743nm、811nm、917nm处的光谱值，分别记为其中，(α，β，γ)为每个像素点的坐标信息，α为横坐标信息，β为纵坐标信息，γ为切面标识信息；
步骤2.4：将步骤2.3中的果实切块Nz的各个切面每个像素点在硬度特征波长处的光谱值分别代入步骤1.6中的定量线性回归方程二中，分别得到每个像素点的对应的硬度信息计算值yz，(α，β，γ)，并根据每个像素点在所在侧面的高光谱图像中的空间坐标，分别形成果实切块Nz的各个侧面的硬度分布图；
步骤2.5：根据步骤2.4获得的果实切块Nz的各个侧面的硬度分布图，以及果实切块Nz在桃果实待测样本N中的空间坐标，采用双三次插值算法分别获得桃果实待测样本N内部硬度的网格状空间立体分布图，实现桃内部果实硬度空间可视化成像。
作为优选，步骤1.2中，每个果实切块至少有两个相对的侧面面积均不少于1cm2。
本发明将桃果实样本切割成小块，分别获取各果实切块的硬度和光谱平均值，并建立关联模型，可以根据桃果实预测待测样本的各果实切块各侧面每个像素点的光谱值预测其每个像素点对应的硬度，根据桃果实待测样本果实切块每个像素点的空间分布坐标来建立硬度在桃果实待测样本中的像素级三维分布模型，实现桃果实内部硬度分布的空间可视化成像。
具体实施方式
下面通过具体实施例对本发明做进一步说明。
实施例：一种桃果实内部硬度空间分布成像的方法，其特征在于，包括以下步骤：
步骤1：建立基于桃果实光谱检测硬度含量的定量线性回归方程：
步骤1.1：采集n个桃果实样本分别记为M1、M2、M3、…、Mn；
步骤1.2：对于每一个桃果实样本Mi，1＜i＜n，去除果皮后，环绕果柄轴线将桃果实样本切开形成若干块扇形切块，并将各扇形切块顺着果柄轴线方向、垂直果柄轴线方向分别进行2-4等分切割，切块时候围绕果核进行，并且等分切割时不计算果核尺寸，最终形成m个果实切块，分别记为Mi，j，1＜j＜m；每个果实切块至少有两个相对的侧面面积均不少于1cm2。
步骤1.3：采集果实切块Mi，j的每个侧面的高光谱图像；
步骤1.4：采用国家标准方法测定每个果实切块Mi，j的硬度，并分别作为果实切块Mi，j的硬度参考值yi，j；
步骤1.5：选取381nm、384nm、388nm、392nm、394nm、400nm、425nm、456nm、478nm、505nm、537nm、564nm、653nm、679nm、688nm、743nm、811nm、917nm波长为桃果实硬度检测的特征波长，并获取果实切块Mi，j的各个侧面的高光谱图像所有像素点在各个硬度特征波长的光谱平均值，分别记为
步骤1.6：采用多元线性回归将步骤1.5中果实切块Mi，j在各个硬度特征波长的光谱平均值与步骤1.4中的果实切块Mi，j的硬度参考值yi，j进行对应关联拟合，采用以下线性回归方程一进行关联拟合：
建立用于桃果实光谱检测硬度预测的定量线性回归方程二，为：
步骤2：获得待测桃内部果实的硬度分布图：
步骤2.1：选取桃果实待测样本N；
步骤2.2：对桃果实待测样本N进行切块，切块方法按照步骤1.2描述方法进行，共获得p个果实切块，分别记为Nz，1＜z＜p；记录下切块Nz在桃果实待测样本N中的空间坐标；
步骤2.3：采集果实切块Nz的各个侧面的高光谱图像，并分别获得图像中每个像素点在硬度特征波长381nm、384nm、388nm、392nm、394nm、400nm、425nm、456nm、478nm、505nm、537nm、564nm、653nm、679nm、688nm、743nm、811nm、917nm处的光谱值，分别记为其中，(α，β，γ)为每个像素点的坐标信息，α为横坐标信息，β为纵坐标信息，γ为切面标识信息；
步骤2.4：将步骤2.3中的果实切块Nz的各个切面每个像素点在硬度特征波长处的光谱值分别代入步骤1.6中的定量线性回归方程二中，分别得到每个像素点的对应的硬度信息计算值yz，(α，β，γ)，并根据每个像素点在所在侧面的高光谱图像中的空间坐标，分别形成果实切块Nz的各个侧面的硬度分布图；
步骤2.5：根据步骤2.4获得的果实切块Nz的各个侧面的硬度分布图，以及果实切块Nz在桃果实待测样本N中的空间坐标，采用双三次插值算法分别获得桃果实待测样本N内部硬度的网格状空间立体分布图，实现桃内部果实硬度空间可视化成像。