CN111707636A - 结合傅里叶变换红外光谱和化学计量学分析鉴别生前死后骨折的方法 - Google Patents
结合傅里叶变换红外光谱和化学计量学分析鉴别生前死后骨折的方法 Download PDFInfo
- Publication number
- CN111707636A CN111707636A CN202010442559.5A CN202010442559A CN111707636A CN 111707636 A CN111707636 A CN 111707636A CN 202010442559 A CN202010442559 A CN 202010442559A CN 111707636 A CN111707636 A CN 111707636A
- Authority
- CN
- China
- Prior art keywords
- fracture
- mortem
- fourier transform
- post
- identifying
- 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
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000005033 Fourier transform infrared spectroscopy Methods 0.000 title claims abstract description 20
- 238000004458 analytical method Methods 0.000 title claims abstract description 16
- 230000009466 transformation Effects 0.000 claims abstract description 7
- 210000000988 bone and bone Anatomy 0.000 claims description 26
- 241000283973 Oryctolagus cuniculus Species 0.000 claims description 14
- 238000010521 absorption reaction Methods 0.000 claims description 11
- 238000006731 degradation reaction Methods 0.000 claims description 10
- 238000001157 Fourier transform infrared spectrum Methods 0.000 claims description 9
- 150000001408 amides Chemical class 0.000 claims description 9
- 238000002329 infrared spectrum Methods 0.000 claims description 9
- 239000000843 powder Substances 0.000 claims description 9
- 238000010171 animal model Methods 0.000 claims description 5
- 238000012935 Averaging Methods 0.000 claims description 2
- 238000000227 grinding Methods 0.000 claims description 2
- 239000000126 substance Substances 0.000 abstract description 9
- 230000003595 spectral effect Effects 0.000 abstract description 8
- 238000012545 processing Methods 0.000 abstract description 3
- 206010017076 Fracture Diseases 0.000 description 103
- 208000010392 Bone Fractures Diseases 0.000 description 89
- 241000282414 Homo sapiens Species 0.000 description 16
- 238000001228 spectrum Methods 0.000 description 14
- 230000008859 change Effects 0.000 description 9
- 238000012549 training Methods 0.000 description 9
- 238000010239 partial least squares discriminant analysis Methods 0.000 description 7
- 230000015556 catabolic process Effects 0.000 description 6
- 238000001514 detection method Methods 0.000 description 5
- 210000002303 tibia Anatomy 0.000 description 5
- 230000002146 bilateral effect Effects 0.000 description 4
- 238000007781 pre-processing Methods 0.000 description 4
- 206010002199 Anaphylactic shock Diseases 0.000 description 3
- 208000000202 Diffuse Axonal Injury Diseases 0.000 description 3
- 208000028979 Skull fracture Diseases 0.000 description 3
- 208000003455 anaphylaxis Diseases 0.000 description 3
- 230000009521 diffuse axonal injury Effects 0.000 description 3
- 230000008439 repair process Effects 0.000 description 3
- 102000013455 Amyloid beta-Peptides Human genes 0.000 description 2
- 108010090849 Amyloid beta-Peptides Proteins 0.000 description 2
- 238000004566 IR spectroscopy Methods 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000013145 classification model Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 230000035876 healing Effects 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 150000002894 organic compounds Chemical class 0.000 description 2
- IOLCXVTUBQKXJR-UHFFFAOYSA-M potassium bromide Chemical compound [K+].[Br-] IOLCXVTUBQKXJR-UHFFFAOYSA-M 0.000 description 2
- 238000002360 preparation method Methods 0.000 description 2
- 238000000513 principal component analysis Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 235000004252 protein component Nutrition 0.000 description 2
- 108090000623 proteins and genes Proteins 0.000 description 2
- 102000004169 proteins and genes Human genes 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 210000001519 tissue Anatomy 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- ZFXQNADNEBRERM-BJDJZHNGSA-N Ala-Ala-Pro-Pro Chemical compound C[C@H](N)C(=O)N[C@@H](C)C(=O)N1CCC[C@H]1C(=O)N1[C@H](C(O)=O)CCC1 ZFXQNADNEBRERM-BJDJZHNGSA-N 0.000 description 1
- 208000002109 Argyria Diseases 0.000 description 1
- 206010003497 Asphyxia Diseases 0.000 description 1
- 102000019034 Chemokines Human genes 0.000 description 1
- 108010012236 Chemokines Proteins 0.000 description 1
- 102000008186 Collagen Human genes 0.000 description 1
- 108010035532 Collagen Proteins 0.000 description 1
- 206010018852 Haematoma Diseases 0.000 description 1
- 208000032843 Hemorrhage Diseases 0.000 description 1
- 206010020649 Hyperkeratosis Diseases 0.000 description 1
- OUYCCCASQSFEME-QMMMGPOBSA-N L-tyrosine Chemical compound OC(=O)[C@@H](N)CC1=CC=C(O)C=C1 OUYCCCASQSFEME-QMMMGPOBSA-N 0.000 description 1
- 241000124008 Mammalia Species 0.000 description 1
- 241001465754 Metazoa Species 0.000 description 1
- 206010037423 Pulmonary oedema Diseases 0.000 description 1
- 206010049418 Sudden Cardiac Death Diseases 0.000 description 1
- 208000004367 Tibial Fractures Diseases 0.000 description 1
- 208000027418 Wounds and injury Diseases 0.000 description 1
- DZHSAHHDTRWUTF-SIQRNXPUSA-N amyloid-beta polypeptide 42 Chemical compound C([C@@H](C(=O)N[C@@H](C)C(=O)N[C@@H](CCC(O)=O)C(=O)N[C@@H](CC(O)=O)C(=O)N[C@H](C(=O)NCC(=O)N[C@@H](CO)C(=O)N[C@@H](CC(N)=O)C(=O)N[C@@H](CCCCN)C(=O)NCC(=O)N[C@@H](C)C(=O)N[C@H](C(=O)N[C@@H]([C@@H](C)CC)C(=O)NCC(=O)N[C@@H](CC(C)C)C(=O)N[C@@H](CCSC)C(=O)N[C@@H](C(C)C)C(=O)NCC(=O)NCC(=O)N[C@@H](C(C)C)C(=O)N[C@@H](C(C)C)C(=O)N[C@@H]([C@@H](C)CC)C(=O)N[C@@H](C)C(O)=O)[C@@H](C)CC)C(C)C)NC(=O)[C@H](CC=1C=CC=CC=1)NC(=O)[C@@H](NC(=O)[C@H](CC(C)C)NC(=O)[C@H](CCCCN)NC(=O)[C@H](CCC(N)=O)NC(=O)[C@H](CC=1N=CNC=1)NC(=O)[C@H](CC=1N=CNC=1)NC(=O)[C@@H](NC(=O)[C@H](CCC(O)=O)NC(=O)[C@H](CC=1C=CC(O)=CC=1)NC(=O)CNC(=O)[C@H](CO)NC(=O)[C@H](CC(O)=O)NC(=O)[C@H](CC=1N=CNC=1)NC(=O)[C@H](CCCNC(N)=N)NC(=O)[C@H](CC=1C=CC=CC=1)NC(=O)[C@H](CCC(O)=O)NC(=O)[C@H](C)NC(=O)[C@@H](N)CC(O)=O)C(C)C)C(C)C)C1=CC=CC=C1 DZHSAHHDTRWUTF-SIQRNXPUSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 208000034158 bleeding Diseases 0.000 description 1
- 230000000740 bleeding effect Effects 0.000 description 1
- 208000029028 brain injury Diseases 0.000 description 1
- 210000005013 brain tissue Anatomy 0.000 description 1
- 238000012569 chemometric method Methods 0.000 description 1
- 229920001436 collagen Polymers 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000002790 cross-validation Methods 0.000 description 1
- 230000006378 damage Effects 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000000228 dispersive infrared spectroscopy Methods 0.000 description 1
- 238000001035 drying Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 125000000524 functional group Chemical group 0.000 description 1
- 239000003102 growth factor Substances 0.000 description 1
- 238000007490 hematoxylin and eosin (H&E) staining Methods 0.000 description 1
- 239000001257 hydrogen Substances 0.000 description 1
- 229910052739 hydrogen Inorganic materials 0.000 description 1
- 238000011532 immunohistochemical staining Methods 0.000 description 1
- 238000000338 in vitro Methods 0.000 description 1
- 230000002757 inflammatory effect Effects 0.000 description 1
- 208000014674 injury Diseases 0.000 description 1
- 238000012067 mathematical method Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000001575 pathological effect Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 235000018102 proteins Nutrition 0.000 description 1
- 238000011552 rat model Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 210000003625 skull Anatomy 0.000 description 1
- 238000004611 spectroscopical analysis Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 206010043827 tibia fracture Diseases 0.000 description 1
- 238000011282 treatment Methods 0.000 description 1
- OUYCCCASQSFEME-UHFFFAOYSA-N tyrosine Natural products OC(=O)C(N)CC1=CC=C(O)C=C1 OUYCCCASQSFEME-UHFFFAOYSA-N 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
- G01N21/3563—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light for analysing solids; Preparation of samples therefor
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/70—Machine learning, data mining or chemometrics
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
- G01N2021/3595—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light using FTIR
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Chemical & Material Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Crystallography & Structural Chemistry (AREA)
- Software Systems (AREA)
- Artificial Intelligence (AREA)
- Theoretical Computer Science (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Immunology (AREA)
- Computing Systems (AREA)
- Medical Informatics (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Databases & Information Systems (AREA)
- Data Mining & Analysis (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Pathology (AREA)
- Evolutionary Computation (AREA)
Abstract
本发明公开了一种结合傅里叶变换红外光谱和化学计量学分析鉴别生前死后骨折的方法：利用傅里叶变换红外光谱采集骨折断面样本的光谱信息，利用化学计量学对选定的光谱信息进行处理，根据均值中心化和标准正态变量变换处理后的光谱数据，利用偏最小二乘法判别分析建立生前死后骨折鉴别模型，利用该模型客观、快速、准确的鉴别生前死后骨折。
Description
技术领域
本发明涉及刑事科学技术和法医领域的骨折时间判断，具体涉及结合傅里叶变换红外光谱(FTIR)技术和化学计量学分析鉴别生前死后骨折的方法。
背景技术
法医领域通常将骨折时间分为生前骨折、濒死期骨折和死后骨折。生前骨折是指确定死亡前发生的骨折。濒死期骨折是死亡前后发生的损伤。死后骨折是指死后发生的骨折。骨折时间的精确估计对于明确骨折的机制、死亡原因、案件性质、侦查方向有着重大的意义。目前对于生前骨折、濒死期骨折和死后骨折之间的鉴别主要采用观察宏观骨折形态、微观骨折形态、骨折周围颜色变化和影像形态的改变等方法，但是这些传统方法存在主观、耗时、准确性低的问题。
傅里叶变换红外光谱(Fourier transform infrared spectrometry，FTIRspectrometry)技术是兴起于20世纪中后期的第三代红外光谱技术，不同于色散型红外分光的原理，其将收集到的经过干涉的原始红外光谱进行傅里叶变换后再输出谱图，具有信息量大、信噪比高的优点。分子在红外光谱上的强度与其振动跃迁的概率相关，而跃迁概率与振动时偶极矩的变化大小有关，偶极矩变化愈大，谱带强度愈大。偶极矩的变化与基团本身固有的偶极矩有关，故基团极性越强，振动时偶极矩变化越大，吸收谱带越强；分子的对称性越高，振动时偶极矩变化越小，吸收谱带越弱。因此，红外光谱技术非常适合对由大量碳氢键、碳氧键等非对称化学键构成的有机化合物的检测。几乎所有的有机化合物在红外光谱上都有吸收峰，根据吸收峰的位置、强度、形态等信息即可推知样本的分子构成。
化学计量学(Chemometrics)又称化学统计学，由瑞典学者Wold在1971年提出，是一门通过统计学或数学方法将对化学体系的测量值与体系的状态之间建立联系的学科。近年来，随着分析化学的发展，各种现代化测量仪器层出不穷，获得了大量关于物质化学组成的数据，应用化学计量学的方法可以从中挖掘出有价值的信息。
Yang等人以弥漫性轴索损伤(diffuse axonal injury，DAI)大鼠模型为样本，将FTIR显微光谱(FTIR Microspectroscopy，FTIRM)分析结果与HE染色、β淀粉样前体蛋白(β-Amyloid Precursor Protein，βAAPP)免疫组织化学染色和银染色相对比确证，发现大鼠脑组织DAI区域的酰胺Ⅱ带的吸收峰明显增强。Lin等人同样利用FTIRM技术将过敏性休克致死者的肺水肿液与机械性窒息、脑损伤和心源性猝死的死者进行对比分析，发现过敏性休克组含有更多的α螺旋蛋白结构，而富含酪氨酸的蛋白较少，并开发了用于鉴定过敏性休克死亡的判别模型。但相对于病理样本，生前死后骨折样本难以获取，且骨折部分的组织随着腐败降解不断发生变化(例如，生前骨折尽管具有周围组织出血、血肿形成和愈伤组织的形成等明显地骨折愈合迹象，但是这些骨折愈合迹象信息会随着腐败降解时间的延长而消失)，影响了生前死后骨折的鉴别模型的建立，这是导致难以将生前骨折与死后骨折鉴别区分的原因之一。目前尚未见到将傅里叶变换红外光谱和化学计量学结合，并应用于分析鉴别生前死后骨折的报道。
发明内容
针对目前难以快速、客观、准确地鉴别生前死后骨折的问题，本发明提供了一种结合傅里叶变换红外光谱和化学计量学分析鉴别生前死后骨折的方法。
为实现上述目的，本发明采用了以下技术方案：
1)提取待分析鉴别生前死后骨折的骨组织样本(一般取自案发现场尸体的骨折部分，例如，颅骨骨折等)，采集该样本的傅里叶变换红外光谱数据；
2)对采集的傅里叶变换红外光谱数据，选取该红外光谱数据中的骨组织指纹区(2000-400cm-1)数据进行校准；
3)将经过校准的骨组织指纹区数据(具体指指纹区内选定的各个波长点对应的吸收值，波长点在指纹区内间隔分布，波长点数为800～850)，输入生前死后骨折鉴别模型，由该模型输出鉴别结果(即鉴别所述骨组织处发生的骨折属于生前骨折还是死后骨折)；其中，所述生前死后骨折鉴别模型是利用生前及死后不同骨折时间的动物模型在死后降解过程中多个时间点处相应的骨组织样本傅里叶变换红外光谱数据(参照步骤1、2获取)，并采用偏最小二乘法判别分析建立的。
优选的，所述骨组织样本提取自骨折断面至断面以下3mm的区域内。
优选的，所述傅里叶变换红外光谱数据的采集方法包括以下步骤：将骨折断面磨为粉状后制成多个压片，对每个压片平行进行多次傅里叶变换红外光谱分析，然后将多次分析的结果平均，得到对应骨组织样本的傅里叶变换红外光谱数据。
优选的，所述骨组织指纹区包括蛋白质成分酰胺Ⅰ和酰胺Ⅱ的吸收峰。
优选的，所述校准包括以下步骤：对选取的红外光谱数据进行均值中心变换和标准正态变量变换。
优选的，所述生前及死后不同骨折时间选自死前7天以内以及死后60天以内；所述死后降解过程对应的时间长度为死后60天以内。
优选的，所述动物模型建立于家兔的胫骨。
本发明的有益效果体现在：
本发明利用提取自骨折动物模型(例如，兔子胫骨骨折)样品，建立可用于有效推断人体骨折时间(生前死后骨折鉴别)的PLS-DA模型，有效克服了人体生前死后骨折样品难以获得的问题，通过建立生前死后骨折鉴别模型，仅需要采集骨组织样品的傅里叶变换红外光谱数据，经自动校准后即可应用鉴别模型，从而鉴别生前死后骨折。相较于传统的鉴别方法，本发明具有样品需求量小、过程简便快速，以及结果客观、准确的优点。
进一步的，本发明可以在不同生前骨折时间(7天以内)、死后骨折时间(0-120天)和死后降解时间(0-120天)条件下鉴别生前死后骨折，预测能力强。
进一步的，本发明选取的指纹区，可以更有效的利用骨组织红外光谱图像的信息(剔除了无效数据)，缩短了数据分析时间。
进一步的，本发明利用采集的平行光谱数据的平均光谱代表一个样品的红外光谱采样数据，增强了样品的代表性。
进一步的，本发明选取骨折断面至断面以下3mm深度的骨组织进行提取，保留了更多的骨折信息，同时有效减少了样品提取量。
附图说明
图1为制作骨折检测样品流程图。
图2为检测样品平均光谱图及对应吸收峰的官能团图。
图3为生前死后骨折PLS-DA鉴别模型图。
图4为家兔预测集随机样品鉴别情况图。
图5为人骨随机样品鉴别情况图。
具体实施方式
下面结合附图和实施例对本发明作进一步详细说明。所述的实施例仅用于解释本发明，而非对本发明保护范围的限制。
(一)建立生前死后骨折鉴别模型的依据
人类的离体生前骨折和死后骨折样本是很难获得的，针对目前难以鉴别生前死后骨折的问题，关键是建立生前死后骨折鉴别模型。采用同为哺乳动物的兔子的胫骨骨折样本进行建模，兔子骨折修复物质变化过程与人骨折修复物质变化过程几乎完全相同，所以利用兔骨所建模型推断人的骨折时间是可行的。针对骨折断面获得原始骨粉，以求保留最大化骨折信息，通过化学计量学发现对于鉴别生前死后骨折贡献最大的是代表蛋白质成分的酰胺Ⅰ和酰胺Ⅱ带吸收峰。这更能说明生前骨折和死后骨折的主要区别在于生前骨折具有修复阶段，期间有大量炎症因子、趋化因子、生长因子和胶原蛋白产生，而这些变化在死后120天内仍然可以被傅里叶变换红外光谱所发现，从而区别于死后骨折。根据以上发现，本发明提出了一种结合傅里叶变换红外光谱和化学计量学分析鉴别生前死后骨折的方法，适合各种不同型号和厂家的FTIR红外光谱仪。
(二)鉴别生前死后骨折的实例
步骤1：建立训练样品集和预测样品集
1.1建立动物实验模型
用102只家兔的双侧胫骨(共204个样本)建立实验模型，其中72只兔子的双侧胫骨，一侧为生前骨折模型，一侧为死后骨折模型，用于建立训练模型，具体分组参见表1。剩余的30只兔子的双侧胫骨建立预测集。此外，收集了4个人类颅骨骨折样品，将遭遇车祸死亡的2例人颅骨骨折样品定性为生前骨折，死后人为造骨折的2例定性为死后骨折。
表1.样品训练集分组
针对表1，在实践中，死亡近期的生前骨折(4天、7天)和濒死期骨折(生前2小时)对案件的性质起着决定性作用，所以将这两个时间段的骨折统一归为生前骨折，与死后骨折进行鉴别。对照组是指在死后降解一段时间时(2小时、15天、60天)进行骨折处理。降解条件：置于恒温箱中，温度25±1℃，湿度50％±5％。
1.2取样和样品傅里叶变换红外光谱采集
参见图1，将骨折断面的骨组织(骨组织选自骨折断面至断面以下3mm)磨成骨粉，将骨粉过200目筛，得原粉，将原粉置于60℃的烘干箱内干燥6小时，将0.9～1.1mg的干燥后的原粉与100mg的溴化钾粉末混合均匀后压片，得到检测样品，用赛默飞Nicolet-5700傅里叶红外光谱仪进行检测(放入红外漫反射附件中，同时收集背景进行基线校正)，扫描次数为32次，分辨率为4cm-1，收集4000-400cm-1范围的红外光谱，每个样品制成3个压片，每个压片平行测定3次(测试温度：常温)，最终9个光谱数据的平均光谱代表样品光谱，截取2000-400cm-1范围内的光谱数据，作为指纹区(参见图2)。
1.3光谱数据预处理
用均值中心化(mean centering)的方法，即通过光谱均值中心变换对选取的指纹图谱进行预处理，预处理中先计算出样品训练集的平均光谱，再用各样品的原始光谱(即1.2得到的样品光谱)减去样品训练集的平均光谱，得到均值中心化光谱数据。用标准正态变量变换(standard normal variate transformation，SNV)对选取的指纹区(均值中心化光谱数据)进行预处理，得到标准正态光谱数据。最终总共得到家兔训练集(144个数据)，家兔预测集(60个数据)，人体待测光谱数据4个。
步骤2：建立训练模型
应用化学计量学软件，将经过预处理的训练集样品的光谱数据(具体指指纹区内等间隔分布的829个波长点对应的吸收值)导入MATLAB中，运用PLS Toolbox 8.6(Eigenvector Research,WA,USA)工具中主成分分析(Principal Component Analysis，PCA)，发现贡献最大区域位于酰胺Ⅰ和酰胺Ⅱ带吸收峰。
将生前骨折时间跨度为2小时-7天，死后骨折时间跨度为0-60天，死后降解时间跨度为0-60天的72只家兔的双侧胫骨(共144个样本)作为训练集用于建立模型，每个分类模型都分别独立进行优化，优化主要受潜变量数(latent variables，LVs)的影响，为最优化各分类模型，针对每个分类模型，利用20折交叉验证(20-fold cross-validations，CV)选取最佳潜变量数为3(LVs从1到20中选取)，然后用偏最小二乘法判别分析(Partial leastsquares discrimination analysis,PLS-DA)建立最优生前死后骨折鉴别PLS-DA模型(图3)。
步骤3：预测样品集验证PLS-DA模型稳健性和准确性
将兔子预测集的60个数据导入PLS-DA模型进行生前死后骨折预测，将预测结果与样本实际骨折时间进行比较，发现成功识别为生前骨折的准确率为83.3％，成功识别为死后骨折的准确率为90％(见图4)，确定PLS-DA模型具有稳健性和准确性。
步骤4：人体骨折时间判定
对本实例中的4个人颅骨样品进行骨折时间判定，预测结果见图5。由图5可以看出，4个样品中有3个被正确识别，有一个生前骨折被识别为死后骨折，但是几乎在分界线边缘，由于人体样本数量有限，识别死后骨折的准确率为100％，识别生前骨折的准确率为50％。以上鉴别人体生前死后骨折的判定结果在准确性远远高于传统方法。
总之，本发明所建立的结合傅里叶变换红外光谱和化学计量学分析鉴别生前死后骨折的PLS-DA模型，具有样品需求量少、样品制备和数据处理简便、鉴别速度快，并且准确性高的优点。
Claims (7)
1.一种鉴别生前死后骨折的方法，其特征在于：包括以下步骤：
1)提取待分析鉴别生前死后骨折的骨组织样本，采集该样本的傅里叶变换红外光谱数据；
2)对采集的傅里叶变换红外光谱数据，选取该红外光谱数据中的骨组织指纹区数据进行校准；
3)将经过校准的骨组织指纹区数据，输入生前死后骨折鉴别模型，由该模型输出鉴别结果；其中，所述生前死后骨折鉴别模型是利用生前及死后不同骨折时间的动物模型在死后降解过程中多个时间点处相应的骨组织样本傅里叶变换红外光谱数据，并采用偏最小二乘法判别分析建立的。
2.根据权利要求1所述一种鉴别生前死后骨折的方法，其特征在于：所述骨组织样本提取自骨折断面至断面以下3mm的区域内。
3.根据权利要求1所述一种鉴别生前死后骨折的方法，其特征在于：所述傅里叶变换红外光谱数据的采集方法包括以下步骤：将骨折断面磨为粉状后制成多个压片，对每个压片平行进行多次傅里叶变换红外光谱分析，然后将多次分析的结果平均，得到对应骨组织样本的傅里叶变换红外光谱数据。
4.根据权利要求1所述一种鉴别生前死后骨折的方法，其特征在于：所述骨组织指纹区包括酰胺Ⅰ和酰胺Ⅱ的吸收峰。
5.根据权利要求1所述一种鉴别生前死后骨折的方法，其特征在于：所述校准包括以下步骤：对选取的红外光谱数据进行均值中心变换和标准正态变量变换。
6.根据权利要求1所述一种鉴别生前死后骨折的方法，其特征在于：所述生前及死后不同骨折时间选自死前7天以内以及死后60天以内；所述死后降解过程对应的时间长度为死后60天以内。
7.根据权利要求1所述一种鉴别生前死后骨折的方法，其特征在于：所述动物模型建立于家兔。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010442559.5A CN111707636B (zh) | 2020-05-22 | 2020-05-22 | 结合傅里叶变换红外光谱和化学计量学分析鉴别生前死后骨折的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010442559.5A CN111707636B (zh) | 2020-05-22 | 2020-05-22 | 结合傅里叶变换红外光谱和化学计量学分析鉴别生前死后骨折的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111707636A true CN111707636A (zh) | 2020-09-25 |
| CN111707636B CN111707636B (zh) | 2021-07-13 |
Family
ID=72537353
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010442559.5A Active CN111707636B (zh) | 2020-05-22 | 2020-05-22 | 结合傅里叶变换红外光谱和化学计量学分析鉴别生前死后骨折的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111707636B (zh) |
Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1242577A1 (en) * | 1999-12-28 | 2002-09-25 | The Board Of Regents Of The University Of Michigan | Process for ex vivo formation of mammalian bone and uses thereof |
| CN101549084A (zh) * | 2009-06-03 | 2009-10-07 | 贵州师范大学 | 吉祥草药材的质量检测方法 |
| CN101657096A (zh) * | 2006-11-14 | 2010-02-24 | Enzo治疗公司 | 骨形成和骨重建的组合物与方法 |
| US20110118227A1 (en) * | 2003-06-18 | 2011-05-19 | White Mountain Pharma, Inc. | Methods for the Treatment of Fibromyalgia and Chronic Fatigue Syndrome |
| CN102128921A (zh) * | 2009-12-18 | 2011-07-20 | 霍夫曼-拉罗奇有限公司 | 用于区分风湿性和非风湿性病症的触发测定方法 |
| US20120015368A1 (en) * | 2010-07-16 | 2012-01-19 | Thomas Jefferson University | Biomarkers for early diagnosis of systemic tissue fibrosis |
| WO2013072504A1 (en) * | 2011-11-16 | 2013-05-23 | Centre National De La Recherche Scientifique (Cnrs) | Crystal structure of the procollagen iii c-propeptide trimer and applications thereof |
| CN109633178A (zh) * | 2019-01-30 | 2019-04-16 | 嘉兴行健生物科技有限公司 | 一种骨钙素-n端肽检测试剂及其检测方法和应用 |
| WO2019209868A1 (en) * | 2018-04-23 | 2019-10-31 | Chris Yu | New apparatus and methods for disease detection |
| CN110726695A (zh) * | 2019-11-19 | 2020-01-24 | 湖北医药学院 | 小鼠脑组织不同部位天麻素近红外光谱特异性识别方法 |
| CN110749575A (zh) * | 2019-10-17 | 2020-02-04 | 夏永刚 | 一种中药多糖二维红外光谱鉴别预测模型及其构建方法和应用 |
-
2020
- 2020-05-22 CN CN202010442559.5A patent/CN111707636B/zh active Active
Patent Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1242577A1 (en) * | 1999-12-28 | 2002-09-25 | The Board Of Regents Of The University Of Michigan | Process for ex vivo formation of mammalian bone and uses thereof |
| US20110118227A1 (en) * | 2003-06-18 | 2011-05-19 | White Mountain Pharma, Inc. | Methods for the Treatment of Fibromyalgia and Chronic Fatigue Syndrome |
| CN101657096A (zh) * | 2006-11-14 | 2010-02-24 | Enzo治疗公司 | 骨形成和骨重建的组合物与方法 |
| CN101549084A (zh) * | 2009-06-03 | 2009-10-07 | 贵州师范大学 | 吉祥草药材的质量检测方法 |
| CN102128921A (zh) * | 2009-12-18 | 2011-07-20 | 霍夫曼-拉罗奇有限公司 | 用于区分风湿性和非风湿性病症的触发测定方法 |
| US20120015368A1 (en) * | 2010-07-16 | 2012-01-19 | Thomas Jefferson University | Biomarkers for early diagnosis of systemic tissue fibrosis |
| WO2013072504A1 (en) * | 2011-11-16 | 2013-05-23 | Centre National De La Recherche Scientifique (Cnrs) | Crystal structure of the procollagen iii c-propeptide trimer and applications thereof |
| WO2019209868A1 (en) * | 2018-04-23 | 2019-10-31 | Chris Yu | New apparatus and methods for disease detection |
| CN109633178A (zh) * | 2019-01-30 | 2019-04-16 | 嘉兴行健生物科技有限公司 | 一种骨钙素-n端肽检测试剂及其检测方法和应用 |
| CN110749575A (zh) * | 2019-10-17 | 2020-02-04 | 夏永刚 | 一种中药多糖二维红外光谱鉴别预测模型及其构建方法和应用 |
| CN110726695A (zh) * | 2019-11-19 | 2020-01-24 | 湖北医药学院 | 小鼠脑组织不同部位天麻素近红外光谱特异性识别方法 |
Non-Patent Citations (1)
| Title |
|---|
| SAMUEL GOURION-ARSIQUAUD: "Use of FTIR Spectroscopic Imaging to Identify Parameters Associated With Fragility Fracture", 《JOURNAL OF BONE AND MINERAL RESEARCH》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111707636B (zh) | 2021-07-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Bhargava | Towards a practical Fourier transform infrared chemical imaging protocol for cancer histopathology | |
| CN111488926B (zh) | 一种基于优化模型的土壤有机质测定方法 | |
| CN108169213A (zh) | 一种激光诱导击穿光谱谱峰元素自动识别方法 | |
| CN106018332A (zh) | 一种柑桔黄龙病的近红外光谱田间检测方法 | |
| CN108844917A (zh) | 一种基于显著性假设检验和偏最小二乘法的近红外光谱数据分析方法 | |
| CN112863665B (zh) | 基于病理及活体组织临床诊断大数据的疾病类型判别方法 | |
| AU2009284671A1 (en) | A method of diagnosing a respiratory disease | |
| CN104965973B (zh) | 一种苹果霉心病多因子无损检测判别模型及其建立方法 | |
| CN106770607B (zh) | 一种利用hs‐imr‐ms鉴别卷烟真伪的方法 | |
| Fan et al. | Direct calibration transfer to principal components via canonical correlation analysis | |
| CN114611582B (zh) | 一种基于近红外光谱技术分析物质浓度的方法及系统 | |
| Fatima et al. | Towards normalization selection of Raman data in the context of protein glycation: Application of validity indices to PCA processed spectra | |
| CN111707636B (zh) | 结合傅里叶变换红外光谱和化学计量学分析鉴别生前死后骨折的方法 | |
| CN111366554B (zh) | 一种基于衰减全反射红外光谱的痕迹遗留时间的测量方法 | |
| JP2009244015A (ja) | トウキの品質鑑定方法 | |
| Lin et al. | An optimization strategy for detection of fertile pigeon egg based on NIR spectroscopy analysis | |
| CN111398250A (zh) | 一种基于光与物质相互作用所产生的分子碎片谱的肿瘤诊断方法 | |
| CN110726695A (zh) | 小鼠脑组织不同部位天麻素近红外光谱特异性识别方法 | |
| ZA202305116B (en) | Establishment of unlabeled quantitative method for detecting collagen in lung tissue | |
| CN106770016A (zh) | 单粒稻谷种子的蛋白质定量分析nir透射光谱测量方法 | |
| CN113341044B (zh) | 一种基于代谢组学标志物鉴定溺死的方法及其应用 | |
| CN110887810A (zh) | 一种基于近红外光谱技术的血塞通滴丸质量一致性评价方法 | |
| Ning et al. | Construction of a ginsenoside content-predicting model based on hyperspectral imaging | |
| Pereira et al. | Evaluating biochemical differences in thyroglobulin from normal and goiter tissues by infrared spectral imaging | |
| CN106568740A (zh) | 一种近红外光谱快速判定茶鲜叶品种的方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |