In open abdominal image-guided liver surgery, sparse measurements of the organ surface can be taken intraoperatively via a laser-range scanning device or a tracked stylus with relatively little impact on surgical workflow. We propose a novel nonrigid registration method which uses sparse surface data to reconstruct a mapping between the preoperative CT volume and the intraoperative patient space. The mapping is generated using a tissue mechanics model subject to boundary conditions consistent with surgical supportive packing during liver resection therapy. Our approach iteratively chooses parameters which define these boundary conditions such that the deformed tissue model best fits the intraoperative surface data. Using two liver phantoms, we gathered a total of five deformation datasets with conditions comparable to open surgery. The proposed nonrigid method achieved a mean target registration error (TRE) of 3.3 mm for targets dispersed throughout the phantom volume, using a limited region of surface data to drive the nonrigid registration algorithm, while rigid registration resulted in a mean TRE of 9.5 mm. In addition, we studied the effect of surface data extent, the inclusion of subsurface data, the trade-offs of using a nonlinear tissue model, robustness to rigid misalignments, and the feasibility in five clinical datasets.

Conventional image-guided neurosurgery relies on preoperative images to provide surgical navigational information and visualization. However, these images are no longer accurate once the skull has been opened and brain shift occurs. To account for changes in the shape of the brain caused by mechanical (e.g., gravity-induced deformations) and physiological effects (e.g., hyperosmotic drug-induced shrinking, or edema-induced swelling), updated images of the brain must be provided to the neuronavigation system in a timely manner for practical use in the operating room. In this paper, a novel preoperative and intraoperative computational processing pipeline for near real-time brain shift correction in the operating room was developed to automate and simplify the processing steps. Preoperatively, a computer model of the patient’s brain with a subsequent atlas of potential deformations due to surgery is generated from diagnostic image volumes. In the case of interim gross changes between diagnosis, and surgery when reimaging is necessary, our preoperative pipeline can be generated within one day of surgery. Intraoperatively, sparse data measuring the cortical brain surface is collected using an optically tracked portable laser range scanner. These data are then used to guide an inverse modeling framework whereby full volumetric brain deformations are reconstructed from precomputed atlas solutions to rapidly match intraoperative cortical surface shift measurements. Once complete, the volumetric displacement field is used to update, i.e., deform, preoperative brain images to their intraoperative shifted state. In this paper, five surgical cases were analyzed with respect to the computational pipeline and workflow timing. With respect to postcortical surface data acquisition, the approximate execution time was 4.5 min. The total update process which included positioning the scanner, data acquisition, inverse model processing, and image deforming was ∼ 11–13 min. In addition, easily implemented hardware, software, and workflow processes were identified for improved performance in the near future.

Purpose Brain shift during neurosurgical procedures must be corrected for in order to reestablish accurate alignment for successful image-guided tumor resection. Sparse-data-driven biomechanical models that predict physiological brain shift by accounting for typical deformation-inducing events such as cerebrospinal fluid drainage, hyperosmotic drugs, swelling, retraction, resection, and tumor cavity collapse are an inexpensive solution. This study evaluated the robustness and accuracy of a biomechanical model-based brain shift correction system to assist with tumor resection surgery in 16 clinical cases. Methods Preoperative computation involved the generation of a patient-specific finite element model of the brain and creation of an atlas of brain deformation solutions calculated using a distribution of boundary and deformation-inducing forcing conditions (e.g., sag, tissue contraction, and tissue swelling). The optimum brain shift solution was determined using an inverse problem approach which linearly combines solutions from the atlas to match the cortical surface deformation data collected intraoperatively. The computed deformations were then used to update the preoperative images for all 16 patients. Results The mean brain shift measured ranged on average from 2.5 to 21.3 mm, and the biomechanical model-based correction system managed to account for the bulk of the brain shift, producing a mean corrected error ranging on average from 0.7 to 4.0 mm. Conclusions Biomechanical models are an inexpensive means to assist intervention via correction for brain deformations that can compromise surgical navigation systems. To our knowledge, this study represents the most comprehensive clinical evaluation of a deformation correction pipeline for image-guided neurosurgery.

Soft-tissue image-guided interventions often require the digitization of organ surfaces for providing correspondence from medical images to the physical patient in the operating room. In this paper, the effect of several inexpensive surface acquisition techniques on target registration error and surface registration error (SRE) for soft tissue is investigated. A systematic approach is provided to compare image-to-physical registrations using three different methods of organ spatial digitization: 1) a tracked laser-range scanner (LRS), 2) a tracked pointer, and 3) a tracked conoscopic holography sensor (called a conoprobe). For each digitization method, surfaces of phantoms and biological tissues were acquired and registered to CT image volume counterparts. A comparison among these alignments demonstrated that registration errors were statistically smaller with the conoprobe than the tracked pointer and LRS (p < 0.01). In all acquisitions, the conoprobe outperformed the LRS and tracked pointer: for example, the arithmetic means of the SRE over all data acquisitions with a porcine liver were 1.73 ± 0.77 mm, 3.25 ± 0.78 mm, and 4.44 ± 1.19 mm for the conoprobe, LRS, and tracked pointer, respectively. In a cadaveric kidney specimen, the arithmetic means of the SRE over all trials of the conoprobe and tracked pointer were 1.50 ± 0.50 mm and 3.51 ± 0.82 mm, respectively. Our results suggest that tissue displacements due to contact force and attempts to maintain contact with tissue, compromise registrations that are dependent on data acquired from a tracked surgical instrument and we provide an alternative method (tracked conoscopic holography) of digitizing surfaces for clinical usage. The tracked conoscopic holography device outperforms LRS acquisitions with respect to registration accuracy.

This article presents a method designed to automatically track cortical vessels in intra-operative microscope video sequences. The main application of this method is the estimation of cortical displacement that occurs during tumor resection procedures. The method works in three steps. First, models of vessels selected in the first frame of the sequence are built. These models are then used to track vessels across frames in the video sequence. Finally, displacements estimated using the vessels are extrapolated to the entire image. The method has been tested retrospectively on images simulating large displacement, tumor resection, and partial occlusion by surgical instruments and on 21 video sequences comprising several thousand frames acquired from three patients. Qualitative results show that the method is accurate, robust to the appearance and disappearance of surgical instruments, and capable of dealing with large differences in images caused by resection. Quantitative results show a mean vessel tracking error (VTE) of 2.4 pixels (0.3 or 0.6 mm, depending on the spatial resolution of the images) and an average target registration error (TRE) of 3.3 pixels (0.4 or 0.8 mm).