一种基于稀疏分解与景深估计的神经元提取方法及装置制造方法及图纸

本发明专利技术公开一种基于稀疏分解与景深估计的神经元提取方法及装置，方法包括：输入原始图像视频数据；利用鲁棒性主成分分析方法分解原始图像视频数据，以得到神经元信号、背景和噪声；计算神经元信号数据的亮暗通道先验，利用两种先验估计神经元信号数据景深对应的传输率；将神经元信号数据以像素为单位点除景深对应的传输率，构建并求解约束性非负矩阵分解的神经元提取框架；输出神经元的提取信息，包括：空间位置及形状大小、动作电位和时序变化。本发明专利技术能够高效地分解背景及噪声等干扰信息，快速估计神经元信号景深以去除数据散射，准确有效地提取神经元的空间位置及形状大小、动作电位和时序变化等特征信息。电位和时序变化等特征信息。电位和时序变化等特征信息。

【技术实现步骤摘要】
一种基于稀疏分解与景深估计的神经元提取方法及装置


[0001]本专利技术涉及神经元提取
，特别是指一种基于稀疏分解与景深估计的神经元提取方法及装置。

技术介绍

[0002]神经元提取是从光学显微仪器所获得的荧光钙成像图像视频数据中，提取神经元信号的空间位置及形状大小、动作电位及时序变化等特征信息。神经元提取作为关键必要环节，用于分析与阐释大尺度或介观尺度神经元细胞群所产生的生命活动及变化机理，已成为脑科学与生命科学等领域的研究热点之一。
[0003]针对钙成像中的神经元提取，现有的神经元提取方法主要有约束性非负矩阵分解和扩展约束性非负矩阵分解等方法。然而，上述方法迭代求解与重构复杂背景信息，需花费大量的运算时间；噪声直接影响神经元的提取性能。同时，这些方法未考虑及解决钙成像的散射问题；在数据矩阵运算中未考虑利用景深信息。因此如何构建一种去背景、去噪声、去散射和景深利用的高效神经元提取方法是本领域技术人员当前亟需解决的问题。

技术实现思路

