拉曼光谱中不同步褪色荧光的消除方法技术

一种拉曼光谱中不同步褪色荧光的消除方法，涉及一种拉曼光谱中消除荧光的方法，该方法是对样本进行连续多次光谱测量，将多次采集的光谱构成多维数组，采用ICA算法进行计算；输入ICA源信号数量，从2开始依次增加，输出分离后的信号；然后求取每个输出信号的信息熵，当出现所有输出的最小熵值时，此熵值即全局最小熵；将全局最小熵值对应的分离信号，作为拉曼光谱信号还原输出，即为消除了荧光的拉曼光谱信号。本发明专利技术可解决拉曼光谱中荧光不同步褪色、ICA中源信号数量判断、ICA输出结果随机性的问题，易于推广使用。

Elimination of asynchronous fading fluorescence in Raman spectroscopy

A method for eliminating asynchronous fading fluorescence in Raman spectra involves a method for eliminating fluorescence in Raman spectra. The method is to measure the samples continuously for many times, and form a multi-dimensional array with ICA algorithm for calculation. The number of input ICA source signals increases in turn from the beginning to output separated signals; and then each output signal is obtained. The information entropy is the global minimum entropy when the minimum entropy of all outputs appears, and the separated signal corresponding to the global minimum entropy is used as the Raman spectrum signal to restore the output, which is to eliminate the Raman spectrum signal of fluorescence. The invention can solve the problems of asynchronous fading of fluorescence in Raman spectra, judgment of the number of source signals in ICA, randomness of ICA output results, and is easy to popularize and use.

