一种术前三维影像与术中透视图像的无标定物配准方法技术

本发明专利技术涉及一种术前三维影像与术中透视图像的无标定物配准方法，包括以下步骤：1)建立一个模拟现实的解剖学和频率相关的系统噪声模型M；2)将CT/MRI数据渲染生成三维体数据，通过为DRR图像添加解剖噪声和频率噪声，生成比较逼真的类似于X光片的影像数据，提取骨骼轮廓；3)手术中，X射线机获取不同角度的两张X‑Ray透视图像；4)根据手术需求生成N个DRR图像，得到每个图像的位姿坐标P；5)将DRR算法生成的影像数据进行处理，然后选择没有骨折部位的感兴趣区域，提取轮廓；6)按同样方法提取真实X‑Ray的轮廓线；7)基于轮廓线和特征点相关的性质为术前三维影像和X‑Ray透视图像精确配准，输出精确配准后的图像。本方法能提高配准精度，使手术导航操作更直观。

【技术实现步骤摘要】
一种术前三维影像与术中透视图像的无标定物配准方法
本专利技术涉及一种术前三维影像与术中透视图像的无标定物配准方法，特别是医疗机器人辅助进行手术导航过程的医学图像配准方法。
技术介绍
目前，在常见医疗机器人辅助手术导航定位过程中，如在股骨颈骨折手术中，传统方法是在股骨或骨盆上植入标定物，以建立机器人坐标系与手术空间坐标系之间的映射关系，并在此基础上进行手术路径规划和定位。该方法能够大大增加股骨骨折复位的精确度，但同样会给患者增加创伤及感染的风险。一种方法是利用2D/3D图像配准算法的图像导航的放射手术方法和系统，其利用了分层和迭代的2D/3D配准算法求出内层面和外层面的变换参数。此配准算法的缺点是准确度较低，速度较慢。另一种方法是使用DRR技术对MRI和X-Ray图像配准的算法和系统，其特点在于先对MRI图像做提取分割，由分割后的图像生成DRR用于配准，缺点是算法需要事先训练，从而获取特征用于分割，受限于训练的集合大小及范围，而且成本较大。形状描述子已经广泛应用于现在的形状匹配算法中，轮廓段的角度信息和弧长弦长的关系表示其边界信息，采样点和质心点之间的距离关系表示其区域信息，所以能够对轮廓段的形状有完全定量的描述，更全面地描述轮廓的形状，显著提高描述能力。本专利技术将形状描述子方法应用于术前三维影像与术中透视图像的无标定物配准过程中，可提高配准精度和速度。
技术实现思路
本专利技术提供了一种术前三维影像与术中透视图像的无标定物配准方法，可应用于外科手术导航系统，使其操作简单，速度快，准确度高。为实现上述目的，本专利技术中术前三维影像与术中透视图像的无标定物配准方法包括以下步骤：1)首先利用人体模型，设计出一个模拟现实的解剖学和频率相关的系统噪声(包括量子噪声)模型M：2)在手术开始前，由CT或MRI获取一系列三维图像，根据手术需求，在冠状位、矢状位或横断位的视角方向上，从0°开始每间隔预设度数根据CT或MRI提供的三维图像生成N张DRR图像，获取对应的位姿参数P；其中，P＝(θx,θy,θz,X,Y,Z)，θx、θy、θz表示旋转方向，X、Y、Z表示坐标系下各个方向下的平移量；3)在手术过程中，C型臂上的X射线机获取不同角度的X-Ray透视图像；4)对DRR图像i应用模型M，添加解剖学噪声和频率相关噪声，透视图像上的每一个像素均做二维高斯加权归一化处理，每一个像素原始灰度值减去高斯加权均值，然后除以高斯加权均方差从而获取归一化后的像素灰度值；5)采用线性直方图匹配算法，校正DRR图像i与X-Ray透射图像的灰度差异，选择灰度值大于投影图像灰度均值的区域进行直方图匹配，消除背景像素对直方图的影响；6)对DRR图像i与X-Ray用OpenCV3.4.0中的canny算子进行边缘检测，再用findContours()分别进行轮廓提取，透射图像进行轮廓和直方图相似度比较，如果两组图像匹配，获取X-Ray透射图像位姿参数，并进入下一步；否则提取DRR图像i+1信息，并返回步骤4)；7)进行轮廓相关性和点配准为术前三维影像和X-Ray透视图像精确配准，输出精确配准后的图像，完成术前三维影像与术中透视图像的无标定物配准；所述步骤2)中，从0°开始每间隔预设度数1°或2°根据CT/MRI提供的三维图像生成N个DRR图像；所述步骤3)中，C型臂上的X射线机获取X-Ray透视图像至少获取两张。所述步骤1)和步骤6)中，所述轮廓匹配相关作为相似度测量的配准算法，轮廓匹配相似度评价方法如下：1)形状描述子：对于一段给定的轮廓段S，对其进行像素点采样，如图3所示，首先采样得到采样点Pi＝(xi,yi)，(i＝1,2,…,N)，N是轮廓段采样点的个数，首先计算该段轮廓的质心点G，然后对每个采样点找到其最远距离点fPi，通过计算采样点Pi到其它所有采样点的距离，最远距离点fPi；函数DS(Pi)计算每个采样点的形状描述子为：式中：di是采样点Pi质心点G的归一化距离，由Pi和fPi之间的距离来归一化；αi是表示由和所形成的角度值，由π归一化；cRi表示弦长chordLeni与弧长radLeni的比值，函数DS(Pi)3个分量都在区间[0,1]内尺度不变；得到每个采样点的形状描述子后，每条轮廓段S的形状描述子SD(S)可以由式(1)组合表示为：式中，SD(S)是3×N维矩阵；其中，每列代表该轮廓段上第i个采样点的形状描述子DS(Pi)；每行代表轮廓段的距离信息，角度信息和弧度信息；得到完整轮廓段的形状描述子之后可以进行模板匹配；通过形状描述子间的相似度来表示两条轮廓段的匹配关系；形状描述子以矩阵表示，二者的矩阵相关系数可以反映其关系得密切程度，就能代表不同轮廓段的匹配程度，不同形状描述子之间的相关系数被称为轮廓的匹配系数；2)评价标准：轮廓匹配过程中，首先要提取DRR图像和X-Ray图像的感兴趣区域轮廓，然后根据50％-IoU和20％-IoU标准得出匹配结果。优选的，必须在DRR图像和X-Ray图像中取一个健康的骨骼边缘轮廓。进一步的，为生成的每张DRR图像添加医学解剖学和频率等噪声，来模拟真实的X-Ray光片。进一步的，对模拟出来的X-Ray图像进行图像预处理，包括亮度、对比度提高，这样可以提高DRR图像股骨边缘轮廓更真实。进一步的，选取感兴趣区域需要有一定医学背景的医生或培训人员进行操作；轮廓边缘点的选取一定要具有代表性；点的匹配如果得出的旋转或者平移矩阵太大，我们认为匹配并不成功。附图说明图1是本专利技术从不同视野角度对术前CT或MRI图像进行数字投影重建示意图；图2是笔形束在CT数据集的第一个后前(PA)切片中进入体素；图3是轮廓段形状描述子示意；图4是为本专利技术流程示意图；图5是分别从X-Ray和DRR图片提取的股骨二值图像轮廓。具体实施示例下面结合附图和实施例对本专利技术进行详细的描述，但本专利技术的的实施方式不限于此。图1～图5所示，本专利技术提供一种术前三维影像与术中透视图像的无标定物配准方法，以股骨为例，其具体包括以下步骤：1)通过实验得出X-Ray图像股骨部分的一般性信噪比。2)在手术开始前，由CT或MRI获取一系列三维图像，根据手术需求，在冠状位、矢状位或横断位的视角方向上；从0°开始每间隔预设度数根据CT或MRI提供的三维图像生成N个DRR图像，如图1，获取对应的位姿参数P；其中，预设度数优选为1°或2°。3)为每张图像添加噪声，以便于提取股骨边缘轮廓；添加噪声函数：DRRp+s＝DRRp+scatter；式中DRRPp+s是添加了线性散布(实验性测量)的主要DRR，scatter是光束散射参数；位姿参数公式：P＝(θx,θy,θz,X,Y,Z)；式中，θx、θy本文档来自技高网...

