一种图像处理方法

权利要求

1.一种图像处理方法，采用视皮层假体的图像处理装置，该装置包括左右微型摄像头、图像预处理器、图像编码器，所述左右微型摄像头将各自同步实时采集的视频信息输出给所述图像预处理器, 图像预处理器将接收的两路所述视频信息进行拼接融合成一幅全视野图像，并基于视觉系统的注意机制,检测图像的显著性区域,进一步完成图像的边缘提取、增强和像素降低的预处理操作，输出给所述图像编码器, 图像编码器将接收的所述经预处理后的全视野图像做进一步的处理与编码，转化为人体视觉皮层能够接收的电信号；

其特征在于：图像处理方法包括图像预处理和图像编码两部分；

一、所述图像预处理按照下述步骤进行：

S1、图像校正

采用模板法对所述左右微型摄像头传输来的实时视频信息进行模板制定、特征点提取 和求解镜头畸变系数；

S2、图像配准

采用比值匹配法对所述左右微型摄像头传输来的实时视频图像进行快速配准；

S3、图像融合

图像融合是为了解决图像匹配结束后所述左右摄像头相邻视角拼接的交界处问题，这 里采用渐入渐出的线性融合方法实现双路视频图像的融合；即：将左右摄像头相邻视角图 像的像素值进行加权平均得到重叠区每一个像素点的值；其中像素点到图像重叠区两边缘 的距离作为融合比重；

S4、显著性区域检测

采用简单线性迭代的超像素图像分割方法与马尔科夫吸收链相结合的方法提取融合 后全视野图像的显著性区域； 即: 采用简单线性迭代的超像素图像分割方法对融合后的 全视野图像进行分割，以图像中超像素点作为节点，连接各节点对图像分割成众多区域，利 用马尔科夫吸收链方法检测出图像中显著性的区域；

S5、边缘提取

采用Canny算子检测算法对显著性区域内的目标图像进行边缘提取；即进行降噪、寻找 图像中的亮度梯度、在图像中跟踪边缘；

S6、图像增强

采用形态学方法中的腐蚀和膨胀进行图像增强；

二、所述图像编码按照下述步骤进行：

C1、借鉴稀疏编码的思想，将预处理后的图像采用图像原子库中的基函数进行表示，得 到每个原子对应的稀疏响应系数；

C2、根据给定响应系数和神经元的特征调谐参数，求解电流刺激方案；即：通过H-H模型 模拟特定神经元的特征调谐特性来确定电流刺激方案，为激活视皮层中的神经元细胞产生 拟定的响应模式提供关键方案与参数；

所述图像编码步骤C1中稀疏响应系数的求解采用如下步骤完成：C11、从自然图像中随机采集5000个2n×2n(n=4，5，6)个图像块， 对每个图像块分别进行去相关性、降维预处理操作，其中所述自然图像来自用于测试稀疏 编码模型的图像库；

C12、利用所述自然图像高阶统计特性与视皮层神经元稀疏响应之间的联系，采用拓扑 独立分量分析的方法对自然图像块组成的集合进行特征提取，得到与视皮层神经元感受野 特性相符的完备原子库；

C13、采用所述原子库，对处理后的图像进行稀疏分解，得到对应的响应系数：C14、设置初始阈值，将响应小于阈值的神经元响应置零，得到响应矩阵；

C15、根据所述响应矩阵重建图像；

C16、计算所述重建图像与原始图像的误差率；

C17、若所述误差大于10%，则返回步骤 C14降低阈值，若小于10%则结束，输出响应系数矩 阵，从而减少在视皮层激活神经元的数量。

2.根据权利要求1所述的图像处理方法，其特征在于：所述左右微型摄像头 均为CMOS图像传感器摄像头，分别对应安装在眼镜架的左右眼镜框内，单个CMOS图像传感 器摄像头的视角范围≤90度，分别通过USB接口与所述图像预处理器进行数据传输。

3.根据权利要求1所述的图像处理方法，其特征在于：所述左右微型摄像头在左右眼镜框内的安装位置与角度均可调，融合后的全视野图像的视角θ计算公式如下：式中，α为左微型摄像头的左极限视界与正前方的夹角，等同于右微型摄像头的右极限 视界与正前方的夹角，d为两个微型摄像头的中心间距，h为拍摄点m到左、右微型摄像头连 线的垂直距离；所述α由单个微型摄像头的视角范围及其与两微型摄像头所呈直线间夹角 共同决定，通过调整α和d，能够扩大总视角θ。

4.根据权利要求1所述的图像处理方法，其特征在于：所述图像预处理步骤S2中的比值匹配法采用如下步骤实现图像的快速配准：S21、读取所述左摄像头所采集图像最右侧一列像素点和每间隔5个像素点对应位置处 的像素点，共取n列像素点， n为大于1的自然数，设每列像素点个数为m，m为大于1的自然 数；

S22、计算所述每列像素点相邻行之间的比值，即得到（m-1）×n个比值，将所述比值结果 直接存放在数组P中；

S23、读取所述右摄像头采集的图像中任意每间隔5个像素点的n列像素点，取k组，k为大 于1的自然数；

S24、计算所述k组像素点每一列相邻行像素点之间的比值，将所述比值结果直接存放在 数组Q中；

S25、根据左摄像头所采集的图像中的特征集模板，即数组P，在右摄像头所采集图像的 特征集合中寻找相应的匹配，得出配准区域。

5.根据权利要求1所述的图像处理方法，其特征在于：所述图像预处理步骤S6后还需要进行像素降 低，目的是为了降低图像的像素值，保证与所述图像编码部分训练原子库所选图像块像素 值吻合，具体步骤如下：n n

S7、根据实际要求将原始图像进行分割，将原始图像平均分割成2 ×2 (n=4、5或6)个图像块；

S8、计算每个分割区域像素点的平均值，然后将该平均值赋给该分割区域中的所有点；

S9、将灰度级设定为5级，即：灰度值为0、60、120、180、255，然后，根据原图像每个像素 点对应的灰度范围为其赋新的灰度值。

6.根据权利要求1所述的图像处理方法，其特征在于：所述图像编码步骤C2中H-H模型的求解，具体步骤如下：C21、建立H-H神经元模型；

C22、采用Runge-Kutta算法和迭代方法求解H-H神经元模型，建立电刺激强度与膜电压 的关系；

C23、模拟不同电刺激方案下神经元的发放情况，包括直流电刺激、脉冲波刺激、双相不 对称矩形脉冲刺激和双相不对称梯度脉冲刺激，从功耗低、脑组织损伤小的角度分析它们的优缺点，最终确定双相不对称梯度脉冲刺激为最佳的电刺激方案。

说明书

一种图像处理方法

技术领域

[0001] 本发明涉及视皮层假体，尤其是涉及视皮层假体的图像处理装置及方法。

背景技术

[0002] 视觉假体是一种通过电刺激视觉神经系统使得盲人重建光明的装置，目前，很多 国家的科技人员在研究采用视觉假体的方法来修复盲人的视觉功能。根据视觉假体的植入 位置不同，可分为视网膜假体、视神经假体和视皮层假体。其中，视皮层假体能够绕过视觉 通路中可能发生病变的视网膜、视神经等，直接刺激于视皮层，适用于更多的盲人，因此，相 对上述另外两种视觉假体具有更加广泛的应用。

[0003] 视皮层假体的设计原理是：首先，采用图像采集模块采集图像信息；然后，图像处 理器对采集到的图像信息进行处理和编码，并输出刺激代码；最后，通过电刺激装置刺激相 应的视觉皮层，使盲人患者产生人工视觉。不同类型的视觉假体因为刺激位置不同，因此所 采用的编码策略和电刺激方案存在特异性。目前现有的视觉假体系统多侧重于在视网膜和 视神经处实施电刺激，例如，欧洲发明专利（公开号EP2155327，公开日为2010年02月24日） 和美国发明专利（公开号为US9180296，公开日为2015年11月10日）均自述是为植入式视觉 假体提供刺激输入的系统，给出的电刺激方案仅适用于视网膜假体；美国发明专利（公开号 WO2014121288，公开日期为2014年08月07日）自述为一种视皮层假体，但仅仅提供了包含针 对大脑皮层矩形沟的植入式器件（柔性电极阵列）和植入方案，未涉及图像编码的方法；中 国发明专利（公开号CN101239008A，公开日2008年8月13日）公开了一种视觉假体的图像处 理装置与方法，自述其装置和方法适用于视网膜假体、视神经假体和视皮层假体，然而对于 对图像的编码运算中未考虑视皮层神经元的编码特性，且未提及图像编码的结果与电极阵 列刺激模块如何对接。