【技术实现步骤摘要】
拉曼光谱中不同步褪色荧光的消除方法
本专利技术涉及一种拉曼光谱中消除荧光的方法，特别是一种拉曼光谱中不同步褪色荧光的消除方法。
技术介绍
荧光是干扰拉曼光谱测量的一个重要原因，在与本申请最接近的现有技术中，一种是利用荧光随照射时间减弱的褪色效应，实现了同步褪色荧光的消除，其假定是体系中引起荧光来源单一，在微小时间段内，荧光强度随时间同步减弱，当找出此荧光褪色微元，通过逆差分就能找出体系所包含的荧光，然后扣除，从而获得不受荧光干扰的拉曼信号。还有一种现有技术是采用离散序列的信息熵作为荧光成分扣减判据，从混合信号中逐步扣减荧光，当信息熵达到最小时，混合信号中的荧光被完全消除。但是这两种技术适用于同步荧光测量对象，如果荧光的来源不同步，且随时间变化也不同步，将无法满足统一荧光微元假定，从而多源荧光信号难以兼顾，实现同时处理，使得荧光干扰无法被完全消除。多个来源的信号辨识是盲源信号分析领域研究的主要问题，目前最成功的方法之一是独立分量分析，即ICA（IndependentComponentAnalysis）。在确定条件下，ICA方法可以识别混合的多个独立源信号。这些确定的条件包括，各个源信号之间需要满足统计独立性，为ICA定义的源信号数量与真实信号源一致等。另外，ICA随机搜索至收敛而得到的输出结果，使得信号在幅度和相位上具有随机性，并不会与实际的真实信号一致。
技术实现思路
本专利技术要解决的技术问题是：提供一种拉曼光谱中不同步褪色荧光的消除方法，以解决现有技术中存在的拉曼光谱中荧光不同步褪色问题、ICA中源信号数量判断问题、ICA输出结果随机性问题。解决上述技术问题的技术方案是：一种拉曼光谱中不同步褪色荧光的消除方法，该方法对样本进行连续多次光谱测量，将多次采集的光谱构成多维数组，采用ICA算法进行计算；输入ICA源信号数量，从2开始依次增加，输出分离后的信号；然后求取每个输出信号的信息熵，当出现所有输出的最小熵值时，此熵值即全局最小熵；将全局最小熵值对应的分离信号，作为拉曼光谱信号还原输出，即为消除了荧光的拉曼光谱信号。本专利技术的进一步技术方案是：该方法包括以下步骤：S1.样本的光谱测量；S2.ICA计算：将连续采集的光谱构成多维数组，输入ICA算法中，依次输入不同的源信号数量，进行ICA计算，分别计算各个输出信号的信号熵，比较不同信号源数量输出的熵值，至最小熵值出现后，停止计算；将不同信号源数量输出的最小熵值作为全局最小熵值；S3.将全局最小熵值对应的拉曼光谱信号还原输出：S31.判断ICA算法中输出全局最小熵值的拉曼光谱信号对应的分量是否为反向，即，如果该分量绝对值最大处为负值，将该分量乘以-1，保证输出的拉曼光谱信号为正。S32.还原全局最小熵值对应的拉曼光谱信号强度；S33.输出调整后的全局最小熵值对应的分量，即为消除了荧光的拉曼光谱信号。本专利技术的进一步技术方案是：所述的步骤S1.样本的光谱测量中，包括以下内容：按照荧光随时间变化的复杂程度，在连续激光照射条件下，按照相同积分时间，连续采集多张光谱。本专利技术的再进一步技术方案是：在步骤S1中，所述的多张光谱为大于5张，少于20张。本专利技术的再进一步技术方案是：所述的步骤S2.ICA计算包括以下具体内容：将连续采集的光谱构成多维数组，输入ICA算法中，输入ICA源信号数进行计算；源信号数量从2开始，计算输出的2个信号的信号熵，然后增加源信号数量，依次计算各个信号的信号熵，比较不同信号源数量输出的熵值，至最小熵值出现后，停止计算；将不同信号源数量输出的最小熵值作为全局最小熵值；全局最小熵值对应的信号为样本的拉曼光谱信号；本专利技术的再进一步技术方案是：所述的步骤S32.还原全局最小熵值对应的拉曼光谱信号强度包括以下内容：S321.求取出现全局最小熵值的拉曼光谱信号对应的分量的一阶导数，计算模量NIC；S322.将直接测量的系列光谱求取均值后计算一阶导数，计算其模量NSamp；S323.将出现全局最小熵值的拉曼光谱信号对应的分量所含各点乘以系数，系数为NSamp/NIC，调整至与实际测量值一致。由于采用上述技术方案，本专利技术之拉曼光谱中不同步褪色荧光的消除方法与现有技术相比，具有以下有益效果：1.可解决拉曼光谱中荧光不同步褪色问题由于拉曼光谱中荧光、拉曼光、高频噪声是仪器测量信号的组成部分，各部分信号在频率上差异显著，而不同步褪色噪声满足了信号的时序独立；拉曼信号系统的这些特点正好满足了ICA方法所要求的信号统计独立性质，本专利技术是对样本进行连续多次光谱测量，将多次采集的光谱构成多维数组，采用ICA算法进行计算；输入ICA源信号数量，从2开始依次增加，输出分离后的信号；然后求取每个输出信号的信息熵，当出现所有输出的最小熵值时，此熵值即全局最小熵；将全局最小熵值对应的分离信号，作为拉曼光谱信号还原输出，即为消除了荧光的拉曼光谱信号。因此，本专利技术以ICA为核心方法来消除荧光的拉曼光谱信号，即可解决了拉曼光谱中不同步褪色荧光的消除问题。2.可解决ICA中源信号数量判断问题当定义的源信号数目与真实源信号数量不一致时，ICA会由于过分辨或欠分辨，导致输出结果错误；本专利技术采用信息熵判断，即依次定义不同的源信号数量，进行ICA计算，分别计算各个输出信号的信息熵，直至出现最小熵值，其中的全局最小熵值的信号对应了实际拉曼光谱所包含的信息。3.可解决ICA输出结果随机性问题ICA所输出对应于真实拉曼信息的信号，存在幅度和相位的不确定性；因为拉曼信号强度恒为正，反向ICA信号即可直接校正；另外，在导数下，荧光强度显著下降，利用一阶或者高阶导数下荧光这一性质，用信号导数的模量，将ICA输出信号校正至真实拉曼信号强度，将ICA信号还原为真实、稳定的拉曼信号输出，因此，本专利技术可解决ICA输出结果随机性问题。下面，结合附图和实施例对本专利技术之拉曼光谱中不同步褪色荧光的消除方法的技术特征作进一步的说明。附图说明图1：实施例一中直接连续测量盐酸氟桂利嗪片剂的系列光谱图，图2：实施例一中信号源定义为3的ICA信号分离效果图，图2-1：实施例一中与拉曼光谱信号对应的第Ⅱ分量图，图3：直接测量的均值信号谱图（含荧光原始谱图）和消除荧光后的拉曼信号谱图对比图。具体实施方式一种拉曼光谱中不同步褪色荧光的消除方法，该方法是对样本进行连续多次光谱测量，将多次采集的光谱构成多维数组，采用ICA算法进行计算；输入ICA源信号的数量，从2开始依次增加，输出分离后的信号；然后求取每个输出信号的信息熵，当出现所有输出的最小熵值时，此熵值即全局最小熵；将全局最小熵值对应的分离信号，作为拉曼光谱信号还原输出，即为消除了荧光的拉曼光谱信号。该方法包括以下步骤：S1.样本的光谱测量：按照荧光随时间变化的复杂程度，在连续激光照射条件下，按照相同积分时间，连续采集多张光谱，光谱数量为大于5张，少于20张；S2.ICA计算：将连续采集的光谱构成多维数组，输入ICA算法中，输入ICA源信号数量进行计算，源信号数量从2开始，计算输出的2个信号的信号熵，然后增加源数量，依次计算各个信号的信号熵，比较不同信号源数量输出的熵值，至最小熵值出现后，停止计算；将不同信号源数量输出的最小熵值作为全局最小熵值；全局最小熵值对应的信号为样本的拉曼光谱信号；本文档来自技高网...

