A framework for correcting brain retraction based on an eXtended Finite Element Method using a laser range scanner

Biomechanical modelling and computer aided simulation of deep brain retraction in neurosurgery

BACKGROUND AND OBJECTIVES

Surgical simulators are widely used to promote faster and safer surgical training. They not only provide a platform for enhancing surgical skills but also minimize risks to the patient's safety, operation theatre usage, and financial expenditure. Retracting the soft brain tissue is an unavoidable procedure during any surgery to access the lesioned tissue deep within the brain. Excessive retraction often results in damaging the brain tissue, thus requiring advanced skills and prior training using virtual platforms. Such surgical simulation platforms require an anatomically correct computational model that can accurately predict the brain deformation in real-time.

METHODS

In this study, we present a 3D finite element brain model reconstructed from MRI dataset. The model incorporates precisely the anatomy and geometrical features of the canine brain. The brain model has been used to formulate and solve a quasi-static boundary value problem for brain deformation during brain retraction. The visco-hyperelastic framework within the theory of non-linear elasticity has been used to set up the boundary value problem. Consequently, the derived non-linear field equations have been solved using finite element solver ABAQUS.

RESULTS

The retraction simulations have been performed for two scenarios: retraction pressure in the brain and forces required to perform the surgery. The brain was retracted by 5 mm and retained at that position for 30 minutes, during which the retraction pressure attenuates to 36% of its peak value. Both the model predictions as well as experimental observations on canine brain indicate that brain retraction up to 30 minutes did not cause any significant risk of induced damage. Also, the peak retraction pressure level indicates that intermittent retraction is a safer procedure as compared to the continuous retraction, for the same extent of retraction.

CONCLUSIONS

The results of the present study indicate the potential of a visco-hyperelastic framework for simulating deep brain retraction effectively. The simulations were able to capture the dominant characteristics of brain tissue undergoing retraction. The developed platform could serve as a basis for the development of a detailed model in the future that can eventually be used for effective preoperative planning and training purposes.

A Surface-Based Spatial Registration Method Based on Sense Three-Dimensional Scanner

OBJECTIVE

The purpose of this study was to investigate the feasibility of a surface-based registration method based on a low-cost, hand-held Sense three-dimensional (3D) scanner in image-guided neurosurgery system.

METHODS

The scanner was calibrated prior and fixed on a tripod before registration. During registration, a part of the head surface was scanned at first and the spatial position of the adapter was recorded. Then the scanner was taken off from the tripod and the entire head surface was scanned by moving the scanner around the patient's head. All the scan points were aligned to the recorded spatial position to form a unique point cloud of the head by the automatic mosaic function of the scanner. The coordinates of the scan points were transformed from the device space to the adapter space by a calibration matrix, and then to the patient space. A 2-step patient-to-image registration method was then performed to register the patient space to the image space.

RESULTS

The experimental results showed that the mean target registration error of 15 targets on the surface of the phantom was 1.61±0.09 mm. In a clinical experiment, the mean target registration error of 7 targets on the patient's head surface was 2.50±0.31 mm, which was sufficient to meet clinical requirements.

CONCLUSIONS

It is feasible to use the Sense 3D scanner for patient-to-image registration, and the low-cost Sense 3D scanner can take the place of the current used scanner in the image-guided neurosurgery system.

A novel method for texture-mapping conoscopic surfaces for minimally invasive image-guided kidney surgery

PURPOSE

Organ-level registration is critical to image-guided therapy in soft tissue. This is especially important in organs such as the kidney which can freely move. We have developed a method for registration that combines three-dimensional locations from a holographic conoscope with an endoscopically obtained textured surface. By combining these data sources clear decisions as to the tissue from which the points arise can be made.

METHODS

By localizing the conoscope's laser dot in the endoscopic space, we register the textured surface to the cloud of conoscopic points. This allows the cloud of points to be filtered for only those arising from the kidney surface. Once a valid cloud is obtained we can use standard surface registration techniques to perform the image-space to physical-space registration. Since our methods use two distinct data sources we test for spatial accuracy and characterize temporal effects in phantoms, ex vivo porcine and human kidneys. In addition we use an industrial robot to provide controlled motion and positioning for characterizing temporal effects.

RESULTS

Our initial surface acquisitions are hand-held. This means that we take approximately 55 s to acquire a surface. At that rate we see no temporal effects due to acquisition synchronization or probe speed. Our surface registrations were able to find applied targets with submillimeter target registration errors.

CONCLUSION

The results showed that the textured surfaces could be reconstructed with submillimetric mean registration errors. While this paper focuses on kidney applications, this method could be applied to any anatomical structures where a line of sight can be created via open or minimally invasive surgical techniques.

In Vivo Investigation of the Effectiveness of a Hyper-viscoelastic Model in Simulating Brain Retraction