[0004] 理论上，视觉假体恢复人工视觉的成像效果与植入微电极阵列的密集程度与针数 成正比；然而，微电极阵列的密集程度与数量又直接影响到微电极阵列在脑内的植入寿命 和图像处理的复杂度。如何采用可能少的电极数目让患者获得尽可能多样的视觉信息，是 视皮层假体研究的关键问题。2016年6月发表的硕士论文《基于初级视皮层的人工视觉假体 总体设计与前端实现》提供了对应的解决方案，借鉴稀疏编码的思想利用稀疏的神经元响 应表征图像中目标的边界信息，然而未利用视觉的注意机制实现对图像显著性区域的检 测，文中给出的图像处理方案仅适用于当目标占据视野较大范围的情况，若图像中的显著 性目标仅占据视野的小部分区域，对整幅图像进行像素降低后，图像中关于显著性目标的 信息会大大降低。

发明内容

[0005] 本发明目的在于提供一种视皮层假体的图像处理装置及方法。

[0006] 为实现上述目的，本发明采取下述技术方案：

[0007] 本发明所述图像处理方法，采用视皮层假体的图像处理装置，该装置包括左右微型摄像头、图像预处理器、图像编码器，所述左右微型摄像头将各自同步实时采集的视频信息输出给所述图像预处理器, 图像预处理器将接收的两路所述视频信息进行拼接融合成一幅全视野图像，并基于视觉系统的注意机制,检测图像的显著性区域,进一步完成图像的边缘提取、增强和像素降低的预处理操作，输出给所述图像编码器, 图像编码器将接收的所述经预处理后的全视野图像做进一步的处理与编码，转化为人体视觉皮层能够接收的电信号；该图像处理方法包括图像预处理和图像编码两部分；

[0008] 一、所述图像预处理按照下述步骤进行：

[0009] S1、图像校正

[0010] 采用模板法对所述左右微型摄像头传输来的实时视频信息进行模板制定、特征点提取 和求解镜头畸变系数；

[0011] S2、图像配准

[0012] 采用比值匹配法对所述左右微型摄像头传输来的实时视频图像进行快速配准；

[0013] S3、图像融合

[0014] 图像融合是为了解决图像匹配结束后所述左右摄像头相邻视角拼接的交界处问题，这 里采用渐入渐出的线性融合方法实现双路视频图像的融合；即：将左右摄像头相邻视角图 像的像素值进行加权平均得到重叠区每一个像素点的值；其中像素点到图像重叠区两边缘 的距离作为融合比重；

[0015] S4、显著性区域检测

[0016] 采用简单线性迭代的超像素图像分割方法与马尔科夫吸收链相结合的方法提取融合 后全视野图像的显著性区域； 即: 采用简单线性迭代的超像素图像分割方法对融合后的 全视野图像进行分割，以图像中超像素点作为节点，连接各节点对图像分割成众多区域，利 用马尔科夫吸收链方法检测出图像中显著性的区域；

[0017] S5、边缘提取

[0018] 采用Canny算子检测算法对显著性区域内的目标图像进行边缘提取；即进行降噪、寻找 图像中的亮度梯度、在图像中跟踪边缘；

[0019] S6、图像增强

[0020] 采用形态学方法中的腐蚀和膨胀进行图像增强；

[0021] 二、所述图像编码按照下述步骤进行：

[0022] C1、借鉴稀疏编码的思想，将预处理后的图像采用图像原子库中的基函数进行表示，得 到每个原子对应的稀疏响应系数；

[0023] C2、根据给定响应系数和神经元的特征调谐参数，求解电流刺激方案；即：通过H-H模型 模拟特定神经元的特征调谐特性来确定电流刺激方案，为激活视皮层中的神经元细胞产生 拟定的响应模式提供关键方案与参数；

[0024] 所述图像编码步骤C1中稀疏响应系数的求解采用如下步骤完成：

[0025] C11、从自然图像中随机采集5000个2n×2n (n=4，5，6)个图像块， 对每个图像块分别进行去相关性、降维预处理操作，其中所述自然图像来自用于测试稀疏 编码模型的图像库；

[0026] C12、利用所述自然图像高阶统计特性与视皮层神经元稀疏响应之间的联系，采用拓扑 独立分量分析的方法对自然图像块组成的集合进行特征提取，得到与视皮层神经元感受野 特性相符的完备原子库；

[0027] C13、采用所述原子库，对处理后的图像进行稀疏分解，得到对应的响应系数：

[0028] C14、设置初始阈值，将响应小于阈值的神经元响应置零，得到响应矩阵；

[0029] C15、根据所述响应矩阵重建图像；

[0030] C16、计算所述重建图像与原始图像的误差率；

[0031] C17、若所述误差大于10%，则返回步骤C14降低阈值，若小于10%则结束，输出响应系数矩 阵，从而减少在视皮层激活神经元的数量。

[0032] 所述左右微型摄像头 均为CMOS图像传感器摄像头，分别对应安装在眼镜架的左右眼镜框内，单个CMOS图像传感 器摄像头的视角范围≤90度，分别通过USB接口与所述图像预处理器进行数据传输。

[0033] 所述左右微型摄像头在左右眼镜框内的安装位置与角度均可调，融合后的全视野图像的视角θ计算公式如下：

[0034]

[0035] 式中，α为左微型摄像头的左极限视界与正前方的夹角，等同于右微型摄像头的右极限 视界与正前方的夹角，d为两个微型摄像头的中心间距，h为拍摄点m到左、右微型摄像头连 线的垂直距离；所述α由单个微型摄像头的视角范围及其与两微型摄像头所呈直线间夹角 共同决定，通过调整α和d，能够扩大总视角θ。

[0036] 所述图像预处理步骤S2中的比值匹配法采用如下步骤实现图像的快速配准：

[0037] S21、读取所述左摄像头所采集图像最右侧一列像素点和每间隔5个像素点对应位置处 的像素点，共取n列像素点， n为大于1的自然数，设每列像素点个数为m，m为大于1的自然 数；

[0038] S22、计算所述每列像素点相邻行之间的比值，即得到（m-1）×n个比值，将所述比值结果 直接存放在数组P中；

[0039] S23、读取所述右摄像头采集的图像中任意每间隔5个像素点的n列像素点，取k组，k为大 于1的自然数；

[0040] S24、计算所述k组像素点每一列相邻行像素点之间的比值，将所述比值结果直接存放在 数组Q中；

[0041] S25、根据左摄像头所采集的图像中的特征集模板，即数组P，在右摄像头所采集图像的特征集合中寻找相应的匹配，得出配准区域。

[0042] 所述图像预处理步骤S6后还需要进行像素降低，目的是为了降低图像的像素值，保证与所述图像编码部分训练原子库所选图像块像素 值吻合，具体步骤如下：

[0043] S7、根据实际要求将原始图像进行分割，将原始图像平均分割成2n×2n(n=4、5或6)个图像块；

[0044] S8、计算每个分割区域像素点的平均值，然后将该平均值赋给该分割区域中的所有点；

[0045] S9、将灰度级设定为5级，即：灰度值为0、60、120、180、255，然后，根据原图像每个像素点对应的灰度范围为其赋新的灰度值。

[0046] 所述图像编码步骤C2中H-H模型的求解，具体步骤如下：

[0047] C21、建立H-H神经元模型；

[0048] C22、采用Runge-Kutta算法和迭代方法求解H-H神经元模型，建立电刺激强度与膜电压 的关系；

[0049] C23、模拟不同电刺激方案下神经元的发放情况，包括直流电刺激、脉冲波刺激、双相不 对称矩形脉冲刺激和双相不对称梯度脉冲刺激，从功耗低、脑组织损伤小的角度分析它们的优缺点，最终确定双相不对称梯度脉冲刺激为最佳的电刺激方案。

[0050] 本发明为视皮层假体的设计提供了核心的技术方案。本发明提供的图像处理装置能够实现将双摄像头采集的图像信息进行实时融合与显著性区域检测，并将显著区域的目标信息转化为用于刺激视皮层的电流信号；图像处理装置所采用的方法能够自动检测全视野图像中的显著性区域，符合视觉系统的注意机制，此外借鉴了稀疏编码的思想，利用尽可能少的神经元发放表征了图像中的显著性区域的目标信息，一方面去除了图像中相对不关键的信息，减少了编码的复杂度；另一方面减少激活神经元的数目，从而减少对组织的损伤。同时，本发明最终选定的双相不对称梯度脉冲刺激方案，在保证诱发神经元发放等量动作电位的前提条件下，能够有效降低刺激电流的强度，从而达到有效减少脑组织损伤的目的。

附图说明

[0051] 图1是本发明所述图像处理装置的结构示意图。

[0052] 图2是本发明所述图像处理装置的使用状态示意图。

[0053] 图3是本发明所述渐入渐出的线性融合原理示意图。

[0054] 图4是本发明所述神经纤维膜并联电导模型的电路原理图。

[0055] 图5是本发明所述H-H模型求解的软件流程图。

[0056] 图6是本发明所述电刺激强度Id取不同值得到的仿真结果图。

