CN103913476A - 一种获取植物类药材粉末衍射文件的方法 - Google Patents
一种获取植物类药材粉末衍射文件的方法 Download PDFInfo
- Publication number
- CN103913476A CN103913476A CN201410105727.6A CN201410105727A CN103913476A CN 103913476 A CN103913476 A CN 103913476A CN 201410105727 A CN201410105727 A CN 201410105727A CN 103913476 A CN103913476 A CN 103913476A
- Authority
- CN
- China
- Prior art keywords
- medicinal material
- file
- plant medicinal
- derivative
- sample
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Landscapes
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本发明公开了一种获取植物类药材粉末衍射文件的方法,先获取植物类药材的X射线衍射图谱文件,再将其转化为仅含衍射角度和相应衍射强度的TXT格式的文本文件,然后采用Savitzky-Golay平滑滤波器对图谱信号进行去噪平滑,再进行B样条插值拟合,对所得拟合曲线函数求二阶导数,获得二阶导数指纹图谱,并采用对称零面积变换法进行寻峰处理,获得所寻衍射峰角度及其衍射强度,再根据布拉格公式计算晶面间距,最后依据国际标准粉末衍射文件格式制作该植物类药材的粉末衍射文件。本发明方法借助信号处理技术,能充分挖掘植物类药材的图谱指纹信息,解决了植物类药材X射线衍射图谱的解析问题,拓展了XRD的应用。
Description
技术领域
本发明属于药材X射线衍射图谱分析技术领域,具体涉及一种获得植物类药材粉末衍射文件的方法。
背景技术
近年来,药材分析领域引入X射线衍射指纹图谱,用于表征和鉴定药材,取得了不错的效果,但也仅限于结晶度好的矿物类和化石类药材。植物类药材由于含有大量有机成分,其XRD指纹图谱均呈弥散性宽峰,利用传统的XRD定性定量分析软件处理,无明显特征峰或根本无峰可寻,无法对其进行准确表征,得不到其粉末衍射文件(PDF卡片),限制了XRD指纹图谱的应用。因此借助数据处理技术,充分挖掘图谱指纹信息,能够拓展XRD的应用。
发明内容
本发明所要解决的技术问题在于克服上述植物类药材弥散性XRD图谱无明显特征峰或根本无峰可寻的缺点,提供一种获取植物类药材粉末衍射文件的方法。
解决上述技术问题所采用的技术方案由下述步骤组成:
1、获取植物类药材的X射线衍射图谱文件
将植物类药材烘干、粉碎、过200目筛,得到植物类药材粉末样品,将样品压片,用X射线衍射仪测试,得到样品的X射线衍射图谱文件。
2、获取植物类药材的二阶导数指纹图谱
将样品的X射线衍射图谱文件转化为仅含衍射角度和相应衍射强度的TXT格式的文本文件,采用Savitzky-Golay平滑滤波器对TXT格式的文本文件对应的图谱信号进行去噪平滑,对平滑后的图谱信号进行B样条插值拟合,并对所得拟合曲线函数求二阶导数,然后对负值部分取绝对值、正值置零,获得样品的二阶导数指纹图谱。
3、寻峰处理
采用对称零面积变换法对样品的二阶导数指纹图谱进行寻峰处理,取阈值因子为0.1~1,获得所寻衍射峰角度及其衍射强度。
4、制作植物类药材的粉末衍射文件
根据步骤3得到的衍射峰角度和布拉格公式2dsinθ=nλ计算相应的晶面间距,依据国际标准粉末衍射文件格式制作该植物类药材的粉末衍射文件。
本发明的获取植物类药材的二阶导数指纹图谱步骤2中,进一步优选将样品的X射线衍射图谱文件转化为仅含衍射角度和相应衍射强度的TXT格式的文本文件,采用基于Savitzky-Golay平滑滤波器的MATLAB内置函数smooth(signal,span,‘sgolay’,degree)对该TXT格式的文本文件对应的图谱信号进行去噪平滑,其中degree取2~11的整数,span取3~45的整数,然后采用MATLAB内置函数spapi(span_B,Theta,Intensity)对平滑后的图谱信号进行B样条插值拟合,其中span_B取2~11的整数,并对所得拟合曲线函数求二阶导数,然后对负值部分取绝对值、正值置零,获得样品的二阶导数指纹图谱。
本发明的寻峰处理步骤3中,采用对称零面积变换法对样品的二阶导数指纹图谱进行寻峰处理,进一步优选取阈值因子为0.2,获得所寻衍射峰角度及其衍射强度。
本发明借助信号处理技术,引入Savitzky-Golay平滑滤波器对XRD图谱进行平滑去噪,以B样条插值拟合计算二阶导数指纹图谱,选择对称零面积变换法进行寻峰处理,能充分挖掘图谱指纹信息,解决了植物类药材X射线衍射图谱的解析问题,拓展了XRD的应用。
附图说明
图1是实施例1获得的药材五味子的MDI格式X射线衍射图谱。
图2是实施例1获得的五味子的平滑后图谱和二阶导数指纹图谱。
图3是实施例1获得的五味子的二阶导数指纹图谱的寻峰结果。
图4是实施例1获得的五味子的PDF卡片。
图5是实施例2获得的黄秦花的PDF卡片。
具体实施方式
下面结合附图和实施例对本发明进一步详细说明,但本发明的保护范围不限于这些实施例。
实施例1
1、获取药材五味子的X射线衍射图谱
将药材五味子(Schisandra chinensls(Turcz)Baill)置于烘干箱内70℃烘干,然后用粉碎机反复粉碎,过200目筛,得到五味子粉末样品,将五味子粉末样品压片,用Y-4Q型X射线衍射仪(由辽宁丹东仪器设备厂提供)进行测试,测试条件为:步宽0.03度,管压35kV,管流20mA,铜靶的Kα辐射,滤波片为镍,连续扫描,时间常数0.5秒,扫描范围10°~80°,得到五味子的MDI格式X射线衍射图谱,见图1。
2、获取药材五味子的二阶导数指纹图谱
去除五味子的MDI格式X射线衍射图谱文件的头信息,将其转化为仅包含衍射角度和衍射强度的TXT格式的文本文件,采用基于Savitzky-Golay平滑滤波器的MATLAB内置函数smooth(signal,span,‘sgolay’,degree)对该TXT格式的文本文件对应的图谱信号进行去噪平滑,其中signal指图谱信号,span=5,指平滑窗口为5个数据点,‘sgolay’指采用Savitzky-Golay平滑滤波器,degree=2,指以二阶多项式进行插值拟合,然后采用MATLAB内置函数spapi(span_B,Theta,Intensity)对平滑后的图谱信号进行B样条插值拟合,其中span_B=5,指样条曲线阶数为5,Theta指衍射角度,Intensity指衍射强度,并对所得拟合曲线函数求二阶导数,然后对所得二阶导数的负值部分取绝对值、正值置零,获得五味子的二阶导数指纹图谱,见图2。
3、寻峰处理
采用对称零面积变换法对五味子的二阶导数指纹图谱进行寻峰处理,具体寻峰算法参照毕云峰等公开的《LIBS/Raman光谱对称零面积变化自动寻峰方法研究》(光谱学与光谱分析,Vol.33,No.22,pp438-443,February,2013)中的高斯线型函数作为对称零面积变换函数,其中高斯线型半高全宽为4、变换函数的窗口宽度为25,取阈值因子为0.2,获得所寻衍射峰角度及其衍射强度,结果见图3。
4、制作五味子的粉末衍射文件
根据布拉格公式2dsinθ=nλ,利用步骤3中所得衍射峰角度信息,计算相应的晶面间距d,然后收集五味子中文名称、拉丁文名称、产地和生存环境等信息,依据国际标准粉末衍射文件格式制作药材五味子的粉末衍射文件,即五味子的PDF卡片。由于步骤3所寻衍射峰较多,这里只对衍射强度较大的前25个峰进行录入,见图4。
实施例2
按照实施例1中方法制作药材五味子用药材黄秦花(Gentiana macrophylla)的粉 末衍射文件,其中步骤1的扫描范围为3°~90°,其他条件和步骤与实施例1相同,得到药材黄秦花的PDF卡片,结果见图5。
为了证明本发明的有益效果,采用传统的XRD分析方法分别对药材五味子和黄秦花进行处理,获得五味子和黄秦花的原始X射线衍射图谱,然后利用JADE6.5图谱分析软件对其进行寻峰处理,平滑参数和扣除背景选择默认值,试验结果见表1。
表1传统XRD分析方法对五味子和黄秦花的寻峰结果
由表1可以看出,五味子没有寻到一个峰,根本无法对药材五味子进行表征;黄秦花仅寻到一个衍射峰,不足以表征药材黄秦花图谱丰富的指纹信息;而从本发明的实施例1和2的寻峰情况可以看出,其衍射峰较多,峰位明确,挖掘出了药材图谱本身所蕴含的大量指纹特征。对比看来,本发明能够有效地解决植物类药材弥散模糊性XRD图谱的解析问题,给出足够丰富的特征峰信息,以制作药材相应的PDF卡片,为药材真伪鉴别和质量控制提供了有价值的数据。
以上列举的仅是本发明的具体实施例子,本领域技术人员在本发明公开内容的基础上可以想到的任何变形都在本发明的保护范围内。
Claims (3)
1.一种获取植物类药材粉末衍射文件的方法,其特征在于它由下述步骤组成:
(1)获取植物类药材的X射线衍射图谱文件
将植物类药材烘干、粉碎、过200目筛,得到植物类药材粉末样品,将样品压片,用X射线衍射仪测试,得到样品的X射线衍射图谱文件;
(2)获取植物类药材的二阶导数指纹图谱
将样品的X射线衍射图谱文件转化为仅含衍射角度和相应衍射强度的TXT格式的文本文件,采用Savitzky-Golay平滑滤波器对该TXT格式的文本文件对应的图谱信号进行去噪平滑,然后对平滑后的图谱信号进行B样条插值拟合,并对所得拟合曲线函数求二阶导数,最后对负值部分取绝对值、正值置零,获得样品的二阶导数指纹图谱;
(3)寻峰处理
采用对称零面积变换法对样品的二阶导数指纹图谱进行寻峰处理,取阈值因子为0.1~1,获得所寻衍射峰角度及其衍射强度;
(4)制作植物类药材的粉末衍射文件
根据步骤(3)得到的衍射峰角度和布拉格公式2dsinθ=nλ计算相应的晶面间距,依据国际标准粉末衍射文件格式制作该植物类药材的粉末衍射文件。
2.根据权利要求1所述的获取植物类药材粉末衍射文件的方法,其特征在于:所述的获取植物类药材的二阶导数指纹图谱步骤(2)中,将样品的X射线衍射图谱文件转化为仅含衍射角度和相应衍射强度的TXT格式的文本文件,采用基于Savitzky-Golay平滑滤波器的MATLAB内置函数smooth(signal,span,‘sgolay’,degree)对该TXT格式的文本文件对应的图谱信号进行去噪平滑,其中degree取2~11的整数,span取3~45的整数,然后采用MATLAB内置函数spapi(span_B,Theta,Intensity)对平滑后的图谱信号进行B样条插值拟合,其中span_B取2~11的整数,并对所得拟合曲线函数求二阶导数,然后对负值部分取绝对值、正值置零,获得样品的二阶导数指纹图谱。
3.根据权利要求1所述的获取植物类药材粉末衍射文件的方法,其特征在于:所述的寻峰处理步骤(3)中,采用对称零面积变换法对样品的二阶导数指纹图谱进行寻峰处理,取阈值因子为0.2,获得所寻衍射峰角度及其衍射强度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410105727.6A CN103913476B (zh) | 2014-03-20 | 2014-03-20 | 一种获取植物类药材粉末衍射文件的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410105727.6A CN103913476B (zh) | 2014-03-20 | 2014-03-20 | 一种获取植物类药材粉末衍射文件的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103913476A true CN103913476A (zh) | 2014-07-09 |
CN103913476B CN103913476B (zh) | 2016-08-17 |
Family
ID=51039316
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410105727.6A Expired - Fee Related CN103913476B (zh) | 2014-03-20 | 2014-03-20 | 一种获取植物类药材粉末衍射文件的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103913476B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106932418A (zh) * | 2017-03-20 | 2017-07-07 | 西南民族大学 | 七叶一枝花或云南重楼的鉴别方法 |
CN111398324A (zh) * | 2020-03-31 | 2020-07-10 | 深圳晶泰科技有限公司 | 一种粉末x射线衍射图谱的衍射峰定标方法 |
WO2021103368A1 (zh) * | 2020-03-31 | 2021-06-03 | 深圳晶泰科技有限公司 | 一种粉末x射线衍射图谱的衍射峰定标方法 |
CN113887658A (zh) * | 2021-10-22 | 2022-01-04 | 深圳晶泰科技有限公司 | 晶体粉末衍射图谱的聚类方法、装置及存储介质 |
CN113945592A (zh) * | 2021-10-14 | 2022-01-18 | 曲靖师范学院 | 一种单晶xrd精修方法及其平台 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101281147A (zh) * | 2007-04-13 | 2008-10-08 | 江苏天瑞信息技术有限公司 | 用于x荧光光谱仪的数字化信噪比增强处理方法 |
CN101788507A (zh) * | 2010-02-03 | 2010-07-28 | 北京矿冶研究总院 | 一种光谱分析方法 |
JP2013145188A (ja) * | 2012-01-16 | 2013-07-25 | Sumitomo Metal Mining Co Ltd | 粒子存在比率算出方法及び粒子結晶サイズ算出方法 |
-
2014
- 2014-03-20 CN CN201410105727.6A patent/CN103913476B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101281147A (zh) * | 2007-04-13 | 2008-10-08 | 江苏天瑞信息技术有限公司 | 用于x荧光光谱仪的数字化信噪比增强处理方法 |
CN101788507A (zh) * | 2010-02-03 | 2010-07-28 | 北京矿冶研究总院 | 一种光谱分析方法 |
JP2013145188A (ja) * | 2012-01-16 | 2013-07-25 | Sumitomo Metal Mining Co Ltd | 粒子存在比率算出方法及び粒子結晶サイズ算出方法 |
Non-Patent Citations (2)
Title |
---|
杨建华: "中药XRD二阶指纹图谱方法学研究", 《中国优秀博硕士学位论文全文数据库 (硕士) 医药卫生科技辑》 * |
谢锟、胡成西: "基于样条法的米兰X射线衍射二阶导数指纹图谱", 《光谱实验室》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106932418A (zh) * | 2017-03-20 | 2017-07-07 | 西南民族大学 | 七叶一枝花或云南重楼的鉴别方法 |
CN111398324A (zh) * | 2020-03-31 | 2020-07-10 | 深圳晶泰科技有限公司 | 一种粉末x射线衍射图谱的衍射峰定标方法 |
WO2021103368A1 (zh) * | 2020-03-31 | 2021-06-03 | 深圳晶泰科技有限公司 | 一种粉末x射线衍射图谱的衍射峰定标方法 |
CN113945592A (zh) * | 2021-10-14 | 2022-01-18 | 曲靖师范学院 | 一种单晶xrd精修方法及其平台 |
CN113945592B (zh) * | 2021-10-14 | 2024-07-19 | 曲靖师范学院 | 一种单晶xrd精修方法及其平台 |
CN113887658A (zh) * | 2021-10-22 | 2022-01-04 | 深圳晶泰科技有限公司 | 晶体粉末衍射图谱的聚类方法、装置及存储介质 |
WO2023065515A1 (zh) * | 2021-10-22 | 2023-04-27 | 深圳晶泰科技有限公司 | 晶体粉末衍射图谱的聚类方法、装置及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN103913476B (zh) | 2016-08-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103913476A (zh) | 一种获取植物类药材粉末衍射文件的方法 | |
Luo et al. | ChinaCropPhen1km: a high-resolution crop phenological dataset for three staple crops in China during 2000–2015 based on leaf area index (LAI) products | |
Zhang et al. | Quantitative analysis of routine chemical constituents in tobacco by near-infrared spectroscopy and support vector machine | |
Shi et al. | Development of a national VNIR soil-spectral library for soil classification and prediction of organic matter concentrations | |
Lao et al. | Predicting the contents of soil salt and major water-soluble ions with fractional-order derivative spectral indices and variable selection | |
Fan et al. | Fourier Transform Infrared (FT-IR) Spectroscopy for discrimination of Rhizoma gastrodiae (Tianma) from different producing areas | |
Guo et al. | Evaluating calibration and spectral variable selection methods for predicting three soil nutrients using Vis-NIR spectroscopy | |
Lee et al. | Discrimination and prediction of cultivation age and parts of Panax ginseng by Fourier-transform infrared spectroscopy combined with multivariate statistical analysis | |
CN101916440B (zh) | 基于数字图像形态学理论的高光谱异常检测方法 | |
Lavanya et al. | Terrain mapping of LandSat8 images using MNF and classifying soil properties using ensemble modelling | |
Niu et al. | A 30 m annual maize phenology dataset from 1985 to 2020 in China | |
Greig et al. | Spectral-temporal modelling of bamboo-dominated forest succession in the Atlantic Forest of Southern Brazil | |
CN105468922A (zh) | 一种现场电磁环境信号实时分析与威胁评估方法 | |
CN105629290A (zh) | 一种数字核脉冲信号墨西哥草帽小波成形方法 | |
CN106568759A (zh) | 一种基于多谱融合的道地药材质量鉴别方法 | |
Liu et al. | Improving potato above ground biomass estimation combining hyperspectral data and harmonic decomposition techniques | |
CN108090462A (zh) | 一种基于盒维数的辐射源指纹特征提取方法 | |
Xiao et al. | Inversion study of cadmium content in soil based on reflection spectroscopy and MSC-ELM model | |
Zhao | Vegetation and climate reconstructions on different time scales in China: a review of Chinese palynological research | |
CN110243885A (zh) | 一种基于时频特征的低产气油井井口含水率预测方法 | |
Jiang et al. | Remotely sensed estimation of total iron content in soil with harmonic analysis and BP neural network | |
Zhou et al. | Spectral response feature bands extracted from near standard soil samples for estimating soil Pb in a mining area | |
CN102135496A (zh) | 基于多尺度回归的红外光谱定量分析方法和装置 | |
CN104076003B (zh) | 一种矿物光谱吸收特征参数的提取方法 | |
CN103500421B (zh) | 一种基于频率特征的耕地复种指数提取方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160817 Termination date: 20170320 |