[发明专利]基于太赫兹时域光谱的峰谷位鉴别转基因稻米的方法

权利要求书

1.基于太赫兹时域光谱的峰谷位鉴别转基因稻米的方法，包括如下的步骤：

1）制作转基因稻米以及待测稻米的样本

将转基因稻米以及待测稻米分别烘干、磨粉、过筛，然后将粉末与聚乙烯粉末以一定比例混匀，并用压片机一一压成薄片，获得多个转基因稻米样本以及待测稻米样本；

2）采集所有转基因稻米样本的太赫兹时域光谱

在0.1-3.5THz光谱频宽区间内，用太赫兹时域光谱系统一一采集转基因稻米样本的太赫兹时域光谱，获得所有转基因稻米样本各自独立的光谱数据组；

3）采用数据处理方法对所获得每个转基因稻米的谱线幅值数据组以及对应时间数据组分别进行平滑处理：得平滑处理后的新谱线幅值数据组以及对应时间数据组；

4）在每个转基因稻米样本的新谱线幅值数据组以及对应时间的数据组中搜索，寻找太赫兹时域光谱的峰值和谷值出现的时间；然后对所有转基因稻米样本的峰值和谷值出现时间进行统计，获得转基因稻米的峰值时间段和谷值时间段，并记录为时间段1、2、3和4；

5）对待测样本按步骤2）和步骤3）方法处理，然后在平滑处理后的数据组中搜索峰值和谷值出现的时间；如果峰值和谷值出现的时间都落在步骤4）的时间段1、2、3和4内，则判断样本为转基因稻米；否则就不是转基因稻米。

2.根据权利要求1所述的基于太赫兹时域光谱的峰谷位鉴别转基因稻米的方法，其特征在于所述步骤4）中寻找太赫兹时域光谱的峰值和谷值出现的时间，是在规定延迟时间区域内搜索转基因稻米谱线数据组中出现第一个波谷的出现位置、第一个波峰的出现位置、第二个波谷的出现位置和第二个波峰的出现位置，具体是在规定延迟时间区域内分别寻找出现最小波谷、最大波峰、最小波谷和最大波峰四个谱线幅值时的对应延迟时间；

所述规定延迟时间区域完整包括每个转基因稻米样本所有谱线出现的时间。

3.根据权利要求2所述的基于太赫兹时域光谱的峰谷位鉴别转基因稻米的方法，其特征在于所述步骤5）中搜索峰值和谷值出现时间的方法是：

在时间段1内寻找谷，当相邻一阶导数乘积为负值，并且前面的一阶导数小于零时，确认此时间段内存在谷，并通过mark1[0]记录；在时间段2内寻找峰，当相邻一阶导数乘积为负值，并且前面的一阶导数大于零时，确认此段内存在峰，并通过mark1[1]记录；用同样方法在时间段3和4内分别寻找谷和峰，分别用mark1[2]和mark1[3]记录；当在规定时间段内峰谷都出现的时候，就判定为转基因稻米，否则就不是。

4.根据权利要求2或3所述的基于太赫兹时域光谱的峰谷位鉴别转基因稻米的方法，其特征在于所述每个样本的粉末质量在0.1-0.3克，压片机压力为15-20MPa。

5.根据权利要求4所述的基于太赫兹时域光谱的峰谷位鉴别转基因稻米的方法，其特征在于所述每个样本的粉末与聚乙烯粉末的混合重量比例为1:18-22。

6.根据权利要求5所述的基于太赫兹时域光谱的峰谷位鉴别转基因稻米的方法，其特征在于所述太赫兹时域光谱系统，在充氮气的情况下采集样本的太赫兹时域光谱；测量环境湿度要求<0.2%。

7.根据权利要求6所述的基于太赫兹时域光谱的峰谷位鉴别转基因稻米的方法，其特征在于所述平滑处理，是对每个转基因稻米的谱线幅值以及对应时间分别以连续排列的h个数值进行加权平均，获得新的脉冲幅值；即作变换t=（x1+x2+……+xh）/h，得平滑后的新数据组；t为新产生的脉冲幅值，h为连续排列的脉冲幅值个数。

8.根据权利要求7所述的基于太赫兹时域光谱的峰谷位鉴别转基因稻米的方法，其特征在于所述转基因稻米样本的数量为10-300个。