The invention relates to a method of Chinese medicine identification based on high-dimensional image database of TCM chromatography and mass spectrometry, which includes the following steps: (1) establishing high-dimensional image database of TCM chromatography and mass spectrometry: (1) acquiring and processing original chromatographic and mass spectrometric data of known TCM samples; (2) generating high-dimensional data and high-dimensional images of known TCM samples; (3) establishing high-dimensional chromatographic and mass spectrometric data of known TCM samples. Image database; _Application of high-dimensional image database of TCM Chromatography and Mass Spectrometry: 1) Acquisition of image data of unknown samples; 2) Recognition of unknown samples; 3) Verification of recognition results of unknown samples. Compared with the prior art, the present invention utilizes the high-dimensional image technology of TCM chromatography and mass spectrometry to fully characterize the spatial information of a large number of compounds in TCM samples, and utilizes the above spatial information to realize matching and recognition between unknown samples and known samples of TCM.
【技术实现步骤摘要】
一种基于中药色谱-质谱高维图像数据库的中药识别方法
本专利技术涉及一种基于中药色谱-质谱高维图像数据库的中药识别方法。
技术介绍
复杂样本的化合物构成具有极端复杂性。中药即是典型的复杂样本,所含成分极其复杂,结构多样、种类繁多,常见类型包括酚类、生物碱类、皂苷类、萜类、黄酮类、内酯类、蒽酮类、有机酸类以及鞣质类等,单一中药即包含数百上千计的次生代谢产物和小分子成分,由多种中药组合的中药复方制剂的成分则更多。相应地,复杂样本中蕴含海量信息。如中药化合物之间的相互关系、不同中药的药性药效差异、同属药材化学成分异同及产地、年份、生长环境对药材质量的影响等科学问题都蕴藏其中。目前对于复杂样本的研究面临两个重要瓶颈:一方面,研究大多采用碎片化、点状的低维数据,如色谱保留时间,m/z值,子离子碎片信息等,这些低维数据忽视也无法体现上述大量化学成分间的关联。高维数据恰是海量信息的有力载体。与低维数据相比,高维数据能够有效地表示样本中各数据点的空间信息从而反映它们的空间关系。因此,获取复杂样本化合物的高维数据才能真正实现从复杂样本中得到、处理、挖掘那些高价值信息。另一方面,实验产生的数据资源庞大却零散,相关研究产生的数据不能整合利用,导致科研工作中人力、物力、时间等投入的成本高,产出却不显著。数据库技术是一种计算机辅助管理、整合数据的方法。将高维数据与数据库技术结合建立高维数据数据库正是解决上述难题的方向。高维数据的获取需要联用仪器来实现。色谱-质谱联用技术将应用范围极广的分离方法-色谱法与灵敏、专属、能提供分子量和结构信息的质谱法结合起来,显然是复杂样本高维数据获取的理 ...
【技术保护点】
1.一种基于中药色谱‑质谱高维图像数据库的中药识别方法,其特征在于,按以下步骤进行:⑴中药色谱‑质谱高维图像数据库的建立:1)获取与处理已知中药样本的原始色谱‑质谱(X‑MS)数据:使用色谱和质谱获取已知中药样本的原始X‑MS数据,将已知中药样本原始X‑MS数据导入峰提取软件(如Progenesis QI)中对色谱‑质谱联用原始X‑MS数据进行数据处理;2)生成已知中药样本的高维数据与图像:获取样本中每个化合物的m/z、t、I、m、z值,产生高维数据矩阵(如m/z‑t‑I矩阵、m‑z‑t‑I矩阵或m‑t‑I矩阵),生成已知中药样本色谱‑质谱联用高维数据;将高维数据导入图像生成软件(如Matlab等)生成图像,使高维数据中的每个离子与构成图像中的点一一对应,每个点拥有自己的坐标信息(t,m/z或m或m与z),每个点的强度由点的大小或/和亮度的强弱表示,高维数据图像中的点与高维数据一一对应;3)建立已知中药样本的色谱‑质谱高维图像数据库:将获得的1类或2类以上已知中药样本高维数据图像作为中药色谱‑质谱高维图像数据库(类别数≧1),每类已知中药样本中的样本数为1个或2个以上;中药色谱‑质谱高 ...
【技术特征摘要】
1.一种基于中药色谱-质谱高维图像数据库的中药识别方法,其特征在于,按以下步骤进行:⑴中药色谱-质谱高维图像数据库的建立:1)获取与处理已知中药样本的原始色谱-质谱(X-MS)数据:使用色谱和质谱获取已知中药样本的原始X-MS数据,将已知中药样本原始X-MS数据导入峰提取软件(如ProgenesisQI)中对色谱-质谱联用原始X-MS数据进行数据处理;2)生成已知中药样本的高维数据与图像:获取样本中每个化合物的m/z、t、I、m、z值,产生高维数据矩阵(如m/z-t-I矩阵、m-z-t-I矩阵或m-t-I矩阵),生成已知中药样本色谱-质谱联用高维数据;将高维数据导入图像生成软件(如Matlab等)生成图像,使高维数据中的每个离子与构成图像中的点一一对应,每个点拥有自己的坐标信息(t,m/z或m或m与z),每个点的强度由点的大小或/和亮度的强弱表示,高维数据图像中的点与高维数据一一对应;3)建立已知中药样本的色谱-质谱高维图像数据库:将获得的1类或2类以上已知中药样本高维数据图像作为中药色谱-质谱高维图像数据库(类别数≧1),每类已知中药样本中的样本数为1个或2个以上;中药色谱-质谱高维图像数据库,包括已知中药样本的样本信息、原始X-MS数据信息、高维数据信息、高维图像数据信息;数据库类型包括但不局限于文件夹数据集、网页数据库、基于商业化工作站或基于用户自研发工作站的数据库;⑵中药色谱-质谱高维图像数据库的应用:1)未知样本高维图像数据的获取:采用与步骤(1)相同或相似的操作参数和条件,按步骤(1)中1)~2)操作,针对待分析的未知样本进行分析,获取未知样本原始X-MS数据和高维数据;利用图像生成软件将X-MS数据得到未知样本的X-MS高维图像;2)未知样本的识别;A、利用机器学习中的图像分割工具(如Matlab2016b自带的分割程序),或聚类工具(如K-Means,DBSCAN或Fanny等中的一种或二种以上),将未知样本X-MS高维图像中的点分割为n个点簇(n≥1整数);点簇指的是在空间上距离接近的点的集合,点簇内点的个数n≥3;每个点簇可以有自己的中心点,点簇的形状可以为任意形状;B、将提取点簇后的未知样本X-MS高维图像与中药X-MS高维图像数据库中的已知中药样本X-MS高维图像逐一进行分别扫描和匹配;扫描时,将两个X-MS高维图像的原点、t轴和m/z(m)轴对齐;扫描时,点簇作为一个整体,移动的范围为0-Tk,Tk为已知中药样本对应的最大分析时间;扫描时,未知样本的每个点簇保留其m/z(或m)轴的位置和几何形状,沿时间轴(t)进行扫描;通过扫描,寻找未知样本点簇与已知中药样本X-MS高维图像中能够在t和m/z(或m)能够准确匹配的共同点;扫描过程中,在未知样本中的一个点簇中的点与已知中药样本X-MS高维图像中的点进行匹配时,每个点允许的t绝对偏移值(ttolerance)为≥T,T等于未知样本X-MS数据采集时色谱仪允许的保留时间平均偏移值(绝对值,可用1个或1个以上标准物质,或某样本中的1个或1个以上化合物的多次重复测定计算)与已知中药样本X-MS数据采集时色谱仪允许的保留时间平均偏移值(绝对值,可用1个或1个以上标准物质,或某样本中的1个或1个以上化合物的多次重复测定计算)之和;扫描过程中,在未知样本中的一个点簇中的点与已知中药样本X-MS高维图像中的点进行匹配时,每个点允许的m/z(或m)绝对测定误差[m/z(或m)tolerance]≥A,A等于未知和已知中药样本X-MS数据采集时质谱仪扫描时允许的质量平均偏差(绝对值,可由仪器所用的校正液多次重复测定)之和;当未知样本点簇内一个点与已知中药样本的某个点满足t偏差和m/z(或m)偏差时,认为该点符合匹配要求;扫描时,点簇沿时间轴(t)扫描的步长≤T,通常情况下,0s<T<10000s;C、当一个点簇移动到已知中药样本X-MS高维图像的t轴的每一个位置时,记录匹配点的个数、每个匹配点的坐标和点簇几何中心点的坐标;D、计算每一个位置时,未知样本一个点簇(i,i≥1整数)与该已知中药样本X-MS高维图像之间的匹配度(Si),匹配度的大小可利用统计工具(如Matlab)计算点簇(i)与已知中药样本X-MS高维图像所匹配的点数、或相似度(如图像相似度计算中的欧氏距离法)、或相关度(如Matlab中的2D-correlationcoefficient)中的一种或二种以上进行计算;由上述三种方法得到的匹配度分别由点数(或点数的函数)、相似度或相关度表示;点簇匹配度大小,与点簇匹配的点数、坐标位置(t,m/z)和强度这四个变量呈线性或非线性相关;计算点数(或点数的函数)、相似度或相关度的基础是基于四个变量的关系变换;可选用不同的匹配度计算方法分别计算点簇和已知中药样本X-MS高维图像的整体匹配度;匹配点的个数指的是点簇符合匹配条件点的个数;基于上述步骤,对未知样本X-MS高维图像中每一个点簇的最大匹配度(Si)进行数学加权处理(如加和、平均或取对数),得到未知样本X-MS高维图像与已知中药样本X-MS高维图像的整体匹配度(Sc);E、重复上述步骤,逐一分析未知样本X-MS高维图像与其它已知中药样本X-MS高维图像之间的匹配度,得到其与每一个已知中药样本的整体匹配度(Sc);F、未知样本的所属类别可不借助阈值或借助阈值进行判定;当不借助阈值时,利用上述步骤,将未知样本与已知中药样本进行匹配,对匹配度从大到小进行排序,若未知样本与某一已知中药样本的匹配度排名越靠前,表明未知样本为与该样本的可能性越大,反之越小;或,当借助阈值时,设定阈值γ,用于判断不同来源未知样本与同类已知中药样本匹配的可信范围;阈值可以根据统计学的方法设定:按步骤(1)中1)~2)操作,采用相同或相近的操作参数和条件,选取2个以上同类别的已知中药样本作为某一类样本的训练样本,进行分析,获取X-MS原始数据;利用图像生成软件(如Matlab2016b)将X-MS原始数据或多维信息文本转化为X-MS高维图像,得到该类样本的训练X-MS高维图像集;利用训练X-MS高维图像集,与同类已知中药样本X-MS高维图像进行匹配,通过统计学的方法(如概率,比率等)发现匹配度分布区间,选定分布区间中匹配度的下限作为该类样本的阈值γ;或,阈值可以利用文献报道或实验观察得到某类样本(n≥2)与已知中药样本匹配度分布区间(按步骤(1)中1)~2)操作采用相同或相近的操作参数和条件所得到的分析结果),选定分布区间中匹配度的下限作为该类样...
【专利技术属性】
技术研发人员:张晓哲,赵楠,程孟春,
申请(专利权)人:中国科学院大连化学物理研究所,
类型:发明
国别省市:辽宁,21
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。