【技术保护点】
1.一种术前三维影像与术中透视图像的无标定物配准方法，其特征在于，包括以下步骤：/n1)首先利用人体模型，设计出一个模拟现实的解剖学和频率相关的系统噪声(包括量子噪声)模型M；/n

【技术特征摘要】
1.一种术前三维影像与术中透视图像的无标定物配准方法，其特征在于，包括以下步骤：
1)首先利用人体模型，设计出一个模拟现实的解剖学和频率相关的系统噪声(包括量子噪声)模型M；



2)在手术开始前，通过CT或MRI获取一系列三维图像，根据手术需求可在冠状位、矢状位或横断位的视角方向上，从0°开始每间隔预设度数根据CT或MRI提供的三维图像生成N张DRR图像，获取对应的位姿参数P：P＝(θx,θy,θz,X,Y,Z)，θx、θy、θz表示旋转方向，X、Y、Z表示坐标系下各个方向下的平移量；
3)在手术过程中，C型臂上的X射线机获取不同角度的X-Ray透视图像；
4)对DRR图像i应用模型M，添加解剖学噪声和频率相关噪声，透视图像上的每一个像素均做二维高斯加权归一化处理，每一个像素原始灰度值减去高斯加权均值，然后除以高斯加权均方差从而获取归一化后的像素灰度值；
添加噪声函数：DRRp+s＝DRRp+scatter；
5)采用线性直方图匹配算法，校正DRR图像i与X-Ray透射图像的灰度差异，选择灰度值大于投影图像灰度均值的区域进行直方图匹配，消除背景像素对直方图的影响；
6)对DRR图像i与X-Ray用OpenCV3.4.0中的canny算子进行边缘检测，再用findContours()分别进行轮廓提取，透射图像进行轮廓和直方图相似度比较，如果两组图像匹配，获取X-Ray透射图像位姿参数，并进入下一步；否则提取DRR图像i+1信息，并返回步骤4)；
7)输出精确配准后的图像，完成术前三维影像与术中透视图像的无标定物配准；
所述步骤2)中，从0°开始每间隔预设度数1°或2°根据CT/MRI提供的三维图像生成N个DRR图像；
所述步骤3)中，C型臂上的X射线机获取X-Ray透视图像至少获取两张。


2.根据权利要求1所述的一种术前三维影像与术中透视图像的无标定物配准方法，其特征在于：必须有一种模拟生成真实...

【专利技术属性】
技术研发人员：刘勇，武利成，李亮，邵言亮，
申请(专利权)人：安徽埃克索医疗机器人有限公司，
类型：发明
国别省市：安徽;34