首页 / 专利库 / 电气元件和设备 / 介质强度恢复 / 单次蛛网膜下腔阻滞术中信息收集处理系统及其处理方法

单次蛛网膜下腔阻滞术中信息收集处理系统及其处理方法

阅读:940发布:2020-11-09

专利汇可以提供单次蛛网膜下腔阻滞术中信息收集处理系统及其处理方法专利检索,专利查询,专利分析的服务。并且本 发明 属于医疗技术领域,公开了一种单次蛛网膜下腔阻滞术中信息收集处理系统及其处理方法,信息收集处理系统包括: 图像采集 模 块 、脑电波采集模块、主控模块、弛豫值测量模块、数据映射模块、数据存储模块、显示模块。本发明通过图像 数据处理 模块大大缩短了T1图像采集时间,因而不易引起用户不适,减少采集图像中的运动伪影;同时通过数据映射模块将用户大脑的MRI或CT用作针对映射的解剖模型。一旦MRI数据已被汇编,那么通过知道用户脑部的确切 位置 和方位以及每个TMS脉冲,就可以将用户脑电 波数 据图像准确地映射到解剖模型;降低 风 险。,下面是单次蛛网膜下腔阻滞术中信息收集处理系统及其处理方法专利的具体信息内容。

1.一种单次蛛网膜下腔阻滞术中信息收集处理方法,其特征在于,所述单次蛛网膜下腔阻滞术中信息收集处理方法包括:
采集单次蛛网膜图像数据信息;结合采集单次蛛网膜图像影响因素构建神经网络模型BP,产生神经网络模型BP的初始权值;
对神经网络模型BP的权值和阈值进行动态改进,获得动态神经网络模型DBP,产生动态神经网络模型DBP的权值和阈值;具体包括:调整神经网络模型BP隐含层与输出层之间的权值wkj;
*
调整wkj的目的是希望输出节点j的新输出opj比当前输出opj更接近目标值tpj,定义:
其中α代表接近度,在每个训练周期保持不变,并随隐含层节点数H的调整而变小,不考虑阈值,则有:
其中wkj和w*kj分别为更新前后的权值,ypk为隐含层输出,Δwkj为wkj的改变量;
根据公式 得到Δwkj的求解方程:
其中,
根据最小平方和误差原则求解方程 得到Δwkj的近似解:
对每一个连接到输出节点j的隐含层节点k,计算k与j之间的权值变化Δwkj,更新权值并计算平方和误差E,然后在k∈[1,H]区间选择一个最优的k,使得E最小;
调整神经网络模型BP输入层与隐含层之间的权值vik;
调整vik的目的是一旦神经网络算法陷入局部极小点,修改权值能够跳出该极小点,判断神经网络算法陷入局部极小点的条件是误差E的变化率ΔE=0,且E>0;
不考虑阈值,隐含层节点k的权值的改变通过以下方程求解:
其中δpj=f-1(ypk+Δypk)-f-1(ypk),M为自然数,则隐含层输出ypk求解公式为:
其中Δypk为ypk的改变量,则有:
根据最小平方和误差原则求解公式 构建的矩阵方程,可以算
出:
计算隐含层与输出层之间权值的动态平均改变量
计算输入层与隐含层之间权值的动态平均改变量
公式中M取10~20之间的自然数,获得神经网络模型BP的动态平均权值,根据神经网络模型BP的动态平均权值获得动态神经网络模型DBP;
运用自适应免疫遗传AIGA算法优化动态神经网络模型DBP,获得预测模型AIGA-DBP,根据预测模型AIGA-DBP计算采集单次蛛网膜图像预测真实数据;
采集单次蛛网膜下腔阻滞术中脑电波数据信息;重建脑电波数据衰减系数;
采集模偏移一定度放置时采集投影数据;重建脑电波数据散射系数;计算脑电波的吸收系数,获得准确脑电波数据图像信息;
重建脑电波数据散射系数包括:
表示Γ-上入射的光子位置 且方向为 时所传播的距离,则:
其中 为 处的辐射度,表示单位立体角内、单位时间内、通过垂直于单位方向矢量的单位面积上的平均功率通量密度,量纲是W/(m2.Sr);K0为引入的光子弹道传播算子,表示脑电波数据在 处的总的衰减系数, 表示脑电波数据的吸收
系数, 表示散射系数;
再:
其中K为引入的光子散射传播算子, 为归一化的散射相位函数,表示光子从方向v′散射到 方向的概率,满足 dΩ′表示单位方向矢量v′对应的立体角微元;
定义m0=K0gin, 则有nn+1=Kmn(n≥0),从而在 处的总的辐射率为:
其中mn表示经n次散射到达 的辐射度分量;当光子在扩散区域传播时,K的谱半径ρ(K)值接近1,当光子在弱散射区域传播时,ρ(K)>>1,在此种情况下,当n→∞时, 快速收敛;
然后,在输出边界Γ+上描述探测器接收到的数据总量gout,即 从而:
其中A是描述光子传输的矩阵,A0、A1和A2分别描述弹道传输、一次散射传输和多次散射传输部分,定义g0=A0gin,g1=A1gin分别表示测量数值中的弹道传输分量和一次散射分量,则知:
入射光方向为 经一次散射后其方向为 则上式中关于 的积分仅在一特定角度上有值,取 系数k的取值由相位函数 确定,同时定义
分别表示光子发生散射后
和散射前的衰减量,则有:
通过弛豫值测量方法测量大脑纵向弛豫值数据信息;
将用户脑电波数据图像准确地映射到解剖模型;
对采集的数据信息进行存储;并显示采集的图像及脑电波数据信息。
2.如权利要求1所述的单次蛛网膜下腔阻滞术中信息收集处理方法,其特征在于,产生神经网络模型BP初始权值的方法为:
将权值初始化为[a,b]之间的随机数,其中a,b为满足以下方程的整数:
其中H为网络隐含层节点数。
3.如权利要求1所述的单次蛛网膜下腔阻滞术中信息收集处理方法,其特征在于,重建脑电波数据衰减系数包括:
在成像中采用空间均匀分布的平行光对术中待处理部位gin进行照射,通过采集模块集成的摄像结构采集无样品遮挡的照射光测得入射光强度;对
左右两边同除以gin并取负对数,则:
采集到360度的测量数据G0后,采用精确高效的滤波反投影重建算法实现逆Radon变换即计算出衰减系数,即μt=FBP(G0)。
4.如权利要求1所述的单次蛛网膜下腔阻滞术中信息收集处理方法,其特征在于,弛豫值测量方法包括:
首先,施加IR射频脉冲至预设数量的层,使得所述预设数量的层内的磁化矢量翻转180度,所述层是指对大脑进行磁共振扫描形成的成像层;
然后,对所述预设数量的层内的磁化矢量恢复曲线采用交错采集模式采集T1图像,所述采用交错采集模式采集T1图像过程中采用了并行成像和部分傅里叶变换加速技术;
对所述预设数量的层内的磁化矢量恢复曲线采用交错采集模式采集T1图像的步骤,包括:
在第一个重复时间内,施加预设的小角度射频脉冲,通过梯度场选中预设数量的层中的第一层,并将所述第一层内的磁化矢量翻转所述小角度;
通过梯度线圈在相位方向发射的梯度场对所述翻转小角度的磁化矢量进行相位编码;
通过所述梯度线圈在读出方向发射的梯度场对所述磁化矢量进行频率编码;
通过模拟数字转换器进行数据采集得到所述第一层的T1图像;
重复执行上述步骤,直至采集完成所述预设数量的层中的每一层的T1图像。
5.如权利要求1所述的单次蛛网膜下腔阻滞术中信息收集处理方法,其特征在于,数据映射模方法包括:
(1)获得大脑的实时图像;
(2)获得脑电波数据图像映射,包括大脑的解剖模型和与大脑有关的NBS功能数据;
(3)从大脑的实时图像中识别大脑的至少一个解剖标志;
(4)识别大脑的解剖模型上的所识别出的解剖标志中的至少一个;
(5)修改功能映射,使得模型中的所识别的至少一个解剖标志在尺寸和方位上对应于大脑的实时图像中的相应的至少一个解剖标志;
(6)根据相应的对准的解剖标志将至少所述NBS功能数据以数字方式重叠在所述大脑的实时图像上。
6.一种实现权利要求1~5任意一项所述单次蛛网膜下腔阻滞术中信息收集处理方法的计算机程序
7.一种实现权利要求1~5任意一项所述单次蛛网膜下腔阻滞术中信息收集处理方法的计算机。
8.一种计算机可读存储介质,包括指令,当其在计算机上运行时,使得计算机执行如权利要求1-5任意一项所述的单次蛛网膜下腔阻滞术中信息收集处理方法。
9.一种实现权利要求1所述单次蛛网膜下腔阻滞术中信息收集处理方法的单次蛛网膜下腔阻滞术中信息收集处理系统,其特征在于,所述单次蛛网膜下腔阻滞术中信息收集处理系统包括:
图像采集模块,与主控模块连接,用于采集单次蛛网膜图像数据信息;
脑电波采集模块,与主控模块连接,用于采集单次蛛网膜下腔阻滞术中脑电波数据信息;
主控模块,与图像采集模块、脑电波采集模块、弛豫值测量模块、数据映射模块、数据存储模块、显示模块连接,用于控制各个模块正常工作;
弛豫值测量模块,与主控模块连接,用于测量大脑纵向弛豫值数据信息;
数据映射模块,与主控模块连接,用于将用户脑电波数据图像准确地映射到解剖模型;
数据存储模块,与主控模块连接,用于对采集的数据信息进行存储;
显示模块,与主控模块连接,用于通过显示器显示采集的图像及脑电波数据信息。
10.一种搭载权利要求1所述单次蛛网膜下腔阻滞术中信息收集处理系统的单次蛛网膜下腔阻滞术中信息收集处理设备。

