人骨遗骸死后间隔评估的任务特异性多模态成像与光谱学：一种探索性决策支持框架

《Diagnostics》：Task-Specific Multimodal Imaging and Spectroscopy for Post-Mortem Interval Assessment in Human Skeletal Remains: An Exploratory Decision-Support Framewor

　　

背景/目的：人骨遗骸的死后间隔（post-mortem interval，PMI）估计具有挑战性，因为骨组织经历的结构、分子和光学变化随时间演变不同，并受埋藏学条件强烈影响。本研究不假设单一多模态分类器在所有PMI区间表现相同，而是旨在确定哪种成像或光谱模态对特定法医决策任务最具信息量，并基于样本水平内部评估推导探索性任务特异性诊断决策支持框架。方法：人股骨骨样本被分为0–2周至>100年的5个PMI类，采用微计算机断层扫描（micro-CT）、高光谱成像（HSI）、手持与显微拉曼光谱（Raman spectroscopy）及NIR-ONE近红外光谱检测；重复采集在物理样本水平聚合。定义四类目标诊断对比：考古类5对比类1–4、早期类1+2对比晚期类4+5、拉曼兼容法医用类1+2对比类4、类1对比类2；另增类4对比类5探索性直接配对分析。采用参数级ROC分析、bootstrap置信区间、约95%特异度操作点及重复分层交叉验证，所有多变量建模预处理均在折叠内完成。结果：诊断性能具强任务依赖性。类5与类1–4区分中，micro-CT Mean2（AUC 0.984）、NIR-ONE 1944 nm反射率（AUC 0.980）、HSI组织水指数（TWI；AUC 0.961）最优；直接C4对C5中micro-CT Mean2完全分离（AUC 1.000），NIR-ONE 1944 nm仍优（AUC 0.963），HSI TWI仅中等（AUC 0.742）。早期类1+2对类4+5中HSI StO2指数最高（AUC 0.940），TWI（AUC 0.894）与NIR 1550–1950 nm斜率（AUC 0.860）次之；拉曼碳酸/磷酸与结晶度参数提供互补分子信息；类1对类2仅为中等；校正五类模型平衡准确率中等且未优于单独NIR-ONE。结论：micro-CT Mean2与高波长NIR-ONE特征对考古兼容C5材料区分最强，HSI光学指数对早晚期对比最有用，拉曼在C1–C4内提供互补分子信息；因C5仅6个物理样本且考古背景与年代混淆，高AUC为探索性假设生成估计；该框架需独立外部验证方可法医学应用。

研究背景方面，人骨遗骸的死后间隔（post-mortem interval，PMI）估计是法医医学、法医人类学与死亡调查核心难题。软组织的指标消失后，骨组织受生物、化学、物理与环境因子影响，发生分解、成岩作用、矿物溶解重结晶、胶原降解、水合变化、微生物活动、土壤成分、pH、温湿度、埋葬与风化等改变，使跨早期法医、晚期法医与考古时间间隔的可靠定年长期无通用方法。以往微计算机断层扫描（micro-CT）反映结构密度变化，拉曼光谱（Raman spectroscopy）反映矿物与基质分子变化，近红外（near-infrared，NIR）光谱反映水合与有机矩阵宽谱带，高光谱成像（hyperspectral imaging，HSI）融合空间与光谱生成设备派生指数；但尚不清楚这些模态诊断价值是否在整个PMI范围均匀。研究人员因此开展回顾性多模态研究，在唯一物理骨样本水平再分析，以任务特异性二元对比替代单一全局多类分类，假设性能随法医问题变化，并构建探索性分层工作流。该研究论文发表于《Diagnostics》。

关键技术方法上，研究人员整合micro-CT、手持Mira与显微Senterra拉曼、NIR-ONE（1550–1950 nm）与HSI四种模态；样本为法医与考古来源股骨骨干皮质骨，主队列107个唯一物理样本（C1=33、C2=47、C3=11、C4=10、C5=6），重复测量在样本水平聚合；统计采用Kruskal–Wallis与Mann–Whitney U、ROC与AUC、20000次分层bootstrap、约95%特异度操作点、5折20次重复分层交叉验证，折叠内拟合缩放/SNV/PCA，线性SVM固定C=1.0，不做嵌套调参；拉曼仅用于C1–C4，因C5荧光过强。

