关键词:
图像分割
DRR
图像配准
点云配准
Huber损失函数
摘要:
随着医学成像技术和计算能力的不断发展,微创手术广泛应用于临床实践,这对于影像引导系统、图像的分割与配准技术、可视化技术以及人机相互系统的要求更高。在术中引导过程中,需要全方位的检测手术介入全过程,给予临床医生对每个细小环节的把控,这需要多种模态的医学图像共同作业,构建解剖结构的时间与空间信息互补。譬如,临床上关于脊椎相关疾病的诊断,特别是脊椎部位,无论是术前病理诊断与分析、术中引导介入过程,还是术后疗效评估,都需要扫描CT或MRI作为辅助手段,其中,CT图像对骨骼结构敏感,成像效果好,在术前诊断分析过程中,临床医生可借助CT图像快速、精准定位椎骨的病灶区域和位置。医生从图像中获取患者病灶模型是一项繁琐而耗时的工作,此外,对于同一个病人,在两个不同的观察者之间可能存在着显著的差异性,最终的图像分析高度依赖于医生的专业知识。因此,一种可校准的、可预测的以及自动化的图像分割方法非常重要,它可以节省大量时间并提供更加精确的分割结果。为了向临床医生提供实时手术环境的分割与配准模型,本文提供新的基于脊椎图像的分割与配准算法,其中包括图像相似度的测量标准。然而,目前基于CT脊柱图像的分割与配准仍然是存在一些缺陷,如分割边界模糊、分割过程中图像细节信息的丢失、配准耗时较长以及迭代效率低下等。因此,本文以下方面做了重要研究:(1)由于医学影像图像边界较为模糊,传统的语义分割网络往往缺乏低层信息、缺乏详细的语义信息,基于U-Net网络,本文提出一种改进的网络结构,该方法在U-Net的图像特征提取中引入多尺度信息和新的残差连接,可以更好地结合高层次和低层次的特征信息,防止图像的细节在上采样过程中被稀释。通过实验分析表明,本文提出的分割模型使得图像的分割边界更加清晰、骨骼检测和定位的准确度将大大提高。(2)不同的X光设备的型号不同,拍摄出的X射线图像的质量会存在较大差异,如边界模糊、特征不明显以及对比度差等情况。因此,在医学图像配准过程中,传统的相似度计算方法难以准确地评估其有效特征,针对这一缺陷,本文提出一种局部归一化互相关加权的相似度方法,通过计算单节脊椎的相似度方法,提高算法的鲁棒性。为了验证这一方案的有效性,利用DRR模拟图像进行测试,在相同的变化参数范围内,本文所提出的方案计算出的相似度具有更加明显的变化。因此,该方法能够较大调高后续配准的精确度以及缩短配准的时间。(3)搜索策略是2D/3D医学图像配准中至关重要的一步,它决定着算法能否在更快的时间内计算出最佳的相似度值。传统的优化算法,如变邻域搜索的算法(VNS)和蚁群优化算法(ACO)具有迭代次数多,搜索步长难以确定等弊端。协方差自适应进化算法(CMA_ES)具有较强的收敛速度和搜索步长自适应变化的特性,但对于较大数据空间,很容易陷入局部最优解中。针对上述问题,本文提出在CMA_ES中引入近端优化策略(PPO),通过在每次迭代中对少量样本的更新以及搜索步长的自适应变化,来提高搜索效率,从而获得全局最佳变换参数矩阵。(4)在基于点特征的配准方法中,最常用的方法是用快速点特征直方图(FPFH)匹配,损失函数选最小二乘法(TLS),它对有噪声和干扰的数据具有很强的鲁棒性,然而,该方法在存在大量离群值的条件下会产生一个较差的估计。为了针对这一缺陷,本文在TLS函数中引入Huber损失函数,可以减少数据中对于异常值的敏感性。实验结果表明,本文提出的方法可以有效降低离群值对于配准结果的影响。