说明书全文

单次蛛网膜下腔阻滞术中信息收集处理系统及其处理方法

技术领域

[0001] 本发明属于医疗技术领域,尤其涉及一种单次蛛网膜下腔阻滞术中信息收集处理系统及其处理方法。

背景技术

[0002] 目前,业内常用的现有技术是这样的:
[0003] 脑和脊髓外面有三层包膜:硬脑(脊)膜、蛛网膜、软脑(脊)膜是对脑脊髓组织起保护作用的。蛛网膜是一层半透明的膜,处于硬脑(脊)膜和软脑 (脊)膜之间。蛛网膜的外面是硬膜下腔,里面是蛛网膜下腔,这两个腔都是临床上很重要的腔隙。蛛网膜跨越脑,被覆于脑的表面,与软脑膜之间有较大的间隙,称为网膜下腔,腔内充满脑脊液。在一定部位,蛛网膜下腔扩展并加深,成为蛛网膜下池。其中最大的是小脑延髓池,它通过正中孔和前侧孔与第四脑室相通:桥池位于脑桥腹侧:脚间池位于脚间凹;交叉池位于视交叉前方。然而,T1值是磁共振成像(MRI)过程中纵向磁化矢量弛豫时间常数,是脑电波数据的固有属性。传统的T1值定量估计非常耗时,T1图像采集过程消耗大量时间,极易引起用户不适且采集图像易产生运动伪影;同时现有外科医生通常使用光学显微镜和/或专用成像装置来察看大脑;长期以来都非常需要以高精度和及时的方式将局部脑功能显示给外科医生,精确的方法需要在操作开始时花大量的时间,这给用户增加了险。反过来,在操作期间需要较少时间来设定的方法不够准确,也会带来相应的风险。
[0004] 综上所述,现有技术存在的问题是:
[0005] 传统的T1值定量估计非常耗时,T1图像采集过程消耗大量时间,极易引起用户不适且采集图像易产生运动伪影;同时现有外科医生通常使用光学显微镜和/或专用成像装置来察看大脑;长期以来都非常需要以高精度和及时的方式将局部脑功能显示给外科医生,精确的方法需要在操作开始时花大量的时间,这给用户增加了风险。反过来,在操作期间需要较少时间来设定的方法不够准确,也会带来相应的风险。
[0006] 现有技术的图像采集获得数据失真度大。
[0007] 现有技术中,成像由于生物样品未经透明化处理,因此存在光子散射问题。常规的成像算法只考虑生物组织对光子的吸收特性,也只能重建光子的吸收系数不能重建光子的散射系数。当进行成像时,散射的影响不可忽略,并且会和光子的吸收特性混合在一起,从而导致传统的成像空间分辨率的降低以及重建结果的不准确。