Intraoperative brain retraction leads to a misalignment between the intraoperative positions of the brain structures and their previous positions, as determined from preoperative images. In vitro swine brain sample uniaxial tests showed that the mechanical response of brain tissue to compression and extension could be described by the hyper-viscoelasticity theory. The brain retraction caused by the mechanical process is a combination of brain tissue compression and extension. In this paper, we first constructed a hyper-viscoelastic framework based on the extended finite element method (XFEM) to simulate intraoperative brain retraction. To explore its effectiveness, we then applied this framework to an in vivo brain retraction simulation. The simulation strictly followed the clinical scenario, in which seven swine were subjected to brain retraction. Our experimental results showed that the hyper-viscoelastic XFEM framework is capable of simulating intraoperative brain retraction and improving the navigation accuracy of an image-guided neurosurgery system (IGNS).

Marker-less tracking of brain surface deformations by non-rigid registration integrating surface and vessel/sulci features

PURPOSE

To compensate for brain shift in image-guided neurosurgery, we propose a new non-rigid registration method that integrates surface and vessel/sulci feature to noninvasively track the brain surface.

METHOD

Textured brain surfaces were acquired using phase-shift three-dimensional (3D) shape measurement, which offers 2D image pixels and their corresponding 3D points directly. Measured brain surfaces were noninvasively tracked using the proposed method by minimizing a new energy function, which is a weighted combination of 3D point corresponding estimation and surface deformation constraints. Initially, the measured surfaces were divided into featured and non-featured parts using a Frangi filter. The corresponding feature/non-feature points between intraoperative brain surfaces were estimated using the closest point algorithm. Subsequently, smoothness and rigidity constraints were introduced in the energy function for a smooth surface deformation and local surface detail conservation, respectively. Our 3D shape measurement accuracy was evaluated using 20 spheres for bias and precision errors. In addition, the proposed method was evaluated based on root mean square error (RMSE) and target registration error (TRE) with five porcine brains for which deformations were produced by gravity and pushing with different displacements in both the vertical and horizontal directions.

RESULTS

The minimum and maximum bias errors were 0.32 and 0.61 mm, respectively. The minimum and maximum precision errors were 0.025 and 0.30 mm, respectively. Quantitative validation with porcine brains showed that the average RMSE and TRE were 0.1 and 0.9 mm, respectively.

CONCLUSION

The proposed method appeared to be advantageous in integrating vessels/sulci feature, robust to changes in deformation magnitude and integrated feature numbers, and feasible in compensating for brain shift deformation in surgeries.

A new markerless patient-to-image registration method using a portable 3D scanner

PURPOSE

Patient-to-image registration is critical to providing surgeons with reliable guidance information in the application of image-guided neurosurgery systems. The conventional point-matching registration method, which is based on skin markers, requires expensive and time-consuming logistic support. Surface-matching registration with facial surface scans is an alternative method, but the registration accuracy is unstable and the error in the more posterior parts of the head is usually large because the scan range is limited. This study proposes a new surface-matching method using a portable 3D scanner to acquire a point cloud of the entire head to perform the patient-to-image registration.

METHODS

A new method for transforming the scan points from the device space into the patient space without calibration and tracking was developed. Five positioning targets were attached on a reference star, and their coordinates in the patient space were measured prior. During registration, the authors moved the scanner around the head to scan its entire surface as well as the positioning targets, and the scanner generated a unique point cloud in the device space. The coordinates of the positioning targets in the device space were automatically detected by the scanner, and a spatial transformation from the device space to the patient space could be calculated by registering them to their coordinates in the patient space that had been measured prior. A three-step registration algorithm was then used to register the patient space to the image space. The authors evaluated their method on a rigid head phantom and an elastic head phantom to verify its practicality and to calculate the target registration error (TRE) in different regions of the head phantoms. The authors also conducted an experiment with a real patient's data to test the feasibility of their method in the clinical environment.

RESULTS

In the phantom experiments, the mean fiducial registration error between the device space and the patient space, the mean surface registration error, and the mean TRE of 15 targets on the surface of each phantom were 0.34 ± 0.01 mm and 0.33 ± 0.02 mm, 1.17 ± 0.02 mm and 1.34 ± 0.10 mm, and 1.06 ± 0.11 mm and 1.48 ± 0.21 mm, respectively. When grouping the targets according to their positions on the head, high accuracy was achieved in all parts of the head, and the TREs were similar across different regions. The authors compared their method with the current surface registration methods that use only a part of the facial surface on the elastic phantom, and the mean TRE of 15 targets was 1.48 ± 0.21 mm and 1.98 ± 0.53 mm, respectively. In a clinical experiment, the mean TRE of seven targets on the patient's head surface was 1.92 ± 0.18 mm, which was sufficient to meet clinical requirements.

CONCLUSIONS

The proposed surface-matching registration method provides sufficient registration accuracy even in the posterior area of the head. The 3D point cloud of the entire head, including the facial surface and the back of the head, can be easily acquired using a portable 3D scanner. The scanner does not need to be calibrated prior or tracked by the optical tracking system during scanning.