Extending dynamic range of detector in non-contact diffuse optical tomography system using deep learning

利用深度学习扩展非接触扩散光学层析成像系统中探测器的动态范围 Extending dynamic range of detector in non-contact diffuse optical tomography system using deep learning 1 2 corresp first-author E-mail： zhaoxingtjnk@nankai.edu.cn E-mail： zhaoxingtjnk@nankai.edu.cn 赵　星（1980-），男，北京人，博士，教授，2007年于南开大学获得博士学位，主要从事光学信息处理和光学系统设计等方面的研究。E-mail：zhaoxingtjnk@nankai.edu.cn http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24007207&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24007215&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24007211&type= 22.01332855 32.00401306 赵　星 （1980-），男，北京人，博士，教授，2007年于南开大学获得博士学位，主要从事光学信息处理和光学系统设计等方面的研究。E-mail： zhaoxingtjnk@nankai.edu.cn zhaoxingtjnk@nankai.edu.cn 1 刘桐君（1995-），男，四川德阳人，硕士研究生，2018年于西北工业大学获学士学位，主要从事扩散光学层析成像方面的研究。E-mail：609345567@qq.com http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24007219&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24007228&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24007224&type= 22.01333618 32.00387573 刘桐君 （1995-），男，四川德阳人，硕士研究生，2018年于西北工业大学获学士学位，主要从事扩散光学层析成像方面的研究。E-mail： 609345567@qq.com 609345567@qq.com 1 1 1 2 南开大学 现代光学研究所，天津 300350 南开大学 现代光学研究所，天津 300350 Institute of Modern Optics， Nankai University， Tianjin 300350， China Institute of Modern Optics， Nankai University， Tianjin 300350， China 天津市微尺度光学信息技术科学重点实验室，天津 300350 天津市微尺度光学信息技术科学重点实验室，天津 300350 Tianjin Key Laboratory of Micro-scale Optical Information Science and Technology， Tianjin 300350， China Tianjin Key Laboratory of Micro-scale Optical Information Science and Technology， Tianjin 300350， China 非接触式扩散光学层析成像（Diffuse Optical Tomography，DOT）技术和传统DOT技术相比具有更高的空间采集密度，系统中探测器的动态范围对实验结果有很大的影响。通过仿真和实验研究了探测器动态范围过小对重建结果的影响，并提出了一种利用深度学习拓展探测器动态范围的方法。在NIRFSAT软件包中设置不同的吸收散射的仿体及入射光场参数，批量生成训练样本，搭建全连接网络进行模型训练，并利用该模型修复仿真和光学实验条件下探测器获得的数据。实验数据修复及重建结果表明，该模型修复了低动态范围探测器获得的数据，将探测器的动态范围从256扩展到10 9 。该方法能够有效地减小探测器动态范围过小带来的重建误差，为使用普通低动态范围探测器进行非接触式DOT提供了一种有效的技术解决方案。 Non-contact diffuse optical tomography （DOT） is a novel technique exhibiting higher spatial detection density compared to traditional DOT. The dynamic range of the detector in non-contact DOT considerably influences the accuracy of the reconstruction of the absorption coefficient. In this study， the effect of low dynamic range of the detector on the reconstruction results was examined by conducting experiments and simulations. A method based on deep learning was proposed to extend the dynamic range of the detector. The training samples were generated in the NIRFAST software by setting different phantom absorption and scattering coefficients and incident light field parameters. A fully connected network was built for model training. The model was used to recover the detection data in the simulations and optical experiments. The data recovery and reconstruction results of the experiments indicate that the model recovers the error data of the detector with a low dynamic range while extending the dynamic range of the detector from 256 to 10 9 . The proposed method can effectively reduce the reconstruction error caused by the low dynamic range of the detector and provide an effective technical solution for non-contact DOT using detectors having a low dynamic range. 扩散光学层析成像 深度学习 非接触系统 数据修复 diffuse optical tomography deep learning non-contact system data recovery 1 引　言 1 扩散光学层析成像（Diffuse Optical Tomography， DOT） ［ 1 1 ］ 是一种能够快速无损探测组织吸收散射系数分布的技术，被广泛应用于脑功能成像 ［ 2 2 - 3 3 2-3 ］ ，乳腺癌诊断 ［ 4 4 - 5 5 4-5 ］ ，关节成像 ［ 6 6 - 7 7 6-7 ］ ，农产品质量检测 ［ 8 8 ］ 等领域。其本质是求解一个包含光子传输模型的逆问题，具有内在的非线性及病态性 ［ 9 9 ］ 。而获取更多的信息用于重建能够有效地提高DOT重建图像的空间分辨率和对比度，并消除伪影。高密度DOT系统是最直接的方案，如Eggebrecht等人 ［ 2 2 ］ 在2014年搭建的接触式高密度DOT系统，该系统对人脑部进行功能成像并和功能核磁共振进行对比，展示了高密度DOT系统的有效性与准确性。虽然该系统具有很高的准确率，但系统组件非常多且复杂。非接触式DOT系统则能有效克服这个问题，其光源和探测器不再以光纤探头的方式接触式地分布于组织表面，光源是直接在空间中传输的激光束，探测器则使用EM-CCD ［ 10 10 - 11 11 10-11 ］ 等面成像装置，特点是能够一次性大面积地采集组织表面的能流分布。Takamasa ［ 12 12 ］ 等人使用纳秒级脉冲激光和时间提取的图像传感器制作了非接触式时域高密度DOT系统，该系统无需使用密集的光纤探头就可以轻松获取高密度测量点。使用这套系统分别对模拟头部吸收系数的仿体和真人的脑部进行了测量，研究了人脑在完成简单计算任务时的活动情况。当然该系统也存在着限制，探测到的部分光信号会被具有较高吸收系数和散射系数的组织或介质衰减，这意味着随着探测范围的加大，光的衰减将进一步加强，数据的动态范围也进一步变大。因此，在已有的非接触式DOT系统中，为了获得足够精确的探测数据，一般采用调整硬件参数的方法，如选择较小尺寸或较小吸收系数的探测目标以减小对探测器动态范围的需求，或选择使用较为昂贵的高动态范围CCD来满足探测需要。虽然这些方法保证了测量精度，但却限制了非接触DOT系统的应用空间及前景。使用算法修复有误差的数据则能够在保证DOT重建精度的同时，降低非接触DOT系统对于探测器动态范围的需求。 深度学习是机器学习的一种，适用于拥有大量数据的非线性问题的建模，已经应用在图像分辨率增强 ［ 13 13 ］ 、模糊图像恢复 ［ 14 14 ］ 、光纤传输后散斑恢复 ［ 15 15 ］ 、高光谱分类 ［ 16 16 ］ 和遥感图像检索 ［ 17 17 ］ 等领域，并展示了非常强的数据恢复能力和建模能力。本文采用深度学习的方法，使用可进行DOT计算和重建的NIRFAST开源程序对不同二维吸收散射特性的介质进行计算，通过数据截断方法模拟低动态范围探测器的信号获取过程，生成大量训练样本。在此基础上搭建了一维全连接网络并完成了模型训练，最后，通过仿真和实验方法对模型的性能进行了分析及验证。该方法能够有效地修复受探测器动态范围限制而获得的误差较大的探测数据，并成功重建出介质中吸收系数的分布特性。 2 基本原理 1 2.1　DOT简介及重建算法 2 DOT是一种利用光通过吸收散射介质表面能流密度分布和光子传输模型计算重建介质内部吸收系数和散射系数分布的技术。其基本的测量过程示意图如 图1 图1 所示。图中的介质具有一定的吸收散射特性，介质上表面的点表示光源的入射位置，下表面的点表示探测位置。当光依次从介质上表面的光源位置入射时，进入介质的光发生吸收及多次散射，最后从介质下表面出射。利用点探测器或面探测器探测从介质下表面出射的扩散光的能流密度，在得到各探测位置的一组数据后，利用DOT求解算法即可计算重建介质内部的吸收系数分布。 DOT求解算法包含正问题及反问题的求解。正问题求解指的是根据已知的光源探测器分布以及介质参数分布，求解光子传输模型以获取介质表面的能流密度分布。反问题求解指的是通过已知光源探测器分布、介质表面的出射光能流密度分布以及光子传输模型求解介质中的参数分布。由此可见，反问题是一个包含正问题的优化求解问题，求解算法流程如 图2 图2 所示。 常见的正问题求解模型有麦克斯韦方程、辐射传输方程和扩散方程等。扩散方程是光源各向同性情况下从辐射传输方程简化而来的，相对辐射传输方程来说具有计算难度低、计算时间快等优点。NIRFAST软件包 ［ 18 18 ］ 是一款被广泛使用的DOT算法开源软件包，所使用的正问题求解模型即为扩散方程。这里以NIRFAST软件包使用的正问题求解模型及反问题算法为例介绍 图2 图2 中的流程。首先，输入介质的尺寸参数和出射表面测得的能流密度 Ф M ，然后设定初始的吸收散射系数，代入软件包中的扩散方程进行正向问题求解，得到计算的能流密度分布 Ф C ，扩散方程一般为： <math display="block" id="M1"><mo>-</mo><mo>∇</mo><mo>⋅</mo><mi>k</mi><mo stretchy="false">(</mo><mi>r</mi><mo stretchy="false">)</mo><mo>∇</mo><mi>Φ</mi><mo stretchy="false">(</mo><mi>r</mi><mo>,</mo><mi>ω</mi><mo stretchy="false">)</mo><mo>+</mo><mfenced separators="|"><mrow><msub><mrow><mi>μ</mi></mrow><mrow><mi mathvariant="normal">a</mi></mrow></msub><mo stretchy="false">(</mo><mi>r</mi><mo stretchy="false">)</mo><mo>+</mo><mfrac><mrow><mi>i</mi><mi>ω</mi></mrow><mrow><msub><mrow><mi>c</mi></mrow><mrow><mi mathvariant="normal">m</mi></mrow></msub><mo stretchy="false">(</mo><mi>r</mi><mo stretchy="false">)</mo></mrow></mfrac></mrow></mfenced><mi>Φ</mi><mo stretchy="false">(</mo><mi>r</mi><mo>,</mo><mi>ω</mi><mo stretchy="false">)</mo><mo>=</mo><msub><mrow><mi>q</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo stretchy="false">(</mo><mi>r</mi><mo>,</mo><mi>ω</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24007255&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24007252&type= 92.03266144 10.49866676 ， 式中： k 是约化散射系数， Ф 是强度 I 对角度的积分，也就是能流密度， μ a 是吸收系数， ω 是光的调制频率，对于稳态DOT来说，调制频率为0， q 0 是源项。求解该方程一般使用数值方法，包括有限元、蒙特卡洛等。 如果计算所得和测量所得的能流密度分布的差值大于设定的阈值，则根据差值使用 式（2） 式（2） 所示的列文伯格-马夸特法确定迭代的参数更新值，进行反向问题求解。 <math display="block" id="M2"><mo stretchy="false">(</mo><msup><mrow><mi mathvariant="bold">J</mi></mrow><mrow><mi mathvariant="normal">T</mi></mrow></msup><msup><mrow><mi mathvariant="bold">J</mi><mo>-</mo><mover><mrow><mi>λ</mi></mrow><mo>¯</mo></mover><mi mathvariant="bold">I</mi><mo stretchy="false">)</mo></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msup><msup><mrow><mi mathvariant="bold">J</mi></mrow><mrow><mi mathvariant="normal">T</mi></mrow></msup><mi>δ</mi><mi>Φ</mi><mo>=</mo><mi>δ</mi><mi>μ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24007261&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24007258&type= 36.91466904 4.23333359 ， 式中： <math id="M3"><mover accent="true"><mi>λ</mi><mo>¯</mo></mover></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24007267&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24007264&type= 1.77800000 3.30200005 是正则系数， I 表示单位矩阵， δФ 代表能流密度测量值和计算值之间的差值， δμ 代表吸收系数及散射系数的更新值， J 是雅阁比矩阵，是正向计算值对于参数的导数矩阵 <math id="M4"><mrow><mrow><mo>∂</mo><msub><mrow><mi>Φ</mi></mrow><mrow><mi mathvariant="normal">C</mi></mrow></msub></mrow><mo>/</mo><mrow><mo>∂</mo><mi>μ</mi></mrow></mrow></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24007274&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24007270&type= 11.26066685 4.14866638 ，可以表示为： <math display="block" id="M5"><mi mathvariant="bold">J</mi><mo>=</mo><mfenced open="[" close="]" separators="|"><mrow><mtable columnalign="left"><mtr><mtd><mfrac><mrow><mo>∂</mo><mi mathvariant="normal">l</mi><mi mathvariant="normal">n</mi><mtext> </mtext><msub><mrow><mi>Φ</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></mrow><mrow><mo>∂</mo><msub><mrow><mi>μ</mi></mrow><mrow><mi mathvariant="normal">a</mi><mn mathvariant="normal">1</mn></mrow></msub></mrow></mfrac><mo>,</mo><mfrac><mrow><mo>∂</mo><mi mathvariant="normal">l</mi><mi mathvariant="normal">n</mi><mtext> </mtext><msub><mrow><mi>Φ</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></mrow><mrow><mo>∂</mo><msub><mrow><mi>μ</mi></mrow><mrow><mi mathvariant="normal">a</mi><mn mathvariant="normal">2</mn></mrow></msub></mrow></mfrac><mo>,</mo><mo>⋯</mo><mfrac><mrow><mo>∂</mo><mi mathvariant="normal">l</mi><mi mathvariant="normal">n</mi><mtext> </mtext><msub><mrow><mi>Φ</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></mrow><mrow><mo>∂</mo><msub><mrow><mi>μ</mi></mrow><mrow><mi mathvariant="normal">a</mi><mi mathvariant="normal">N</mi><mi mathvariant="normal">N</mi></mrow></msub></mrow></mfrac><mo>;</mo><mfrac><mrow><mo>∂</mo><mi mathvariant="normal">l</mi><mi mathvariant="normal">n</mi><mtext> </mtext><msub><mrow><mi>Φ</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></mrow><mrow><mo>∂</mo><msub><mrow><mi>k</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></mrow></mfrac><mo>,</mo><mfrac><mrow><mo>∂</mo><mi mathvariant="normal">l</mi><mi mathvariant="normal">n</mi><mtext> </mtext><msub><mrow><mi>Φ</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></mrow><mrow><mo>∂</mo><msub><mrow><mi>k</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></mrow></mfrac><mo>,</mo><mo>⋯</mo><mfrac><mrow><mo>∂</mo><mi mathvariant="normal">l</mi><mi mathvariant="normal">n</mi><mtext> </mtext><msub><mrow><mi>Φ</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></mrow><mrow><mo>∂</mo><msub><mrow><mi>k</mi></mrow><mrow><mi>N</mi></mrow></msub></mrow></mfrac></mtd></mtr><mtr><mtd><mfrac><mrow><mo>∂</mo><mi mathvariant="normal">l</mi><mi mathvariant="normal">n</mi><mtext> </mtext><msub><mrow><mi>Φ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></mrow><mrow><mo>∂</mo><msub><mrow><mi>μ</mi></mrow><mrow><mi mathvariant="normal">a</mi><mn mathvariant="normal">1</mn></mrow></msub></mrow></mfrac><mo>,</mo><mfrac><mrow><mo>∂</mo><mi mathvariant="normal">l</mi><mi mathvariant="normal">n</mi><mtext> </mtext><msub><mrow><mi>Φ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></mrow><mrow><mo>∂</mo><msub><mrow><mi>μ</mi></mrow><mrow><mi mathvariant="normal">a</mi><mn mathvariant="normal">2</mn></mrow></msub></mrow></mfrac><mo>,</mo><mo>⋯</mo><mfrac><mrow><mo>∂</mo><mi mathvariant="normal">l</mi><mi mathvariant="normal">n</mi><mtext> </mtext><msub><mrow><mi>Φ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></mrow><mrow><mo>∂</mo><msub><mrow><mi>μ</mi></mrow><mrow><mi mathvariant="normal">a</mi><mi mathvariant="normal">N</mi></mrow></msub></mrow></mfrac><mo>;</mo><mfrac><mrow><mo>∂</mo><mi mathvariant="normal">l</mi><mi mathvariant="normal">n</mi><mtext> </mtext><msub><mrow><mi>Φ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></mrow><mrow><mo>∂</mo><msub><mrow><mi>k</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></mrow></mfrac><mo>,</mo><mfrac><mrow><mo>∂</mo><mi mathvariant="normal">l</mi><mi mathvariant="normal">n</mi><mtext> </mtext><msub><mrow><mi>Φ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></mrow><mrow><mo>∂</mo><msub><mrow><mi>k</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></mrow></mfrac><mo>,</mo><mo>⋯</mo><mfrac><mrow><mo>∂</mo><mi mathvariant="normal">l</mi><mi mathvariant="normal">n</mi><mtext> </mtext><msub><mrow><mi>Φ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></mrow><mrow><mo>∂</mo><msub><mrow><mi>k</mi></mrow><mrow><mi>N</mi></mrow></msub></mrow></mfrac></mtd></mtr><mtr><mtd><mtable><mtr><mtd><mtext>    </mtext><mtable><mtr><mtd><mo>⋮</mo><mtext>          </mtext></mtd><mtd><mo>⋯</mo></mtd><mtd><mtext> </mtext><mo>⋮</mo></mtd><mtd><mtext>    </mtext><mo>⋯</mo></mtd></mtr></mtable><mtext>      </mtext></mtd><mtd><mo>⋯</mo></mtd><mtd><mtext>        </mtext><mo>⋮</mo></mtd><mtd><mtext>   </mtext><mo>⋯</mo></mtd></mtr></mtable><mtext>   </mtext><mo>⋮</mo></mtd></mtr><mtr><mtd><mfrac><mrow><mo>∂</mo><mi mathvariant="normal">l</mi><mi mathvariant="normal">n</mi><mtext> </mtext><msup><mrow><mi>Φ</mi></mrow><mrow><mi mathvariant="normal">d</mi><mi mathvariant="normal">s</mi></mrow></msup></mrow><mrow><mo>∂</mo><msub><mrow><mi>μ</mi></mrow><mrow><mi mathvariant="normal">a</mi><mn mathvariant="normal">1</mn></mrow></msub></mrow></mfrac><mo>,</mo><mfrac><mrow><mo>∂</mo><mi mathvariant="normal">l</mi><mi mathvariant="normal">n</mi><mtext> </mtext><msup><mrow><mi>Φ</mi></mrow><mrow><mi mathvariant="normal">d</mi><mi mathvariant="normal">s</mi></mrow></msup></mrow><mrow><mo>∂</mo><msub><mrow><mi>μ</mi></mrow><mrow><mi mathvariant="normal">a</mi><mn mathvariant="normal">2</mn></mrow></msub></mrow></mfrac><mo>,</mo><mo>⋯</mo><mfrac><mrow><mo>∂</mo><mi mathvariant="normal">l</mi><mi mathvariant="normal">n</mi><mtext> </mtext><msup><mrow><mi>Φ</mi></mrow><mrow><mi mathvariant="normal">d</mi><mi mathvariant="normal">s</mi></mrow></msup></mrow><mrow><mo>∂</mo><msub><mrow><mi>μ</mi></mrow><mrow><mi mathvariant="normal">a</mi><mi mathvariant="normal">N</mi></mrow></msub></mrow></mfrac><mo>;</mo><mfrac><mrow><mo>∂</mo><mi mathvariant="normal">l</mi><mi mathvariant="normal">n</mi><mtext> </mtext><msup><mrow><mi>Φ</mi></mrow><mrow><mi mathvariant="normal">d</mi><mi mathvariant="normal">s</mi></mrow></msup></mrow><mrow><mo>∂</mo><msub><mrow><mi>k</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></mrow></mfrac><mo>,</mo><mfrac><mrow><mo>∂</mo><mi mathvariant="normal">l</mi><mi mathvariant="normal">n</mi><mtext> </mtext><msup><mrow><mi>Φ</mi></mrow><mrow><mi mathvariant="normal">d</mi><mi mathvariant="normal">s</mi></mrow></msup></mrow><mrow><mo>∂</mo><msub><mrow><mi>k</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></mrow></mfrac><mo>,</mo><mo>⋯</mo><mfrac><mrow><mo>∂</mo><mi mathvariant="normal">l</mi><mi mathvariant="normal">n</mi><mtext> </mtext><msup><mrow><mi>Φ</mi></mrow><mrow><mi mathvariant="normal">d</mi><mi mathvariant="normal">s</mi></mrow></msup></mrow><mrow><mo>∂</mo><msub><mrow><mi>k</mi></mrow><mrow><mi>N</mi></mrow></msub></mrow></mfrac></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24007281&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24007276&type= 100.16066742 37.33800125 ， 式中： μ a 和 k 分别代表吸收系数和约化散射系数， ds 为光源位置数量与探测点位置数量的乘积， N 表示网格节点数。根据 式（2） 式（2） 所得的参数更新值重新进行正向求解，直到能流密度计算值和测量值之间的误差小于阈值。此时的吸收和散射系数分布即为计算重建结果。从 式（3） 式（3） 中可以看到，在实际求解过程中，采用了能流密度数值的对数，这要求DOT系统探测器必须具有较大的动态范围，才能利用软件算法实现高精度的计算重建。 2.2　深度学习网络模型及损失函数 2 深度学习是机器学习的一种，它使用大量样本进行模型训练，建立起输入样本集和目标样本集之间的模型。在本文中参与训练的数据是一维数据，且输入数据和目标数据长度相同。为了实现数据的训练，选择了全连接网络作为深度学习网络模型。理论上全连接网络对于复杂问题的建模能力略好于卷积神经网络 ［ 19 19 ］ ，且能够灵活地调整网络中每一层的大小，能够更容易地实现下采样和上采样。使用的网络模型示意图如 图3 图3 所示。 网络的输入为一个1×40的向量，经过3次下采样后变为1×10的向量，之后再经过3次上采样变回1×40的向量。训练的目标是修复由于探测器动态范围受限而测量误差较大的数据，所以在实际训练中，使用NIRFAST开源程序对不同二维吸收散射特性的介质进行透射能流密度的计算，通过数据截断方法模拟低动态范围探测器的信号获取过程，生成大量训练样本。考虑到这些数据动态范围大的特点，为了使所有数据在训练过程中具有相近的权重，同时为了使训练结果适用于DOT重建算法，采用 式（4） 式（4） 所示的对数处理方法，对由NIRFAST计算得到的透射能流密度和数据截断后的透射能流密度进行预处理，分别得到网络的目标数据和输入数据，即： <math display="block" id="M6"><mi>y</mi><mo>=</mo><mfenced open="|" close="|" separators="|"><mrow><mi mathvariant="normal">l</mi><mi mathvariant="normal">n</mi><mo stretchy="false">(</mo><msub><mrow><mi>Φ</mi></mrow><mrow><mi>i</mi></mrow></msub><mo stretchy="false">)</mo></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24007299&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24007296&type= 20.06599998 4.91066647 ， 式中： y 为网络目标数据和输入数据， Ф 为经过NIRFAST计算出的能流密度及其截断处理后的值。网络的损失函数被设置为有权重的均方差（Mean Square Error， MSE），即有： <math display="block" id="M7"><mi>L</mi><mo>=</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><msub><mrow><mi>y</mi></mrow><mrow><mi mathvariant="normal">t</mi><mi mathvariant="normal">u</mi><mi mathvariant="normal">r</mi><mi mathvariant="normal">e</mi></mrow></msub></mrow></mfrac><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mi>M</mi></mrow></mfrac><munder><mo movablelimits="false">∑</mo><mrow><mi>M</mi></mrow></munder><mrow><mo stretchy="false">(</mo><msub><mrow><mi>y</mi></mrow><mrow><mi mathvariant="normal">t</mi><mi mathvariant="normal">u</mi><mi mathvariant="normal">r</mi><mi mathvariant="normal">e</mi></mrow></msub><mo>-</mo><msub><mrow><mi>y</mi></mrow><mrow><mi mathvariant="normal">t</mi><mi mathvariant="normal">r</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi></mrow></msub><msup><mrow><mo stretchy="false">)</mo></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24007304&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=24007302&type= 43.18000031 8.80533314 ， 式中： L 即为损失值， y ture 即训练的目标值， M 为数据长度， y train 为经过网络模型输出的数据。该损失函数可以使训练过程中绝对值差异较大的数据具备一致的模型训练收敛性。 3 仿真及实验结果 1 3.1　深度学习网络模型训练 2 在模型训练中，吸收散射介质为65 mm×15 mm的矩形仿体，参考了实际应用中的DOT系统的光源密度和探测器密度 ［ 2 2 - 3 3 2-3 ， 6 6 ］ ，把光源位置数设置为4个，探测点数量设置为10个，分别等间距地分布于介质上下表面，总探测次数等于光源与探测器位置数量的乘积，即40次。训练样本中仿体背景的吸收系数分布于0.3~0.5 mm - 1 之间，异质体的吸收系数分布于0.4~1 mm - 1 之间。通过设置不同的背景吸收系数、异质体的大小、位置以及吸收系数等参数，计算出的能流密度分布在10 - 5 ~10 - 25 内，经过 式（4） 式（4） 预处理后，共产生了含有4 000个样本的训练集、100个样本的验证集以及7个样本的测试集。 为了模拟低动态范围探测器（动态范围为10 3 ）的信号获取过程，所有能流密度探测数据小于最大值千分之一的数据都会被人为截断，作为一维网络的输入值，利用上述网络模型和损失函数开展训练。在实际训练中，为了展示不同样本量的训练集对于网络预测能力的影响，分别使用了1 000个样本和4 000个样本构成训练集进行网络训练，并采用同样的验证集进行检测验证，训练过程中训练集和验证集的损失值下降如 图4 图4 所示。 从 图4 图4 中可以发现，网络训练过程中，不论是1 000样本量还是4 000样本量，训练集和验证集的误差均呈现下降趋势，这反映了模型训练的有效性。 3.2　仿真验证结果 2 为对模型性能进行验证和分析，在网络模型训练完成后，把没有参与训练的7个测试样本输入到网络中，可获得如 图5 图5 所示的其中3个样本的测试结果。 图5 图5 展示了按照 式（4） 式（4） 预处理前样本能流密度的完整数据、截断数据和修复数据。从图中可以看到，网络模型修复后的数据与截断前的完整数据非常接近，这说明训练好的网络修复了被截断的数据。为了进一步验证网络模型对低动态范围探测器测量数据的修复效果，使用 图2 图2 所示流程进行了DOT重建，重建效果如 图6 图6 和 图7 图7 所示。 图6 图6 为利用 图5 图5 中所示的被截断处理的数据重建得到的吸收系数分布结果。可以发现，不论仿体和异质体的参数如何设置，模拟低动态范围探测器采集的数据重建结果具有相同的特点，即吸收系数分布中会出现一些高吸收系数（10 mm -1 左右）的异常区域，这些高吸收区域的分布呈现空间对称性。该结果说明使用低动态范围探测器采集吸收散射光信号所产生的数据误差严重影响了DOT的重建效果，对DOT计算重建来说是不容忽视的。 图7 图7 为分别利用经过 图5 图5 中所示完整的数据和网络模型修复后的数据重建得到的吸收系数分布结果。虽然 图7 图7 （b）和 图7 图7 （c）中吸收系数重建结果在对比度和空间分布上和 图7 图7 （a）中的目标值重建结果有偏差，但与 图6 图6 中使用截断数据重建得到的吸收系数分布结果相比，使用经过网络模型修复后的预测数据能够成功进行重建。为了进一步分析重建结果，使用均方根误差（Root Mean Square Error， RMSE）和皮尔森相关系数（Pierson Correlation Coefficient， PCC）对 图7 图7 中的吸收系数分布进行量化比较，均方根误差越小，皮尔森相关系数越接近于1，这表示由网络模型修复数据重建的吸收系数分布越接近于截断前完整数据的重建结果。数据如 表1 表1 所示。 2 测试样本号 2 训练集样本量1 000的网络 2 训练集样本量4 000的网络 RMSE PCC RMSE PCC 1 0.018 3 0.628 6 0.017 2 0.701 7 2 0.047 0 0.616 4 0.042 7 0.695 1 3 0.054 8 0.654 2 0.051 0 0.714 8 从 表1 表1 中可以发现，使用更多训练样本训练出的模型具有更好的数据修复性能，而且修复的数据重建出的吸收系数和截断前完整数据重建出的吸收系数是强相关的。 由此可见，使用深度学习方法对模拟动态范围为10 3 的探测器采集的有误差的DOT能流密度数据进行修复，能够将探测器的动态范围扩展至10 20 ，经由DOT重建算法对数据计算后，可以获得较好的吸收系数分布重建结果。模拟结果说明，该方法可以有效扩展探测器的动态范围，减小低动态范围探测器采集信号对DOT重建结果精度的影响。 3.3　实验结果 2 在仿真分析的基础上，进一步对利用深度学习方法扩展非接触DOT系统中探测器的动态范围进行实验验证。首先，参考文献报道 ［ 10 10 ］ 搭建了非接触式稳态DOT系统，系统示意图和实物图分别如 图8 图8 （a）和 图8 图8 （b）所示。整个实验系统由三部分组成，分别是在一维方向移动的光源，具有一定吸收散射系数的仿体，以及对仿体后表面进行成像的探测器。 光源部分由光纤输出的激光器（B & W TEK INC. BRM-785-0.55-100-0.22-SMA）、准直透镜（Ocean Insight 74-VIS）、夹持器和电动平移台（Thorlabs LTS300/M）组成。探测器部分由镜头和CCD（DAHENG MV-EM510M/C）组成，其中CCD的灰阶为8位，动态范围为256。实验中仿体是采用环氧树脂作为基质，掺杂印度墨水和氧化锌粉末，并加入固化剂进行固化得到的长方体固体，仿体上钻孔填充一定浓度的印度墨水溶液作为异质体，其横截面结构尺寸及吸收系数分布如 图8 图8 （c）所示。仿体溶液配置过程中，使用吸收光谱法可以确定掺杂印度墨水后的环氧树脂溶液吸收系数，通过调整印度墨水的掺杂浓度，确定溶液吸收系数为0.1 mm - 1 。根据氧化锌粉末的参数以及Mie散射理论 ［ 20 20 ］ 可计算获得溶液散射系数的理论值，通过调整氧化锌粉末的掺杂浓度，近似确定溶液散射系数为1 mm - 1 。异质体溶液配置过程中，采用与仿体相同的方法确定其吸收系数为0.3 mm - 1 。 在实验测量过程中，光源及探测器的数量和分布情况同模拟仿真中一致。激光束依次从仿体前表面预设的位置入射，经仿体吸收散射后从后表面出射并被CCD相机探测。为获得仿体后表面出射光的能流密度分布，通过实验标定法建立了仿体后表面能流密度和CCD像素灰度值的关系 ［ 21 21 ］ 。按照实验中仿体设定的吸收系数及散射系数进行计算，可知能流密度的理论动态范围大约为10 10 ，这远大于实验中CCD的动态范围。根据3.2节中仿真模拟得到的结论，在这种情况下想要获得无重建误差的结果，需要对采集到的数据进行修复。 利用3.1中训练得到的网络模型，对实验采集到的能流密度数据进行修复，结果如 图9 图9 所示。从 图9 图9 可以发现，实际测量值的动态范围大约为200，修复后数据的动态范围被扩展到10 9 左右。这说明深度学习方法修复了低动态范围探测器采集的能流密度数据，并扩展了探测器的动态范围。使用 图9 图9 中的数据进行了吸收系数重建，得到的吸收系数分布如 图10 图10 所示。 从 图10 图10 （c）可以发现，使用实验中探测得到的能流密度进行DOT重建得到的吸收系数分布结果与 图8 图8 （c）中的预设参数相比有较大的误差，而且吸收系数的空间分布与 图6 图6 中仿真截断后数据的重建结果相似。这从实验上验证了探测器动态范围过小会导致探测到的能流密度的误差较大，从而严重影响DOT的重建精度。 图10 图10 （a）和 图10 图10 （b）中的结果则说明了尽管使用不同样本量的训练集训练的网络模型对于数据的修复能力有差异，但是都可以显著提升DOT重建吸收系数分布结果的精度。在此基础上计算了 图10 图10 （a）和 图10 图10 （b）中的结果与 图8 图8 （c）中设定值的RMSE，分别为0.545和0.542。虽然重建的吸收系数分布和 图8 图8 （c）的理论预设值之间存在一定的误差，但考虑到实验中运动精度、器件噪声、仿体固化等不确定性因素的影响，该误差是在可接受范围内的。 实验结果表明，深度学习方法对实验中动态范围为256的探测器采集的有误差的DOT能流密度数据进行修复后，能够将探测器的动态范围扩展至10 9 ，经由DOT重建算法对数据计算后，吸收系数分布的重建结果准确。这从实验上证明深度学习方法可以有效扩展探测器的动态范围，减小低动态范围探测器采集信号对DOT重建结果精度的影响。 4 结　论 1 本文提出了一种利用深度学习扩展非接触扩散光学层析成像系统中探测器动态范围的方法。首先，使用NIRFAST软件包生成了一系列模拟低动态范围采集的DOT数据。其次，使用这些数据完成了一维全连接网络模型的训练。之后分别通过模拟仿真和实验对该模型扩展非接触DOT系统中探测器动态范围的能力进行了验证。在模拟中，该方法修复了模拟动态范围为10 3 的探测器采集到的含有误差的测试集数据，显著减小了数据误差对DOT重建结果精度的影响，理论上将探测器的动态范围扩展到10 20 ；在实验中，制备了实验仿体并搭建了一套非接触DOT系统，使用动态范围为256的CCD相机完成了DOT系统的实验数据采集。利用本文提出的方法对采集到的数据修复后，可重建得到与预设仿体吸收系数分布基本一致的重建结果，将DOT系统中探测器的动态范围从256扩展到10 9 。虽然模拟和实验中通过网络模型修复的数据重建出的吸收系数结果的精度还有待提高，但该方法能够有效修复低动态范围探测器所采集的有误差的DOT数据，减小它对DOT重建结果精度的影响，为利用非接触DOT系统探测高吸收散射介质或大面积探测提供了一个可行的解决方案，使得非接触式DOT系统有了更广阔的应用空间。 目前，该方法只对单层的仿体进行了验证。如果进一步修改DOT求解模型，使之能适用于多层不同介质组成的仿体，该方法还可以应用到更复杂的多层仿体或生物组织上。 参考文献 ARRIDGE S R ， HEBDEN J C . Optical imaging in medicine： II. Modelling and reconstruction ［J］. Physics in Medicine and Biology ， 1997 ， 42 （ 5 ）： 841 - 853 . doi: 10.1088/0031-9155/42/5/008 http://dx.doi.org/10.1088/0031-9155/42/5/008 Optical imaging in medicine： II. Modelling and reconstruction EGGEBRECHT A T ， FERRADAL S L ， ROBICHAUX-VIEHOEVER A ， et al . Mapping distributed brain function and networks with diffuse optical tomography ［J］. Nature Photonics ， 2014 ， 8 （ 6 ）： 448 - 454 . doi: 10.1038/nphoton.2014.107 http://dx.doi.org/10.1038/nphoton.2014.107 Mapping distributed brain function and networks with diffuse optical tomography ALFONSO ， GALDERISI ， SABRINA ， et al . Long-term continuous monitoring of the preterm brain with diffuse optical tomography and electroencephalography： a technical note on cap manufacturing ［J］. Neurophotonics ， 2016 ， 3 （ 4 ）： 045009 . doi: 10.1117/1.nph.3.4.045009 http://dx.doi.org/10.1117/1.nph.3.4.045009 Long-term continuous monitoring of the preterm brain with diffuse optical tomography and electroencephalography： a technical note on cap manufacturing NTZIACHRISTOS V . Concurrent Diffuse Optical Tomography ， Spectroscopy and Magnetic Resonance Imaging of Breast Cancer ［D］. Philadelphia ： University of Pennsylvania . 2000 . doi: 10.1364/bosd.2000.wb4 http://dx.doi.org/10.1364/bosd.2000.wb4 Concurrent Diffuse Optical Tomography， Spectroscopy and Magnetic Resonance Imaging of Breast Cancer JIANG H B ， IFTIMIA N V ， XU Y ， et al . Near-infrared optical imaging of the breast with model-based reconstruction ［J］. Academic Radiology ， 2002 ， 9 （ 2 ）： 186 - 194 . doi: 10.1016/s1076-6332(03)80169-1 http://dx.doi.org/10.1016/s1076-6332(03)80169-1 Near-infrared optical imaging of the breast with model-based reconstruction YUAN Z ， ZHANG Q Z ， SOBEL E S ， et al . Tomographic x-ray-guided three-dimensional diffuse optical tomography of osteoarthritis in the finger joints ［J］. Journal of Biomedical Optics ， 2008 ， 13 （ 4 ）： 044006 . doi: 10.1117/1.2965547 http://dx.doi.org/10.1117/1.2965547 Tomographic x-ray-guided three-dimensional diffuse optical tomography of osteoarthritis in the finger joints XU Y ， IFTIMIA N ， JIANG H B ， et al . Imaging of in vitro and in vivo bones and joints with continuous-wave diffuse optical tomography ［J］. Optics Express ， 2001 ， 8 （ 7 ）： 447 - 451 . doi: 10.1364/oe.8.000447 http://dx.doi.org/10.1364/oe.8.000447 Imaging of in vitro and in vivo bones and joints with continuous-wave diffuse optical tomography 李江涛 ， 胡文雁 ， 赵龙莲 ， 等 . 频域近红外光学成像法的苹果内部病变检测精度 ［J］. 光谱学与光谱分析 ， 2019 ， 39 （ 9 ）： 2836 - 2841 . doi: 10.1016/j.postharvbio.2019.03.014 http://dx.doi.org/10.1016/j.postharvbio.2019.03.014 频域近红外光学成像法的苹果内部病变检测精度 LI J T ， HU W Y ， ZHAO L L ， et al . Study of the accuracy of apple internal lesion detection based on frequency domain diffuse optical tomography ［J］. Spectroscopy and Spectral Analysis ， 2019 ， 39 （ 9 ）： 2836 - 2841 . （in Chinese） . doi: 10.1016/j.postharvbio.2019.03.014 http://dx.doi.org/10.1016/j.postharvbio.2019.03.014 Study of the accuracy of apple internal lesion detection based on frequency domain diffuse optical tomography YOKO HOSHI M D ， YAMADA Y . Overview of diffuse optical tomography and its clinical applications ［J］. Journal of Biomedical Optics ， 2016 ， 21 （ 9 ）： 091312 . doi: 10.1117/1.jbo.21.9.091312 http://dx.doi.org/10.1117/1.jbo.21.9.091312 Overview of diffuse optical tomography and its clinical applications ZHOU X Q ， FAN Y ， HOU Q ， et al . Spatial-frequency-compression scheme for diffuse optical tomography with dense sampling dataset ［J］. Applied Optics ， 2013 ， 52 （ 9 ）： 1779 - 1792 . doi: 10.1364/ao.52.001779 http://dx.doi.org/10.1364/ao.52.001779 Spatial-frequency-compression scheme for diffuse optical tomography with dense sampling dataset 范颖 . 基于空间光测量的光学层析成像系统研究 ［D］. 天津 ： 天津大学 ， 2012 . 基于空间光测量的光学层析成像系统研究 FAN Y . Research on Continues Wave Non-Contact Optical Tomography System ［D］. Tianjin ： Tianjin University ， 2012 . （in Chinese） Research on Continues Wave Non-Contact Optical Tomography System ANDO T ， NAKAMURA T ， FUJII T ， et al . Non-contact acquisition of brain function using a time-extracted compact camera ［J］. Scientific Reports ， 2019 ， 9 ： 17854 . doi: 10.1038/s41598-019-54458-7 http://dx.doi.org/10.1038/s41598-019-54458-7 Non-contact acquisition of brain function using a time-extracted compact camera RIVENSON Y ， GÖRÖCS Z ， GÜNAYDIN H ， et al . Deep learning microscopy ［J］. Optica ， 2017 ， 4 （ 11 ）： 1437 - 1443 . doi: 10.1364/optica.4.001437 http://dx.doi.org/10.1364/optica.4.001437 Deep learning microscopy LI S ， DENG M ， LEE J ， et al . Imaging through glass diffusers using densely connected convolutional networks ［J］. Optica ， 2018 ， 5 （ 7 ）： 803 - 813 . doi: 10.1364/optica.5.000803 http://dx.doi.org/10.1364/optica.5.000803 Imaging through glass diffusers using densely connected convolutional networks BORHANI N ， KAKKAVA E ， MOSER C ， et al . Learning to see through multimode fibers ［J］. Optica ， 2018 ， 5 （ 8 ）： 960 - 966 . doi: 10.1364/optica.5.000960 http://dx.doi.org/10.1364/optica.5.000960 Learning to see through multimode fibers 闫敬文 ， 陈宏达 ， 刘蕾 . 高光谱图像分类的研究进展 ［J］. 光学 精密工程 ， 2019 ， 27 （ 3 ）： 680 - 693 . doi: 10.3788/ope.20192703.0680 http://dx.doi.org/10.3788/ope.20192703.0680 高光谱图像分类的研究进展 YAN J W ， CHEN H D ， LIU L . Overview of hyperspectral image classification ［J］. Opt. Precision Eng. ， 2019 ， 27 （ 3 ）： 680 - 693 . （in Chinese） . doi: 10.3788/ope.20192703.0680 http://dx.doi.org/10.3788/ope.20192703.0680 Overview of hyperspectral image classification 李宇 ， 刘雪莹 ， 张洪群 ， 等 . 基于卷积神经网络的光学遥感图像检索 ［J］. 光学 精密工程 ， 2018 ， 26 （ 1 ）： 200 - 207 . doi: 10.3788/ope.20182601.0200 http://dx.doi.org/10.3788/ope.20182601.0200 基于卷积神经网络的光学遥感图像检索 LI Y ， LIU X Y ， ZHANG H Q ， et al . Optical remote sensing image retrieval based on convolutional neural networks ［J］. Opt. Precision Eng. ， 2018 ， 26 （ 1 ）： 200 - 207 . （in Chinese） . doi: 10.3788/ope.20182601.0200 http://dx.doi.org/10.3788/ope.20182601.0200 Optical remote sensing image retrieval based on convolutional neural networks DEHGHANI H ， EAMES M E ， YALAVARTHY P K ， et al . Near infrared optical tomography using NIRFAST： Algorithm for numerical model and image reconstruction ［J］. Communications in Numerical Methods in Engineering ， 2009 ， 25 （ 6 ）： 711 - 732 . doi: 10.1002/cnm.1162 http://dx.doi.org/10.1002/cnm.1162 Near infrared optical tomography using NIRFAST： Algorithm for numerical model and image reconstruction KRIZHEVSKY A ， SUTSKEVER I ， HINTON G E . ImageNet classification with deep convolutional neural networks ［J］. Communications of the ACM ， 2017 ， 60 （ 6 ）： 84 - 90 . doi: 10.1145/3065386 http://dx.doi.org/10.1145/3065386 ImageNet classification with deep convolutional neural networks WANG L V ， WU H I . Biomedical Optics ［M］. Hoboken ： John Wiley & Sons， Inc. ， 2009 . doi: 10.1002/tdm_license_1.1 http://dx.doi.org/10.1002/tdm_license_1.1 Biomedical Optics 秦莉 ， 董丽丽 ， 许文海 ， 等 . CCD图像灰度与照度的转换标定方法 ［J］. 仪器仪表学报 ， 2015 ， 36 （ 3 ）： 639 - 644 . CCD图像灰度与照度的转换标定方法 QIN L ， DONG L L ， XU W H ， et al . Method for conversion calibration between CCD image gray value and illumination ［J］. Chinese Journal of Scientific Instrument ， 2015 ， 36 （ 3 ）： 639 - 644 . （in Chinese） Method for conversion calibration between CCD image gray value and illumination 10.37188/OPE.2021.0219 O439 TP391.4 A 1004-924X（2021）11-2529-10 赵星,刘桐君,陈平等.利用深度学习扩展非接触扩散光学层析成像系统中探测器的动态范围[J].光学精密工程,2021,29(11):2529-2538. 赵星,刘桐君,陈平等.利用深度学习扩展非接触扩散光学层析成像系统中探测器的动态范围[J].光学精密工程,2021,29(11):2529-2538. ZHAO Xing,LIU Tong-jun,CHEN Ping,et al.Extending dynamic range of detector in non-contact diffuse optical tomography system using deep learning[J].Optics and Precision Engineering,2021,29(11):2529-2538. ZHAO Xing,LIU Tong-jun,CHEN Ping,et al.Extending dynamic range of detector in non-contact diffuse optical tomography system using deep learning[J].Optics and Precision Engineering,2021,29(11):2529-2538. 国家重点研发计划资助项目（No.2018YFB0504400）；国家自然科学基金资助项目（No.62075106）；天津市科技支撑计划资助项目（No.19YFZCSY00250）；天津市自然科学基金资助项目（No.19JCZDJC36600） 2021-04-23 2021-05-13 2021-11-15 2022-01-25 光学精密工程 11 29