发明内容

[0008] 针对现有技术存在的问题,本发明提供了一种单次蛛网膜下腔阻滞术中信息收集处理系统及其处理方法。
[0009] 本发明是这样实现的,一种单次蛛网膜下腔阻滞术中信息收集处理方法,包括:
[0010] 采集单次蛛网膜图像数据信息;结合采集单次蛛网膜图像影响因素构建神经网络模型BP,产生神经网络模型BP的初始权值;
[0011] 对神经网络模型BP的权值和阈值进行动态改进,获得动态神经网络模型 DBP,产生动态神经网络模型DBP的权值和阈值;具体包括:调整神经网络模型BP隐含层与输出层之间的权值wkj;
[0012] 调整wkj的目的是希望输出节点j的新输出 比当前输出opj更接近目标值 tpj,定义:
[0013]
[0014] 其中α代表接近度,在每个训练周期保持不变,并随隐含层节点数H的调整而变小,不考虑阈值,则有:
[0015]
[0016] 其中wkj和w*kj分别为更新前后的权值,ypk为隐含层输出,Δwkj为wkj的改变量;
[0017] 根据公式 得到Δwkj的求解方程:
[0018]
[0019] 其中,
[0020]
[0021] 根据最小平方和误差原则求解方程 得到Δwkj的近似解:
[0022]
[0023] 对每一个连接到输出节点j的隐含层节点k,计算k与j之间的权值变化Δ wkj,更新权值并计算平方和误差E,然后在k∈[1,H]区间选择一个最优的k,使得E最小;
[0024] 调整神经网络模型BP输入层与隐含层之间的权值vik;
[0025] 调整vik的目的是一旦神经网络算法陷入局部极小点,修改权值能够跳出该极小点,判断神经网络算法陷入局部极小点的条件是误差E的变化率ΔE=0,且 E>0;
[0026] 不考虑阈值,隐含层节点k的权值的改变通过以下方程求解:
[0027]
[0028] 其中δpj=f-1(ypk+Δypk)-f-1(ypk),M为自然数,则隐含层输出ypk求解公式为:
[0029]
[0030] 其中Δypk为ypk的改变量,则有:
[0031]
[0032] 根据最小平方和误差原则求解公式 构建的矩阵方程,可以算出:
[0033]
[0034] 计算隐含层与输出层之间权值的动态平均改变量
[0035]
[0036] 计算输入层与隐含层之间权值的动态平均改变量
[0037]
[0038] 公式中M取10~20之间的自然数,获得神经网络模型BP的动态平均权值,根据神经网络模型BP的动态平均权值获得动态神经网络模型DBP;
[0039] 运用自适应免疫遗传AIGA算法优化动态神经网络模型DBP,获得预测模型AIGA-DBP,根据预测模型AIGA-DBP计算采集单次蛛网膜图像预测真实数据;
[0040] 采集单次蛛网膜下腔阻滞术中脑电波数据信息;重建脑电波数据衰减系数;
[0041] 采集模偏移一定度放置时采集投影数据;重建脑电波数据散射系数;计算脑电波的吸收系数,获得准确脑电波数据图像信息;
[0042] 重建脑电波数据散射系数包括:
[0043] 表示Γ-上入射的光子到位置 且方向为 时所传播的距离,则:
[0044]
[0045] 其中 为 处的辐射度,表示单位立体角内、单位时间内、通过垂直于单位方向2
矢量 的单位面积上的平均功率通量密度,量纲是W/(m .Sr);K0为引入的光子弹道传播算子, 表示脑电波数据在 处的总的衰减系数, 表示脑电波数据的吸
收系数, 表示散射系数;
[0046] 再:
[0047]
[0048] 其中K为引入的光子散射传播算子, 为归一化的散射相位函数,表示光子从方向v′散射到 方向的概率,满足 dΩ′表示单位方向矢量v′对应的立体角微元;定义m0=K0gin, 则有nn+1=Kmn(n≥0),从而在 处的总的辐射率为:
[0049]
[0050] 其中mn表示经n次散射到达 的辐射度分量;当光子在扩散区域传播时,K 的谱半径ρ(K)值接近1,当光子在弱散射区域传播时,ρ(K)>>1,在此种情况下,当n→∞时,快速收敛;
[0051] 然后,在输出边界Γ+上描述探测器接收到的数据总量gout,即 从而:
[0052]
[0053] 其中A是描述光子传输的矩阵,A0、A1和A2分别描述弹道传输、一次散射传输和多次散射传输部分,定义g0=A0gin,g1=A1gin分别表示测量数值中的弹道传输分量和一次散射分量,则知:
[0054]
[0055]
[0056] 入射光方向为 经一次散射后其方向为 则上式中关于 的积分仅在一特定角度上有值,取 系数k的取值由相位函数 确定,同时
定义 分别表示光子发生散
射后和散射前的衰减量,则有:
[0057]
[0058] 通过弛豫值测量方法测量大脑纵向弛豫值数据信息;
[0059] 将用户脑电波数据图像准确地映射到解剖模型;
[0060] 对采集的数据信息进行存储;并显示采集的图像及脑电波数据信息。
[0061] 进一步,产生神经网络模型BP初始权值的方法为:
[0062] 将权值初始化为[a,b]之间的随机数,其中a,b为满足以下方程的整数:
[0063]
[0064] 其中H为网络隐含层节点数。
[0065] 进一步,重建脑电波数据衰减系数包括:
[0066] 在成像中采用空间均匀分布的平行光对术中待处理部位gin进行照射,通过采集模块集成的摄像结构采集无样品遮挡的照射光测得入射光强度;对左右两边同除以gin并取负对数,则:
[0067]
[0068] 采集到360度的测量数据G0后,采用精确高效的滤波反投影重建算法实现逆Radon变换即计算出衰减系数,即μt=FBP(G0)。
[0069] 进一步,弛豫值测量方法包括:
[0070] 首先,施加IR射频脉冲至预设数量的层,使得所述预设数量的层内的磁化矢量翻转180度,所述层是指对大脑进行磁共振扫描形成的成像层;
[0071] 然后,对所述预设数量的层内的磁化矢量恢复曲线采用交错采集模式采集 T1图像,所述采用交错采集模式采集T1图像过程中采用了并行成像和部分傅里叶变换加速技术;
[0072] 对所述预设数量的层内的磁化矢量恢复曲线采用交错采集模式采集T1图像的步骤,包括:
[0073] 在第一个重复时间内,施加预设的小角度射频脉冲,通过梯度场选中预设数量的层中的第一层,并将所述第一层内的磁化矢量翻转所述小角度;
[0074] 通过梯度线圈在相位方向发射的梯度场对所述翻转小角度的磁化矢量进行相位编码;
[0075] 通过所述梯度线圈在读出方向发射的梯度场对所述磁化矢量进行频率编码;
[0076] 通过模拟数字转换器进行数据采集得到所述第一层的T1图像;
[0077] 重复执行上述步骤,直至采集完成所述预设数量的层中的每一层的T1图像。
[0078] 进一步,数据映射模方法包括:
[0079] (1)获得大脑的实时图像;
[0080] (2)获得脑电波数据图像映射,包括大脑的解剖模型和与大脑有关的NBS 功能数据;
[0081] (3)从大脑的实时图像中识别大脑的至少一个解剖标志;
[0082] (4)识别大脑的解剖模型上的所识别出的解剖标志中的至少一个;
[0083] (5)修改功能映射,使得模型中的所识别的至少一个解剖标志在尺寸和方位上对应于大脑的实时图像中的相应的至少一个解剖标志;
[0084] (6)根据相应的对准的解剖标志将至少所述NBS功能数据以数字方式重叠在所述大脑的实时图像上。
[0085] 本发明另一目的在于提供一种实现所述单次蛛网膜下腔阻滞术中信息收集处理的计算机程序
[0086] 本发明另一目的在于提供一种实现所述单次蛛网膜下腔阻滞术中信息收集处理方法的计算机。
[0087] 本发明另一目的在于提供一种计算机可读存储介质,包括指令,当其在计算机上运行时,使得计算机执行所述的单次蛛网膜下腔阻滞术中信息收集处理方法。
[0088] 本发明另一目的在于提供一种实现所述单次蛛网膜下腔阻滞术中信息收集处理方法的单次蛛网膜下腔阻滞术中信息收集处理系统,所述单次蛛网膜下腔阻滞术中信息收集处理系统包括:
[0089] 图像采集模块,与主控模块连接,用于采集单次蛛网膜图像数据信息;
[0090] 脑电波采集模块,与主控模块连接,用于采集单次蛛网膜下腔阻滞术中脑电波数据信息;
[0091] 主控模块,与图像采集模块、脑电波采集模块、弛豫值测量模块、数据映射模块、数据存储模块、显示模块连接,用于控制各个模块正常工作;
[0092] 弛豫值测量模块,与主控模块连接,用于测量大脑纵向弛豫值数据信息;
[0093] 数据映射模块,与主控模块连接,用于将用户脑电波数据图像准确地映射到解剖模型;
[0094] 数据存储模块,与主控模块连接,用于对采集的数据信息进行存储;
[0095] 显示模块,与主控模块连接,用于通过显示器显示采集的图像及脑电波数据信息。
[0096] 本发明另一目的在于提供一种搭载所述单次蛛网膜下腔阻滞术中信息收集处理系统的单次蛛网膜下腔阻滞术中信息收集处理设备。
[0097] 本发明的优点及积极效果为:
[0098] 本发明通过图像数据处理模块能够将识别为有异常信号的蛛网膜下腔位置显示出来,供医务人员参考,为医务人员的诊断提供精确的资料,以辅助对蛛网膜下腔出血进行准确判定,能够降低蛛网膜下腔出血的误诊/本发明提供的弛豫值测量模块由于交错采集模式采集T1图像过程中采用了并行成像和部分傅里叶变换加速技术,只需在一次IR激励内即可完成一幅图谱所需的所有数据的采集,大大缩短了T1图像采集时间,因而不易引起用户不适,减少采集图像中的运动伪影;同时通过数据映射模块将用户大脑的MRI或CT用作针对映射的解剖模型。一旦MRI数据已被汇编,那么通过知道用户脑部的确切位置和方位以及每个TMS脉冲,就可以将用户脑电波数据图像准确地映射到解剖模型;降低风险。
[0099] 本发明通过采集单次蛛网膜图像数据信息;结合采集单次蛛网膜图像影响因素构建神经网络模型BP,产生神经网络模型BP的初始权值;
[0100] 对神经网络模型BP的权值和阈值进行动态改进,获得动态神经网络模型 DBP,产生动态神经网络模型DBP的权值和阈值;具体包括:调整神经网络模型BP隐含层与输出层之间的权值wkj;
[0101] 调整wkj的目的是希望输出节点j的新输出 比当前输出opj更接近目标值 tpj,定义:
[0102]
[0103] 其中α代表接近度,在每个训练周期保持不变,并随隐含层节点数H的调整而变小,不考虑阈值,则有:
[0104] 本发明图像采集获得数据失真度小、清晰度高。
[0105] 本发明针对成像存在弱散射的问题,结合OPT成像数据采集的自身特点,构建相应的数学模型,并通过额外测量一组与入射光保持一定倾斜角度的数据的方法,对弹道传输分量和一次散射传输分量进行分离,进而实现生物组织吸收系数和散射系数的三维重建,从而既可以有效解决OPT成像存在的散射问题;提高了OPT成像质量,同时更丰富了OPT技术提供的信息量,使得OPT技术可以从吸收系数和散射系数两个角度描述生物的组织结构特性。附图说明
[0106] 图1是本发明实施提供的单次蛛网膜下腔阻滞术中信息收集处理方法流程图
[0107] 图2是本发明实施提供的单次蛛网膜下腔阻滞术中信息收集处理系统结构框图
[0108] 图中:1、图像采集模块;2、参数配置模块;3、中央控制模块;4、图像数据处理模块;5、麻醉模块;6、定位模块;7、显示模块。