[0057] 图7是本发明所述刺激周期T=50ms占空比n=0.5时不同电刺激强度Id下得到的神经元仿真图。

[0058] 图8是本发明所述设定脉冲的刺激强度为20μA/cm2，脉冲刺激周期T=50ms，改变脉冲的占空比n的值得到的神经元仿真图。

[0059] 图9是本发明所述刺激强度Id与产生动作电位所需最小占空比n的关系图。

[0060] 图10是本发明所述双相不对称矩形脉冲刺激神经元时的神经元电位变化图。

[0061] 图11是本发明所述双相不对称梯度脉冲一个周期的波形图。

[0062] 图12是本发明所述改进方案刺激神经元电位变化图。

[0063] 图13是本发明所述图像配准算法流程图。

[0064] 图14是本发明所述左、右微型摄像头在眼睛框内的安装位置与角度关系图。

具体实施方式

[0065] 如图1、2所示，本发明所述视皮层假体的图像处理装置，包括分别对应安装在眼镜 架1左右眼镜框内的左右微型摄像头2、3，和安装在左右眼镜框连接架上的图像预处理器4， 以及图像编码器5和无线传输模块。左右微型摄像头2、3的安装位置根据双眼融合后全视野 图像的视角范围来调整。保证全视野图像的视角范围为120 150度，且连续可调。左右微型 ~摄像头2、3将各自同步实时采集的视频信息输出给图像预处理器4，图像预处理器4将接收 的两路所述视频信息进行预处理并拼接融合成一幅全视野图像输出给图像编码器5，图像 编码器5将接收的全视野图像做进一步的处理与编码，转化为人体视觉系统接受的电信号。 左右微型摄像头2、3均选择CMOS图像传感器摄像头。

[0066] 选择CMOS图像传感器摄像头理由如下：

[0067] 摄像头根据感光器件的不同分为：CMOS和CCD两种类型。CCD的功耗远远大于CMOS；而 CMOS往往通透性一般，对实物的色彩还原能力偏弱，成像效果没有CCD好。然而本发明中更 关心功耗问题，对于成像质量要求不高，故本发明选择CMOS摄像头进行图像采集。

[0068] 考虑到若摄像头的视角范围大于90度成像会产生难以忍受的畸变，不利于后续的 分析，因此限制左右CMOS图像传感器摄像头的视角范围≤90度，分别通过USB接口与图像预 处理器4进行数据传输。

[0069] 通过调整左、右微型摄像头2、3在眼睛框内的安装位置和角度（如图14 所示），可 使得输出的全视野图像视角优于人眼视角，计算公式如下：

[0070]

[0071] 式中，α为左微型摄像头2的左极限视界与正前方的夹角，等同于右微型摄像头3的右极限视界与正前方的夹角，d为两个微型摄像头2、3的中心间距，h为拍摄点m到左、右微型摄像头2、3连线的垂直距离；所述α由单个微型摄像头的视角范围以及单微型摄像头与两微型摄像头所呈直线间的夹角共同决定，通过调整α和d，能够扩大总视角θ。

[0072] 所述视皮层假体的图像处理方法包括图像预处理和图像编码两部分；

[0073] 一、所述图像预处理按照下述步骤进行：

[0074] S1、图像校正

[0075] 采用模板法对所述左、右微型摄像头2、3传输来的实时视频信息进行模板制定、特征点提取和求解镜头畸变系数；

[0076] S2、图像配准

[0077] 如图13所示，采用比值匹配法实现图像的配准，具体步骤如下：

[0078] S21、取所述左摄像头2所采集的图像的最右侧一列像素点和每间隔5个像素点对应位置处的像素点，共取n组像素点，n为大于1的自然数；

[0079] S22、计算所述n组像素点的比值，将所述比值结果直接存放在数组m中，m为大于1的自然数；

[0080] S23、在所述右摄像头3采集的图像中同样每间隔5个像素点的两列各自取出m＋n个像素点，计算其比值，然后将该比值结果存入数组k中，k为大于1的自然数；

[0081] S24、利用左摄像头2所采集的图像中的比值模板在右摄像头3所采集的图像中寻找相应的匹配，得出结果；

[0082] S3、图像融合

[0083] 图像融合是为了解决图像匹配结束后左、右摄像头2、3相邻视角拼接的交界处问题，这里采用渐入渐出的线性融合方法实现双路视频图像的融合；即：将左、右摄像头2、3相邻视角图像的像素值进行加权平均得到重叠区每一个像素点的值；其中像素点到图像重叠区两边缘的距离作为融合比重；其原理图如图3所示。

