CN111366554B - 一种基于衰减全反射红外光谱的痕迹遗留时间的测量方法 - Google Patents
一种基于衰减全反射红外光谱的痕迹遗留时间的测量方法 Download PDFInfo
- Publication number
- CN111366554B CN111366554B CN202010452130.4A CN202010452130A CN111366554B CN 111366554 B CN111366554 B CN 111366554B CN 202010452130 A CN202010452130 A CN 202010452130A CN 111366554 B CN111366554 B CN 111366554B
- Authority
- CN
- China
- Prior art keywords
- trace
- sample
- remaining time
- infrared spectrum
- total reflection
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 38
- 238000005102 attenuated total reflection Methods 0.000 title claims abstract description 26
- 238000004566 IR spectroscopy Methods 0.000 title claims abstract description 23
- 238000002329 infrared spectrum Methods 0.000 claims abstract description 60
- 239000000463 material Substances 0.000 claims abstract description 51
- 238000001514 detection method Methods 0.000 claims abstract description 43
- 238000013528 artificial neural network Methods 0.000 claims abstract description 39
- 210000004369 blood Anatomy 0.000 claims description 43
- 239000008280 blood Substances 0.000 claims description 43
- 238000012549 training Methods 0.000 claims description 33
- 238000012795 verification Methods 0.000 claims description 23
- 210000004243 sweat Anatomy 0.000 claims description 5
- 238000005259 measurement Methods 0.000 claims description 4
- 238000010200 validation analysis Methods 0.000 claims 1
- 238000002360 preparation method Methods 0.000 abstract description 4
- 238000004088 simulation Methods 0.000 abstract description 3
- 238000011840 criminal investigation Methods 0.000 abstract description 2
- 239000000523 sample Substances 0.000 description 70
- 150000001408 amides Chemical class 0.000 description 20
- 238000011160 research Methods 0.000 description 17
- 238000005516 engineering process Methods 0.000 description 16
- 238000010521 absorption reaction Methods 0.000 description 14
- 238000005070 sampling Methods 0.000 description 10
- 239000000758 substrate Substances 0.000 description 6
- 238000000338 in vitro Methods 0.000 description 5
- 238000007689 inspection Methods 0.000 description 5
- 230000003595 spectral effect Effects 0.000 description 5
- 238000004611 spectroscopical analysis Methods 0.000 description 5
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N Silicium dioxide Chemical compound O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 238000011835 investigation Methods 0.000 description 4
- 210000002569 neuron Anatomy 0.000 description 4
- 239000000741 silica gel Substances 0.000 description 4
- 229910002027 silica gel Inorganic materials 0.000 description 4
- 238000003860 storage Methods 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 3
- 210000003743 erythrocyte Anatomy 0.000 description 3
- 238000005457 optimization Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 229920004933 Terylene® Polymers 0.000 description 2
- 238000000862 absorption spectrum Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000000701 chemical imaging Methods 0.000 description 2
- 239000003814 drug Substances 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 125000000524 functional group Chemical group 0.000 description 2
- 229920000728 polyester Polymers 0.000 description 2
- 239000005020 polyethylene terephthalate Substances 0.000 description 2
- 238000004154 testing of material Methods 0.000 description 2
- 238000001931 thermography Methods 0.000 description 2
- 210000001519 tissue Anatomy 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 238000004435 EPR spectroscopy Methods 0.000 description 1
- 238000005033 Fourier transform infrared spectroscopy Methods 0.000 description 1
- 235000010363 Vitex negundo Nutrition 0.000 description 1
- 244000248021 Vitex negundo Species 0.000 description 1
- 125000003368 amide group Chemical group 0.000 description 1
- 239000003146 anticoagulant agent Substances 0.000 description 1
- 229940127219 anticoagulant drug Drugs 0.000 description 1
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 1
- 238000001210 attenuated total reflectance infrared spectroscopy Methods 0.000 description 1
- 239000012496 blank sample Substances 0.000 description 1
- 210000004556 brain Anatomy 0.000 description 1
- 239000000969 carrier Substances 0.000 description 1
- 239000000919 ceramic Substances 0.000 description 1
- 210000003710 cerebral cortex Anatomy 0.000 description 1
- 238000002939 conjugate gradient method Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000008602 contraction Effects 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000011496 digital image analysis Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000001962 electrophoresis Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000003203 everyday effect Effects 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000004817 gas chromatography Methods 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 238000011478 gradient descent method Methods 0.000 description 1
- 238000004128 high performance liquid chromatography Methods 0.000 description 1
- 238000002129 infrared reflectance spectroscopy Methods 0.000 description 1
- 210000000231 kidney cortex Anatomy 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 210000004185 liver Anatomy 0.000 description 1
- 210000004072 lung Anatomy 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 210000003205 muscle Anatomy 0.000 description 1
- 230000001537 neural effect Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 239000001301 oxygen Substances 0.000 description 1
- 229910052760 oxygen Inorganic materials 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 108090000623 proteins and genes Proteins 0.000 description 1
- 238000000611 regression analysis Methods 0.000 description 1
- 210000002027 skeletal muscle Anatomy 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
- 210000000952 spleen Anatomy 0.000 description 1
- 210000000225 synapse Anatomy 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
- 239000004753 textile Substances 0.000 description 1
- 230000002861 ventricular Effects 0.000 description 1
- 239000002023 wood Substances 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/3577—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light for analysing liquids, e.g. polluted water
-
- 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
- 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/55—Specular reflectivity
- G01N21/552—Attenuated total reflection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
-
- 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)
- General Health & Medical Sciences (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Biochemistry (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Analytical Chemistry (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Artificial Intelligence (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- Data Mining & Analysis (AREA)
- Evolutionary Computation (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Biomedical Technology (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
一种基于衰减全反射红外光谱的痕迹遗留时间的测量方法,采集检材痕迹的红外光谱数据,制作样本痕迹,分类样本,采集样本数据集,使用神经网络建立样本痕迹的红外光谱数据与遗留时间的对应关系模型并验证精度,将检材痕迹的红外光谱数据录入符合精度要求的模型,模拟检材痕迹的遗留时间,评估模拟结果,必要时扩充样本数据集;该方法重现性好,制样简便,结果稳定可靠,检测精度高,满足刑侦领域中进行快速无损鉴定的需要。
Description
技术领域
本发明属于刑事科学技术领域,具体涉及一种基于衰减全反射红外光谱的痕迹遗留时间的测量方法。
背景技术
在勘验检查案件现场时,除对常规的指纹、足迹、DNA等痕迹提取检验以及对血痕的基本形态、色泽、厚度、分布形态等记录检验外,对现场痕迹(包括血痕、笔迹、汗液指纹等)遗留时间的正确鉴别判断有助于客观反映案发时刻现场人员的活动状况,也会成为证实案发时间的有力依据。
受现场条件的限制以及基层科研水平的制约,在勘验检查案件现场时,对观测判断现场痕迹的遗留时间并无深入的研究进展,对尸体组织的变化规律研究集中于法医学检验,鲜见针对案件现场离体血痕的规律性探索;而针对血液的研究方向多为医用,由变化规律进行现场血痕遗留时间判定的方法比较少,对现场血痕的规律进行研究的方法多需要复杂的实验流程和昂贵的器材,且多数判断方法存在较大误差。
对痕迹的有损检验方法包括电泳技术、高效液相色谱技术、气相色谱法、氧电极技术、基因分析技术等,由于痕迹的量微小,无损耗检验痕迹遗留时间的研究方法更有利于痕迹的保存,也能够重复进行制样采样。对痕迹遗留时间的无损检验方法包括立体显微镜、扫描电子显微镜、原子力显微镜、近红外反射光谱技术、光纤光谱技术、高光谱照相法、电子顺磁共振技术等,然而,现有的无损检验多集中于实验研究。
邓晓军等(血痕红外热成像变化与血痕经过时间关系研究[J].法医检验,2015,(8):72~74)采用红外热像技术结合计算机图像分析技术,对玻璃、纺织物、A4纸、瓷砖、木块等五类载体上的血痕温度进行观测、运算,记录和统计了15分钟内五类客体上血痕自制样时起每分钟的温度情况,通过回归分析,给出了较好的回归方程,但其考察时间过短,进一步研究应用的空间较小。
研究表明,衰减全反射红外光谱法(ATR)能够揭示实验样本的微小变化所造成的红外谱图的变化,可以满足无损快速检验痕迹的需要。
黄平等人(应用FT红外光谱光谱技术推断死亡时间[J].中国法医学杂志,2011,(2):104~109)运用KBr压片法和ATR法分别对人体和大鼠的左心室肌、肺脏、肝脏、脾脏、肾脏皮质、骨骼肌、大脑皮质进行傅里叶红外光谱检测,并选择光谱数据中的主要吸收峰为研究指标,指出了人尸体离体组织用衰减全反射红外光谱的变化规律,但在死亡时间和器官离体时间方面的判断较模糊,有一定的局限性。
神经网络是“人工神经网络”的简称,神经网络是一种信息处理工具,是模仿大脑神经突触联接结构进行信息处理的数学模型,通过两个节点间的信号会被赋予加权值,建立起的神经网络经过对大量信号输入的学习训练后,调整加权实现对信息的分类,向该学习训练过的神经网络输入新的信号后,神经网络按照已确定的权重值对信息实现自动归类,是大量信息并行处理和大规模并行计算的基础。神经网络与光谱学的结合在当前的各个领域应用广泛,车牌的自动识别、白酒年份鉴别、大数据分析等领域都有应用。
张宝菊等人(基于BP神经网络的人体血液中红细胞浓度无创检测[J].光谱学与光谱分析,2012,(9):2508~2511)对基于BP神经网络对人体血液中红细胞浓度的无创检测方法进行了研究,其以新鲜血液的光谱数据作为建模数据,对新鲜血液的红细胞浓度进行预测。
中国专利申请CN110196233A、硕士论文《基于高光谱成像技术的血迹检测研究》中公开了一种基于高光谱成像的血迹时间预测方法,其以高光谱血迹三维图像立方体为数据支撑,提取血迹样本的光谱曲线,再经过预处理矩阵化,运用非线性迭代结合非线性拟合的方法,构建了血迹时间预测相关的时间-光谱域模型,然而,其模型构建过程复杂,并且,其样本为在25℃恒温条件下保存的实验室采集的血迹(干性样本),且只采集1h、24h、30天的样本数据,以此建立模型,其可靠程度存疑,以医院血库采集的血液样本(液体样本,且低温保存、含有抗凝剂)进行测试检验,采样时间不合理及样本与检材保存环境的差异会显著影响实验结论的准确性。
传统的刑事科学技术的勘验、鉴定过程在反映人与人、人与物、物与物的空间关系问题上较为成熟,却只能粗略实现时间关系上的推断,对现场痕迹的遗留时间进行分析鉴定一直是刑事科学技术工作的难点,目前,将衰减全反射红外光谱法与神经网络结合用于判断痕迹遗留时间的研究未见报道。
发明内容
本发明的目的在于提供一种基于衰减全反射红外光谱(ATR-FTIR)的痕迹遗留时间的测量方法,结合传统刑事科学技术研究、光谱学和计算机科学,通过用痕迹样本的红外光谱数据对神经网络进行训练,进而对痕迹的遗留时间进行判断,痕迹样本的红外吸收特征谱带明显,重现性好,制样简便,结果稳定可靠,检测精度高,测量天数值与真实天数值间的判定系数R2>0.8,满足刑侦领域中进行快速无损鉴定的需要。
为了达到上述目的,本发明提供如下技术方案:
一种基于衰减全反射红外光谱的痕迹遗留时间的测量方法,包括以下步骤:
1)采集检材痕迹的红外光谱数据
检材痕迹来自于案发现场,检材痕迹包括但不限于血痕、笔迹、汗液指纹等,设置固定的波数范围和分辨率,使用衰减全反射红外光谱法采集检材痕迹的红外光谱数据,
2)模拟检材痕迹制作样本痕迹,分类样本
使用与检材痕迹载体相同的材料为载体制作样本痕迹,将样本痕迹按照7:3比例,随机划分为训练用样本痕迹和验证用样本痕迹,
3)采集样本数据集
模拟检材痕迹的遗留环境,使用衰减全反射红外光谱法,要求波数范围和分辨率与前述相同,以均匀的时间间隔,采集不同遗留时间的样本痕迹的红外光谱数据,与对应的遗留时间归纳为样本数据集,记录样本痕迹的遗留时间最大值,保留样本数据集中红外光谱数据与遗留时间的对应关系,将样本数据集按照训练用和验证用,划分为训练集和验证集,
4)使用神经网络建立样本痕迹的红外光谱数据与遗留时间的对应关系模型并验证精度,以训练集中的红外光谱数据和遗留时间为变量,训练神经网络,建立样本痕迹的红外光谱数据与遗留时间的对应关系模型,使用验证集验证模型精度,如判定系数R2大于0.8,则该神经网络可用,否则重复“训练神经网络”至“使用验证集验证模型精度”的流程,直至判定系数R2大于0.8,
5)将检材痕迹的红外光谱数据录入符合精度要求的模型,模拟检材痕迹的遗留时间,向可用模型输入检材痕迹的红外光谱数据,得到模拟的检材痕迹的遗留时间,
6)评估检材痕迹遗留时间
比较模拟的检材痕迹的遗留时间和样本痕迹的遗留时间最大值,如比值小于等于0.9,则评估模拟的检材痕迹的遗留时间即为测量的检材痕迹遗留时间,测量完成,如比值大于0.9,则扩充样本数据集,延长样本痕迹的遗留时间,重复步骤3)至步骤5),比较扩充样本数据集后模拟的检材痕迹的遗留时间和样本痕迹的遗留时间最大值,直至比值小于等于0.9,则评估模拟的检材痕迹的遗留时间即为测量的检材痕迹遗留时间,测量完成。
优选地,步骤3)中,所述不同遗留时间痕迹样本的采样时间间隔24小时,红外光谱的采集波数范围为:4000cm-1-600cm-1。
本发明的技术方案还可以如下表述:
一种基于衰减全反射红外光谱的痕迹遗留时间的测量方法,包括以下步骤:
a.采集红外光谱数据集
采集待测案发现场的检材痕迹的红外光谱数据;
制作模拟案发现场的不同遗留时间的各样本痕迹,采集其红外光谱,获取红外吸收光谱数据,得到红外光谱数据集;
其中,所述样本痕迹包括训练集样本痕迹和验证集样本痕迹,红外光谱数据集包括相应的训练集样本痕迹的红外光谱数据集和验证集样本痕迹的红外光谱数据集;
b.构建神经网络
利用数据处理工具,以各样本痕迹的遗留时间和红外光谱数据为变量构建神经网络架构,设定层数,选择训练类型和误差类型,隐含层和输出层的传递函数为TANSIG,设置隐含层的数量和神经元的数量;
根据所选择的训练类型设定神经元数量及训练参数,利用训练集样本痕迹的红外光谱数据集训练构建的神经网络架构,之后再输入验证集痕迹样本的红外光谱数据,输出验证集样本痕迹的痕迹遗留时间模拟值,并与对应验证集样本痕迹遗留时间的真值进行比较,得到判定系数R2,若判定系数R2>0.8,则完成构建,否则对神经网络再次训练,直到判定系数R2>0.8,得到神经网络;
c.检测待测样本
将采集的待测案发现场的样本痕迹的红外光谱数据输入神经网络,神经网络输出待测痕迹样本的遗留时间。
进一步,步骤a)中,所述不同遗留时间样本痕迹的采样时间间隔24小时。
优选地,步骤a)中,所提取的红外光谱数据为采集波数范围内全段或部分波数段内所有吸收峰的出峰位置及峰强数据。
进一步,步骤a)中,在各痕迹样本的红外光谱中,提取自酰胺I、酰胺II、酰胺III、酰胺A、PO2-反对称伸缩振动至PO2对称伸缩振动对应波数区间段内的吸收峰数据作为红外光谱数据集。
优选地,步骤a)中,红外光谱的采集波数范围为:4000cm-1-600cm-1。
又,步骤c)中,待检测痕迹样本的红外吸收光谱数据的采集波数范围为1800cm-1-1300cm-1。
优选地,步骤b)中,所述神经网络的训练类型为量化共轭梯度法、L-M优化算法或有动量的梯度下降法。
又,步骤b)中,得到的神经网络判定系数R2>0.9。
优选地,所述痕迹样本为血痕、笔迹、汗斑或指印。
本发明对离体痕迹的相关性质随时间环境的变化规律进行探索,在传统的现场勘验基础上,拓展现场勘验采痕取证、检验鉴定思路,从微量物证角度出发,对在案件现场遗留的痕迹(如血迹、汗液痕迹或笔迹)进行检验分析和归纳总结,对痕迹离体或产生后随时间和环境因素的变化规律进行了研究。
本发明使用衰减全反射红外光谱结合神经网络方法对痕迹的遗留时间进行预测,使用衰减全反射红外光谱技术方法对痕迹进行红外光谱采集,能够清晰反映痕迹的基本特征。对同一份基板上遗留痕迹后的不同时间采集到的红外光谱,能够反映出基板上痕迹的红外性质随时间变化的基本规律,且这个规律不因遗留痕迹个体的变化发生转移。
本发明中,酰胺I与酰胺II的吸收峰是血痕在红外光谱图中峰强度最高的两处峰值,二者之间线性关系良好,在波数为1536cm-1处的酰胺I峰强度数值相对于波数为1650cm-1处的酰胺I峰强度数值存在良好的线性关系,一方面说明衰减全反射红外光谱法采集到的血痕红外光谱数据具备相对稳定可用的条件;另一方面说明该两处峰值之间的相对关系受到的样本不均匀、实验操作差异等因素的影响较小,有利于针对重点峰位进行细化研究。
与现有技术相比,本发明具有如下有益效果:
本发明将传统刑事科学技术研究与光谱学、计算机科学进行结合,对刑事科学技术痕迹遗留时间分析鉴定进行探索,以样本痕迹的遗留时间和红外光谱数据为变量对神经网络进行训练,建立样本痕迹的红外光谱数据与遗留时间的对应关系模型,通过验证筛选出较高精度的模型,能够实现对痕量检材的无损快速检测,使用时投入的人工成本和后期维护、耗材成本可控,具备普及条件,填补刑事科学技术领域中痕迹样本遗留时间检测的技术空白。
本发明采用全波段谱图数据训练神经网络,测量精度较高,通过调整重点波数范围可能进一步提高测量精度,痕迹的红外吸收特征谱带明显,重现性好,制样简便,结果稳定可靠。
附图说明
图1-7依次为本发明实施例中样品D1YP1-1、D2YP1-1、D3YP1-1、D4YP1-1、D5YP1-1、D6YP1-1、D7YP1-1的红外数据谱图。
具体实施方式
以下结合具体实施例对本发明作进一步说明。
实施例 本发明基于衰减全反射红外光谱的痕迹遗留时间的测量方法,包括以下步骤:
1.采集检材痕迹的红外光谱数据
采集一份以涤纶基板层析硅胶为载体的血痕的红外光谱数据,其采集条件如下:
实验仪器:布鲁克V70傅里叶变换红外光谱仪(带显微镜附件和ATR附件);
扫描范围:4000cm-1-600cm-1;
空间分辨率:25μm-1;
测试方法:衰减全反射红外光谱法(红外显微镜);
扫描次数:32次;
温度:21℃;
湿度:40%。
使用带有ATR附件的红外显微镜,以4000cm-1-600cm-1的波数段为扫描范围,如表1,其中,“……”表示前后波数区间内未全部列出的部分红外光谱数据:
表1
波数 | 吸收峰强度 | 波数 | 吸收峰强度 | 波数 | 吸收峰强度 |
4001.49232 | 0.00033 | 2985.20969 | 0.00054 | 1000.85519 | 0.03265 |
3999.56389 | 0.0009 | 2983.28126 | 0.00029 | 998.92676 | 0.03175 |
3997.63546 | 0.00132 | …… | …… | 996.99833 | 0.03082 |
3995.70703 | 0.0021 | 2015.20938 | 0.00059 | 995.0699 | 0.02997 |
3993.7786 | 0.00269 | 2013.28095 | 0.00056 | 993.14147 | 0.02955 |
3991.85017 | 0.00241 | 2011.35252 | 0.00052 | 991.21304 | 0.02964 |
3989.92174 | 0.00166 | 2009.42409 | 0.00054 | 989.28461 | 0.02996 |
3987.99331 | 0.00109 | 2007.49566 | 0.00044 | 987.35618 | 0.02996 |
3986.06488 | 0.0009 | 2005.56723 | 0.00037 | 985.42775 | 0.02932 |
3984.13645 | 0.00104 | 2003.6388 | 0.00058 | 983.49932 | 0.02887 |
…… | …… | 2001.71037 | 0.00085 | …… | …… |
3010.27928 | 0.00123 | 1999.78194 | 0.0008 | 622.8829 | 0.00852 |
3008.35085 | 0.00097 | 1997.85351 | 0.00051 | 620.95447 | 0.00642 |
3006.42242 | 0.00098 | 1995.92508 | 0.00017 | 619.02604 | 0.00503 |
3004.49399 | 0.00138 | 1993.99665 | 0 | 617.09761 | 0.006 |
3002.56556 | 0.00166 | 1992.06822 | 0.00027 | 615.16918 | 0.00765 |
3000.63713 | 0.00143 | 1990.13979 | 0.00053 | 613.24075 | 0.00557 |
2998.7087 | 0.00125 | 1988.21136 | 0.00051 | 611.31232 | 0.0018 |
2996.78027 | 0.00062 | …… | …… | 609.38389 | 0 |
2994.85184 | 0 | 1010.49734 | 0.03978 | 607.45546 | 0.00402 |
2992.92341 | 0.00008 | 1008.56891 | 0.03851 | 605.52703 | 0.00985 |
2990.99498 | 0.00091 | 1006.64048 | 0.03707 | 603.5986 | 0.01127 |
2989.06655 | 0.00168 | 1004.71205 | 0.03539 | 601.67017 | 0.00616 |
2987.13812 | 0.00143 | 1002.78362 | 0.03381 | 599.74174 | 0 |
2.模拟检材痕迹制作样本痕迹,分类样本
以涤纶基板层析硅胶为载体,制成样本痕迹;
采样过程:收集到10名志愿者提供的新鲜静脉血样本,每人取150μL血液滴于涤纶基板的层析硅胶片上,制成样本血痕,静置于温度为21℃,湿度为40%的室内保存,以其中7份血痕样本为训练集样本,以其中3份样本为验证集样本。
3.采集样本数据集
将10份实验样本置于实验环境内,采血后开始计时,采血后24小时为1d采样时间,以此类推至7d。分别于采血后1d、2d、3d、4d、5d、6d、7d进行采样。采样时将实验样本置于实验仪器载物台上,在红外显微镜下标记3个采样点,每天对该3个采样点的红外光谱进行重复采样,另制得涤纶基板的层析硅胶片表面空白样本一份。记录样本痕迹的遗留时间最大值7。
以OPUS软件记录采集样本数据,对采集的样品数据进行编号,以D1YP1-1为例进行说明,D为天数,YP1-1代表1号样品上的第一个点上的数据,共得有效数据210份,按照前述样本分类将其划分为147个训练集血痕样本和63个验证集血痕样本。
分别采集147个训练集血痕样本和63个验证集血痕样本的红外光谱,采集每个样本在波数段4000cm-1-600cm-1间的1765个吸收峰的出峰位置的峰强值,得到红外光谱数据集,包括训练集血痕样本的红外光谱数据集(简称训练集,下同)和验证集血痕样本的红外光谱数据集(简称验证集,下同)。
采集到的210份血痕样本的红外光谱图在酰胺I、酰胺II、酰胺A等重要的官能团波数段上表现清晰,以1号样品上第一个点为例,其从第1天到第7天的红外谱图见附图1-7。
将各样品的数据酰胺I、酰胺I、酰胺II、酰胺III、酰胺A、PO2-反对称伸缩振动、PO2对称伸缩振动对应波数区间的最强峰值抽取出来进行考量,参见表2。
表2
样品名 | 酰胺Ⅰ1650 | 酰胺Ⅱ1536 | 酰胺A3289 | 酰胺Ⅲ1315 | PO2-反对称伸缩振动1080 | PO2反对称收缩振动1238 |
D1YP1-1 | 0.08089 | 0.05406 | 0.01598 | 0.00294 | 0.1304 | 0.03477 |
D1YP1-2 | 0.18241 | 0.11802 | 0.04344 | 0.00515 | 0.22214 | 0.04731 |
D1YP1-3 | 0.09355 | 0.05604 | 0.01979 | 0.00134 | 0.33532 | 0.02809 |
D1YP2-1 | 0.01106 | 0.01362 | 0.00363 | 0.00477 | 0.01466 | 0.00187 |
D1YP2-2 | 0.02197 | 0.01493 | 0.0051 | 0.00118 | 0.12493 | 0.01393 |
D1YP2-3 | 0.16845 | 0.1106 | 0.03129 | 0.00438 | 0.28947 | 0.02753 |
D1YP3-1 | 0.22003 | 0.15852 | 0.04593 | 0.00799 | 0.12556 | 0.00569 |
D1YP3-2 | 0.01013 | 0.01072 | 0.00367 | 0.0045 | 0.00742 | 0.02028 |
D1YP3-3 | 0.05871 | 0.0387 | 0.01103 | 0.00228 | 0.16031 | 0.02467 |
D1YP4-1 | 0.03004 | 0.02278 | 0.00688 | 0.00226 | 0.10585 | 0.03093 |
D1YP4-2 | 0.15073 | 0.0979 | 0.03321 | 0.00707 | 0.21574 | 0.05183 |
D1YP4-3 | 0.01859 | 0.0143 | 0.00263 | 0.00338 | 0.02942 | 0.01941 |
D1YP5-1 | 0.04717 | 0.03486 | 0.0053 | 0.00296 | 0.02087 | 0.01301 |
D1YP5-2 | 0.11054 | 0.06811 | 0.02179 | 0.00345 | 0.19702 | 0.04177 |
D1YP5-3 | 0.00666 | 0.00716 | 0.00214 | 0.00145 | 0.07659 | 0.02617 |
D1YP6-1 | 0.21577 | 0.17509 | 0.04426 | 0.031 | 0.02088 | 0.05368 |
D1YP6-2 | 0.05492 | 0.03784 | 0.00929 | 0.00122 | 0.14774 | 0.02496 |
D1YP6-3 | 0.02757 | 0.01927 | 0.00344 | 0.00109 | 0.11637 | 0.01164 |
D1YP7-1 | 0.04081 | 0.02379 | 0.0077 | 0.00024 | 0.22219 | 0.03846 |
D1YP7-2 | 0.02192 | 0.01451 | 0.00316 | 0.00036 | 0.16243 | 0.02503 |
D1YP7-3 | 0.21577 | 0.17509 | 0.04426 | 0.031 | 0.02088 | 0.05368 |
D1YP8-1 | 0.05492 | 0.03784 | 0.00929 | 0.00122 | 0.14774 | 0.02496 |
D1YP8-2 | 0.02757 | 0.01927 | 0.00344 | 0.00109 | 0.11637 | 0.01164 |
D1YP8-3 | 0.04081 | 0.02379 | 0.0077 | 0.00024 | 0.22219 | 0.03846 |
D1YP9-1 | 0.02192 | 0.01451 | 0.00316 | 0.00036 | 0.16243 | 0.02503 |
D1YP9-2 | 0.01081 | 0.0122 | 0.00287 | 0.00316 | 0.03168 | 0.00046 |
D1YP9-3 | 0.03723 | 0.02209 | 0.00674 | 0.00015 | 0.25687 | 0.02738 |
D1YP10-1 | 0.01301 | 0.01057 | 0.00311 | 0.00138 | 0.08874 | 0.02734 |
D1YP10-2 | 0.03412 | 0.02175 | 0.00469 | 0.0003 | 0.22982 | 0.04224 |
D1YP10-3 | 0.01707 | 0.01495 | 0.00387 | 0.00238 | 0.09031 | 0.03215 |
D2YP1-1 | 0.14099 | 0.09589 | 0.03454 | 0.00119 | 0.22462 | 0.01499 |
D2YP1-2 | 0.105 | 0.07966 | 0.02657 | 0.00343 | 0.08786 | 0.03049 |
D2YP1-3 | 0.0638 | 0.04679 | 0.01451 | -0.00004 | 0.16789 | 0.01567 |
D2YP2-1 | 0.19096 | 0.13358 | 0.04617 | 0.00614 | 0.16077 | 0.03129 |
D2YP2-2 | 0.03523 | 0.0257 | 0.01158 | 0.00142 | 0.16833 | 0.03542 |
D2YP2-3 | 0.20601 | 0.17759 | 0.05072 | 0.03391 | 0.04814 | 0.05525 |
D2YP3-1 | 0.11147 | 0.0817 | 0.02686 | 0.004 | 0.10574 | 0.02959 |
D2YP3-2 | 0.01037 | 0.01042 | 0.00702 | 0.00076 | 0.02013 | 0.01037 |
D2YP3-3 | 0.06237 | 0.04308 | 0.01542 | 0.00189 | 0.15756 | 0.02826 |
D2YP4-1 | 0.04722 | 0.03337 | 0.01009 | 0.0007 | 0.18717 | 0.01154 |
D2YP4-2 | 0.07786 | 0.05543 | 0.02085 | 0.00163 | 0.0989 | 0.02065 |
D2YP4-3 | 0.19836 | 0.13065 | 0.0455 | 0.00797 | 0.25106 | 0.04816 |
D2YP5-1 | 0.01942 | 0.01451 | 0.00462 | 0.00095 | 0.0405 | 0.00752 |
D2YP5-2 | 0.04742 | 0.03801 | 0.00796 | 0.00388 | 0.02295 | 0.01432 |
D2YP5-3 | 0.05091 | 0.03441 | 0.00644 | 0.00081 | 0.1479 | 0.01042 |
D2YP6-1 | 0.06275 | 0.04736 | 0.01061 | 0.00232 | 0.05311 | 0.00309 |
D2YP6-2 | 0.06432 | 0.04094 | 0.01093 | 0.00142 | 0.26904 | 0.01901 |
D2YP6-3 | 0.0504 | 0.03368 | 0.00932 | 0.00013 | 0.22259 | 0.02952 |
D2YP7-1 | 0.04535 | 0.02795 | 0.00799 | -0.00026 | 0.32813 | 0.04209 |
D2YP7-2 | 0.02254 | 0.01488 | 0.00604 | 0.00111 | 0.15686 | 0.03195 |
D2YP7-3 | 0.06275 | 0.04736 | 0.01061 | 0.00232 | 0.05311 | 0.00309 |
D2YP8-1 | 0.06432 | 0.04094 | 0.01093 | 0.00142 | 0.26904 | 0.01901 |
D2YP8-2 | 0.0504 | 0.03368 | 0.00932 | 0.00013 | 0.22259 | 0.02952 |
D2YP8-3 | 0.04535 | 0.02795 | 0.00799 | -0.00026 | 0.32813 | 0.04209 |
D2YP9-1 | 0.02254 | 0.01488 | 0.00604 | 0.00111 | 0.15686 | 0.03195 |
D2YP9-2 | 0.08489 | 0.05132 | 0.02128 | 0.00063 | 0.32017 | 0.03696 |
D2YP9-3 | 0.03501 | 0.02431 | 0.00596 | 0.00036 | 0.23966 | 0.02576 |
D2YP10-1 | 0.02242 | 0.0158 | 0.00656 | 0.00043 | 0.1319 | 0.01342 |
D2YP10-2 | 0.03179 | 0.02284 | 0.00673 | -0.00019 | 0.22954 | 0.02373 |
D2YP10-3 | 0.03256 | 0.02177 | 0.00857 | 0.00043 | 0.2301 | 0.03148 |
D3YP1-1 | 0.05401 | 0.03719 | 0.01132 | 0.00306 | 0.10351 | 0.01695 |
D3YP1-2 | 0.14379 | 0.0919 | 0.03807 | 0.00447 | 0.25191 | 0.0479 |
D3YP1-3 | 0.04599 | 0.0274 | 0.00919 | 0.00414 | 0.15445 | 0.03734 |
D3YP2-1 | 0.02212 | 0.01639 | 0.00735 | 0.0018 | 0.07849 | 0.00626 |
D3YP2-2 | 0.02939 | 0.01584 | 0.00744 | 0.00247 | 0.21956 | 0.04261 |
D3YP2-3 | 0.22217 | 0.18965 | 0.0513 | 0.05108 | 0.06586 | 0.07621 |
D3YP3-1 | 0.12576 | 0.1041 | 0.03022 | 0.01658 | 0.13519 | 0.02832 |
D3YP3-2 | 0.00588 | 0.01866 | 0.00316 | 0.01938 | 0.03328 | 0.04826 |
D3YP3-3 | 0.09794 | 0.07345 | 0.02738 | 0.01488 | 0.24178 | 0.05386 |
D3YP4-1 | 0.04119 | 0.03907 | 0.01143 | 0.01663 | 0.14282 | 0.04577 |
D3YP4-2 | 0.05268 | 0.04416 | 0.01483 | 0.01628 | 0.1113 | 0.04119 |
D3YP4-3 | 0.18374 | 0.11984 | 0.0397 | 0.00955 | 0.24753 | 0.03944 |
D3YP5-1 | 0.00709 | 0.00769 | 0.00223 | 0.00676 | 0.01087 | 0.0131 |
D3YP5-2 | 0.03435 | 0.02669 | 0.00659 | 0.00767 | 0.02317 | 0.02075 |
D3YP5-3 | 0.05818 | 0.03813 | 0.01041 | 0.00124 | 0.15935 | 0.00698 |
D3YP6-1 | 0.04345 | 0.03471 | 0.01156 | 0.00864 | 0.03397 | 0.02176 |
D3YP6-2 | 0.03527 | 0.02128 | 0.00805 | 0.00347 | 0.13724 | 0.03263 |
D3YP6-3 | 0.06817 | 0.04079 | 0.01063 | 0.00319 | 0.28102 | 0.04666 |
D3YP7-1 | 0.02951 | 0.01602 | 0.00588 | 0.00167 | 0.22565 | 0.03591 |
D3YP7-2 | 0.0391 | 0.02164 | 0.00749 | 0.00138 | 0.25287 | 0.04196 |
D3YP7-3 | 0.04345 | 0.03471 | 0.01156 | 0.00864 | 0.03397 | 0.02176 |
D3YP8-1 | 0.03527 | 0.02128 | 0.00805 | 0.00347 | 0.13724 | 0.03263 |
D3YP8-2 | 0.06817 | 0.04079 | 0.01063 | 0.00319 | 0.28102 | 0.04666 |
D3YP8-3 | 0.02951 | 0.01602 | 0.00588 | 0.00167 | 0.22565 | 0.03591 |
D3YP9-1 | 0.0391 | 0.02164 | 0.00749 | 0.00138 | 0.25287 | 0.04196 |
D3YP9-2 | 0.07131 | 0.04436 | 0.01572 | 0.00092 | 0.25308 | 0.02301 |
D3YP9-3 | 0.03165 | 0.01882 | 0.00376 | 0.00022 | 0.35489 | 0.01367 |
D3YP10-1 | 0.0065 | 0.0048 | 0.00289 | 0.00258 | 0.04835 | 0.01818 |
D3YP10-2 | 0.04318 | 0.02143 | 0.00506 | 0.00068 | 0.37513 | 0.03453 |
D3YP10-3 | 0.03607 | 0.02148 | 0.00874 | 0.00118 | 0.18703 | 0.03476 |
D4YP1-1 | 0.0239 | 0.01951 | 0.00763 | 0.00095 | 0.04035 | 0.00524 |
D4YP1-2 | 0.115 | 0.07917 | 0.02549 | 0.00156 | 0.14182 | 0.01966 |
D4YP1-3 | 0.05883 | 0.03798 | 0.00868 | 0.00196 | 0.23941 | 0.01007 |
D4YP2-1 | 0.02724 | 0.02068 | 0.00667 | 0.00144 | 0.0758 | 0.00737 |
D4YP2-2 | 0.01509 | 0.01295 | 0.0052 | 0.00099 | 0.0759 | 0.0091 |
D4YP2-3 | 0.19176 | 0.12844 | 0.0407 | 0.00264 | 0.20994 | 0.02328 |
D4YP3-1 | 0.08068 | 0.06284 | 0.01677 | 0.00244 | 0.07657 | 0.00753 |
D4YP3-2 | 0.00944 | 0.01029 | 0.00571 | 0.00194 | 0.00702 | 0.00914 |
D4YP3-3 | 0.10097 | 0.06426 | 0.0192 | 0.00163 | 0.25192 | 0.02816 |
D4YP4-1 | 0.05803 | 0.03907 | 0.01143 | 0.00139 | 0.1496 | 0.0248 |
D4YP4-2 | 0.03347 | 0.02179 | 0.00472 | 0.00103 | 0.25085 | 0.0182 |
D4YP4-3 | 0.13959 | 0.09677 | 0.01815 | 0.00256 | 0.40862 | 0.0082 |
D4YP5-1 | 0.16087 | 0.10958 | 0.02664 | 0.00249 | 0.28953 | 0.02241 |
D4YP5-2 | 0.16171 | 0.09837 | 0.0315 | 0.00211 | 0.27407 | 0.02358 |
D4YP5-3 | 0.14848 | 0.0941 | 0.02838 | 0.00165 | 0.25105 | 0.03131 |
D4YP6-1 | 0.15228 | 0.10442 | 0.0311 | 0.00279 | 0.16577 | 0.0209 |
D4YP6-2 | 0.05589 | 0.03289 | 0.00838 | 0.00094 | 0.20004 | 0.02251 |
D4YP6-3 | 0.08265 | 0.05007 | 0.00961 | 0.00114 | 0.2776 | 0.02019 |
D4YP7-1 | 0.0587 | 0.03248 | 0.00937 | 0.00066 | 0.3017 | 0.02693 |
D4YP7-2 | 0.05713 | 0.03143 | 0.00761 | 0.00034 | 0.32705 | 0.03402 |
D4YP7-3 | 0.15228 | 0.10442 | 0.0311 | 0.00279 | 0.16577 | 0.0209 |
D4YP8-1 | 0.05589 | 0.03289 | 0.00838 | 0.00094 | 0.20004 | 0.02251 |
D4YP8-2 | 0.08265 | 0.05007 | 0.00961 | 0.00114 | 0.2776 | 0.02019 |
D4YP8-3 | 0.0587 | 0.03248 | 0.00937 | 0.00066 | 0.3017 | 0.02693 |
D4YP9-1 | 0.05713 | 0.03143 | 0.00761 | 0.00034 | 0.32705 | 0.03402 |
D4YP9-2 | 0.03859 | 0.02051 | 0.00686 | 0.00005 | 0.25089 | 0.02736 |
D4YP9-3 | 0.02302 | 0.01483 | 0.004 | 0.00075 | 0.19234 | 0.01555 |
D4YP10-1 | 0.02623 | 0.01584 | 0.0035 | 0.00028 | 0.23632 | 0.0252 |
D4YP10-2 | 0.03019 | 0.01852 | 0.00428 | 0.00083 | 0.18727 | 0.03251 |
D4YP10-3 | 0.01042 | 0.00876 | 0.00286 | 0.00044 | 0.07556 | 0.01226 |
D5YP1-1 | 0.07723 | 0.05588 | 0.02856 | 0.0126 | 0.12958 | 0.02034 |
D5YP1-2 | 0.09295 | 0.06794 | 0.03664 | 0.0183 | 0.2128 | 0.04516 |
D5YP1-3 | 0.05697 | 0.04745 | 0.02115 | 0.0275 | 0.17975 | 0.03629 |
D5YP2-1 | 0.07045 | 0.04795 | 0.02713 | 0.00916 | 0.18077 | 0.02758 |
D5YP2-2 | 0.02231 | 0.02252 | 0.01208 | 0.02171 | 0.02524 | 0.04046 |
D5YP2-3 | 0.15494 | 0.13254 | 0.06704 | 0.05114 | 0.22358 | 0.09948 |
D5YP3-1 | 0.21006 | 0.15889 | 0.07658 | 0.037 | 0.22867 | 0.06468 |
D5YP3-2 | 0.18802 | 0.138 | 0.06494 | 0.04148 | 0.21518 | 0.07512 |
D5YP3-3 | 0.11732 | 0.08058 | 0.04961 | 0.02443 | 0.2776 | 0.06904 |
D5YP4-1 | 0.06137 | 0.04362 | 0.02141 | 0.01637 | 0.17511 | 0.03459 |
D5YP4-2 | 0.02648 | 0.02423 | 0.01115 | 0.02183 | 0.09718 | 0.04703 |
D5YP4-3 | 0.20292 | 0.14723 | 0.07778 | 0.04303 | 0.24488 | 0.08241 |
D5YP5-1 | 0.12585 | 0.11091 | 0.04334 | 0.0667 | 0.11672 | 0.06435 |
D5YP5-2 | 0.05202 | 0.04499 | 0.01392 | 0.02668 | 0.02703 | 0.03691 |
D5YP5-3 | 0.0627 | 0.0495 | 0.02202 | 0.0209 | 0.16673 | 0.02881 |
D5YP6-1 | 0.05943 | 0.04892 | 0.02641 | 0.01796 | 0.06576 | 0.02534 |
D5YP6-2 | 0.08646 | 0.07716 | 0.02567 | 0.07512 | 0.19755 | 0.09069 |
D5YP6-3 | 0.1498 | 0.1089 | 0.0587 | 0.03587 | 0.39988 | 0.07066 |
D5YP7-1 | 0.04234 | 0.02448 | 0.02187 | 0.0054 | 0.32279 | 0.04489 |
D5YP7-2 | 0.02233 | 0.01456 | 0.01063 | 0.00569 | 0.20359 | 0.03308 |
D5YP7-3 | 0.05943 | 0.04892 | 0.02641 | 0.01796 | 0.06576 | 0.02534 |
D5YP8-1 | 0.08646 | 0.07716 | 0.02567 | 0.07512 | 0.19755 | 0.09069 |
D5YP8-2 | 0.1498 | 0.1089 | 0.0587 | 0.03587 | 0.39988 | 0.07066 |
D5YP8-3 | 0.04234 | 0.02448 | 0.02187 | 0.0054 | 0.32279 | 0.04489 |
D5YP9-1 | 0.02233 | 0.01456 | 0.01063 | 0.00569 | 0.20359 | 0.03308 |
D5YP9-2 | 0.02638 | 0.01781 | 0.01764 | 0.01051 | 0.1764 | 0.04667 |
D5YP9-3 | 0.04781 | 0.0369 | 0.023 | 0.02459 | 0.42332 | 0.06618 |
D5YP10-1 | 0.02744 | 0.01904 | 0.01882 | 0.01153 | 0.1601 | 0.03803 |
D5YP10-2 | 0.02656 | 0.01776 | 0.01678 | 0.00713 | 0.17441 | 0.03989 |
D5YP10-3 | 0.03258 | 0.02194 | 0.02225 | 0.00513 | 0.19853 | 0.02527 |
D6YP1-1 | 0.05348 | 0.04235 | 0.01289 | 0.00366 | 0.08298 | 0.01314 |
D6YP1-2 | 0.10796 | 0.07701 | 0.02608 | 0.00651 | 0.20752 | 0.04086 |
D6YP1-3 | 0.06804 | 0.05061 | 0.01653 | 0.00418 | 0.16127 | 0.0297 |
D6YP2-1 | 0.07973 | 0.0575 | 0.02295 | 0.00193 | 0.1816 | 0.01905 |
D6YP2-2 | 0.02395 | 0.02332 | 0.01026 | 0.00751 | 0.086 | 0.03351 |
D6YP2-3 | 0.1078 | 0.08041 | 0.02045 | 0.00931 | 0.25514 | 0.07536 |
D6YP3-1 | 0.16883 | 0.12768 | 0.04446 | 0.01012 | 0.18353 | 0.05078 |
D6YP3-2 | 0.03811 | 0.0336 | 0.01102 | 0.00838 | 0.06324 | 0.03079 |
D6YP3-3 | 0.13154 | 0.09113 | 0.03044 | 0.00547 | 0.35945 | 0.05266 |
D6YP4-1 | 0.10439 | 0.07674 | 0.02551 | 0.00228 | 0.22212 | 0.02632 |
D6YP4-2 | 0.02966 | 0.02931 | 0.01139 | 0.00866 | 0.08187 | 0.03972 |
D6YP4-3 | 0.11724 | 0.08606 | 0.02434 | 0.00825 | 0.13049 | 0.03288 |
D6YP5-1 | 0.05708 | 0.04879 | 0.01381 | 0.01517 | 0.01199 | 0.03031 |
D6YP5-2 | 0.03981 | 0.03213 | 0.00595 | 0.00563 | 0.024 | 0.01281 |
D6YP5-3 | 0.00736 | 0.00744 | 0.00228 | 0.00274 | 0.14909 | 0.02556 |
D6YP6-1 | 0.03718 | 0.03072 | 0.00835 | 0.00292 | 0.0294 | 0.0094 |
D6YP6-2 | 0.04269 | 0.02563 | 0.0129 | 0.00165 | 0.1597 | 0.02212 |
D6YP6-3 | 0.08828 | 0.05934 | 0.01547 | 0.00215 | 0.38835 | 0.02412 |
D6YP7-1 | 0.05209 | 0.03461 | 0.01104 | 0.00215 | 0.28096 | 0.03988 |
D6YP7-2 | 0.05158 | 0.03207 | 0.01013 | 0.00264 | 0.30859 | 0.04272 |
D6YP7-3 | 0.03718 | 0.03072 | 0.00835 | 0.00292 | 0.0294 | 0.0094 |
D6YP8-1 | 0.04269 | 0.02563 | 0.0129 | 0.00165 | 0.1597 | 0.02212 |
D6YP8-2 | 0.08828 | 0.05934 | 0.01547 | 0.00215 | 0.38835 | 0.02412 |
D6YP8-3 | 0.05209 | 0.03461 | 0.01104 | 0.00215 | 0.28096 | 0.03988 |
D6YP9-1 | 0.05158 | 0.03207 | 0.01013 | 0.00264 | 0.30859 | 0.04272 |
D6YP9-2 | 0.04418 | 0.02546 | 0.00922 | 0.00147 | 0.20001 | 0.03113 |
D6YP9-3 | 0.03638 | 0.02299 | 0.00666 | 0.00039 | 0.31614 | 0.02562 |
D6YP10-1 | 0.02969 | 0.04283 | 0.00505 | 0.04396 | 0.14356 | 0.06491 |
D6YP10-2 | 0.02871 | 0.02411 | 0.00782 | 0.00572 | 0.15739 | 0.04436 |
D6YP10-3 | 0.01984 | 0.01458 | 0.00628 | 0.00093 | 0.17616 | 0.02171 |
D7YP1-1 | 0.02692 | 0.02187 | 0.00822 | 0.00754 | 0.02555 | 0.02032 |
D7YP1-2 | 0.02645 | 0.02363 | 0.0097 | 0.0074 | 0.05327 | 0.02282 |
D7YP1-3 | 0.07286 | 0.05166 | 0.0131 | 0.00878 | 0.08689 | 0.02909 |
D7YP2-1 | 0.03814 | 0.02962 | 0.01204 | 0.00759 | 0.09519 | 0.03854 |
D7YP2-2 | 0.20016 | 0.17207 | 0.04393 | 0.04629 | 0.04212 | 0.06953 |
D7YP2-3 | 0.03352 | 0.02614 | 0.00862 | 0.008 | 0.1238 | 0.0439 |
D7YP3-1 | 0.04508 | 0.03253 | 0.01041 | 0.00595 | 0.14172 | 0.0209 |
D7YP3-2 | 0.15413 | 0.09291 | 0.03021 | 0.01185 | 0.34282 | 0.06181 |
D7YP3-3 | 0.05645 | 0.04456 | 0.01448 | 0.0099 | 0.0804 | 0.04408 |
D7YP4-1 | 0.07666 | 0.05692 | 0.01524 | 0.05481 | 0.21794 | 0.08633 |
D7YP4-2 | 0.14207 | 0.10034 | 0.0303 | 0.01027 | 0.15727 | 0.04022 |
D7YP4-3 | 0.1328 | 0.12076 | 0.04006 | 0.04428 | 0.17335 | 0.07534 |
D7YP5-1 | 0.08788 | 0.06016 | 0.01497 | 0.00542 | 0.09872 | 0.02238 |
D7YP5-2 | 0.01664 | 0.01167 | 0.00278 | 0.00168 | 0.16591 | 0.02152 |
D7YP5-3 | 0.13792 | 0.1058 | 0.03979 | 0.01358 | 0.08637 | 0.03336 |
D7YP6-1 | 0.06704 | 0.04478 | 0.01362 | 0.00367 | 0.19577 | 0.03014 |
D7YP6-2 | 0.10989 | 0.07425 | 0.02263 | 0.00739 | 0.29172 | 0.04205 |
D7YP6-3 | 0.07005 | 0.04272 | 0.01391 | 0.00214 | 0.30298 | 0.04212 |
D7YP7-1 | 0.18129 | 0.22104 | 0.21365 | 0.36247 | 0.54057 | 0.38038 |
D7YP7-2 | 0.13792 | 0.1058 | 0.03979 | 0.01358 | 0.08637 | 0.03336 |
D7YP7-3 | 0.06704 | 0.04478 | 0.01362 | 0.00367 | 0.19577 | 0.03014 |
D7YP8-1 | 0.10989 | 0.07425 | 0.02263 | 0.00739 | 0.29172 | 0.04205 |
D7YP8-2 | 0.07005 | 0.04272 | 0.01391 | 0.00214 | 0.30298 | 0.04212 |
D7YP8-3 | 0.18129 | 0.22104 | 0.21365 | 0.36247 | 0.54057 | 0.38038 |
D7YP9-1 | 0.03862 | 0.02573 | 0.01052 | 0.00089 | 0.18869 | 0.02811 |
D7YP9-2 | 0.02662 | 0.01576 | 0.00472 | 0.00131 | 0.15501 | 0.0223 |
D7YP9-3 | 0.00826 | 0.00773 | 0.0048 | 0.00282 | 0.04177 | 0.01449 |
D7YP10-1 | 0.05086 | 0.02783 | 0.01095 | 0.0013 | 0.34475 | 0.03973 |
D7YP10-2 | 0.03712 | 0.01818 | 0.00634 | 0.00014 | 0.29957 | 0.0246 |
D7YP10-3 | 0.01442 | 0.01096 | 0.00466 | 0.00461 | 0.13774 | 0.0326 |
结合表2,由采样的210份样本数据可见,所有样品在波数1080cm-1左右、波数1650cm-1和波数1535cm-1左右均出现了较强的吸收峰,在波数为1315cm-1、1392cm-1、1452cm-1、2854cm-1、2873cm-1、2935cm-1、2960cm-1、3012cm-1、3289cm-1左右处均有吸收峰出现。绝大部分样本在波数1650cm-1和波数1535cm-1左右的两个吸收峰峰高相近,绝大部分样本在波数1650cm-1处的峰强度大于在波数1535cm-1处的峰强度。绝大部分样本在波数1080cm-1左右处的吸收峰强度最大,所有样本在波数3289cm-1处有中等强度的吸收峰,该处的峰强度小于波数1080cm-1、波数1650cm-1和波数1535cm-1处的吸收峰并大于波数4000cm-1-1300cm-1的官能团区的其他峰值。
说明衰减全反射红外光谱法采集到的血痕红外光谱数据具备相对稳定可用的条件;另一方面,波数为1536cm-1处的峰强度数值相对于波数为1650cm-1处的峰强度数值符合y=0.7744x-0.002942(R2=0.9281)的线性方程,说明酰胺I和酰胺II两处峰值之间的相对关系受到的样本不均匀、实验操作差异等因素的影响较小,有利于针对重点峰位进行细化研究。
4.使用神经网络建立样本痕迹的红外光谱数据与遗留时间的对应关系模型并验证精度,
利用MATLAB数据分析软件中的神经网络工具箱功能,以痕迹样本的遗留时间和红外光谱数据为变量构建BP神经网络,训练类型为TRAINLM即L-M优化算法;误差类型选择MSE;层数选2层(1个隐含层和1个输出层),隐含层和输出层的传递函数选用系统默认的TANSIG,因L-M优化算法需占用大量内存资源,过大的神经元数对计算机要求较高,以神经元个数分别设定为1,3,5构建三个神经网络network_all_trainlm_n1、network_all_trainlm_n3、network_all_trainlm_n5。
利用上述训练集对已建立结构的三个神经网络分别进行训练,训练参数中设置time为180(network_all_trainlm_n5的训练时间为300),max_fail数值设置为20,其余为默认参数,以变量AAall为输入数据inputs,以TrainDay为目标数据Targets,对神经网络进行训练,待神经网络按照预设参数训练完成后,以APall为输入数据inputs对其对应天数进行预测,预测天数分别保存为变量network_all_trainlm_n1_outputs、network_all_trainlm_n3_outputs、network_all_trainlm_n5_outputs。
利用上述验证集进行验证,将验证集中的红外光谱数据录入上述训练得到的三组模型,三组模型反馈其模拟的验证集中对应遗留天数。
利用数据处理工具比较三组模型模拟的验证集中对应遗留天数与真实的遗留天数,计算其判定系数R2,分别为0.7071、0.8435、0.5677,因此,使用神经网络network_all_trainlm_n3建立的模型可用。
本发明中的判定系数R2又叫决定系数,其主要作用是评估回归模型对因变量产生变化的解释程度,是判断模型精确程度的指标,一般认为,判定系数R2大于0.8表示模型拟合度好,其计算公式如下:
R2=回归平方和/总平方和,其中回归平方和=总平方和-残差平方和。
5.将检材痕迹的红外光谱数据录入符合精度要求的模型,得到模拟检材痕迹的遗留时间
将检材痕迹的红外光谱数据录入上述神经网络network_all_trainlm_n3建立的模型,模型反馈模拟的检材痕迹遗留时间D=5.5。
6.评估检材痕迹遗留时间
模型反馈模拟的检材痕迹遗留天数D=5.5,其与样本痕迹的遗留时间最大值7的比值≤0.9,因此,检材痕迹的遗留天数评估为5.5。
Claims (5)
1.一种基于衰减全反射红外光谱的痕迹遗留时间的测量方法,其特征在于,包括以下步骤:
1)使用衰减全反射红外光谱法采集检材痕迹的红外光谱数据,所述检材痕迹来自于案发现场,检材痕迹为血痕、笔迹、汗斑或指印;
2)模拟检材痕迹制作样本痕迹,分类样本,
3)采集样本数据集,所述数据集包括不同遗留时间样本痕迹的红外光谱数据与对应的遗留时间,
4)使用神经网络建立样本痕迹的红外光谱数据与遗留时间的对应关系模型并验证精度,
5)将检材痕迹的红外光谱数据录入符合精度要求的模型,得到模拟检材痕迹的遗留时间,
6)评估检材痕迹遗留时间,比较模拟的检材痕迹的遗留时间和样本痕迹的遗留时间最大值,如比值≤0.9,则评估模拟的检材痕迹的遗留时间即为测量的检材痕迹遗留时间,测量完成,如比值>0.9,则扩充样本数据集,延长样本痕迹的遗留时间,重复步骤3)至步骤5),比较扩充样本数据集后模拟的检材痕迹的遗留时间和样本痕迹的遗留时间最大值,直至比值≤0.9,则评估模拟的检材痕迹的遗留时间即为测量的检材痕迹遗留时间,测量完成。
2.根据权利要求1所述基于衰减全反射红外光谱的痕迹遗留时间的测量方法,其特征在于,步骤2)中,使用与检材痕迹载体相同的材料为载体制作样本痕迹,将样本痕迹按照7:3比例,随机划分为训练用样本痕迹和验证用样本痕迹。
3.根据权利要求1所述基于衰减全反射红外光谱的痕迹遗留时间的测量方法,其特征在于,步骤3)中,模拟检材痕迹的遗留环境,使用衰减全反射红外光谱法,要求波数范围和分辨率与步骤1)相同,以均匀的时间间隔,采集不同遗留时间的样本痕迹的红外光谱数据,与对应的遗留时间归纳为样本数据集,记录样本痕迹的遗留时间最大值,保留样本数据集中红外光谱数据与遗留时间的对应关系,将样本数据集按照训练用和验证用,划分为训练集和验证集。
4.根据权利要求3所述基于衰减全反射红外光谱的痕迹遗留时间的测量方法,其特征在于,步骤4)中,以训练集中的红外光谱数据和遗留时间为变量,训练神经网络,建立样本痕迹的红外光谱数据与遗留时间的对应关系模型,使用验证集验证模型精度,如判定系数R2>0.8,则该模型可用,否则重复“训练神经网络”至“使用验证集验证模型精度”的流程,直至判定系数R2>0.8。
5.根据权利要求1所述基于衰减全反射红外光谱的痕迹遗留时间的测量方法,其特征在于,步骤5)中,向可用模型输入检材痕迹的红外光谱数据,得到模拟的检材痕迹的遗留时间。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010452130.4A CN111366554B (zh) | 2020-05-26 | 2020-05-26 | 一种基于衰减全反射红外光谱的痕迹遗留时间的测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010452130.4A CN111366554B (zh) | 2020-05-26 | 2020-05-26 | 一种基于衰减全反射红外光谱的痕迹遗留时间的测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111366554A CN111366554A (zh) | 2020-07-03 |
CN111366554B true CN111366554B (zh) | 2020-08-11 |
Family
ID=71207539
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010452130.4A Active CN111366554B (zh) | 2020-05-26 | 2020-05-26 | 一种基于衰减全反射红外光谱的痕迹遗留时间的测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111366554B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112098362B (zh) * | 2020-09-15 | 2024-02-27 | 南京工程学院 | 一种基于近红外频谱特征的松质骨离体时间评估方法及系统 |
CN113293217B (zh) * | 2021-05-27 | 2022-02-18 | 山西医科大学 | 一种鉴定血痕遗留时间的方法及其mRNA组合 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090036781A1 (en) * | 2007-07-31 | 2009-02-05 | Kei Utsugi | External condition control device based on measurement of brain functions |
US20190285550A1 (en) * | 2015-07-30 | 2019-09-19 | The Research Foundation For The State University Of New York | Spectroscopic methods for body fluid age determination |
-
2020
- 2020-05-26 CN CN202010452130.4A patent/CN111366554B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090036781A1 (en) * | 2007-07-31 | 2009-02-05 | Kei Utsugi | External condition control device based on measurement of brain functions |
US20190285550A1 (en) * | 2015-07-30 | 2019-09-19 | The Research Foundation For The State University Of New York | Spectroscopic methods for body fluid age determination |
Also Published As
Publication number | Publication date |
---|---|
CN111366554A (zh) | 2020-07-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Baker et al. | Investigating FTIR based histopathology for the diagnosis of prostate cancer | |
CN111366554B (zh) | 一种基于衰减全反射红外光谱的痕迹遗留时间的测量方法 | |
CN109799195B (zh) | 一种激光诱导击穿光谱的高精度定量分析方法 | |
CN104198432B (zh) | 应用近红外光谱技术无损鉴别钢结构防火涂料品牌的方法 | |
CN109064553B (zh) | 基于近红外光谱分析的实木板材节子形态反演方法 | |
CN106124508A (zh) | 一种文物鉴定系统及鉴定方法 | |
CN111553877B (zh) | 基于太赫兹陶瓷基复合材料叶片损伤识别和寿命评估方法 | |
CN113008805B (zh) | 基于高光谱成像深度分析的白芷饮片质量预测方法 | |
CN109669023A (zh) | 一种基于多传感器融合的土壤属性预测方法 | |
CN106841083A (zh) | 基于近红外光谱技术的芝麻油品质检测方法 | |
CN102841070B (zh) | 由二维相关红外同步光谱识别原油种类的方法 | |
CN110068544A (zh) | 物质识别网络模型训练方法及太赫兹光谱物质识别方法 | |
Thickett et al. | Using non-invasive non-destructive techniques to monitor cultural heritage objects | |
CN110208238A (zh) | 一种基于svm模型结合图像对肺癌组织的精确定位方法 | |
CN110082310A (zh) | 一种橡胶树叶片氮素含量的近红外波段高光谱诊断方法 | |
CN110441254A (zh) | 一种用于识别塑料的近红外光频梳光谱仪 | |
CN110070004B (zh) | 一种应用于深度学习的近地高光谱数据扩展方法 | |
CN108548794A (zh) | 一种生物制品近红外模型转移方法 | |
CN110567894A (zh) | 一种含水溶液pH值快速检测方法和检测装置 | |
CN108489928B (zh) | 一种短波红外消光光谱纺织品纤维成分检测方法 | |
CN112716447A (zh) | 一种基于拉曼检测光谱数据深度学习的口腔癌分类系统 | |
CN111220575A (zh) | 一种基于太赫兹近场光谱的细胞检测方法 | |
JPH02271211A (ja) | 塗装鮮映性評価方法 | |
CN115541621A (zh) | 一种基于微波扫频技术和Stacking集成模型的茶叶籽仁含水率检测方法及系统 | |
CN109659031A (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 |