具体实施方式

[0109] 为了使本发明的目的、技术方案及优点更加清楚明白,以下结合实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
[0110] 下面结合附图及具体实施例对本发明的应用原理作进一步描述。
[0111] 如图1所示,本发明提供的一种单次蛛网膜下腔阻滞术中信息收集处理系统及其处理方法包括以下步骤:
[0112] S101,通过图像采集模块采集单次蛛网膜图像数据信息;
[0113] S102,通过脑电波采集模块采集单次蛛网膜下腔阻滞术中脑电波数据信息;
[0114] S103,主控模块通过弛豫值测量模块测量大脑纵向弛豫值数据信息;
[0115] S104,通过数据映射模块将用户脑电波数据图像准确地映射到解剖模型;
[0116] S105,通过数据存储模块对采集的数据信息进行存储;并通过显示模块显示采集的图像及脑电波数据信息。
[0117] 如图2所示,本发明提供的单次蛛网膜下腔阻滞术中信息收集处理系统包括:图像采集模块1、脑电波采集模块2、主控模块3、弛豫值测量模块4、数据映射模块5、数据存储模块6、显示模块7。
[0118] 图像采集模块1,与主控模块3连接,用于采集单次蛛网膜图像数据信息;
[0119] 脑电波采集模块2,与主控模块3连接,用于采集单次蛛网膜下腔阻滞术中脑电波数据信息;
[0120] 主控模块3,与图像采集模块1、脑电波采集模块2、弛豫值测量模块4、数据映射模块5、数据存储模块6、显示模块7连接,用于控制各个模块正常工作;
[0121] 弛豫值测量模块4,与主控模块3连接,用于测量大脑纵向弛豫值数据信息;
[0122] 数据映射模块5,与主控模块3连接,用于将用户脑电波数据图像准确地映射到解剖模型;
[0123] 数据存储模块6,与主控模块3连接,用于对采集的数据信息进行存储;
[0124] 显示模块7,与主控模块3连接,用于通过显示器显示采集的图像及脑电波数据信息。
[0125] 本发明提供的弛豫值测量模块4测量方法如下:
[0126] 首先,施加IR射频脉冲至预设数量的层,使得所述预设数量的层内的磁化矢量翻转180度,所述层是指对大脑进行磁共振扫描形成的成像层;
[0127] 然后,对所述预设数量的层内的磁化矢量恢复曲线采用交错采集模式采集 T1图像,所述采用交错采集模式采集T1图像过程中采用了并行成像和部分傅里叶变换加速技术。
[0128] 本发明提供的对所述预设数量的层内的磁化矢量恢复曲线采用交错采集模式采集T1图像的步骤,包括:
[0129] 在第一个重复时间内,施加预设的小角度射频脉冲,通过梯度场选中预设数量的层中的第一层,并将所述第一层内的磁化矢量翻转所述小角度;
[0130] 通过梯度线圈在相位方向发射的梯度场对所述翻转小角度的磁化矢量进行相位编码;
[0131] 通过所述梯度线圈在读出方向发射的梯度场对所述磁化矢量进行频率编码;
[0132] 通过模拟数字转换器进行数据采集得到所述第一层的T1图像;
[0133] 重复执行上述步骤,直至采集完成所述预设数量的层中的每一层的T1图像。
[0134] 本发明提供的数据映射模块5映射方法如下:
[0135] (1)获得大脑的实时图像;
[0136] (2)获得脑电波数据图像映射,包括大脑的解剖模型和与大脑有关的NBS 功能数据;
[0137] (3)从大脑的实时图像中识别大脑的至少一个解剖标志;
[0138] (4)识别大脑的解剖模型上的所识别出的解剖标志中的至少一个;
[0139] (5)修改功能映射,使得模型中的所识别的至少一个解剖标志在尺寸和方位上对应于大脑的实时图像中的相应的至少一个解剖标志;
[0140] (6)根据相应的对准的解剖标志将至少所述NBS功能数据以数字方式重叠在所述大脑的实时图像上。
[0141] 下面结合具体分析对本发明作进一步描述。
[0142] 本发明实施例提供的单次蛛网膜下腔阻滞术中信息收集处理方法,包括:
[0143] 采集单次蛛网膜图像数据信息;结合采集单次蛛网膜图像影响因素构建神经网络模型BP,产生神经网络模型BP的初始权值;
[0144] 对神经网络模型BP的权值和阈值进行动态改进,获得动态神经网络模型 DBP,产生动态神经网络模型DBP的权值和阈值;具体包括:调整神经网络模型BP隐含层与输出层之间的权值wkj;
[0145] 调整wkj的目的是希望输出节点j的新输 比当前输出opj更接近目标值tpj,定义:
[0146]
[0147] 其中α代表接近度,在每个训练周期保持不变,并随隐含层节点数H的调整而变小,不考虑阈值,则有:
[0148]
[0149] 其中wkj和w*kj分别为更新前后的权值,ypk为隐含层输出,Δwkj为wkj的改变量;
[0150] 根据公式 得到Δwkj的求解方程:
[0151]
[0152] 其中,
[0153]
[0154] 根据最小平方和误差原则求解方程 得到Δwkj的近似解:
[0155]
[0156] 对每一个连接到输出节点j的隐含层节点k,计算k与j之间的权值变化Δ wkj,更新权值并计算平方和误差E,然后在k∈[1,H]区间选择一个最优的k,使得E最小;
[0157] 调整神经网络模型BP输入层与隐含层之间的权值vik;
[0158] 调整vik的目的是一旦神经网络算法陷入局部极小点,修改权值能够跳出该极小点,判断神经网络算法陷入局部极小点的条件是误差E的变化率ΔE=0,且 E>0;
[0159] 不考虑阈值,隐含层节点k的权值的改变通过以下方程求解:
[0160]
[0161] 其中δpj=f-1(ypk+Δypk)-f-1(ypk),M为自然数,则隐含层输出ypk求解公式为:
[0162]
[0163] 其中Δypk为ypk的改变量,则有:
[0164]
[0165] 根据最小平方和误差原则求解公式 构建的矩阵方程,可以算出:
[0166]
[0167] 计算隐含层与输出层之间权值的动态平均改变量
[0168]
[0169] 计算输入层与隐含层之间权值的动态平均改变量
[0170]
[0171] 公式中M取10~20之间的自然数,获得神经网络模型BP的动态平均权值,根据神经网络模型BP的动态平均权值获得动态神经网络模型DBP;
[0172] 运用自适应免疫遗传AIGA算法优化动态神经网络模型DBP,获得预测模型AIGA-DBP,根据预测模型AIGA-DBP计算采集单次蛛网膜图像预测真实数据;
[0173] 采集单次蛛网膜下腔阻滞术中脑电波数据信息;重建脑电波数据衰减系数;
[0174] 采集模块偏移一定角度放置时采集投影数据;重建脑电波数据散射系数;计算脑电波的吸收系数,获得准确脑电波数据图像信息;
[0175] 重建脑电波数据散射系数包括:
[0176] 表示Γ-上入射的光子到位置 且方向为 时所传播的距离,则:
[0177]
[0178] 其中 为 处的辐射度,表示单位立体角内、单位时间内、通过垂直于单位方向矢量 的单位面积上的平均功率通量密度,量纲是W/(m2.Sr);K0为引入的光子弹道传播算子, 表示脑电波数据在 处的总的衰减系数, 表示脑电波数据的吸收系数, 表示散射系数;
[0179] 再:
[0180]
[0181] 其中K为引入的光子散射传播算子, 为归一化的散射相位函数,表示光子从方向v′散射到 方向的概率,满足 dΩ′表示单位方向矢量v′对应的立体角微元;定义m0=K0gin, 则有nn+1=Kmn(n≥0),从而在 处的总的辐射率为:
[0182]
[0183] 其中mn表示经n次散射到达 的辐射度分量;当光子在扩散区域传播时,K 的谱半径ρ(K)值接近1,当光子在弱散射区域传播时,ρ(K)>>1,在此种情况下,当n→∞时,快速收敛;
[0184] 然后,在输出边界Γ+上描述探测器接收到的数据总量gout,即 从而:
[0185]
[0186] 其中A是描述光子传输的矩阵,A0、A1和A2分别描述弹道传输、一次散射传输和多次散射传输部分,定义g0=A0gin,g1=A1gin分别表示测量数值中的弹道传输分量和一次散射分量,则知:
[0187]
[0188]
[0189] 入射光方向为 经一次散射后其方向为 则上式中关于 的积分仅在一特定角度上有值,取 系数k的取值由相位函数 确定,同时
定义 分别表示光子发生散
射后和散射前的衰减量,则有:
[0190]
[0191] 通过弛豫值测量方法测量大脑纵向弛豫值数据信息;
[0192] 将用户脑电波数据图像准确地映射到解剖模型;
[0193] 对采集的数据信息进行存储;并显示采集的图像及脑电波数据信息。
[0194] 进一步,产生神经网络模型BP初始权值的方法为:
[0195] 将权值初始化为[a,b]之间的随机数,其中a,b为满足以下方程的整数:
[0196]
[0197] 其中H为网络隐含层节点数。
[0198] 重建脑电波数据衰减系数包括:
[0199] 在成像中采用空间均匀分布的平行光对术中待处理部位gin进行照射,通过采集模块集成的摄像结构采集无样品遮挡的照射光测得入射光强度;对左右两边同除以gin并取负对数,则:
[0200]
[0201] 采集到360度的测量数据G0后,采用精确高效的滤波反投影重建算法实现逆Radon变换即计算出衰减系数,即μt=FBP(G0)。
[0202] 重复执行上述步骤,直至采集完成所述预设数量的层中的每一层的T1图像。
[0203] 在上述实施例中,可以全部或部分地通过软件硬件固件或者其任意组合来实现。当使用全部或部分地以计算机程序产品的形式实现,所述计算机程序产品包括一个或多个计算机指令。在计算机上加载或执行所述计算机程序指令时,全部或部分地产生按照本发明实施例所述的流程或功能。所述计算机可以是通用计算机、专用计算机、计算机网络、或者其他可编程装置。所述计算机指令可以存储在计算机可读存储介质中,或者从一个计算机可读存储介质向另一个计算机可读存储介质传输,例如,所述计算机指令可以从一个网站站点、计算机、服务器数据中心通过有线(例如同轴电缆、光纤、数字用户线(DSL)或无线(例如红外、无线、微波等)方式向另一个网站站点、计算机、服务器或数据中心进行传输)。所述计算机可读取存储介质可以是计算机能够存取的任何可用介质或者是包含一个或多个可用介质集成的服务器、数据中心等数据存储设备。所述可用介质可以是磁性介质,(例如,软盘硬盘、磁带)、光介质(例如,DVD)、或者半导体介质(例如固态硬盘Solid State Disk(SSD))等。
[0204] 以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