一种智能分频调制的多通道近红外光谱脑功能成像方法
【专利摘要】本发明公开了一种智能分频调制的多通道近红外光谱脑功能成像方法，涉及医学成像技术领域。所述方法包括频率段划分、初始频率生成、分空间频率调制和其他波长光源频率调制的步骤，成像时，所有光源按照调制频率同时点亮，可以显著提高并行探测效率，并提高时间分辨率；本发明根据空间距离和通道连接等参数进行分空间频率调制，考虑了光子信号强度随探测距离的增加呈指数衰减的特点，降低了光源之间频谱泄露的干扰，可以提高信噪比；本发明提出的成像方法，由于进行分空间频率调制，实现了同时多个光源的频率并行调制，并且各个频率分量之间无泄露，在多通道fNIRS成像中有无可比拟的应用优势。
【专利说明】
一种智能分频调制的多通道近红外光谱脑功能成像方法
技术领域
[0001] 本发明涉及医学成像技术领域，具体涉及一种智能分频调制的多通道近红外光谱 脑功能成像方法。
【背景技术】
[0002] 近年来，在脑功能成像方法中，近红外光谱（Funct ional near-inf rared spectroscopy，fNIRS)已发展为脑功能研究和临床诊断中不可或缺的新一代成像模态。 fNIRS利用波长在650-1000nm范围的近红外光谱，检测大脑代谢引起的光学吸收特性的改 变，计算出该区域脱氧血红蛋白（HbR)、携氧血红蛋白（HbO 2)的浓度相对变化量。相较于其 他脑功能成像方法，如脑电图、功能核磁共振成像等，fNIRS具有性价比高、时间分辨率高、 空间分辨率适中、功能参数丰富、对运动不敏感等优势。
[0003] fNIRS成像系统中，在头部相应部位放置光源探头，并在距光源探头一定距离处放 置探测器。光源发射的光(通常为双波长或三波长)经波分复用后，由光纤传输到达置于头 部的光源探头，随后穿透头皮层、颅骨等入射到脑组织，在历经一系列吸收、散射后，仍会有 一部分光子到达头皮层表面，这些光信号可以被探测器接收到，然后根据Beer-Lambert定 律，计算出HbO 2、HbR的浓度相对变化量。光源探头放置处为光源入射点，但一个光源入射点 处入射的光包含不同波长(通常为双波长或三波长），因此一个光源入射点处实际有2个(双 波长)或3个(三波长)重合的光源，为进行区分，称重合的多波长光源为光源组，而单波长光 源即简称为光源，探测器放置处为探测点，相邻的光源和探测器可组成一个探测通道，从探 测点检测的信息便可反映对应探测通道的信息。不同的光源探测器（source-detector)排 布，可组成不同的探测通道。多个光源和多个探测器的排列，可以形成多通道系统，从而获 得更多探测通道对应的大脑信号。
[0004] 目前，多通道fNIRS成像系统中，等间距探头拓扑排布被广泛应用。在等间距拓扑 结构的fNIRS成像中，主要有两种方式实现多通道光谱信号的检测。
[0005] 第一种是分时方法，在每个特定时间只点亮一个光源，通过轮询实现信号探测。在 多通道多光源的系统中，完成一次轮询需要顺序遍历所有的光源，所需的时间较长，系统的 时间分辨率的提高受到限制。
[0006] 第二种方法为频率并行调制，给不同探头(不同排布位置）的、不同波长的光源设 定不同的调制频率(通常6-12KHZ之间，方波调制），所有光源同时点亮。由于同一时刻，探测 器接收到的来自不同光源的光信号频率不相同，因而可以通过频谱提取(例如快速傅里叶 变换或带通滤波器)把来自不同光源的信号区分开来，实现同时采集多个探测通道的信号。 频率调制还可显著提高信噪比，大幅度降低来自电子电路固有的低频率Ι/f噪声以及环境 光干扰。例如，美国Techon公司的CW6产品将32个光源在6.4到12.8KHz范围均匀分配了调制 频率(激光调制频率的间隔是200Hz的整数倍）。然而，已有方法中，同时被调制的光源不宜 过多，否则各个频率分量之间会有泄露，造成探测通道间的串扰。而对于可覆盖全脑的 fNIRS成像系统，其光源数目较多，受限于噪声、数模转换器(AD)性能和成本等因素，同时频 率调制时，不同光源间的频率间隔较小，频率泄露使得相邻探测通道之间存在串扰，降低信 号采集的准确度和信噪比。
【发明内容】
[0007] 针对上述不足，本发明提出了一种智能频率调制方案。注意到，光子信号强度随探 测距离的增加呈指数衰减，大约探测距离每增加1厘米，信号强度就降低一个数量级。如果 空间距离较远的光源同时以相同的频率点亮，由于空间距离较远，彼此之间的干扰可以忽 略，因而可以分空间调制频率，实现并行探测信号。为进一步降低彼此间频谱泄露干扰，不 同空间的调制频率略有差异，距离越近偏移越大。
[0008] 本发明提供一种智能分频调制的多通道近红外光谱脑功能成像方法，主要包括如 下步骤：
[0009] 第一步，频率段划分。根据所采用光源波长个数(2个或3个)和频率范围，将频率范 围均匀划分为与波长个数相等的频率段，使得一个波长对应一个频率段，并计算频率段之 间的频率间隔；
[0010] 第二步，初始频率生成。选择其中一个波长对应的频率段，根据光源入射点总个 数，将该频率段范围内的频率均匀划分，作为初始频率;所选波长对应的光源的频率只能从 生成的初始频率中选取；
[0011] 第三步，分空间频率调制。在选中的波长中，依次遍历本波长所有光源，完成该波 长下，所有光源的频率调制；
[0012] 第四步，其他波长光源频率调制。对于相同光源入射点处其他波长的光源，由于第 三步中已经调制其中一个波长下的光源频率，其他波长光源只需加上或减去第一步中计算 的频率间隔即可得出其对应频率，从而完成所有光源的频率调制。
[0013] 第五步，所有光源同时点亮，其发光频率为第四步计算出的调制频率，随后进行成 像。
[0014] 本发明的优点是：
[0015] (1)本发明提出的智能分频调制的多通道近红外光谱脑功能成像方法，根据空间 距离和通道连接等参数进行分空间频率调制，考虑了光子信号强度随探测距离的增加呈指 数衰减的特点，降低了光源之间频谱泄露的干扰，可以提高信噪比；
[0016] (2)本发明提出的成像方法中，所有波长下的所有光源同时点亮，可以显著提高并 tx探测效率，并提尚时间分辨率；
[0017] (3)本发明提出的成像方法，由于进行分空间频率调制，实现了同时多个光源的频 率并行调制，并且各个频率分量之间无泄露，在多通道fNIRS成像中有无可比拟的应用优 势。
【附图说明】
[0018] 图1.光源探测器排布图。
[0019] 图2.干扰距离的定义和计算。
[0020]图3A.光子信号强度随距离的增加呈指数衰减曲线。
[0021]图3B.表示光源5"对其他光源的干扰强度随干扰距离呈指数衰减曲线。
[0022] 图4.确定光源频率。
[0023] 图5.本发明提供的成像方法流程图。
[0024]图6.fNIRS成像系统硬件控制示意图。
【具体实施方式】
[0025] 下面结合附图对本发明提出的智能分频调制的多通道近红外光谱脑功能成像方 法做进一步的说明。
[0026] 以调制频率范围在6-12KHZ之间的三波长（例如785nm，808nm，850nm)fNIRS成像仪 器为例，包含72个光源，每个光源入射点处有3个不同波长的光源重合，因此有24个可见光 源分布在24个不同的光源入射点上，如图1所示，空心四边形为光源探测器，实心黑色圆代 表光源入射点，光源入射点和光源探测器之间的连线代表探测通道。为方便叙述，选择 785nm波长下的24个光源为例，将光源分别表不为Si，S2,......，S24。
[0027] 采用本发明提供的智能分频调制的多通道近红外光谱脑成像方法，如图5所示流 程，具体实现步骤如下：
[0028] 第一步，频率段划分。根据光源波长个数(2个或3个)和光源频率范围，将频率均匀 划分为R个频率段，R与光源波长个数相等;使得一个波长对应一个频率段，并计算频率段之 间的频率间隔。
[0029] 本例的fNIRS成像系统调制频率在6-12KHZ之间，且为三波长，则在6-12KHZ将频率 均匀划分为3个频率段：6-8KHZ、8-10KHZ、10-12KH Z。使得785nm的光源的频率段对应为6-8KHz，808nm的光源的频率段对应为8-10KHz，850nm的光源的频率段对应为10-12KHz，频率 段之间的频率间隔为2KHz。因此，在同一个光源入射点处，只需确定785nm光源的频率段，将 频率偏移2KHz，得到808nm光源的频率段;偏移4KHz，可得到850nm光源的频率段。
[0030] 第二步，在选定的频率段内，生成初始频率。选择光源中一个波长对应的频率段， 根据光源入射点总个数，根据公式
[0031]
[0032] 计算得出频率间隔μ，其中ε为正数，N为选定频率段内光源总个数，F2为选定频率 段内的最大频率，F 1为选定频率段内的最小频率，并且Ν+ε的取值应使μ为整数，将该频率段 范围内的频率均匀划分，作为初始频率，所选波长的其他光源的频率只能从生成的初始频 率中选取。
[0033] 本实施例选择785nm波长下的24个光源，将光源分别表不为Si，S2，......，S24，如图1 所示，其对应的频率段为6-8KHz，光源频率的初始间隔设置为，
[0034]
[0035] 其中，N为785nm光源总个数，这里N = 24，ε为正数，为使光源频率易于调节，Ν+ε的 取值应使得μ为整数，这里取ε = 16，得μ = 0.05ΚΗζ。因此，6-8ΚΗΖ的频率段内以μ为间隔生成 40个不同的频率，如6KHz，6.05KHz，6. IOKHz……，7.95ΚΗζ，这40个频率就是初始频率。所有 785nm光源的频率只能从生成的40个初始频率中选取。
[0036]第三步，分空间频率调制，在选中的波长中，依次遍历本波长所有光源，完成该波 长下，所有光源的频率调制，具体为：
[0037]计算光源3"对其它光源的干扰距离。在选中的波长中，首先，对于第一个光源，根 据空间距离、探测通道连接等参数，计算该光源与其他光源之间的最小偏移频率，由此可在 初始频率中选择第一个光源的发光频率;然后根据该方法，获得第二个光源的发光频率;同 样的，依次遍历本波长所有光源，完成该波长下，所有光源的频率调制。
[0038] 本例中，分空间调制波长为785nm的所有光源，即光源5:到光源&4,按光源编号从 小到大的顺序进行频率调制。现以光源Sn为例，为说明分空间调制频率的计算方法，本发明 提出干扰距离的概念，如图2所示。在一个确定的光源探测器排布中，若某一频率段（以 785nm波长下的频率段为例）中，光源总数为N，对于任意光源S n，其中KnSN，光源Sn对其他 光源的干扰距离定义为:光源3"对某一频率已经确定的光源S 1的信号干扰，表现为对光源51 周围M(通常为四）个探测器检测到的信号的干扰，计算光源Sn与光源S1周围M个光源探测器 之间的距离，M个距离中的最小距离即为光源5"对光源S 1的干扰距离，记为dn(i)，干扰距离 与空间距离和探测通道连接等有关，需要注意的是，若某一光源还未确定频率，则认为光源 Sn对该光源没有干扰，干扰距离为无穷大。光源Sn对光源探测器排布中其他N-I个光源的干 扰距离分别记为dn(l)，d n(2)，......，dn(i)，......，dn(N-l)。
[0039] 易得，光源Sn对本例中其他23个光源的干扰距离分别为dn( I)，dn(2)，……，dn (1) ......，dn( 23)。然后根据干扰距离和最小频率偏移之间的关系：
[0040] Δ fn(i)=-Blg(Adn(i))
[00411其中，η = 1，2，……，N，i = 1，2，……，N-I川为选定频率段下光源总个数，A和B是不 为零的常量，dn(i)表示光源3"对光源S1的干扰距离，Afn(i)表示光源3"对光源最小频 率偏移。
[0042]如图3A和3B所示，光子信号强度随距离的增加呈指数衰减，光源Sn对其他光源的 干扰强度随干扰距离呈指数衰减，干扰距离越小，最小频率偏移越大。计算光源3"对其他光 源的最小频率偏移：Afn(I)，Af n(2)，……，Afn(i)，……Afn(23)。
[0043] 如图4所示，在一个确定的光源探测器排布中，若光源（以785nm波长的光源为例） 总数为N，发光频率分别为，f2，……，A，……，ft，对于任意光源S n，可根据干扰距离的定义 分别计算出光源3"对光源探测器排布中其他N-I个光源的最小频率偏移△ fn(l)，Δ fn (2) ，……，Δ fn(i)，……Δ fn(N-l).如在6-8kHz(785nm波长对应的频率段）的频率范围内， 其频率分布如频率轴所示。由于光源S n对光源S1的最小频率偏移为△ fn(l)，光源S1的频率 记为fi，则在频率轴上可以形成以心为圆心，Af n(I)为半径的圆，当光源5"取该圆外的频率 时，对光源31的信号基本无干扰。因此，在确定光源5"的频率时，可以在6-8KHz的频率轴上分 别作其他23个以自身光源频率为圆心，Af n(l)，Afn(2)，……，Afn(i)，……Afn(23)为半 径的圆。意味着所有圆以内的频率不可取，而这些圆之外的频率均符合条件，可选择符合条 件的任一初始频率作为光源S n的发光频率。圆之外的任一频率fn均满足：
[0044] fn-fi| ^ Afn(i)
[0045] 可在初始频率中为光源5"寻找符合要求的发光频率。遍历该波长所有光源，完成 光源频率的智能调制。
[0046]第四步，其他波长光源频率调制。在相同光源入射点处，785nm的光源频率已经根 据第三步确定，将频率偏移2KHz得到808nm光源的频率，偏移4KHz得到850nm光源的频率，遍 历所有光源入射点，从而完成所有光源的频率调制，使得距离越近的光源频率偏移越大，降 低了彼此间频谱泄露干扰。
[0047]第五步，所有光源同时点亮，其发光频率为第四步计算出的调制频率，随后进行成 像。
[0048]上述的智能频率调制方法依fNIRS成像系统实现。该fNIRS成像系统主要由以下部 分组成:光源、探测器、信号滤波及A/D转换模块、底层信号处理模块、高速传输接口和计算 机PC，如图6所示。光源可以采用激光二极管或LED(发光二极管）。探测器可以采用PMT(光电 倍增管）、APD(雪崩二极管)或ro(光电二极管）等。探测器采集到的光源发出的信号经信号 滤波及A/D转换模块转换后，使用底层信号处理模块进行底层信号处理，包括数据解调等， 并使用高速传输接口进行数据传输。所述的高速传输接口可采用千兆网、USB3.0等。所述的 底层信号处理模块采用FPGA(现场可编程门阵列）或DSP(数字信号处理技术)芯片。然后利 用计算机PC或直接利用底层FPGA、底层DSP等底层信号处理模块对数据进行解调重建，从而 得到有效的脑功能血流参数。用户的操作指令也可由计算机发送至高速传输接口，经FPGA、 DSP芯片等底层信号处理后驱动光源，实现对光源的实时控制。
[0049]本发明提出的智能分频调制的多通道近红外光谱脑功能成像方法，根据空间距离 和通道连接等参数进行分空间调制频率，考虑了光子信号强度随探测距离的增加呈指数衰 减的特点，降低了光源之间频谱泄露的干扰，提高了系统的信噪比。系统中所有光源可以同 时点壳，能显者提尚并彳丁探测效率，并提尚系统的时间分辨率。由于进彳丁分空间频率调制， 实现了同时多个光源的频率并行调制，并且各个频率分量之间无泄露，在多通道fNIRS系统 中有无可比拟的应用优势。
【主权项】
1. 一种智能分频调制的多通道近红外光谱脑功能成像方法，其特征在于:包括如下步 骤， 第一步，频率段划分； 根据所采用光源入射点处的波长个数和频率范围，将频率范围均匀划分为与波长个数 相等数量的频率段，使得一个波长对应一个频率段，并计算频率段之间的频率间隔； 第二步，初始频率生成； 选定其中一个波长对应的频率段，根据光源入射点总个数，将该频率段范围内的频率 均匀划分，作为初始频率;所选波长对应的光源的频率只能从生成的初始频率中选取； 第三步，分空间频率调制； 在选中的波长中，依次遍历本波长下所有光源，为每一个光源选择发光频率，完成所有 光源的频率调制； 第四步，其他波长光源频率调制； 对于相同光源入射点处其他波长的光源，由于第三步中已经调制其中一个波长下的光 源频率，其他波长光源只需加上或减去第一步中计算的频率间隔即可得出其对应频率，从 而完成所有光源的频率调制； 第五步，所有光源按照各自的调制频率同时点亮，进行成像。2. 根据权利要求1所述的一种智能分频调制的多通道近红外光谱脑功能成像方法，其 特征在于:所述的最小偏移频率根据干扰距离得到，具体为， Δ fn(i) =-Blg(Adn(i)) 其中，n=l，2,……，N，i = l，2,……，N-1;N为选定频率段下光源总个数，A和B是不为零 的常量，dn(i)表示光源5"对光源Si的干扰距离，△ fn(i)表示光源5"对光源Si的最小频率偏 移。3. -种实现权利要求1所述的智能分频调制的多通道近红外光谱脑功能成像方法的 fNIRS成像系统，其特征在于:fNIRS成像系统包括光源、探测器、信号滤波及A/D转换模块、 底层信号处理模块、高速传输接口和计算机，探测器采集到的光源发出的信号经信号滤波 及A/D转换模块转换后，使用底层信号处理模块进行底层信号处理，并使用高速传输接口进 行数据传输给计算机;然后利用计算机或直接利用底层信号处理模块对数据进行解调重 建，从而得到脑功能血流参数。4. 根据权利要求3所述的一种实现权利要求1所述的智能分频调制的多通道近红外光 谱脑功能成像方法的fNIRS成像系统，其特征在于:光源采用激光二极管或发光二极管;探 测器采用光电倍增管、雪崩二极管或光电二极管;所述的底层信号处理模块采用现场可编 程门阵列FPGA或数字信号处理技术DSP芯片。
【文档编号】G02F2/00GK106073716SQ201610490466
【公开日】2016年11月9日
【申请日】2016年6月28日
【发明人】汪恭正
【申请人】丹阳慧创医疗设备有限公司