研究结果如下。数据集结构与多模态重叠：107个唯一样本跨五类，模态可用数不同，多模态仅取完整数据交集，五类融合模型排除10例（9.3%），早晚期对比87例完整三模态；所有分析均在物理样本水平，无外部验证队列。考古兼容类5与类1–4区分：micro-CT Mean2（AUC 0.984）、NIR-ONE 1944 nm（AUC 0.980）、HSI TWI（AUC 0.961）最强；约95%特异度下Mean2敏感1.000；探索性C4对C5中Mean2完全分离（AUC 1.000，C5仅5例可算），NIR-ONE 1944 nm（AUC 0.963）仍优，HSI TWI仅0.742，说明池化高AUC不等于C4–C5边界同样强。早对晚PMI区分：HSI StO₂（AUC 0.940）最高，TWI（0.894）与NIR斜率R1950–R1550（0.860）补充；该对比排除过渡类C3，不能外推连续PMI人群。HSI光学指数与灌注/THI关系：NIR perfusion、StO₂、组织血红蛋白指数（THI）、TWI均为设备派生光学指数而非生理量；Perfusion/(Perfusion+THI)有方向趋势但弱于StO₂与TWI。拉曼分子参数（C1–C4）：Mira结晶度（1/FWHM958）AUC 0.854、碳酸/磷酸A1070/A958 0.849，Senterra同比0.849；类1对类2仅中等（Senterra Amide III 0.787）。NIR-ONE签名：1944 nm反射率区分C5最强，1550–1950 nm斜率助早晚期；单点非孤立而是高波长区。跨方法比较：无单一模态全域最优；C1+C2对C4内HSI平衡准确率0.828高于Mira拉曼0.764与NIR-PCA10 0.726。次级全局五类建模：NIR-ONE SNV+PCA平衡准确率约0.483，HSI+micro-CT+NIR-ONE约0.442，融合未胜单模态。任务特异性决策支持框架：未知来源骨先保守处理与非破坏HSI/NIR筛查；考古兼容判断用NIR-ONE 1944 nm与micro-CT Mean2，TWI补充；早晚期用HSI StO₂/TWI；拉曼做C1–C4分子刻画；类1对类2概率化；结果不一致或质量差则判不确定；框架为探索性而非验证分类器。

讨论部分总结：研究人员强调诊断价值强任务依赖，不可由单一全局五类模型替代。C5不仅代表年代，还混杂考古与埋藏学背景，10–100年区间缺失使连续阈值无法建立。HSI指数在干骨中不是生理氧合/灌注，而是表面光学与矩阵改变派生量。micro-CT Mean2与高波长NIR捕获老骨密度/结构/水相关信号；HSI StO₂/TWI对早晚期更敏感；拉曼受荧光限制只能法医用C1–C4。多模态融合未自动增效，价值在按问题选模。局限包括回顾性、C5仅6例、类别不平衡、模态协议不统一、供体年龄/性别/病变与角度元数据缺失、环境埋藏元数据不足、操作点未外部校准、无校准与跨设备重复、类1对类2中等。未来需前瞻多中心、平衡队列、统一协议、埋藏学元数据、10–100年区间样本、预设参数与独立外测。

结论部分翻译：在本内部评估队列中，多模态成像与光谱学对PMI评估的诊断价值具有强任务依赖性。micro-CT Mean2与高波长NIR-ONE反射率对可用考古兼容C5样本区别于池化C1–C4参考组尤其具信息量；在附加探索性C4对C5分析中，micro-CT Mean2在可用样本中完全观测分离，NIR-ONE 1944 nm反射率保持强区分，而HSI派生TWI对该具体边界仅中等表现。鉴于C5样本极少且年代与考古及埋藏学背景混淆，这些发现应视为假设生成而非经验证年代PMI性能。对刻意分离的早对晚对比，设备派生HSI StO₂与TWI指数最具信息量，拉曼矿物与基质特征在C1–C4内提供互补分子信息；排除过渡C3可放大表观分离，限制向连续或未选PMI人群外推；类1对类2区分中等，校正五类多模态融合未稳定优于最强单模态。综上，结果支持任务特异性模态选择与选择性分析升级，而非全模态一律使用；所提框架为内部分析衍生的假设生成探索性决策支持概念，非验证序贯分类器、法医指南或法律裁决规则，需前瞻多中心、更大平衡队列、埋藏学元数据、10–100年区间代表、统一协议与独立外测后方可进入法医办案。