CN107252305B - 基于全并行锁相光子计数检测模式的nirs脑功能成像系统 - Google Patents

基于全并行锁相光子计数检测模式的nirs脑功能成像系统 Download PDF

Info

Publication number
CN107252305B
CN107252305B CN201710466642.4A CN201710466642A CN107252305B CN 107252305 B CN107252305 B CN 107252305B CN 201710466642 A CN201710466642 A CN 201710466642A CN 107252305 B CN107252305 B CN 107252305B
Authority
CN
China
Prior art keywords
detection
source
channel
phase
light source
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
Application number
CN201710466642.4A
Other languages
English (en)
Other versions
CN107252305A (zh
Inventor
高峰
贺捷
丁雪梅
丁号
赵会娟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nielsen Human Factors (Xiamen) Technology Co.,Ltd.
Original Assignee
Tianjin University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Tianjin University filed Critical Tianjin University
Priority to CN201710466642.4A priority Critical patent/CN107252305B/zh
Publication of CN107252305A publication Critical patent/CN107252305A/zh
Application granted granted Critical
Publication of CN107252305B publication Critical patent/CN107252305B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/0059—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
    • A61B5/0075—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence by spectroscopy, i.e. measuring spectra, e.g. Raman spectroscopy, infrared absorption spectroscopy
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/40—Detecting, measuring or recording for evaluating the nervous system
    • A61B5/4058—Detecting, measuring or recording for evaluating the nervous system for evaluating the central nervous system
    • A61B5/4064—Evaluating the brain
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/40—Detecting, measuring or recording for evaluating the nervous system
    • A61B5/4076—Diagnosing or monitoring particular conditions of the nervous system
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/40—Detecting, measuring or recording for evaluating the nervous system
    • A61B5/4076—Diagnosing or monitoring particular conditions of the nervous system
    • A61B5/4088—Diagnosing of monitoring cognitive diseases, e.g. Alzheimer, prion diseases or dementia
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/40—Detecting, measuring or recording for evaluating the nervous system
    • A61B5/4076—Diagnosing or monitoring particular conditions of the nervous system
    • A61B5/4094—Diagnosing or monitoring seizure diseases, e.g. epilepsy
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Neurology (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Public Health (AREA)
  • Animal Behavior & Ethology (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Pathology (AREA)
  • General Health & Medical Sciences (AREA)
  • Biophysics (AREA)
  • Veterinary Medicine (AREA)
  • Neurosurgery (AREA)
  • Physiology (AREA)
  • Psychiatry (AREA)
  • Signal Processing (AREA)
  • Psychology (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Hospice & Palliative Care (AREA)
  • Developmental Disabilities (AREA)
  • Child & Adolescent Psychology (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)
  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)

Abstract

本发明公开了一种基于全并行锁相光子计数检测模式的NIRS脑功能成像系统，包括可调制LED光源单元、源‑探光纤布配阵列单元，探测单元、FPGA以及计算机控制和数据处理单元；可调制LED光源单元包括LED光源模块、多通道方波信号发生器模块和20束源光纤，所述LED光源模块包括60个三种波长的LED和LED光源驱动电路；多通道方波信号发生器模块输出60路不同频率的方波信号；源‑探光纤布配阵列单元包括用于配置20个光源点和12个探测点的源‑探光纤布配贴片、多根源光纤和探测光纤，探测单元包括12根探测光纤、12个PMT光子计数器和12通道×60路可变选通锁相光子计数检测模块。本发明所提出的成像系统集高灵敏度、大动态范围与高时间分辨率等诸多优势于一身。

Description

基于全并行锁相光子计数检测模式的NIRS脑功能成像系统
技术领域
本发明属于近红外脑功能成像技术领域，具体涉及基于不同频率方波调制LED编码激励的多通道、多波长锁相光子计数全并行检测模式以及基于此模式建立的NIRS脑功能成像系统，同时提出了面向脑枕部视觉功能区探测的源-探布配方案。
背景技术
后脑枕叶部(简称枕部)作为主要的视觉区是视觉刺激脑功能(简称视觉脑功能)研究的重点区域之一[1]。对脑枕部视觉功能区的在体成像研究具有非常重要的意义：在科学研究领域，视觉脑功能成像是研究大脑认知机制的重要一环，有助于探索视觉注意和信息表达等众多未知神经机制的奥秘[1]；临床医学领域，视觉脑功能成像可更有效地用于对诸如多动症、自闭症、阿尔兹海默症、癫痫等神经系统的疾病的病因诊断和疗效分析等[2,3]；应用技术领域，基于视觉脑功能探测的脑机接口技术具有较其他激励方式更高的可靠性(通信精度和分辨精度)[4,5]。
近红外光谱(NIRS)成像方法使用两个以上波长的近红外光(650-900nm)，将源点和探测点在预定脑功能区内成阵列布配，由于生物组织在近红外波段的弱吸收、强散射特征，入射光可穿透头皮、头骨而达到脑皮层，并获得经颅吸收和散射共同作用后的扩散光(或漫反射光)。由于大脑的神经活动会导致周围区域氧合血红蛋白(HbO)和脱氧血红蛋白(HbR)浓度发生变化，而在近红外光测量和治疗的光学窗口，由于反映脑组织代谢和血流动力学特征的HbO和HbR正是近红外光波段内的主要吸收体，因此由应激前后漫反射光强度变化的空间分辨测量可还原出脑皮层HbO和HbR浓度变化(△CHbO和△CHbR)的空间分布，从而为脑科学基础研究和脑疾病诊断提供大脑皮层(负责高级的思维过程)应激响应的血氧代谢功能信息[6-10]。NIRS成像方法因其无创、深度探测、高动态性等诸多优点，是目前可用于脑功能在体测量的首选光学方法，同时它也是目前唯一可用于床边婴儿脑功能成像的方法[6]。
以NIRS为代表的脑功能测量技术将逐步发展成熟并进入临床测试应用阶段，成为现有脑成像模态(fMRI、EEG、MEG等)[11]的一个极为有益的补充。现有的NIRS系统有时域(Time-Domain，TD)[12]、频域(Frequency-Domain，FD)[13]和连续光(Continuous-Wave，CW)[14]三种主要的测量模式。其中TD-NIRS测量系统虽然具有测量信息相对完整的优势，但是其系统价格昂贵、测量时间长，难以在小型实验室中进行普及和测量神经快变信号。FD-NIRS系统一般需要200MHz以上的高频调制以实现相移测量所需的合理信噪比，实现难度较大，单频测量提供的信息有限,而多频测量系统性价比与时域测量相比不占优势。因此，CW系统是脑功能光学成像的主流技术。目前，国内外的科学工作者在CW-DOT技术系统和基础研究领域做了大量的工作，也取得了诸多成果。NIRS方法当前普遍采用的模拟测量方式，更适合于头骨较薄的前额和头侧部的检测。头骨较厚的成人脑枕部检测对系统的灵敏度提出了更高的要求。即便采用高灵敏度的光电倍增管(PMT)和雪崩光电二极管(APD)仍无法可靠地实现成人脑枕部测量所要求的高灵敏度和大动态范围，也无法获取具有较高信噪比的测量数据。
[参考文献]
[1]唐孝威等，《脑科学导论》，浙江大学出版社，2006.
[2]Ann-Christine Ehlis et al,“Application of functional near-infraredspectroscopy in psychiatry,”NeuroImage 85,478-88,2014.
[3]G.Stothart et al,“Early Visual Evoked Potentials and MismatchNegativity in Alzheimer's Disease and Mild Cognitive Impairment,”J.Alzhem.Disease 44,397,2015.
[4]S.-K.Gao et al,“Visual and auditory brain-computer interfaces”,IEEE Trans.Biomed.Eng.61,1436,2014.
[5]Y.Tomita et al,“Biomodal BCI using simultaneously NIRS and EEG”IEEE Trans.Biomed.Eng.61,1274,2014.
[6]S.Lloyd-Fox et al,"Illuminating the developing brain:The past,present and future of functional near infrared spectroscopy,"Neurosci.Biobehav.Rev.34,269-284,2010.
[7]F.Scholkmann et al,"A review on continuous wave functional near-infrared spectroscopy and imaging instrumentation and methodology,"NeuroImage85,6-27,2014.
[8]M.Ferrari and V.Quaresima,"A brief review on the history of humanfunctional near-infrared spctroscopy(fNIRS)development and fields ofapplication,"NeuroImage 63,921-35,2012.
[9]Y.Hoshi,"Towards the next generation of near-infraredspectroscopy,"Phil.Trans.R.Soc.A 369,4425-39,2011.
[10]D.A.Boas et al,"Twenty years of functional near-infraredspectroscopy:introduction for the special issue,"NeuroImage 85,1-5,2014.
[11]Hiroshi Kawaguchi,Tatsuya Koyama,Eiji Okada,Virtual Head Phantomfor Evaluation of Near Infrared Topography,Lasers and Electro-Optics Society,2006；
[12]Hideaki Koizumi,Atsushi Maki,Tsuyoshi Yamamoto,Hiroki Sato,YukariYamamoto,Hideo Kawaguchi,Non-invasive brain-function imaging by opticaltopography,Trends in Analytical Chemistry,Vol.24,No.2,2005；
[13]M.S.Patterson,B.Chance and B.C.Wilson,Time resolved reflectanceand transmittance for the noninvasive measurement of tissue opticalproperties,Appl.Opt.1989,28:2331–2336；
[14]J.S.Dam,C.B.Pedersen,T.Dalgaard et.al,Fiber optic probe fornoninvasive real-time determination of tissue optical properties at multiplewavelengths,Appl.Opt.40,2000:1155-1164；
发明内容
针对现有NIRS技术的不足，本发明旨在整合光子计数的超高灵敏度和锁相检测技术的多通道并行测量优势，发展一种基于不同频率方波调制LED编码激励的多通道、多波长锁相光子计数全并行检测模式并建立基于此模式的NIRS脑功能成像系统，同时提出了面向脑枕部视觉功能区的源-探布配方案，有效实现枕部视觉脑功能数据获取所要求的高灵敏度、大动态范围与高时间分辨率的统一。
为了解决上述技术问题，本发明提出一种基于全并行锁相光子计数检测模式的NIRS脑功能成像系统，包括源-探光纤布配阵列单元、可调制LED光源单元、探测单元、现场可编程门阵列(FPGA)和计算机控制及数据处理单元；所述源-探光纤布配阵列单元包括源-探光纤布配贴片、源光纤和探测光纤，源-探光纤布配贴片采用源-探点空间交叉布置的方式，源-探光纤布配贴片与20个光源点和12个探测点的位置对应，用于布配源光纤和探测光纤的位置，以实现光源从不同光源点入射和探测器从不同的探测点接收出射光；所述可调制LED光源单元包括LED光源模块、多通道方波信号发生器模块和20束源光纤，所述LED光源模块包括60个LED和LED光源驱动电路，所述60个LED包括20个波长为660nm的LED、20个波长为780nm的LED和20个波长为830nm的LED；所述多通道方波信号发生器模块由FPGA实现，用以输出频率为从6.2kHz到18kHz、且以0.2kHz等间隔分布的60路方波调制信号对60个LED分别进行调制；每一束源光纤由分别接通到660nm,780nm和830nm波长LED的三根单模光纤通过共用一个塑料护套耦合而成；共有20束这样的源光纤将60路光源分成20束传导到20个光源点；所述探测单元包括12根探测光纤，12个PMT光子计数器和12通道×60路可变选通锁相光子计数检测模块，其中，12根探测光纤用于将所述12个探测点位置处的扩散光分别传导到所述12个PMT光子计数器；所述12个PMT光子计数器用于将接收到的光信号转化为电脉冲信号；所述12通道×60路可变选通锁相光子计数检测模块由FPGA实现，其中，每通道对应所述12个探测点中的一个，每路对应所述60个LED的调制光源中的一个，从而使所述FPGA同时获取12个探测点处分别对应60路光源的出射光强信息，实行全并行快速测量；所述计算机控制及数据处理单元向多通道方波信号发生器模块发送频率控制字控制其产生选定频率的方波调制信号；同时，向12通道×60路可变选通锁相光子计数检测模块发送门宽控制字控制单次测量时间(累加时间)；等到测量过程结束后通过串口读取12个探测点处分别对应60路光源的出射光强信息；最后，利用扩散光学层析重建算法构建反映脑皮层光学参数变化量的空间分布图。
进一步讲，所述源-探光纤布配阵列单元中采用源-探点空间交叉布置方式布配贴片，用于确定20个光源点和12个探测点的位置；20个源点的布置方案是：按照4行5列以19mm距离等间距呈网格状布置；12个探测点的布置方案是：在每4个相邻源点的几何中心布置一个探测点；针对每个源点，将与该源点距离分别为L1＝13.4mm、L2＝30mm、L3＝40mm和L4＝46.5mm的所有探测点与该源点之间形成的通道依次定义为第一近邻通道、第二近邻通道、第三近邻通道和第四近邻通道；其中包括有44个第二近邻通道，在所有第二近邻通道的中点处均设有一采样点。
与现有技术相比，本发明的有益效果是：
本发明将锁相检测技术的背景噪声抑制和多波长-多通道并行测量能力与光子计数检测的高灵敏度和大动态范围性能优势相结合，发展了一种基于不同频率方波调制LED编码激励的多通道、多波长锁相光子计数全并行检测模式CW-NIRS成像系统,其满足脑枕部测量所要求的时间分辨率、动态范围和灵敏度要求，可以实现对深埋在背景和系统噪声中的枕部视觉功能相关NIRS信号的有效采集。假设所获得信号的信噪比需保持在大于40dB，以800nm波长光为例，计算可得光子计数测量能达到的最大动态范围为92-112dB(假设探测量子效率为1.0)，即可保证脑枕部NIRS系统较高的测量灵敏度和大的动态范围。该系统将多通道、多波长数字锁相光子计数检测技术和优化的源-探布配模式相结合，有效地克服了一般模拟光电检测技术在源-探测器距离较大时测量信噪比不高的局限性，增大了灵敏度、测量动态范围和时间分辨率，同时系统中所用到的可调制LED光源的驱动电路为自己设计制作，提高了系统的性价比，该系统在脑认知功能研究以及多动症、自闭症、阿尔兹海默症、癫痫等神经系统疾病的诊断和疗效分析等领域有着广阔的应用前景。
本发明采用多波长多通道并行激励-测量模式，测量时间分辨率只取决于光子计数的单次累加时间，便于提高系统的时间分辨率，实现多波长出射光信号的快速并行测量。本发明采用适合脑枕部视觉功能区探测的源-探布配方案以及满足测量动态范围的交叠采样NIRS模式，大大提高了图像重建结果的空间分辨率。本发明采用FPGA产生方波信号直接调制LED光源，避免了D/A转换导致的采样信号失真，而且减小设计的难度，降低了设计的成本。
附图说明
图1.本发明中NIRS脑功能成像系统的原理框图；
图2.本发明中LED光源驱动电路图；
图3.本发明中源-探测光纤布配示意图；
图4.本发明中单通道-60路可变选通锁相光子计数工作原理图。
图5.本发明中成像系统工作流程图。
具体实施方式
下面结合附图和具体实施例对本发明技术方案作进一步详细描述，所描述的具体实施例仅对本发明进行解释说明，并不用以限制本发明。
本发明提出的一种基于全并行锁相光子计数检测模式的NIRS脑功能成像系统，包括源-探光纤布配阵列单元、可调制LED光源单元、探测单元、FPGA和计算机控制及数据处理单元。其原理框图如图1所示，其中：
(一)可调制LED光源部分：包括可调制LED光源模块、多通道方波信号调制模块、20束源光纤及电源模块。
(1)可调制LED光源模块：此部分包括三种波长的半导体发光二极管LED(波长分别为660nm、780nm、830nm，每种各20个，共60个)及其对应的LED光源驱动电路。常用的LED驱动方式有恒压驱动和恒流驱动。为了保证系统的稳定性和测量结果的可重复性，要求LED的发光功率稳定。而LED是温度敏感元件，伏安特性存在负温度系数，而发光功率与正向电流成正比，恒压驱动方式不能保证其光功率的稳定性。本发明中，建议采用恒流的驱动方式，为了实现对3种不同波长LED的分别控制，采用了三块原理相同的驱动电路板，每一块电路板单独控制同一种波长的20路LED使其正常工作，其电路原理如图2(此电路原理图为4路LED驱动电路图，2 0路LED驱动电路为在此4路驱动电路基础上的扩展)。该LED驱动电路使用的器件包括运算放大器、精密并联基准电压源、场效应管、模拟开关以及不同阻值的电容和电阻。该电路使用的运算放大器LM321是一款低功耗单运算放大器，压摆率为0.4V/μs，满足该驱动电路切换速度的需要。配合该运放使用的是低功耗、低成本、低温度系数的2.5V精密并联基准电压源LM336，其起到了基准电压源缓冲和电流开关控制的双重作用。模拟开关为ADI公司生产的2.5Ω导通电阻、四通道单刀双掷、19ns开关时间的高速兼容TTL和CMOS电平的模拟开关ADG734。电源模块与驱动电路板接通时，多通道方波信号调制模块产生的方波调制信号通过控制器端输入驱动电路，在方波信号的调制下，4通道单刀双掷模拟开关实现切换，控制基准电压源与运算放大器同相放大电路相连或断开。当方波信号为高电平时，基准电压源与运算放大器同相放大电路接通，由于开环增益作用，使得位于运算放大器反馈环路内的场效应管导通，此时，运算放大器正负极电压相等，处于“虚短”状态，这将强迫电路中电流设置电阻51R上的电压刚好等于基准电压(最大为2.5V，可通过电位器调节)，此时LED所处电路导通产生近红外光波信号；当方波信号为低电平时，基准电压源与电流设置电阻51R断开，并将此电阻接地，进而使LED所处电路断开连接，实现了LED输出与调制信号频率相同的方波光信号。同时，通过调节电位器可以调节LED输出功率使其达到系统要求(1mW-10mW)。
(2)多通道方波信号发生器模块由所述FPGA实现，用于输出不同频率(6.2kHz到18kHz、且以0.2kHz等间隔分布的)60路方波调制信号分别对各个LED光源进行调制，每一束源光纤由分别接通到660nm,780nm和830nm波长LED的三根单模光纤通过共用一个塑料护套耦合而成；共用20束这样的源光纤将60路光源分成20束传导到20个光源点；该60路方波调制信号也作为12通道×60路可变选通锁相光子计数检测模块中锁相解调的参考信号。由于系统对频率精确度要求较高，多通道方波信号发生器模块中60种不同频率方波信号的产生是基于相位累加器的原理，由于不同频率对应不同的频率控制字K，根据频率控制字K的变化，产生所需要频率的方波信号。
(3)20束源光纤，每束源光纤为一端是三根单模光纤，另一端通过采取三根纤芯共用一个塑料护套的形式将三根单模光纤耦合为一束光纤。每3种经过调制的LED光源产生的不同波长(660nm、780nm、830nm)的近红外方波光信号经由源光纤输出，作用到源-探光纤布配贴片的源点上。
(4)电源模块：用于整个光源部分的供电(其输出电压为5V)。
(二)源-探光纤布配阵列部分：包括源-探光纤布配贴片、20束源光纤和12根探测光纤。源-探光纤布配阵列是采用面向脑枕部视觉功能区探测的源-探布配方案(如图3)，它采用源-探点空间交叉布置的方式，排布20个源点和12个探测点在被测组织体表面的位置用于布配源光纤和探测光纤的位置，20个源点的布置方案是按照4行5列以19mm距离等间距呈网格状布置，12个探测点的布置方案是：在每4个相邻源点的几何中心布置一个探测点，从而实现光源从不同源点入射和探测器从不同的探测点接收出射光。
20束源光纤一端与光源系统LED相连，另一端与源-探光纤布配贴片的20个不同的源点相连，用于传导照射到组织体表面不同源位置的入射光；12根探测光纤的一端与源-探光纤布配贴片的12个探测点相连，另一端分别与探测部分的12个PMT光子计数器对应相连，用于传导由组织体表面不同探测位置出射的漫反射光。针对每个源点，将与该源点距离分别为L1＝13.4mm、L2＝30mm、L3＝40mm和L4＝46.5mm的所有探测点与该源点之间形成的通道依次定义为第一近邻通道、第二近邻通道、第三近邻通道和第四近邻通道；共计包括有44个第二近邻通道，在所有第二近邻通道的中点处均设有一采样点。第一近邻通道(NN1)信号主要对头皮层敏感，将用于头皮层干扰信号的抑制；第二近邻采样通道(NN2)探测深度可达脑皮层，将用于获取44个采样点(见图3)功能成像数据。对于NIRS测量，假设探测NN1通道测量时达到最大可探测光强(满足单光子计数条件)，则依据光子计数模式最小92dB的动态范围，最大探测距离可达到>50mm。因此NIRS测量将采用所有源-探距离小于60mm的采样通道(NN1-NN4)数据，其中NN1通道的短距采样数据仅用于抑制头皮层等背景干扰。
(三)探测单元部分：包括12根探测光纤，12个PMT光子计数器和12通道×60路可变选通锁相光子计数检测模块及电源模块。
(1)PMT光子计数器：PMT具有单光子计数的功能，其工作在单光子计数状态，12个PMT分别与12根探测光纤对应连接，将接收到的出射光转换为电脉冲信号，便于下一步的处理。
(2)探测光纤：共有12根探测光纤，12根探测光纤接收由组织体表面从所述12个探测点位置处出射的扩散光并将其传导到所述12个PMT光子计数器中。
(3)12通道×60路可变选通锁相光子计数检测模块：其中的每通道对应所述12个探测点中的一个探测点，每路对应所述60个LED的调制光源中的一个，12通道×60路可变选通锁相光子计数检测模块将采集到的信号进行数字锁相解调，将来自不同频率调制下不同波长的混合光源信号进行分离提取，得到每种调制频率下光源对应的出射光强信息，该模块包括数字参考信号接收模块和数字锁相解调模块，针对单个PMT锁相光子计数通道，锁相光子计数的实现过程如图4，由于每个PMT锁相光子计数通道接收到的为60种频率调制下的混合光源信号，为解调混合调制信号，每个锁相光子计数通道需并行设计60路数字相敏检测器(PSD)子模块，其接收多通道方波信号发生器模块产生的60种不同频率的同相和正交数字参考信号并各自分配给60路数字相敏检测器(PSD)子模块，并将其锁存，分别与60个调制频率对应。数字锁相检测器每检测到一个PMT电脉冲信号的上升沿，就对计数门宽信号和参考信号进行判断；若门宽信号和参考信号都为高电平，则累加器加1；若门宽信号为高电平而参考信号为低电平，则累加器减1。单次累加时间(门宽信号高电平持续时长)内，根据Vs×Vr＝|I|+|Q|(Vr＝1)，可得与该参考信号频率相同的调制光与脑组织作用后出射的光强信息。基于FPGA的并行工作特性，12通道×60路可变选通锁相光子计数检测模块可以同时获取12个探测点处分别对应60路光源的出射光强信息，即可实现全并行快速测量。
(4)电源模块：用于PMT光子计数器的供电(其输出电压为5V)。
(四)计算机控制和数据处理部分：采用RS-232串口通过LabVIEW可视化编程软件实现计算机与FPGA的双向通信：所述计算机控制及数据处理单元向多通道方波信号发生器模块发送频率控制字控制其产生60中不同选定频率的方波调制信号；同时，向12通道×60路可变选通锁相光子计数检测模块发送门宽控制字控制单次测量时间(累加器累加时间)；等到测量过程结束后通过串口读取暂存在FPGA内部的12个探测点处分别对应60路光源的出射光强信息。最后，利用扩散光学层析重建算法结合所测得的组织体边界光流量构建反映脑皮层光学参数变化量空间分布的图像。
测量实例：在人脑处于静息和任务状态时候时，分别用本发明系统测得各源探组合下出射光强信息并且用之重建反映脑皮层光学参数变化量的空间分布图的过程如图5所示：
(1)通信测试：测试计算机和FPGA之间的串口双向通信是否正常；
(2)调制信号生成：计算机通过串口向FPGA发送控制字，FPGA内部实现的多通道方波信号发生器生成以0.2kHz等间隔的6.2～18kHz共60路方波调制信号；另外，这60路方波调制信号也用作检测模块中的参考信号。
(3)对LED进行方波调制：对20个660nm的LED、20个波长为780nm的LED和20个波长为830nm的共60个LED分别用60路不同频率的调制信号进行调制。
(4)源光纤传导：20束源光纤将60路光源分成20组混合光源传导到源-探布配阵列上的20个光源点处。
(5)探测光纤传导：从源-探布配阵列上的12个探测点处接收的组织体出射的扩散光经由12根探测光纤分别传导到12个PMT光子计数器中。
(6)光电转换：PMT光子计数器对接收的微弱光信号进行光电转换，逐级放大，电路整形最终输出规则的方波电脉冲信号。
(7)方波电脉冲上升沿检测以及信号解调：FPGA利用“消抖法“对PMT光子计数器输出的电脉冲信号进行上升沿精准检测。在FPGA内部实现的12通道×60路可变选通锁相光子计数检测模块对混合调制光进行解调，得到12个探测点处分别对应60个光源的出射光强信息。
(8)数据传输：12通道×60路可变选通锁相光子计数检测模块中所测得的数据通过串口传输到计算机中储存。
(9)数据处理：对处于静息和任务状态的人脑分别重复以上过程测得各源探组合下的出射光强信息，结合扩散光学层析重建算法对人脑的光学参数变化量进行重建，从而重建反映脑皮层光学参数变化量的空间分布图。
尽管上面结合附图对本发明进行了描述，但是本发明并不局限于上述的具体实施方式，上述的具体实施方式仅仅是示意性的，而不是限制性的，本领域的普通技术人员在本发明的启示下，在不脱离本发明宗旨的情况下，还可以做出很多变形，这些均属于本发明的保护之内。

Claims (2)

1.一种基于全并行锁相光子计数检测模式的NIRS脑功能成像系统，包括源-探光纤布配阵列单元、可调制LED光源单元、探测单元、现场可编程门阵列(FPGA)和计算机控制及数据处理单元；
所述源-探光纤布配阵列单元包括源-探光纤布配贴片、源光纤和探测光纤，源-探光纤布配贴片采用源-探点空间交叉布置的方式，源-探光纤布配贴片与20个光源点和12个探测点的位置对应，用于布配源光纤和探测光纤的位置，以实现光源从不同光源点入射和探测器从不同的探测点接收出射光；
所述可调制LED光源单元包括LED光源模块、多通道方波信号发生器模块和20束源光纤，所述LED光源模块包括60个LED和LED光源驱动电路，所述60个LED包括20个波长为660nm的LED、20个波长为780nm的LED和20个波长为830nm的LED；所述多通道方波信号发生器模块由FPGA实现，用以输出频率为从6.2kHz到18kHz、且以0.2kHz等间隔分布的60路方波调制信号对60个LED分别进行调制；每一束源光纤由分别接通到660nm,780nm和830nm波长LED的三根单模光纤通过共用一个塑料护套耦合而成；共有20束这样的源光纤将60路光源分成20束传导到20个光源点；
所述探测单元包括12根探测光纤，12个PMT光子计数器和12通道×60路可变选通锁相光子计数检测模块，其中，12根探测光纤用于将所述12个探测点位置处的扩散光分别传导到所述12个PMT光子计数器；所述12个PMT光子计数器用于将接收到的光信号转化为电脉冲信号；所述12通道×60路可变选通锁相光子计数检测模块由FPGA实现，其中，每通道对应所述12个探测点中的一个，每路对应所述60个LED的调制光源中的一个，从而使所述FPGA同时获取12个探测点处分别对应60路光源的出射光强信息，实行全并行快速测量；
所述计算机控制及数据处理单元向多通道方波信号发生器模块发送频率控制字控制其产生选定频率的方波调制信号；同时，向12通道×60路可变选通锁相光子计数检测模块发送门宽控制字控制单次测量时间；等到测量过程结束后通过串口读取12个探测点处分别对应60路光源的出射光强信息，最后，利用扩散光学层析重建算法构建反映脑皮层光学参数变化量的空间分布图；其特征在于，
所述12通道×60路可变选通锁相光子计数检测模块将采集到的信号进行数字锁相解调，将来自不同频率调制下不同波长的混合光源信号进行分离提取，得到每种调制频率下光源对应的出射光强信息，该模块包括数字参考信号接收模块和数字锁相解调模块；每个PMT锁相光子计数通道接收到的为60种频率调制下的混合光源信号，为解调混合调制信号，每个锁相光子计数通道需并行设计60路数字相敏检测器(PSD)子模块，其接收多通道方波信号发生器模块产生的60种不同频率的同相和正交数字参考信号并各自分配给60路数字相敏检测器(PSD)子模块，并将其锁存，分别与60个调制频率对应；数字锁相检测器每检测到一个PMT电脉冲信号的上升沿，就对计数门宽信号和参考信号进行判断；若门宽信号和参考信号都为高电平，则累加器加1；若门宽信号为高电平而参考信号为低电平，则累加器减1；单次累加时间内，根据Vs×Vr＝|I|+|Q|(Vr＝1)，可得与该参考信号频率相同的调制光与脑组织作用后出射的光强信息，基于FPGA的并行工作特性，12通道×60路可变选通锁相光子计数检测模块可以同时获取12个探测点处分别对应60路光源的出射光强信息，即可实现全并行快速测量。
2.根据权利要求1所述基于全并行锁相光子计数检测模式的NIRS脑功能成像系统，其特征在于，所述源-探光纤布配阵列单元中采用源-探点空间交叉布置方式布配贴片，用于确定20个光源点和12个探测点的位置；
20个源点的布置方案是：按照4行5列以19mm距离等间距呈网格状布置；
12个探测点的布置方案是：在每4个相邻源点的几何中心布置一个探测点；
针对每个源点，将与该源点距离分别为L1＝13.4mm、L2＝30mm、L3＝40mm和L4＝46.5mm的所有探测点与该源点之间形成的通道依次定义为第一近邻通道、第二近邻通道、第三近邻通道和第四近邻通道；其中包括有44个第二近邻通道，在所有第二近邻通道的中点处均设有一采样点。
CN201710466642.4A 2017-06-18 2017-06-18 基于全并行锁相光子计数检测模式的nirs脑功能成像系统 Active CN107252305B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710466642.4A CN107252305B (zh) 2017-06-18 2017-06-18 基于全并行锁相光子计数检测模式的nirs脑功能成像系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710466642.4A CN107252305B (zh) 2017-06-18 2017-06-18 基于全并行锁相光子计数检测模式的nirs脑功能成像系统

Publications (2)

Publication Number Publication Date
CN107252305A CN107252305A (zh) 2017-10-17
CN107252305B true CN107252305B (zh) 2020-11-10

Family

ID=60024041

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710466642.4A Active CN107252305B (zh) 2017-06-18 2017-06-18 基于全并行锁相光子计数检测模式的nirs脑功能成像系统

Country Status (1)

Country Link
CN (1) CN107252305B (zh)

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108113652B (zh) * 2017-12-19 2020-10-23 天津大学 用于近红外脑功能成像仪的强度自动调节光源系统及方法
CN108042131A (zh) * 2017-12-26 2018-05-18 武汉智普天创科技有限公司 脑电和近红外光谱联合采集脑信号的便携式帽子
CN108334195B (zh) * 2018-01-17 2019-10-18 西安交通大学 基于调制的生物运动视觉感知的脑-机接口方法
CN109596525B (zh) * 2018-10-23 2021-08-03 浙江亨达光学有限公司 一种检测组织活性的实时测量方法和仪器
CN109770849A (zh) * 2018-12-04 2019-05-21 天津大学 一种基于锁相光子计数和压缩感知的单像素sfdi系统
CN110123281A (zh) * 2019-05-27 2019-08-16 天津大学 基于锁相光子计数技术的并行激励扩散光学层析成像装置
CN111466874B (zh) * 2020-03-12 2021-08-10 西安电子科技大学 一种基于方波调制的扩散光学层析成像系统
CN111134626B (zh) * 2020-01-13 2021-08-24 西安电子科技大学 一种无光纤平板式乳腺扩散光层析成像系统
CN111481170A (zh) * 2020-03-12 2020-08-04 西安电子科技大学 基于手持式的扩散光学探头及输出电压校正的方法
CN112987919B (zh) * 2021-02-07 2023-11-03 江苏集萃脑机融合智能技术研究所有限公司 一种基于间接测量飞行时间技术的脑机接口系统以及实现方法
CN112629683B (zh) * 2021-03-09 2021-07-13 成都成电光信科技股份有限公司 一种单光子计数装置及其方法
CN115426265B (zh) * 2022-11-02 2023-04-18 之江实验室 一种多模态网络下交换资源分配优化方法及装置、介质

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0497021A1 (en) * 1991-01-31 1992-08-05 Hamamatsu Photonics K.K. Oximeter with monitor
JP2006247253A (ja) * 2005-03-14 2006-09-21 Shimadzu Corp 近赤外脳機能イメージング装置
CN106037657A (zh) * 2016-06-28 2016-10-26 丹阳慧创医疗设备有限公司 一种时空频多重耦合的高密度近红外光谱脑功能成像方法

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN100409806C (zh) * 2005-08-09 2008-08-13 中国科学院自动化研究所 正六边形高分辨率近红外光谱脑功能成像头盔
CN101832815B (zh) * 2010-04-30 2011-09-14 天津大学 基于数字锁相解复用的多通道单光子计数测量系统
CN101849821B (zh) * 2010-06-13 2012-07-04 华中科技大学 一种光纤近红外光谱检测仪
CN102327111A (zh) * 2011-09-23 2012-01-25 天津析像光电科技有限公司 一种多通道近红外脑功能成像并行检测系统
CN102599888A (zh) * 2012-03-02 2012-07-25 天津大学 基于数字锁相检测技术的光学拓扑成像系统及方法
CN103156620B (zh) * 2013-03-27 2015-01-07 中国科学院自动化研究所 一种多通道并行近红外光谱成像系统

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0497021A1 (en) * 1991-01-31 1992-08-05 Hamamatsu Photonics K.K. Oximeter with monitor
JP2006247253A (ja) * 2005-03-14 2006-09-21 Shimadzu Corp 近赤外脳機能イメージング装置
CN106037657A (zh) * 2016-06-28 2016-10-26 丹阳慧创医疗设备有限公司 一种时空频多重耦合的高密度近红外光谱脑功能成像方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于多通道数字锁相系统的脑功能光学成像相关研究;孟伟;《中国优秀硕士学位论文全文数据库 医药卫生科技辑》;20150515(第5期);论文摘要、论文正文第12-15、20、30、44-46、50-52页 *

Also Published As

Publication number Publication date
CN107252305A (zh) 2017-10-17

Similar Documents

Publication Publication Date Title
CN107252305B (zh) 基于全并行锁相光子计数检测模式的nirs脑功能成像系统
Koizumi et al. Higher-order brain function analysis by trans-cranial dynamic near-infrared spectroscopy imaging
CA2319456C (en) Imaging and characterization of brain tissue
CN103156620B (zh) 一种多通道并行近红外光谱成像系统
US20080154126A1 (en) High Performance Imaging System for Diffuse Optical Tomography and Associated Method of Use
CN102327111A (zh) 一种多通道近红外脑功能成像并行检测系统
Gratton et al. The event-related optical signal: a new tool for studying brain function
CN102495027B (zh) 混沌介质光学参数多波长空间分辨锁相测量方法与装置
CN104146714A (zh) 一种组织体局部血氧饱和度变化量拓扑成像装置及方法
CN102599888A (zh) 基于数字锁相检测技术的光学拓扑成像系统及方法
Zhang et al. A fast neuronal signal-sensitive continuous-wave near-infrared imaging system
CN110123281A (zh) 基于锁相光子计数技术的并行激励扩散光学层析成像装置
Chiarelli et al. Characterization of a fiber-less, multichannel optical probe for continuous wave functional near-infrared spectroscopy based on silicon photomultipliers detectors: in-vivo assessment of primary sensorimotor response
US6907279B2 (en) Optical system for measuring metabolism in a body
CN210408435U (zh) 一种近红外脑氧检测系统
CN116784841A (zh) 用于脑机接口的新型光纤漫射时间分辨探针系统
CN109998559A (zh) 一种近红外脑氧检测系统及方法
JP2007111461A (ja) 生体の光計測装置
Torricelli et al. Functional brain imaging by multi-wavelength time-resolved near infrared spectroscopy
Chiarelli et al. Wearable, fiber-less, multi-channel system for continuous wave functional near infrared spectroscopy based on silicon photomultipliers detectors and lock-in amplification
Wojtkiewicz et al. Towards optical tomography of an adult human head
Jelzow In vivo quantification of absorption changes in the human brain by time-domain diffuse near-infrared spectroscopy
CN108113652B (zh) 用于近红外脑功能成像仪的强度自动调节光源系统及方法
Zhang et al. Experimental comparison of using continuous-wave and frequency-domain diffuse optical imaging systems to detect heterogeneities
Wabnitz et al. Time-domain diffuse optical imaging of tissue by non-contact scanning

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
TR01 Transfer of patent right

Effective date of registration: 20231219

Address after: Unit 70, 1st Floor, Building A19, No. 2036 Wengjiao West Road, Haicang District, Xiamen City, Fujian Province, 361028

Patentee after: Nielsen Human Factors (Xiamen) Technology Co.,Ltd.

Address before: 300072 Tianjin City, Nankai District Wei Jin Road No. 92

Patentee before: Tianjin University

TR01 Transfer of patent right