Mohammadi Amrollah, Ahmadian Alireza, Azar Amir Darbandi, Sheykh Ahmad Darban, Amiri Faramarz, Alirezaie Javad
Department of Biomedical Engineering, School of Medicine, Tehran University of Medical Sciences and Research Centre for Biomedical Technology and Robotics (RCBTR), Tehran, Iran.
Rajaei Cardiovascular, Medical and Research Center, Tehran, Iran.
Int J Comput Assist Radiol Surg. 2015 Nov;10(11):1753-64. doi: 10.1007/s11548-015-1216-z. Epub 2015 May 10.
Combination of various intraoperative imaging modalities potentially can reduce error of brain shift estimation during neurosurgical operations. In the present work, a new combination of surface imaging and Doppler US images is proposed to calculate the displacements of cortical surface and deformation of internal vessels in order to estimate the targeted brain shift using a Finite Element Model (FEM). Registration error in each step and the overall performance of the method are evaluated.
The preoperative steps include constructing a FEM from MR images and extracting vascular tree from MR Angiography (MRA). As the first intraoperative step, after the craniotomy and with the dura opened, a designed checkerboard pattern is projected on the cortex surface and projected landmarks are scanned and captured by a stereo camera (Int J Imaging Syst Technol 23(4):294-303, 2013. doi: 10.1002/ima.22064 ). This 3D point cloud should be registered to boundary nodes of FEM in the region of interest. For this purpose, we developed a new non-rigid registration method, called finite element drift that is more compatible with the underlying nature of deformed object. The presented algorithm outperforms other methods such as coherent point drift when the deformation is local or non-coherent. After registration, the acquired displacement vectors are used as boundary conditions for FE model. As the second step, by tracking a 2D Doppler ultrasound probe swept on the parenchyma, a 3D image of deformed vascular tree is constructed. Elastic registration of this vascular point cloud to the corresponding preoperative data results the second series of displacement vector applicable to closest internal nodes of FEM. After running FE analysis, the displacement of all nodes is calculated. The brain shift is then estimated as displacement of nodes in boundary of a deep target, e.g., a tumor. We used intraoperative MR (iMR) images as the references for measuring the performance of the brain shift estimator. In the present study, two set of tests were performed using: (a) a deformable brain phantom with surface data and (b) an alive brain of an approximately big dog with surface data and US Doppler images. In our designed phantom, small tubes connected to an inflatable balloon were considered as displaceable targets and in the animal model, the target was modeled by a cyst which was created by an injection.
In the phantom study, the registration error for the surface points before FE analysis and for the target points after running FE model were <0.76 and 1.4 mm, respectively. In a real condition of operating room for animal model, the registration error was about 1 mm for the surface, 1.9 mm for the vascular tree and 1.55 mm for the target points.
The proposed projected surface imaging in conjunction with the Doppler US data combined in a powerful biomechanical model can result an acceptable performance in calculation of deformation during surgical navigation. However, the projected landmark method is sensitive to ambient light and surface conditions and the Doppler ultrasound suffers from noise and 3D image construction problems, the combination of these two methods applied on a FEM has an eligible performance.
多种术中成像方式的结合可能会减少神经外科手术期间脑移位估计的误差。在本研究中,提出了一种表面成像和多普勒超声图像的新组合,以计算皮质表面的位移和内部血管的变形,从而使用有限元模型(FEM)估计目标脑移位。评估了每个步骤中的配准误差和该方法的整体性能。
术前步骤包括从磁共振图像构建有限元模型,并从磁共振血管造影(MRA)中提取血管树。作为术中的第一步,开颅并打开硬脑膜后,将设计好的棋盘图案投影到皮质表面,并用立体相机扫描并捕获投影的地标点(《国际影像系统与技术杂志》23(4):294 - 303, 2013. doi: 10.1002/ima.22064)。这个三维点云应与感兴趣区域内有限元模型的边界节点配准。为此,我们开发了一种新的非刚性配准方法,称为有限元漂移,它与变形对象的潜在性质更兼容。当变形是局部的或非相干的时,所提出的算法优于其他方法,如相干点漂移。配准后,获取的位移向量用作有限元模型的边界条件。作为第二步,通过跟踪在脑实质内扫过的二维多普勒超声探头,构建变形血管树的三维图像。将该血管点云与相应的术前数据进行弹性配准,得到适用于有限元模型最近内部节点的第二组位移向量。运行有限元分析后,计算所有节点的位移。然后将脑移位估计为深部目标(如肿瘤)边界处节点的位移。我们使用术中磁共振(iMR)图像作为测量脑移位估计器性能的参考。在本研究中,使用以下两种情况进行了两组测试:(a)具有表面数据的可变形脑模型;(b)一只体型近似的活犬脑,具有表面数据和超声多普勒图像。在我们设计的模型中,连接到可充气气球的小管子被视为可移位目标,在动物模型中,目标由注射形成的囊肿模拟。
在模型研究中,有限元分析前表面点的配准误差和运行有限元模型后目标点的配准误差分别<0.76毫米和1.4毫米。在动物模型手术室的实际情况下,表面的配准误差约为1毫米,血管树的配准误差为1.9毫米,目标点的配准误差为1.55毫米。
所提出的投影表面成像与多普勒超声数据相结合,并结合强大的生物力学模型,在手术导航期间的变形计算中可产生可接受的性能。然而,投影地标方法对环境光和表面条件敏感,多普勒超声存在噪声和三维图像构建问题,这两种方法在有限元模型上的结合具有合格的性能。