[0084] 图3中T 表示融合过渡带，W 表示权重。Wl 表示左视图图像的权重系数，Wr 表示右视图图像的权重系数，其中W1和Wr满足如下关系：

[0085] （1）

[0086] S4、显著性区域检测

[0087] 采用简单线性迭代（SLIC）的超像素图像分割方法与马尔科夫吸收链相结合的方法提取融合后全视野图像的显著性区域； 具体实现步骤如下：

[0088] 首先，使用 SLIC分割方法进行图像分割，并将图中各超像素点作为节点，定义位于边界上的节点为吸收状态点，其余节点为临时状态点，并使得所有边界上的吸收状态节点保持不相连，所有临时状态节点为相连，计算相邻的节点i和j的边 的权重值 ：

[0089] （2）

[0090] 式中，x表示对应节点处的像素值， 是控制权值强度的常量。定义一个相似矩阵A，其中，元素为：

[0091] (3)

[0092] (4)

[0093]   (5)

[0094] N(i)表示与i节点连接的节点，最终可以得到一个转换矩阵P，而矩阵P中各元素是从给定的一系列状态S= 中，从状态 转移到状态 的概率。假定前t个节点为临时状态，后r个状态为吸收状态，定义矩阵 为在任意临时状态节点间的转移概率， 为从任一临时状态移动到任一吸收状态

[0095] 的概率，则矩阵P的规范形式为：

[0096] (6)

[0097] 这里I是 的标准矩阵。结合(5)(6)可以推导出矩阵Q，进而在吸收链中推导出一个基础矩阵 ，然后对每个临时状态，可以计算出它转移到吸收状态的吸收时间：

[0098] (7)

[0099] 这里c是一个所有元素为1的t维列向量。然后通过标准化吸收时间向量y，可以得到显著性映射S：

[0100] S  ( i )=   i=1, 2,..., t,      (8)

[0101] 由于节点转移到吸收状态时，转移次数多，花费时间长的区域即为显著性区域。故 S(i)值大的地方，生成的8位深度灰度图像中区域越亮，从而检测出图中显著性区域。

[0102] S5、边缘提取

[0103] 对显著性区域内的目标图像转化为灰度图像病进行二值化处理，以增强灰度对比，便于边缘检测，然后进行边缘提取。常用的边缘提取算法有：Canny算子边缘检测，Sobel算子边缘检测，考虑到Canny算子进行的边缘检测是基于灰度的，且处理方法更加符合人眼的视觉机制，故本发明最终选择采用Canny算子检测算法进行边缘提取；核心步骤即进行降噪、寻找图像中的亮度梯度、在图像中跟踪边缘；

[0104] S6、图像增强

[0105] 采用形态学方法中的腐蚀和膨胀进行图像增强；

[0106] ① 胀操作的原理

[0107] 用3×3的结构元素，扫描二值图像的每一个像素，用结构元素与其覆盖的二值图像做“与”运算，如果都为0，结构图像的该像素为0，否则为1。处理后的结果：使二值图像扩大一圈。

[0108] ②腐蚀操作的原理

[0109] 用3×3的结构元素，扫描二值图像的每一个像素，用结构元素与其覆盖的二值图像做“与”运算，如果都为1，结构图像的该像素为1，否则为0。处理后的结果：使二值图像缩小一圈。

[0110]  S7、像素降低

[0111] 由于目前视皮层假体中电极阵列制作水平以及电极植入技术的限制，仅能在皮层植入的有限数量的电极阵列，这就要求拟编码图像具有尽可能低的像素，目的是与电极阵列的针数相匹配。因此，在保证图像质量的前提下，对输入的图像进行像素降低就显得尤为必要。像素处理的方法有很多，本发明根据图像后续编码传输要求设计了有针对性的像素降低方法。具体步骤如下：

[0112] S71、根据实际的要求将输入的图像进行分割；分割的要求就是最终的像素大小，假如图像的最后像素是2n×2n(n=4、5或6)；将原始图像平均分割成2n×2n(n=4、5或6)个图像块；

[0113] S72、计算每个分割区域像素点的平均值，然后将该平均值赋给该分割区域中的所有点；

[0114] S73、将灰度级设定为5级，即：灰度值为0、60、120、180、255，然后，根据原图像每个像素点对应的灰度范围为其赋予新的灰度值。

[0115] 二、所述图像编码按照下述进步骤行：

