电子学报 ›› 2022, Vol. 50 ›› Issue (8): 1985-1991.DOI: 10.12263/DZXB.20210944

• 学术论文 • 上一篇    下一篇

基于蒙特卡罗法硬膜血肿厚度计算的数值仿真研究

范鑫燕1, 罗海军1,2(), 李妍妍1, 向洋1, 罗霞1, 覃睿1, 郭盼1   

  1. 1.重庆师范大学光电功能材料重庆市重点实验室, 重庆 401331
    2.重庆国家应用数学中心, 重庆 401331
  • 收稿日期:2021-07-17 修回日期:2021-11-07 出版日期:2022-08-25
    • 通讯作者:
    • 罗海军
    • 作者简介:
    • 范鑫燕 女,1996年9月出生于四川省绵阳市.重庆师范大学物理与电子工程学院研究生,主要研究方向为近红外功能成像及系统设计.E-mail: y_4099@126.com
      罗海军 男,1983 年 10 出生于四川省资阳市 . 博士、重庆师范大学教授、硕士生导师,主要研究方向是磁感应成像、近红外成像和电阻抗成像. E-mail: luohaijun@cqnu.edu.cn
      郭 盼 女,1987年8月出生,湖北天门人.博士、重庆师范大学副教授、硕士生导师,主要研究方向为电磁测量技术、便携式磁共振技术.E-mail: guopan@cqnu.edu.cn
    • 基金资助:
    • 国家自然科学基金 (51507023); 重庆市自然科学基金 (面上项目) (cstc2020jcyj-msxmX0726); 重庆市教委重点项目 (KJZD-K202100506); 重庆师范大学研究生科研创新项目 (YKC21046); 国家级大学生创新创业训练计划项目 (20210637024)

Numerical Simulation Study on Thickness Calculation of Dural Hematoma Based on Monte Carlo Method

FAN Xin-yan1, LUO Hai-jun1,2(), LI Yan-yan1, XIANG Yang1, LUO Xia1, QIN Rui1, GUO Pan1   

  1. 1.Chongqing Key Laboratory of Optoelectronic Functional Materials,Chongqing Normal University,Chongqing 401331,China
    2.Chongqing National Center for Applied Mathematics,Chongqing 401331,China
  • Received:2021-07-17 Revised:2021-11-07 Online:2022-08-25 Published:2022-09-08
    • Corresponding author:
    • LUO Hai-jun

摘要:

为实现对颅脑硬膜血肿的快速检测与评估,本研究基于近红外光谱技术(Near Infra-Red Spectroscopy,NIRS)与蒙特卡罗算法建立适用于正问题的5层颅脑组织进行仿真,采用控制变量法以10组血肿厚度模型的出射光子数量为基础,定义源-检测器灵敏度(Source-Detector Sensitivity,SDS)变量参数,构造基于幂律衰减模型和指数律衰减模型的正向函数矩阵.对正问题获取的数据进行数学模型重构即逆问题的实现,以6组血肿厚度进行逆向理论计算,对比参照值,血肿厚度小于0.7 cm,两种函数模型平均绝对误差均小于3.6%;血肿厚度为0.75 cm,幂律衰减模型平均绝对误差比指数律衰减模型小于4.3%,误差小于6.388 6%.结果显示该方法对准确判断是否含有血肿以及血肿厚度预测具有可行性,检测灵敏度与检测距离具有对数相关性,幂律衰减模型比指数律衰减模型预测平均值更接近参照值,构建效果更好.

关键词: 近红外光谱技术, 蒙特卡罗, 控制变量法, 灵敏度, 函数矩阵, 逆问题

Abstract:

To realize the rapid detection and evaluation of dural hematoma, this study is based on the near infra-red spectroscopy(NIRS) and Monte Carlo algorithm to establish a simulation of 5 layers of brain tissue suitable for the positive problem. Source-detector sensitivity(SDS) parameters are defined to construct a forward function matrix based on the power-law attenuation model and the exponential model. These parameters are given by the number of transmitted photons from 10 groups of hematoma thickness models via the controlled variable method. The mathematical model reconstruction of the data obtained from the forward problem is the realization of the inverse problem. The 6 groups of hematoma thicknesses are used for reverse theoretical calculations. Compared with the reference value, the thickness of hematoma is less than 0.7 cm, and the average absolute error of the two function models is less than 3.6%. The thickness is equal to 0.75 cm, the average absolute error of the power-law attenuation model is less than 4.3%, and the error is 6.3886% less than the exponential law attenuation model. The results show that the method is feasible for accurately determining whether hematoma is contained and predicting the thickness of hematoma. The detection sensitivity and detection distance have a logarithmic correlation. The power-law attenuation model is closer to the reference value than the exponential law attenuation model, and the build effect is better.

Key words: near infra-red spectroscopy, Monte Carlo, control variable, sensitivity, function matrix, inverse problem

中图分类号: