基于MC模型和Nelder-Mead单纯形算法的时域组织光谱学

Time⁃Resolved Spectroscopy Based on Monte Carlo Model and Nelder⁃Mead Simplex Algorithm

生物组织光学参数的无创在体测量是近红外光谱学(NIRS)研究的基础课题之一。在已发展的NIRS方案中,时间分辨测量具备优良的同时反演吸收系数和散射系数的能力,而且近年来该技术的性价比显著提高,获得了更多关注,但其针对分层组织的参数反演在准确性、稳定性和反演速度等方面尚存在一定的局限性。为此,本文提出了一种时域蒙特卡罗模型支持的Nelder-Mead单纯形算法,该算法利用双源探距下的时间分辨漫反射光信号,以循环迭代方式设置不同的基向量,经线性变换后无须求导便可启发式搜索目标函数的最优值,从而实现了高信效度的分层光学参数反演。模拟实验与仿体实验均表明:在组织光学参数反演方面,所提算法的精度优于传统算法,而且该算法具有良好的噪声鲁棒性和临床适应性,为生物组织光学参数的在体测量提供了新方法。

Objective Changes in optical parameters can reflect the physiological status of biological tissue and constitute a fundamental and important topic in the field of near-infrared spectroscopy.Compared with the continuous-wave and frequency-domain measurement methods,the time-domain measurement method has the best performance in distinguishing and separating absorption and scattering coefficients,particularly in single-point measurement scenarios.Consequently,the time-domain measurement method is more commonly used to measure changes in optical parameters,also known as time-resolved spectroscopy.Currently,the biological tissue model used for time-resolved spectroscopy inversion schemes often assumes that the biological tissue is a single-layer biological tissue model,which hypothesizes that the optical properties are identical throughout the tissue.Although the single-layer biological tissue model simplifies the complexity of light propagation models and inversion algorithms,it is not very suitable for representing the structure of most human biological tissues;biological tissues at different depths exhibit a layered structure owing to variations in structure and function.Considering the high computational complexity and marginal improvement in accuracy associated with multilayer models,recent research has increasingly focused on the double-layer biological tissue model.Currently,the commonly used double-layer biological tissue model parameter inversion methods face several challenges.First,they require a substantial amount of experimental data from multiple sources and detector separation,resulting in extended overall measurement times.Second,considerable time and effort are required to construct precise databases.Third,the need for iterative differentiation leads to prolonged computation time and lower accuracy.Finally,these methods struggle to handle complex scenarios,such as those in which both layers of tissue parameters are entirely unknown.To address these issues,this study introduces the Nelder-Mead simplex algorithm for the first time into the framework of Monte Carlo model-based tissue optical parameter inversion,developing a time-domain Monte Carlo supported Nelder–Mead simplex(MC-NMS)inversion algorithm.Methods The Monte Carlo model can use simulations to customize the optical parameters based on the structural characteristics of layered tissues,which are often employed as transport models for double-layer biological tissue model parameter inversion.This study introduces the Nelder–Mead simplex algorithm into Monte Carlo-based tissue optical parameter inversion for the first time.By utilizing only two source and detector separation time-domain diffuse reflectance data,the need to construct extensive databases in advance is eliminated.Initially,we empirically set the parameter values to find a feasible solution within the feasible region.The variables were incrementally adjusted through a cyclic iterative approach involving the establishment of different base vectors.Through matrix linear transformations,the optimal value of the objective function was determined using a heuristic search method that obviated the need for differentiation.Ultimately,this approach achieves a high-fidelity inversion of the optical parameters in layered tissues under complex conditions.Results and Discussions The numerical simulation experiments demonstrate that in the single-layer tissue optical parameter inversion application scenario,the proposed MC-NMS method yields significantly superior results compared to the TDIA and SDIA methods(Fig.3).Additionally,when the source-detector separation is set to 3 mm,the MC-NMS method yields the best results.(Fig.4).For double-layer biological tissue model optical parameter inversion,the results reveal that changing the upper-layer tissue thickness also requires different optimal source–detector separations(Fig.5,Fig.6).Moreover,the inversion errors obtained using the MC-NMS method are lower than those obtained using the TDIA and SDIA methods.Experimental validation was conducted using a multi-wavelength and multi-source-detector separation time-resolved measurement system based on time-correlated singlephoton counting technology developed by our research group.Data were computed using the 25%rising edge to 20%falling edge time channels of the measured time-point spread function curves.The results indicate that the MC-NMS method achieves inversion errors of 11.64%and 0.89%in the single-layer biological tissue model forμa andμ's,respectively.In the double-layer biological tissue model,the inversion errors forμa1,μa2,andμ's are 28%,20%,and 14%,respectively,representing an overall improvement in the parameter inversion accuracy of approximately 17%compared to the TDIA method.Conclusions The results of the liquid phantom experiments for the single-layer biological tissue and double-layer biological tissue models consistently demonstrate that the MC-NMS method outperforms other approaches.Notably,when theμa is relatively larger,the inversion errors are smaller.The results indicate that the MC-NMS can not only maintain high inversion accuracy in the straightforward application scenario of the single-layer biological tissue model but also ensure the accuracy of complex scenarios,such as those in which all optical parameters of the double-layer biological tissue model layers are unknown.The MC-NMS can be an effective method for clinical applications of time-resolved tissue oxygen measurement and imaging instruments,providing a pragmatic tool for enhancing accuracy in clinical applications.