[0116] C1、从自然图像中随机采集5000个2n×2n(n=4、5或6)个图像块，对每个图像块分别进行去相关性、降维预处理操作，其中所述自然图像来自于广泛用于测试稀疏编码模型的图像库；

[0117] C2、利用所述自然图像高阶统计特性与视皮层神经元稀疏响应之间的联系，采用拓扑独立分量分析（Independent Component Analysis，ICA）的方法对自然图像块组成的集合进行特征提取，得到与视皮层神经元感受特性相符的完备原子库；

[0118] C3、采用所述原子库，对处理后的图像进行稀疏分解，得到对应的响应系数：

[0119] 假设输入的刺激图像为I（x，y），则该函数可以由基函数线性叠加表示，如公式9所示：

[0120] （9）

[0121] 其中， 是基函数的系数，表示的是神经元响应系数，该系数是发明中需要求解得到的重要参数； 是基函数；x表示图像中横坐标；y表示图像中纵坐标。

[0122] C4、设置阈值，将响应小于阈值的神经元响应置零，得到新的响应矩阵ϕ’；

[0123] C5、根据新的所述响应矩阵ϕ’采用公式2重建图像，记为I’；

[0124] C6、采用公式10计算所述重建图像与原始图像的误差率E；

[0125] （10）

[0126] C7、若所述误差大于10%，则返回步骤S4降低阈值，若小于10%则结束，输出响应系数矩阵，从而达到减少需要在视皮层激活神经元的数量；

[0127] C8、根据给定响应系数和神经元的特征调谐参数，求解电流刺激方案；即：通过H-H模型模拟特定神经元的特征调谐特性，来确定适当的电流刺激方案，以激活视皮层中的神经元细胞产生拟定的响应模式，从而使患者感知到对应的图像模式。

[0128] 利用H-H模型（神经元模型）确定电流刺激方案的方法如下：

[0129] 1、H-H模型的建模过程

[0130] 在H-H 模型中将单位膜面积的神经纤维等效成如图4所示的并联电路模型。其中，分别为K+，Na+和漏电离子的通透电导率， 为对应的平衡电压。

[0131] 根据图4的电路模型，有如下的四动力学变量方程，公式如下所示：

[0132] （11）

[0133] （12）

[0134] （13）

[0135] （14）

[0136] 其中， 是单位膜面积的膜电压，n是膜内分子的概率，m是膜外分子的概率，h为未激活分子的概率。 分别为对应变量的系数，针对特定的神经元其值是固定的，通过不同的 值可以求解这些系数的值，对应关系如下：

[0137] （15）

[0138] （16）

[0139] （17）

[0140] （18）

[0141] （19）

[0142] （20）

[0143] 2、H-H模型的求解

[0144] H-H模型的求解过程如下：以 h=0.025m 为步长，计算离子通道的电导率、膜的动作电位以及他们在t+h时候发生的变化，得出离子电流的增减率从而推算下一步时的膜电位。对于这四个微分方程的具体的求解思想是利用中间步点的值来代替高阶导数，按步长进行迭代计算，具体的公式如下：

[0145] （21）

[0146] （22）

[0147] （23）

[0148] （24）

[0149] （25）

[0150] 其中，在常见的标准算法公式中，取 ，。

[0151] 在Runge-Kutta算法中，计算 的值，从而可以得到下一步t+h时的膜电压和各离子通道电导值的大小。

[0152] 使用Runge-Kutta算法和迭代方法求解H-H模型，具体的软件流程如图5所示。

[0153] 3、最佳电刺激方案的确定

[0154] 模拟不同电流刺激方案下的神经元发放情况，包括直流电刺激、脉冲波刺激、双相不对称矩形脉冲刺激和双相不对称梯度脉冲刺激，确定最佳电刺激方案。仿真分析如下：

[0155] （1）直流电刺激下神经元动作电位仿真分析

[0156] 通过改变电刺激强度Id,观察在不同的刺激强度下神经元的动作电位发放情况。不同电刺激强度Id下的仿真结果如图6所示。由图6可以得出如下结果：

[0157] ①当刺激电流Id≤4.9μA/cm2 时，神经元没有动作电位产生。

[0158] ②当刺激电流5.0μA/cm2≤Id≤6.1μA/cm2时，神经元会产生有限个动作电位，随着刺激电流的强度增大，产生的动作电位个数也越来越多。

[0159] ③当刺激电流6.1μA/cm2≤Id＜120μA/cm2时，神经元会产生周期性的动作电位，随着电刺激强度的增加，神经元产生动作电位的峰峰间隔变小，动作电位个数相应增加。