【技术保护点】
1.一种拉曼光谱中不同步褪色荧光的消除方法，其特征在于：该方法对样本进行连续多次光谱测量，将多次采集的光谱构成多维数组，采用ICA算法进行计算；输入ICA源信号数量，从2开始依次增加，输出分离后的信号；然后求取每个输出信号的信息熵，当出现所有输出的最小熵值时，此熵值即全局最小熵；将全局最小熵值对应的分离信号，作为拉曼光谱信号还原输出，即为消除了荧光的拉曼光谱信号。

【技术特征摘要】
1.一种拉曼光谱中不同步褪色荧光的消除方法，其特征在于：该方法对样本进行连续多次光谱测量，将多次采集的光谱构成多维数组，采用ICA算法进行计算；输入ICA源信号数量，从2开始依次增加，输出分离后的信号；然后求取每个输出信号的信息熵，当出现所有输出的最小熵值时，此熵值即全局最小熵；将全局最小熵值对应的分离信号，作为拉曼光谱信号还原输出，即为消除了荧光的拉曼光谱信号。2.根据权利要求1所述的拉曼光谱中不同步褪色荧光的消除方法，其特征在于：该方法包括以下步骤：S1.样本的光谱测量；S2.ICA计算：将连续采集的光谱构成多维数组，输入ICA算法中，依次输入不同的源信号数量，进行ICA计算，分别计算各个输出信号的信号熵，比较不同信号源数量输出的熵值，至最小熵值出现后，停止计算；将不同信号源数量输出的最小熵值作为全局最小熵值；S3.将全局最小熵值对应的拉曼光谱信号还原输出：S31.判断ICA算法中输出全局最小熵值的拉曼光谱信号对应的分量是否为反向，即，如果该分量绝对值最大处为负值，将该分量乘以-1，保证输出的拉曼光谱信号为正。3.S32.还原全局最小熵值对应的拉曼光谱信号强度；S33.输出调整后的全局最小熵值对应的分量，即为消除了荧光的拉曼光谱信号。4.根据权利要求2所述的拉曼光谱中不同步褪色荧光的消除方法，其特...

【专利技术属性】
技术研发人员：姚志湘，姚桔，粟晖，
申请(专利权)人：广西科技大学，
类型：发明
国别省市：广西,45