视神经纤维自动识别与分析关键技术研究

论文摘要

目前,在有关人和动物的神经科学研究中,神经纤维形态分析被广泛应用于正常或病理状态下的神经组织研究或诊断中。神经纤维及轴突的形态分析,对于对正确理解神经发育、损伤及退化程度,掌握神经纤维形态特征规律及分布具有重要意义。神经纤维形态分析能否客观反映神经科学现象,关键取决于兴趣区域(如轴突或髓鞘)的选择提取是否科学、准确。在神经纤维形态分析研究之初,这项工作依赖手动方式提取兴趣区域,即由有经验的人员从显微图像中勾画出兴趣区域的轮廓及边界。随着研究的深入,手动方式的弊端明显体现:第一,手动方式费时费力,分割结果的优劣很大程度上取决于操作者的经验知识,客观性和可重复性较差;第二,对于数量庞大的神经纤维,少量采样可以缓解科研人员的繁重劳动,但是获取数据的科学性和可靠性却受影响。因而,神经纤维形态的自动识别与分析已成为图像处理和神经科学领域科研工作者共同关注的研究热点。从上世纪80年代,相关的神经纤维识别算法和商业化软件不断被推出和改进,使神经纤维的识别效果不断提高。精确的神经纤维形态识别和分析是建立在准确的图像分割基础之上。由于神经系统组织成分复杂多变,给神经纤维识别带来诸多困难,例如细小神经纤维的识别、复杂背景下特征的提取、不规则神经纤维的提取等。这些困难在神经纤维分布密集的神经(如视神经)中表现尤为突出。本课题研究的首要目标是开发出适合视神经纤维分割的有效方法。在纤维识别中,我们充分利用神经纤维的解剖结构特点,先获取纤维髓鞘的中心线,即粗略锁定神经纤维;在此基础上通过中心线的内移、外移及曲线演化,进而获得轴突和髓鞘的精细轮廓。这一技术思想也是本研究的创新主体。在上文所述创新思想的驱动之下,本研究采用改进后的Watersnake进行神经纤维分割。在实验中分两步进行:第一步,由Watershed生成连续且拓扑保持的神经纤维的初始骨架,并剔除假性纤维和边界位置的不完整纤维;第二步,由Snake将初始骨架引导到轴突和髓鞘的精细位置上。本算法提取的神经纤维保持了位置的准确和外形的光滑。在算法的性能评价实验中,随机抽取纤维图像,将自动分割结果与采用PHOTOSHOP软件手动分割结果进行对照分析。在纤维识别率分析中,自动识别正确率95.7%,假阳性率3.3%。在轴突直径频数分布对照分析中,自动识别与手动识别方式下的轴突直径频数分布直方图相吻合,波峰位置未有偏移。从目前的实验结果来看,本研究的识别率稍高于Romero等人的研究(94.8%),而假阳性率方面不及前者(1.0%)。研究中发现,假阳性纤维主要来源于拥挤纤维簇之间形成的闭合区域。在Romero等人的研究中,其研究对象是坐骨神经纤维,而坐骨神经纤维的分布比视神经纤维的分布稀疏,这也许是其假阳性率低于本研究的原因之一。本研究中的图像来源于分辨率较高的透射电镜底片,Watersnake技术对分辨率低的图像(如光镜图像)是否同样适用,在后期研究中我们将继续探索和改进,拓宽该技术的应用范围。

论文目录

• 摘要
• ABSTRACT
• 第一章 绪论
• 1.1.研究背景和意义
• 1.2.神经纤维图像分析的研究现状和发展趋势
• 1.2.1.图像获取方式
• 1.2.2.神经纤维图像处理与分析
• 1.3.本文的研究内容和安排
• 1.4.本章小结
• 第二章 标本制备及图像获取
• 2.1.概述
• 2.1.1.神经轴突、纤维
• 2.1.2.视神经解剖基础
• 2.2.材料
• 2.3.方法
• 2.4.本章小结
• 第三章 医学显微图像的滤波降噪算法
• 3.1.均值滤波和中值滤波
• 3.2.邻域比较法
• 3.3.与邻域平均值比较设阈值置换法
• 3.4.基于灰阶梯度的自适应非线性滤波器
• 3.4.1.梯度计算
• 3.4.2.梯度自适应滤波
• 3.5.结果
• 3.6.本章小结
• 第四章 灰度图像二值化
• 4.1.阈值分割
• 4.1.1.阈值分割原理
• 4.2.基于边缘的分割
• 4.2.1.坎尼算子
• 4.2.2.拉普拉斯算子
• 4.3.区域生长
• 4.3.1.生长准则和过程
• 4.4.霍夫（Hough）变换
• 4.4.1.基本原理
• 4.4.2.椭圆检测
• 4.5.形态学图像处理
• 4.5.1.结构元素的选取
• 4.5.2.神经纤维图像形态学处理
• 4.6.本章小结
• 第五章 距离变换
• 5.1.概述
• 5.2.距离
• 5.2.1.数字图像中常用的距离
• 5.2.2.点到点集的距离
• 5.3.距离计算方法
• 5.3.1.Chamfer距离
• 5.4.本章小结
• 第六章 基于WATERSNAKE的纤维识别
• 6.1.概述
• 6.2.分水岭（watershed）变换
• 6.2.1.变换原理
• 6.2.2.模拟漫水算法
• 6.3.变形模型snake
• 6.3.1.参数变形模型的能量最小公式
• 6.3.2.参数变形模型的实现
• 6.3.3.梯度矢量流（GVF）变形模型
• 6.4.神经纤维分割
• 6.5.本章小结
• 第七章 形态分析及评价
• 7.1.形态数据测量
• 7.2.形状分析
• 7.2.1.圆形度
• 7.3.实验结果与性能评价
• 7.3.1.纤维检测结果分析
• 7.3.2.纤维直径结果评价与分析
• 7.3.3.大鼠视觉神经纤维形态分布谱分析
• 7.4.本章小结
• 第八章 总结和展望
• 8.1.研究总结
• 8.2.研究展望
• 参考文献
• 在读期间发表和录用的论文
• 致谢