一种引导式虚拟内窥镜导航方法技术

本发明专利技术公开了一种引导式虚拟内窥镜导航方法。该方法首先将断层图像分割成待检查被测腔体体素和背景体素，然后采用快速3D距离变换计算每个被测腔体体素到最近背景体素的距离，计算距离梯度作为内在拉力，最后将虚拟摄像机拖拽至接近被测腔体的中心位置。在内在拉力作用下，引导摄像机既能沿着接近指定方向前进，又可以尽可能远离被测腔体内壁，从而达到快速、自动浏览待检查被测腔体，同时尽可能观察到被测腔体所有内壁，视野清晰、开阔。该方法不需要预先提取待检查被测腔体的中心线，从而避免了中心线提取过程中的拓扑结构、噪声干扰等难题。

【技术实现步骤摘要】
一种引导式虚拟内窥镜导航方法
本专利技术属于图像处理与应用领域，涉及一种引导式自动虚拟内窥镜导航方法，确保成像视点逼近中心线，使得观察视野尽量开阔，且实现快速自动导航，并获得平滑的导航轨迹。
技术介绍
传统内窥镜技术通常借助于光学、视频等技术观察被测腔体内表面。虚拟内窥镜从由一系列断层图像组成的体数据场中，重建出类似于传统光学内窥镜所观察到的被测腔体内窥透视视图。优点是：(1)非侵入；(2)可实现任意角度、位置的观察，不存在视觉“死区”；(3)通过使用半透明体绘制技术，还可观察被测腔体内壁的组织；(4)数据能重用，可长期保存。为了有效观察被测腔体内部，避免景深变化带来的透视变形或视觉死区的存在，需要令观察点尽量贴近被测腔体的中心位置。目前的做法是预先提取导航路径，即被测腔体的中心线。主要的中心线提取方法有拓扑细化、最短路径法、距离编码法等。但是，目前的中心线提取方法存在不少问题。例如，拓扑细化法，采用迭代计算，计算量非常较大，造成算法非常耗时；最短路径法、距离编码法，容易造成提取的中心线容易贴近腔壁，影响观察视野；基于距离场的边界距离算法，时间效率也不够高。以上方法都利用了中心线的中心性、连通性、单一性、鲁棒性等特性。上海交通大学改进了当前距离映射提取中心线的方法（WaveletActiveMediaTechnologyandInformationProcessing(ICCWAMTIP),19-21Dec,2014,Chengdu,China），使用边界距离场的平方代表体素的相对位置，以便于后续的开平方计算。经过手工标记、拓扑细化、距离映射和水平集匹配获得中心线。此方法保持了连通性，但是对于表面的形状变化很敏感。华南理工大学提出基于图像分割的路径规划法（ControlConference,July262007-June31，2007,Hunan,China，先利用阈值分割算法和边缘检测方法获得一张切片图像中主体的一部分，随后自适应地选择该区域的种子和增长区域的阈值，所有种子都被分别限定在一个四边形中存入一系列堆栈；这些种子作为导航的关键点，经过分割后便得到能够组成路径的所有关键点，经过平滑后则获得最终规划路径。中国科学院提出包含距离的中心线算法（IEEEInternationalSymposiumonBiomedicalImaging:NanotoMacro,2004，261-264Vol.1），并且利用像素编码自适应地去除伪分支。此方法能够保持拓扑结构，保留内部点和边界间的距离信息，但计算代价较大。BrunoM.等人（IEEEInternationalSymposiumonBiomedicalImaging:NanotoMacro,2004，1541-1544Vol.2）认为图片中的每一个物体都有独自的连接性强度定义，也有专属的种子体素集；每一个物体都可以当做是一个体素集合，此集合中的体素都是真包含于物体的，并且体素与此物体中某个种子点间的连接性高于体素与种子点集中其他种子的连接性；以此理念提出M模糊分割的方法。用面心立方晶格表达三维CT图；利用距离图找到最优体素用B样条曲线表达路径。此方法生成连通图，降低了导航的时间消耗，但复杂度较高并且需要二次分割。
技术实现思路
为解决现有方法的不足，本专利技术提供一种引导式虚拟内窥镜导航方法，避开中心线提取过程，从而避开了中心线不稳定、存在毛刺、拓扑复杂等缺点，而是把距离梯度作为虚拟内在拉力，迫使虚拟摄像机自动贴近中心线向内探索，提高虚拟内窥镜系统的成像效率，能简化操作、节省检查时间，帮助实现无遗漏全面检查。为达到上述目的，本专利技术采用的技术方案是：一种引导式虚拟内窥镜导航方法，其特征在于：将断层图像分割成被测腔体体素和背景体素，用快速3D距离变换计算每个被测腔体体素到最近背景体素的距离；随后利用距离梯度方向和距离值引入虚拟拉力，进行导航方向的校正，修正虚拟摄像机当前的运动方向，使得摄像机尽量远离被测腔体内壁，从而使摄像机接近中心线。进一步，在校正导航方向的过程中，为实现自动导航和获得光滑的导航路径，需要预先计算导航轨迹控制点，可以边计算边前进；采用分段三次Bezier曲线拼接方式，生成具有至少C1连续的光滑导航轨迹，从而防止摄像机产生剧烈“抖动”。更具体地，预先计算后2个单位时刻的2点坐标，将坐标集合分组构成一段光滑Bezier曲线段，在实现一小段组内平滑后，实时更新分组进行组间平滑，最后插值补点令路径具有至少C1连续，以便于摄像机沿平滑中心路径顺利前进展示被测腔体的内部结构。本专利技术首先将断层图像分割成待检查被测腔体体素和背景体素，然后采用快速3D距离变换计算每个被测腔体体素到最近背景体素的距离，计算距离梯度作为内在拉力，最后将虚拟摄像机拖拽至接近被测腔体的中心位置。在内在拉力作用下，引导摄像机既能沿着接近指定方向前进，又可以尽可能远离被测腔体内壁，从而达到快速、自动浏览待检查被测腔体，同时尽可能观察到被测腔体所有内壁，视野清晰、开阔。该方法不需要预先提取待检查被测腔体的中心线，从而避免了中心线提取过程中的拓扑结构、噪声干扰等难题。附图说明图1是内窥镜导航系统方法流程图图2本专利技术方向校正示意图图3本专利技术三次Bezier平滑与插值示意图。具体实施方式以下结合附图对本专利技术的实施例作进一步的详细说明。图1给出内窥镜导航系统方法流程图。从图1可以看出，工作站将医学图像通过装载模块加载入混合场景管理模块，混合场景管理模块将待调整的信息发送给二维处理模块，信息经过常规处理后反馈回混合场景管理模块显示。同样，混合场景管理模块将数据发送给三维处理模块，数据经过处理返回混合场景管理模块显示。系统导航模块与自三维成像模块通信，使虚拟摄影机实时获得当前图像，送入混合场景管理模块显示。图像保存模块接受来自二维处理模块、三维成像模块和系统导航模块的数据存储，实现这些数据的存盘处理。一种引导式虚拟内窥镜导航方法包括步骤如下：1.速度方向校正设被测腔体体素集合为F(Foreground)，背景集合为B(background)，采用快速3D距离变换，计算每个体素v(x,y,z)到背景体素B的最近距离：显然，当时，。然后，计算每个体素的3D距离梯度：(2)根据距离值D及距离梯度方向G，修正虚拟摄像机当前的运动方向V，使得虚拟摄像机朝着距离最大方向运动。如图2所示，我们首先计算当前位置P0处的距离梯度方向G，并把当前的运动方向V投影到当前切平面，得到归一化的切线方向T。最后，根据该归一化的切线方向T，采用公式(3)修正摄像机的运动方向，从而得到修正后的虚拟摄像机运动方向V’：(3)(4)其中，权重。公式(4)是对修正后的运动方向向量归一化。虚拟摄像机离被测腔体中心线越远（即被测腔体内壁越近）时，也就是距离D越小，此时应该越小，虚拟摄像机的前进方向越需要偏向中心线方向修正，故要求系数函数f(D)是单调递增函数，即，同时要求。在其中，f(D)定义如下基于高斯的单调递增函数：(5)其中，σ为距离作用范围参数，其作用是控制运动方向的修正量。导航轨迹控制点的计算为实现自动导航和获得光滑的导航路径，需要预先计算接下来i个单位时刻虚拟摄像机的导航轨迹控制点坐标，以便在导航轨迹控制点坐标之间插出光滑的导航本文档来自技高网...

【技术保护点】
一种引导式虚拟内窥镜导航方法，其特征在于：将断层图像分割成被测腔体体素和背景体素，用快速3D距离变换计算每个被测腔体体素到最近背景体素的距离；随后利用距离梯度方向和距离值引入虚拟拉力，进行导航方向的校正，修正虚拟摄像机当前的运动方向，使得摄像机尽量远离被测腔体内壁，从而使摄像机接近中心线。

【技术特征摘要】
1.一种引导式虚拟内窥镜导航方法，其特征在于：将断层图像分割成被测腔体体素和背景体素，用快速3D距离变换计算每个被测腔体体素到最近背景体素的距离；随后利用距离梯度方向和距离值引入虚拟拉力，进行导航方向的校正，修正虚拟摄像机当前的运动方向，使得摄像机尽量远离被测腔体内壁，从而使摄像机接近中心线；在校正导航方向的过程中，为实现自动导航和获得光滑的导航路径，预先计算导航轨迹控制点，边计算边前进；设被测腔体体素集合为F，背景集合为B，采用快速3D距离变换，计算每个体素v(x,y,z)到背景体素B的最近距离：当v(x,y,z)∈B时，D(x,y,z)＝0；然后，计算每个体素的3D距离梯度：根据距离值D及距离梯度方向G，修正虚拟摄像机当前的运动方向V，使得虚拟摄像机朝着距离最大方向运动；首先计算当前位置P0处的距离梯度方向G，并把当前的运动方向V投影到当前切平面，得到归一化的切线方向T；最后，根据该归一化的切线方向T，采用公式(3)修正摄像机的运动方向，从而得到修正后的虚拟摄像机运动方向V’：v′＝αT+(1-α)G(3)其中，权重α＝f(D)，公式(4)是对修正后的运动方向向量归一化；系数函数f(D)是单调递增函数，即α∝D，同时要求0≤α≤1；其中，f(D)定义如下基于高斯的单调递增函数：2.根据权利要求1所述的引导式虚拟内窥镜导航方法，其特征在于：采用分段三次Bez...

【专利技术属性】
技术研发人员：袁非牛，
申请(专利权)人：袁非牛，
类型：发明
国别省市：江西;36