[0160] ④当刺激电流Id≥120μA/cm2时，神经元动作电位呈现过阻尼震荡现象，这说明神经元放电的时间间隔过小。这种情况对视皮层进行电刺激时会产生明显的疼痛感。

[0161] （2）脉冲波刺激下神经元动作电位仿真分析

[0162] 脉冲波电刺激常用于治疗某些疾病，相较于直流电刺激对细胞的伤害较小。使用脉冲波刺激时，主要考虑两个参数：脉冲的刺激强度和脉冲的占空比。仿真结果如图7所示。

[0163] 图7中a、b、c、d分别表示不同强度的电刺激仿真结果（其中占空比为n为0.5），由图可以得出如下结果：

[0164] ①当刺激电流Id≤3.5μA/cm2 时，神经元没有动作电位产生，最小刺激强度比直流电刺激的最小刺激强度低。

[0165] ②当刺激电流在3.5μA/cm2＜Id＜150μA/cm2范围内时，神经元会产生有限个动作电位，随着刺激电流的强度增大，产生的动作电位个数逐渐增多。但是随着刺激电流强度增大到一定程度之后动作电位的个数不会在改变，一个刺激周期（T=50ms）内，动作电位的最大个数是3。

[0166] ③当刺激电流Id≥120μA/cm2时，随着脉冲刺激强度增加，神经元产生的动作电位的最大值衰减为1，这个和直流电刺激情况一致。

[0167] 图8是设定脉冲的刺激强度为20μA/cm2，脉冲刺激周期T=50ms，改变脉冲的占空比n的值得到的神经元仿真图。从图8可知：

[0168] ①当n≤0.0022时，没有动作电位产生，如图8（a）所示。

[0169] ②当0.0023＜n≤1时，神经元会有动作电位产生，且会随着q值的增大动作电位的个数也会增加，一个刺激周期（T=50ms）内动作电位最多产生4个，如图8（d）所示。

[0170] 通过仿真实验发现，在相同的刺激周期，不同刺激强度的情况下，要使神经元有动作电位产生所需的最小占空比不同，即需要的最小刺激脉冲宽度不同，上述仿真中得到了不同的电刺激强度下对应的最小占空比值。

[0171] 图9表示的是刺激周期（T）为50ms时电刺激强度Id与神经元放电最小的占空比关系。由图可以看出随着脉冲刺激强度的增加，神经元产生动作电位所需的最小脉冲宽度减小，最后趋于稳定。

[0172] （3）双相不对称矩形脉冲刺激下动作电位仿真分析

[0173] 双相不对称矩形脉冲是由阴极相、阳极相和相间间隔三部分组成的矩形脉冲，先使用阴极脉冲刺激，所需要的刺激脉冲强度比单独使用正脉冲刺激需要的脉冲强度更小，并且双相不对称矩形脉冲能够避免在电极-电解液界面直流电荷的积累产生的细胞毒性，因此能够减少组织损伤。

[0174] 双相不对称矩形脉冲刺激的参数主要包括：阴极相强度、阴极相刺激时间、阴极与阳极的强度比和阴极与阳极时间间隔。

[0175] 本发明重点分析了使用双相不对称矩形脉冲刺激神经元时，阳极相的最小刺激强度。其中参数设置为：刺激周期T=50ms，阴极相的刺激时间是1ms，阴极刺激强度（Id1）是3.6μA/cm2，阳极刺激强度（Id2）改变。具体仿真结果如图10所示。

[0176] 从图 10的仿真结果可知，采用双相不对称矩形脉冲刺激与使用矩形脉冲刺激神经元所需要的最小刺激强度大幅度降低，从而大幅度地降低了刺激的能耗。

[0177] （4）双相不对称梯度脉冲刺激下动作电位仿真分析

[0178] 为了进一步减少电刺激的峰值电压，尽最大可能的降低电刺激对脑组织的损伤，本发明在双相脉冲刺激基础上进行改进，将正脉冲波形改变成分段的阶梯波形，提出了一种双相不对称梯度脉冲方案。刺激波形如图11所示。其中，梯度脉冲的刺激强度最大的是一段持续1ms的0.45μA/cm2的矩形脉冲，剩下的部分都是刺激强度低于0.45μA/cm2的矩形脉冲。通过分析可以知道，如果单独使用为持续1ms强度0.45μA/cm2的矩形脉冲刺激神经元不会发放动作电位。由图12可以看出，该刺激方案能够使神经元产生发放。从而确定了双相不对称梯度脉冲波为最佳电刺激方案。