[0004]本专利技术实施例提供了一种基于稀疏分解与景深估计的神经元提取方法及装置。所述技术方案如下：
[0005]一方面，提供了一种基于稀疏分解与景深估计的神经元提取方法，该方法由电子设备实现，该方法包括：
[0006]S1、输入原始图像视频数据；
[0007]S2、利用鲁棒性主成分分析方法分解输入的所述原始图像视频数据，以得到神经元信号、背景和噪声；
[0008]S3、计算神经元信号数据的亮暗通道先验，利用两种先验估计神经元信号数据景深对应的传输率；
[0009]S4、将所述神经元信号数据以像素为单位点除景深对应的传输率，构建约束性非负矩阵分解的神经元提取框架；
[0010]S5、求解约束性非负矩阵分解的神经元提取框架，输出神经元的提取信息，所述提取信息包括：空间位置及形状大小、动作电位和时序变化。
[0011]可选地，所述S2，具体包括：
[0012]对输入数据V中的神经元信号S和背景L进行稀疏分解，对应的目标函数表示为：
[0013]（1）
[0014]其中，和分别为核范数与L0范式，为平衡神经元信号稀疏性与背景低秩性的权重参数；
[0015]利用增广拉格朗日乘子法交替迭代求解S与L的目标函数：
[0016]固定神经元信号S，求解背景L：
[0017]（2）
[0018]其中，和分别为奇异值阈值操作和软阈值操作，A和B是奇异值阈值操作的正交矩阵，是奇异值阈值操作的阈值参数，k和分别为迭代次数和平衡权重参数，x为像素点坐标，Y
k
为残差及噪声；
[0019]固定当前背景L，求解神经元信号S：
[0020]（3）
[0021]其中，为平衡权重参数；
[0022]更新残差及噪声Y：。
[0023]可选地，所述S3，具体包括：
[0024]将神经元信号分解为一系列列向量，然后将这些列向量分别转变为一帧帧神经元信号图像，其中S
i
为第i个神经元信号图像对应的矩阵，N为神经元信号数据的图像帧数；
[0025]计算每帧神经元信号图像的亮通道先验S
bcp
：
[0026]（4）
[0027]计算每帧神经元信号图像的暗通道先验S
dcp
：
[0028]（5）
[0029]其中，S
i
含有r个通道，x为S的当前像素点坐标，是x为中心的邻近区域，y是邻近区域的当前像素点坐标；
[0030]在每帧神经元信号图像的亮暗通道先验基础上，估计每帧神经元信号图像的景深对应传输率M：
[0031]（6）
[0032]其中，除法运算均为像素点除操作，加入权重系数w，以避免传输率的过估计。
[0033]可选地，所述S4，具体包括：
[0034]将得到的神经元信号图像S和估计的对应景深的传输率M进行像素点除运算：
[0035]，i = 1, ..., N (7)
[0036]其中，为极小常数10
‑6，以避免分母趋近于0；
[0037]将点除运算后的全部结果表示为一系列列向量并组成一个矩阵，作为约束性非负矩阵分解的神经元提取框架的初始输入数据。
[0038]可选地，所述S5，具体包括：
[0039]S51、固定神经元的动作电位及时序变化，利用快速分层交替最小二乘算法，求解得到神经元的空间位置及形状大小；
[0040]S52、固定步骤S51得到的神经元的空间位置及形状大小，利用在线有效集算法，求解得到神经元的动作电位及时序变化；
[0041]S53、重复步骤S51
‑
S52直到交替迭代结束，输出神经元的提取信息，所述提取信息包括：空间位置及形状大小、动作电位和时序变化。
[0042]可选地，所述S51，具体包括：
[0043]固定所有神经元的当前动作电位及时序变化，求解所有神经元空间位置及形状大小的目标函数：
[0044](8)
[0045]其中，P为所有神经元的空间位置及形状大小，具有稀疏性和空间局部性，C为所有神经元的动作电位及时序变化，上标符号^为当前参数估计值；
[0046]采用快速分层交替最小二乘算法求解P的目标函数, 对C中每个神经元均进行如下运算操作：
[0047](9)
[0048]其中，, 上标符号T为转置运算，k = 1, ..., K，K为C中所有神经元数目，迭代上述运算直至迭代结束。
[0049]可选地，所述S52，具体包括：
[0050]固定步骤S51得到的所有神经元的当前空间位置及形状大小，推断所有神经元的动作电位C及时序变化S的目标函数构建为：
[0051](10)
[0052]其中，，，t
k
为神经元时间动态变化中所产生的冲击次数，具有稀疏性，G
(k)
为二阶回归系数，利用二阶回归过程对每个神经元的时间动态活动性进行建模；
[0053]采用在线有效集算法求解神经元动作电位C及时序变化T的目标函数，对每个神经元进行如下迭代运算操作：
[0054]（11）
[0055]其中，池化变量，池化长度，为G的自相关系数，更新权重，q
k
初始设置为，迭代此运算以求解得到神经元动作电位C及时序变化T。
[0056]另一方面，提供了一种基于稀疏分解与景深估计的神经元提取装置，所述装置包括：
[0057]输入模块，用于输入原始图像视频数据；
[0058]分解模块，用于利用鲁棒性主成分分析方法分解输入的所述原始图像视频数据，以得到神经元信号、背景和噪声；
[0059]计算模块，用于计算神经元信号数据的亮暗通道先验，利用两种先验估计神经元信号数据景深对应的传输率；
[0060]构建模块，用于将所述神经元信号数据以本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于稀疏分解与景深估计的神经元提取方法，其特征在于，所述方法包括：S1、输入原始图像视频数据；S2、利用鲁棒性主成分分析方法分解输入的所述原始图像视频数据，以得到神经元信号、背景和噪声；S3、计算神经元信号数据的亮暗通道先验，利用两种先验估计神经元信号数据景深对应的传输率；S4、将所述神经元信号数据以像素为单位点除景深对应的传输率，构建约束性非负矩阵分解的神经元提取框架；S5、求解约束性非负矩阵分解的神经元提取框架，输出神经元的提取信息，所述提取信息包括：空间位置及形状大小、动作电位和时序变化。2.根据权利要求1所述的方法，其特征在于，所述S2，具体包括：对输入数据V中的神经元信号S和背景L进行稀疏分解，对应的目标函数表示为：（1）其中，和分别为核范数与L0范式，为平衡神经元信号稀疏性与背景低秩性的权重参数；利用增广拉格朗日乘子法交替迭代求解S与L的目标函数：固定神经元信号S，求解背景L：（2）其中，和分别为奇异值阈值操作和软阈值操作，A和B是奇异值阈值操作的正交矩阵，是奇异值阈值操作的阈值参数，k和分别为迭代次数和平衡权重参数，x为像素点坐标，Y
k
为残差及噪声；固定当前背景L，求解神经元信号S：（3）其中，为平衡权重参数；更新残差及噪声Y：。3.根据权利要求1所述的方法，其特征在于，所述S3，具体包括：将神经元信号分解为一系列列向量，然后将这些列向量分别转变为一帧帧神经元信号图像，其中S
i
为第i个神经元信号图像对应的矩阵，N为神经元信号数据的图像帧数；
计算每帧神经元信号图像的亮通道先验S
bcp
：（4）计算每帧神经元信号图像的暗通道先验S
dcp
：（5）其中，S
i
含有r个通道，x为S的当前像素点坐标，是x为中心的邻近区域，y是邻近区域的当前像素点坐标；在每帧神经元信号图像的亮暗通道先验基础上，估计每帧神经元信号图像的景深对应传输率M：（6）其中，除法运算均为像素点除操作，加入权重系数w，以避免传输率的过估计。4.根据权利要求1所述的方法，其特征在于，所述S4，具体包括：将得到的神经元信号图像S和估计的对应景深的传输率M进行像素点除运算：，i = 1, ..., N
ꢀꢀ
(7)其中，为极小常数10
‑6，以避免分母趋近于0；将点除运算后的全部结果表示为一系列列向量并组成一个矩阵，作为约束性非负矩阵分解的神经元提取框架的初始输入数据。5.根据权利要求1所述的方法，其特征在于，所述S5，具体包括：S51、固定神经元的动作电位及时序变化，利用快速分层交替最小二乘算法，求解得到神经元的空间位置及形状大小；S52、固定步骤S51得到的神经元的空间位置及形状大小，利用在线有效集...

【专利技术属性】
技术研发人员：庄培显，李擎，李江昀，张天翔，徐银梅，洪然，姜天轶，英子瑄，钟祚栋，
申请(专利权)人：北京科技大学，
类型：发明
国别省市：