A computational model for tracking subsurface tissue deformation during stereotactic neurosurgery

Model-based image updating in deep brain stimulation with assimilation of deep brain sparse data

BACKGROUND

Accuracy of electrode placement for deep brain stimulation (DBS) is critical to achieving desired surgical outcomes and impacts the efficacy of treating neurodegenerative diseases. Intraoperative brain shift degrades the accuracy of surgical navigation based on preoperative images.

PURPOSE

We extended a model-based image updating scheme to address intraoperative brain shift in DBS surgery and improved its accuracy in deep brain.

METHODS

We evaluated 10 patients, retrospectively, who underwent bilateral DBS surgery and classified them into groups of large and small deformation based on a 2 mm subsurface movement threshold and brain shift index of 5%. In each case, sparse brain deformation data were used to estimate whole brain displacements and deform preoperative CT (preCT) to generate updated CT (uCT). Accuracy of uCT was assessed using target registration errors (TREs) at the Anterior Commissure (AC), Posterior Commissure (PC), and four calcification points in the sub-ventricular area by comparing their locations in uCT with their ground truth counterparts in postoperative CT (postCT).

RESULTS

In the large deformation group, TREs were reduced from 2.5 mm in preCT to 1.2 mm in uCT (53% compensation); in the small deformation group, errors were reduced from 1.25 to 0.74 mm (41%). Average reduction of TREs at AC, PC and pineal gland were significant, statistically (p ⩽ 0.01).

CONCLUSIONS

With more rigorous validation of model results, this study confirms the feasibility of improving the accuracy of model-based image updating in compensating for intraoperative brain shift during DBS procedures by assimilating deep brain sparse data.

In vivo measurement of human brain material properties under quasi-static loading

Computational modelling of the brain requires accurate representation of the tissues concerned. Mechanical testing has numerous challenges, in particular for low strain rates, like neurosurgery, where redistribution of fluid is biomechanically important. A finite-element (FE) model was generated in FEBio, incorporating a spring element/fluid-structure interaction representation of the pia-arachnoid complex (PAC). The model was loaded to represent gravity in prone and supine positions. Material parameter identification and sensitivity analysis were performed using statistical software, comparing the FE results to human in vivo measurements. Results for the brain Ogden parameters µ, α and k yielded values of 670 Pa, -19 and 148 kPa, supporting values reported in the literature. Values of the order of 1.2 MPa and 7.7 kPa were obtained for stiffness of the pia mater and out-of-plane tensile stiffness of the PAC, respectively. Positional brain shift was found to be non-rigid and largely driven by redistribution of fluid within the tissue. To the best of our knowledge, this is the first study using in vivo human data and gravitational loading in order to estimate the material properties of intracranial tissues. This model could now be applied to reduce the impact of positional brain shift in stereotactic neurosurgery.

A Minireview on Brain Models Simulating Geometrical, Physical, and Biochemical Properties of the Human Brain

There is a growing body of evidences that brain surrogates will be of great interest for researchers and physicians in the medical field. They are currently mainly used for education and training purposes or to verify the appropriate functionality of medical devices. Depending on the purpose, a variety of materials have been used with specific and accurate mechanical and biophysical properties, More recently they have been used to assess the biocompatibility of implantable devices, but they are still not validated to study the migration of leaching components from devices. This minireview shows the large diversity of approaches and uses of brain phantoms, which converge punctually. All these phantoms are complementary to numeric models, which benefit, reciprocally, of their respective advances. It also suggests avenues of research for the analysis of leaching components from implantable devices.

Viscoelastic biomechanical models to predict inward brain-shift using public benchmark data

Brain-shift during neurosurgery compromises the accuracy of tracking the boundaries of the tumor to be resected. Although several studies have used various finite element models (FEMs) to predict inward brain-shift, evaluation of their accuracy and efficiency based on public benchmark data has been limited. This study evaluates several FEMs proposed in the literature (various boundary conditions, mesh sizes, and material properties) by using intraoperative imaging data (the public REtroSpective Evaluation of Cerebral Tumors [RESECT] database). Four patients with low-grade gliomas were identified as having inward brain-shifts. We computed the accuracy (using target registration error) of several FEM-based brain-shift predictions and compared our findings. Since information on head orientation during craniotomy is not included in this database, we tested various plausible angles of head rotation. We analyzed the effects of brain tissue viscoelastic properties, mesh size, craniotomy position, CSF drainage level, and rigidity of meninges and then quantitatively evaluated the trade-off between accuracy and central processing unit time in predicting inward brain-shift across all models with second-order tetrahedral FEMs. The mean initial target registration error (TRE) was 5.78 ± 3.78 mm with rigid registration. FEM prediction (edge-length, 5 mm) with non-rigid meninges led to a mean TRE correction of 1.84 ± 0.83 mm assuming heterogeneous material. Results show that, for the low-grade glioma patients in the study, including non-rigid modeling of the meninges was significant statistically. In contrast including heterogeneity was not significant. To estimate the optimal head orientation and CSF drainage, an angle step of 5° and an CSF height step of 5 mm were enough leading to <0.26 mm TRE fluctuation.

A hybrid, image-based and biomechanics-based registration approach to markerless intraoperative nodule localization during video-assisted thoracoscopic surgery

The resection of small, low-dense or deep lung nodules during video-assisted thoracoscopic surgery (VATS) is surgically challenging. Nodule localization methods in clinical practice typically rely on the preoperative placement of markers, which may lead to clinical complications. We propose a markerless lung nodule localization framework for VATS based on a hybrid method combining intraoperative cone-beam CT (CBCT) imaging, free-form deformation image registration, and a poroelastic lung model with allowance for air evacuation. The difficult problem of estimating intraoperative lung deformations is decomposed into two more tractable sub-problems: (i) estimating the deformation due the change of patient pose from preoperative CT (supine) to intraoperative CBCT (lateral decubitus); and (ii) estimating the pneumothorax deformation, i.e. a collapse of the lung within the thoracic cage. We were able to demonstrate the feasibility of our localization framework with a retrospective validation study on 5 VATS clinical cases. Average initial errors in the range of 22 to 38 mm were reduced to the range of 4 to 14 mm, corresponding to an error correction in the range of 63 to 85%. To our knowledge, this is the first markerless lung deformation compensation method dedicated to VATS and validated on actual clinical data.

Model-Based Image Updating for Brain Shift in Deep Brain Stimulation Electrode Placement Surgery

OBJECTIVE

The efficacy of deep brain stimulation (DBS) depends on accurate placement of electrodes. Although stereotactic frames enable co-registration of image-based surgical planning and the operative field, the accuracy of electrode placement can be degraded by intra-operative brain shift. In this study, we adapted a biomechanical model to estimate whole brain displacements from which we deformed preoperative CT (preCT) to generate an updated CT (uCT) that compensates for brain shift.

METHODS

We drove the deformation model using displacement data derived from deformation in the frontal cortical surface that occurred during the DBS intervention. We evaluated 15 patients, retrospectively, who underwent bilateral DBS surgery, and assessed the accuracy of uCT in terms of target registration error (TRE) relative to a CT acquired post-placement (postCT). We further divided subjects into large (Group L) and small (Group S) deformation groups based on a TRE threshold of 1.6mm. Anterior commissure (AC), posterior commissure (PC) and pineal gland (PG) were identified on preCT and postCT and used to quantify TREs in preCT and uCT.

RESULTS

In the group of large brain deformation, average TREs for uCT were 1.11 ± 0.13 and 1.07 ± 0.38 mm at AC and PC, respectively, compared to 1.85 ± 0.17 and 0.92 ± 0.52 mm for preCT. The model updating approach improved AC localization but did not alter TREs at PC.

CONCLUSION

This preliminary study suggests that our image updating method may compensate for brain shift around surgical targets of importance during DBS surgery, although further investigation is warranted before conclusive evidence will be available.

SIGNIFICANCE

With further development and evaluation, our model-based image updating method using intraoperative sparse data may compensate for brain shift in DBS surgery efficiently, and have utility in updating targeting coordinates.

Image Updating for Brain Shift Compensation During Resection

BACKGROUND

In open-cranial neurosurgery, preoperative magnetic resonance (pMR) images are typically coregistered for intraoperative guidance. Their accuracy can be significantly degraded by intraoperative brain deformation, especially when resection is involved.

OBJECTIVE

To produce model updated MR (uMR) images to compensate for brain shift that occurred during resection, and evaluate the performance of the image-updating process in terms of accuracy and computational efficiency.

METHODS

In 14 resection cases, intraoperative stereovision image pairs were acquired after dural opening and during resection to generate displacement maps of the surgical field. These data were assimilated by a biomechanical model to create uMR volumes of the evolving surgical field. A tracked stylus provided independent measurements of feature locations to quantify target registration errors (TREs) in the original coregistered pMR and uMR as surgery progressed.

RESULTS

Updated MR TREs were 1.66 ± 0.27 and 1.92 ± 0.49 mm in the 14 cases after dural opening and after partial resection, respectively, compared to 8.48 ± 3.74 and 8.77 ± 4.61 mm for pMR, respectively. The overall computational time for generating uMRs after partial resection was less than 10 min.

CONCLUSION

We have developed an image-updating system to compensate for brain deformation during resection using a computational model with data assimilation of displacements measured with intraoperative stereovision imaging that maintains TREs less than 2 mm on average.

In vivo modeling of interstitial pressure in a porcine model: approximation of poroelastic properties and effects of enhanced anatomical structure modeling

The purpose of this investigation is to test whether a poroelastic model with enhanced structure can capture in vivo interstitial pressure dynamics in a brain undergoing mock surgical loads. Using interstitial pressure data from a porcine study, we use an inverse model to reconstruct material properties in an effort to capture these in vivo brain tissue dynamics. Four distinct models for the reconstruction of parameters are investigated (full anatomical condition description, condition without dural septa description, condition without ventricle boundary description, and the conventional fully saturated model). These models are systematic in their development to isolate the influence of three model characteristics: the dural septa, the treatment of the ventricles, and the treatment of the brain as a saturated media. This study demonstrates that to capture appropriate pressure compartmentalization, interstitial pressure gradients, pressure transient effects, and deformations within the brain, the proposed boundary conditions and structural enhancement coupled with a heterogeneous description invoking partial saturation are needed in a biphasic poroelastic model. These findings suggest that with enhanced anatomical modeling and appropriate model assumptions, poroelastic models can be used to capture quite complex brain deformations and interstitial pressure dynamics.

Biphasic analysis of rat brain slices under creep indentation shows nonlinear tension-compression behavior

Biphasic theory can provide a mechanistic description of deformation and transport phenomena in soft tissues, and has been used to model surgery and drug delivery in the brain for decades. Knowledge of corresponding mechanical properties of the brain is needed to accurately predict tissue deformation and flow transport in these applications. Previously in our group, creep indentation tests were conducted for multiple anatomical regions in acute rat brain tissue slices. In the current study, a biphasic finite element model of creep indentation was developed with which to compare these data. Considering the soft tissue structure of brain, the solid matrix was assumed to be composed of a neo-Hookean ground matrix reinforced by continuously distributed fibers that exhibits tension-compression nonlinearity during deformation. By fixing Poisson's ratio of the ground matrix, Young's modulus, fiber modulus and hydraulic permeability were estimated. Hydraulic permeability was found to be nearly independent of the properties of the solid matrix. Estimated modulus (40 Pa to 1.1 kPa for the ground matrix, 3.2-18.2 kPa for fibers) and hydraulic permeability (1.2-5.5×10^-13m⁴/N s) fell within an acceptable range compared with those in previous studies. Instantaneous indentation depth was dominated by tension provided by fibers, while the tissue response at equilibrium was sensitive to Poisson's ratio. Results of sensitivity analysis also point to the necessity of considering tension-compression nonlinearity in the solid phase when the biphasic material undergoes large creep deformation.

Retrospective study comparing model-based deformation correction to intraoperative magnetic resonance imaging for image-guided neurosurgery

Brain shift during tumor resection compromises the spatial validity of registered preoperative imaging data that is critical to image-guided procedures. One current clinical solution to mitigate the effects is to reimage using intraoperative magnetic resonance (iMR) imaging. Although iMR has demonstrated benefits in accounting for preoperative-to-intraoperative tissue changes, its cost and encumbrance have limited its widespread adoption. While iMR will likely continue to be employed for challenging cases, a cost-effective model-based brain shift compensation strategy is desirable as a complementary technology for standard resections. We performed a retrospective study of [Formula: see text] tumor resection cases, comparing iMR measurements with intraoperative brain shift compensation predicted by our model-based strategy, driven by sparse intraoperative cortical surface data. For quantitative assessment, homologous subsurface targets near the tumors were selected on preoperative MR and iMR images. Once rigidly registered, intraoperative shift measurements were determined and subsequently compared to model-predicted counterparts as estimated by the brain shift correction framework. When considering moderate and high shift ([Formula: see text], [Formula: see text] measurements per case), the alignment error due to brain shift reduced from [Formula: see text] to [Formula: see text], representing [Formula: see text] correction. These first steps toward validation are promising for model-based strategies.

A numerical framework for interstitial fluid pressure imaging in poroelastic MRE

A numerical framework for interstitial fluid pressure imaging (IFPI) in biphasic materials is investigated based on three-dimensional nonlinear finite element poroelastic inversion. The objective is to reconstruct the time-harmonic pore-pressure field from tissue excitation in addition to the elastic parameters commonly associated with magnetic resonance elastography (MRE). The unknown pressure boundary conditions (PBCs) are estimated using the available full-volume displacement data from MRE. A subzone-based nonlinear inversion (NLI) technique is then used to update mechanical and hydrodynamical properties, given the appropriate subzone PBCs, by solving a pressure forward problem (PFP). The algorithm was evaluated on a single-inclusion phantom in which the elastic property and hydraulic conductivity images were recovered. Pressure field and material property estimates had spatial distributions reflecting their true counterparts in the phantom geometry with RMS errors around 20% for cases with 5% noise, but degraded significantly in both spatial distribution and property values for noise levels > 10%. When both shear moduli and hydraulic conductivity were estimated along with the pressure field, property value error rates were as high as 58%, 85% and 32% for the three quantities, respectively, and their spatial distributions were more distorted. Opportunities for improving the algorithm are discussed.

A systematic evaluation of intraoperative white matter tract shift in pediatric epilepsy surgery using high-field MRI and probabilistic high angular resolution diffusion imaging tractography

OBJECTIVE Characterization of intraoperative white matter tract (WMT) shift has the potential to compensate for neuronavigation inaccuracies using preoperative brain imaging. This study aimed to quantify and characterize intraoperative WMT shift from the global hemispheric to the regional tract-based scale and to investigate the impact of intraoperative factors (IOFs). METHODS High angular resolution diffusion imaging (HARDI) diffusion-weighted data were acquired over 5 consecutive perioperative time points (MR₁ to MR₅) in 16 epilepsy patients (8 male; mean age 9.8 years, range 3.8-15.8 years) using diagnostic and intraoperative 3-T MRI scanners. MR₁ was the preoperative planning scan. MR₂ was the first intraoperative scan acquired with the patient's head fixed in the surgical position. MR₃ was the second intraoperative scan acquired following craniotomy and durotomy, prior to lesion resection. MR₄ was the last intraoperative scan acquired following lesion resection, prior to wound closure. MR₅ was a postoperative scan acquired at the 3-month follow-up visit. Ten association WMT/WMT segments and 1 projection WMT were generated via a probabilistic tractography algorithm from each MRI scan. Image registration was performed through pairwise MRI alignments using the skull segmentation. The MR₁ and MR₂ pairing represented the first surgical stage. The MR₂ and MR₃ pairing represented the second surgical stage. The MR₃ and MR₄ (or MR₅) pairing represented the third surgical stage. The WMT shift was quantified by measuring displacements between a pair of WMT centerlines. Linear mixed-effects regression analyses were carried out for 6 IOFs: head rotation, craniotomy size, durotomy size, resected lesion volume, presence of brain edema, and CSF loss via ventricular penetration. RESULTS The average WMT shift in the operative hemisphere was 2.37 mm (range 1.92-3.03 mm) during the first surgical stage, 2.19 mm (range 1.90-3.65 mm) during the second surgical stage, and 2.92 mm (range 2.19-4.32 mm) during the third surgical stage. Greater WMT shift occurred in the operative than the nonoperative hemisphere, in the WMTs adjacent to the surgical lesion rather than those remote to it, and in the superficial rather than the deep segment of the pyramidal tract. Durotomy size and resection size were significant, independent IOFs affecting WMT shift. The presence of brain edema was a marginally significant IOF. Craniotomy size, degree of head rotation, and ventricular penetration were not significant IOFs affecting WMT shift. CONCLUSIONS WMT shift occurs noticeably in tracts adjacent to the surgical lesions, and those motor tracts superficially placed in the operative hemisphere. Intraoperative probabilistic HARDI tractography following craniotomy, durotomy, and lesion resection may compensate for intraoperative WMT shift and improve neuronavigation accuracy.

Functional Magnetic Resonance Imaging (fMRI), Pre-intraoperative Tractography in Neurosurgery: The Experience of Sant' Andrea Rome University Hospital

BACKGROUND

The goal of neurosurgery for cerebral intraparenchymal neoplasms of the eloquent areas is maximal resection with the preservation of normal functions, and minimizing operative risk and postoperative morbidity. Currently, modern technological advances in neuroradiological tools, neuronavigation, and intraoperative magnetic resonance imaging (MRI) have produced great improvements in postoperative morbidity after the surgery of cerebral eloquent areas. The integration of preoperative functional MRI (fMRI), intraoperative MRI (volumetric and diffusion tensor imaging [DTI]), and neuronavigation, defined as "functional neuronavigation" has improved the intraoperative detection of the eloquent areas.

METHODS

We reviewed 142 patients operated between 2004 and 2010 for intraparenchymal neoplasms involving or close to one or more major white matter tracts (corticospinal tract [CST], arcuate fasciculus [AF], optic radiation). All the patients underwent neurosurgery in a BrainSUITE equipped with a 1.5 T MR scanner and were preoperatively studied with fMRI and DTI for tractography for surgical planning. The patients underwent MRI and DTI during surgery after dural opening, after the gross total resection close to the white matter tracts, and at the end of the procedure. We evaluated the impact of fMRI on surgical planning and on the selection of the entry point on the cortical surface. We also evaluated the impact of preoperative and intraoperative DTI, in order to modify the surgical approach, to define the borders of resection, and to correlate this modality with subcortical neurophysiological monitoring. We evaluated the impact of the preoperative fMRI by intraoperative neurophysiological monitoring, performing "neuronavigational" brain mapping, following its data to localize the previously elicited areas after brain shift correction by intraoperative MRI.

RESULTS

The mean age of the 142 patients (89 M/53 F) was 59.1 years and the lesion involved the CST in 66 patients (57 %), the language pathways in 24 (21 %), and the optic radiations in 25 (22 %). The integration of tractographic data into the volumetric dataset for neuronavigation was technically possible in all cases. In all patients intraoperative DTI demonstrated a shift of the bundle position caused by the surgical procedure; its dislocation was both outward and inward in the range of +6 mm and -2 mm.

CONCLUSION

We found a high concordance between fMRI/DTI and intraoperative brain mapping; their combination improves the sensitivity of each technique, reducing pitfalls and so defining "functional neuronavigation", increasing the definition of eloquent areas and also reducing the time of surgery.

Intraoperative image updating for brain shift following dural opening

OBJECTIVE Preoperative magnetic resonance images (pMR) are typically coregistered to provide intraoperative navigation, the accuracy of which can be significantly compromised by brain deformation. In this study, the authors generated updated MR images (uMR) in the operating room (OR) to compensate for brain shift due to dural opening, and evaluated the accuracy and computational efficiency of the process. METHODS In 20 open cranial neurosurgical cases, a pair of intraoperative stereovision (iSV) images was acquired after dural opening to reconstruct a 3D profile of the exposed cortical surface. The iSV surface was registered with pMR to detect cortical displacements that were assimilated by a biomechanical model to estimate whole-brain nonrigid deformation and produce uMR in the OR. The uMR views were displayed on a commercial navigation system and compared side by side with the corresponding coregistered pMR. A tracked stylus was used to acquire coordinate locations of features on the cortical surface that served as independent positions for calculating target registration errors (TREs) for the coregistered uMR and pMR image volumes. RESULTS The uMR views were visually more accurate and well aligned with the iSV surface in terms of both geometry and texture compared with pMR where misalignment was evident. The average misfit between model estimates and measured displacements was 1.80 ± 0.35 mm, compared with the average initial misfit of 7.10 ± 2.78 mm between iSV and pMR, and the average TRE was 1.60 ± 0.43 mm across the 20 patients in the uMR image volume, compared with 7.31 ± 2.82 mm on average in the pMR cases. The iSV also proved to be accurate with an average error of 1.20 ± 0.37 mm. The overall computational time required to generate the uMR views was 7-8 minutes. CONCLUSIONS This study compensated for brain deformation caused by intraoperative dural opening using computational model-based assimilation of iSV cortical surface displacements. The uMR proved to be more accurate in terms of model-data misfit and TRE in the 20 patient cases evaluated relative to pMR. The computational time was acceptable (7-8 minutes) and the process caused minimal interruption of surgical workflow.

Brain shift in neuronavigation of brain tumors: A review

PURPOSE

Neuronavigation based on preoperative imaging data is a ubiquitous tool for image guidance in neurosurgery. However, it is rendered unreliable when brain shift invalidates the patient-to-image registration. Many investigators have tried to explain, quantify, and compensate for this phenomenon to allow extended use of neuronavigation systems for the duration of surgery. The purpose of this paper is to present an overview of the work that has been done investigating brain shift.

METHODS

A review of the literature dealing with the explanation, quantification and compensation of brain shift is presented. The review is based on a systematic search using relevant keywords and phrases in PubMed. The review is organized based on a developed taxonomy that classifies brain shift as occurring due to physical, surgical or biological factors.

RESULTS

This paper gives an overview of the work investigating, quantifying, and compensating for brain shift in neuronavigation while describing the successes, setbacks, and additional needs in the field. An analysis of the literature demonstrates a high variability in the methods used to quantify brain shift as well as a wide range in the measured magnitude of the brain shift, depending on the specifics of the intervention. The analysis indicates the need for additional research to be done in quantifying independent effects of brain shift in order for some of the state of the art compensation methods to become useful.

CONCLUSION

This review allows for a thorough understanding of the work investigating brain shift and introduces the needs for future avenues of investigation of the phenomenon.

Learning intervention-induced deformations for non-rigid MR-CT registration and electrode localization in epilepsy patients

This paper describes a framework for learning a statistical model of non-rigid deformations induced by interventional procedures. We make use of this learned model to perform constrained non-rigid registration of pre-procedural and post-procedural imaging. We demonstrate results applying this framework to non-rigidly register post-surgical computed tomography (CT) brain images to pre-surgical magnetic resonance images (MRIs) of epilepsy patients who had intra-cranial electroencephalography electrodes surgically implanted. Deformations caused by this surgical procedure, imaging artifacts caused by the electrodes, and the use of multi-modal imaging data make non-rigid registration challenging. Our results show that the use of our proposed framework to constrain the non-rigid registration process results in significantly improved and more robust registration performance compared to using standard rigid and non-rigid registration methods.

Clinical evaluation of a model-updated image-guidance approach to brain shift compensation: experience in 16 cases

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.

Suitability of poroelastic and viscoelastic mechanical models for high and low frequency MR elastography

PURPOSE

Descriptions of the structure of brain tissue as a porous cellular matrix support application of a poroelastic (PE) mechanical model which includes both solid and fluid phases. However, the majority of brain magnetic resonance elastography (MRE) studies use a single phase viscoelastic (VE) model to describe brain tissue behavior, in part due to availability of relatively simple direct inversion strategies for mechanical property estimation. A notable exception is low frequency intrinsic actuation MRE, where PE mechanical properties are imaged with a nonlinear inversion algorithm.

METHODS

This paper investigates the effect of model choice at each end of the spectrum of in vivo human brain actuation frequencies. Repeat MRE examinations of the brains of healthy volunteers were used to compare image quality and repeatability for each inversion model for both 50 Hz externally produced motion and ≈1 Hz intrinsic motions. Additionally, realistic simulated MRE data were generated with both VE and PE finite element solvers to investigate the effect of inappropriate model choice for ideal VE and PE materials.

RESULTS

In vivo, MRE data revealed that VE inversions appear more representative of anatomical structure and quantitatively repeatable for 50 Hz induced motions, whereas PE inversion produces better results at 1 Hz. Reasonable VE approximations of PE materials can be derived by equating the equivalent wave velocities for the two models, provided that the timescale of fluid equilibration is not similar to the period of actuation. An approximation of the equilibration time for human brain reveals that this condition is violated at 1 Hz but not at 50 Hz. Additionally, simulation experiments when using the "wrong" model for the inversion demonstrated reasonable shear modulus reconstructions at 50 Hz, whereas cross-model inversions at 1 Hz were poor quality. Attenuation parameters were sensitive to changes in the forward model at both frequencies, however, no spatial information was recovered because the mechanisms of VE and PE attenuation are different.

CONCLUSIONS

VE inversions are simpler with fewer unknown properties and may be sufficient to capture the mechanical behavior of PE brain tissue at higher actuation frequencies. However, accurate modeling of the fluid phase is required to produce useful mechanical property images at the lower frequencies of intrinsic brain motions.

Computational Modeling for Enhancing Soft Tissue Image Guided Surgery: An Application in Neurosurgery

With the recent advances in computing, the opportunities to translate computational models to more integrated roles in patient treatment are expanding at an exciting rate. One area of considerable development has been directed towards correcting soft tissue deformation within image guided neurosurgery applications. This review captures the efforts that have been undertaken towards enhancing neuronavigation by the integration of soft tissue biomechanical models, imaging and sensing technologies, and algorithmic developments. In addition, the review speaks to the evolving role of modeling frameworks within surgery and concludes with some future directions beyond neurosurgical applications.

Intraoperative fiducial-less patient registration using volumetric 3D ultrasound: a prospective series of 32 neurosurgical cases

OBJECT

Fiducial-based registration (FBR) is used widely for patient registration in image-guided neurosurgery. The authors of this study have developed an automatic fiducial-less registration (FLR) technique to find the patient-to-image transformation by directly registering 3D ultrasound (3DUS) with MR images without incorporating prior information. The purpose of the study was to evaluate the performance of the FLR technique when used prospectively in the operating room and to compare it with conventional FBR.

METHODS

In 32 surgical patients who underwent conventional FBR, preoperative T1-weighted MR images (pMR) with attached fiducial markers were acquired prior to surgery. After craniotomy but before dural opening, a set of 3DUS images of the brain volume was acquired. A 2-step registration process was executed immediately after image acquisition: 1) the cortical surfaces from pMR and 3DUS were segmented, and a multistart sum-of-squared-intensity-difference registration was executed to find an initial alignment between down-sampled binary pMR and 3DUS volumes; and 2) the alignment was further refined by a mutual information-based registration between full-resolution grayscale pMR and 3DUS images, and a patient-to-image transformation was subsequently extracted.

RESULTS

To assess the accuracy of the FLR technique, the following were quantified: 1) the fiducial distance error (FDE); and 2) the target registration error (TRE) at anterior commissure and posterior commissure locations; these were compared with conventional FBR. The results showed that although the average FDE (6.42 ± 2.05 mm) was higher than the fiducial registration error (FRE) from FBR (3.42 ± 1.37 mm), the overall TRE of FLR (2.51 ± 0.93 mm) was lower than that of FBR (5.48 ± 1.81 mm). The results agreed with the intent of the 2 registration techniques: FBR is designed to minimize the FRE, whereas FLR is designed to optimize feature alignment and hence minimize TRE. The overall computational cost of FLR was approximately 4-5 minutes and minimal user interaction was required.

CONCLUSIONS

Because the FLR method directly registers 3DUS with MR by matching internal image features, it proved to be more accurate than FBR in terms of TRE in the 32 patients evaluated in this study. The overall efficiency of FLR in terms of the time and personnel involved is also improved relative to FBR in the operating room, and the method does not require additional image scans immediately prior to surgery. The performance of FLR and these results suggest potential for broad clinical application.

Stereovision to MR image registration for cortical surface displacement mapping to enhance image-guided neurosurgery

PURPOSE

A surface registration method is presented to align intraoperative stereovision (iSV) with preoperative magnetic resonance (pMR) images, which utilizes both geometry and texture information to extract tissue displacements as part of the overall process of compensating for intraoperative brain deformation in order to maintain accurate neuronavigational image guidance during surgery.

METHODS

A sum-of-squared-difference rigid image registration was first executed to detect lateral shift of the cortical surface and was followed by a mutual-information-based block matching method to detect local nonrigid deformation caused by distention or collapse of the cortical surface. Ten (N = 10) surgical cases were evaluated in which an independent point measurement of a dominant cortical surface feature location was recorded with a tracked stylus in each case and compared to its surface-registered counterpart. The full three-dimensional (3D) displacement field was also extracted to drive a biomechanical brain deformation model, the results of which were reconciled with the reconstructed iSV surface as another form of evaluation.

RESULTS

Differences between the tracked stylus coordinates of cortical surface features and their surface-registered locations were 1.94 ± 0.59 mm on average across the ten cases. When the complete displacement map derived from surface registration was utilized, the resulting images generated from mechanical model updates were consistent in terms of both geometry (1-2 mm of model misfit) and texture, and were generated with less than 10 min of computational time. Analysis of the surface-registered 3D displacements indicate that the magnitude of motion ranged from 4.03 to 9.79 mm in the ten patient cases, and the amount of lateral shift was not related statistically to the direction of gravity (p = 0.73 ≫ 0.05) or the craniotomy size (p = 0.48 ≫ 0.05) at the beginning of surgery.

CONCLUSIONS

The iSV-pMR surface registration method utilizes texture and geometry information to extract both global lateral shift and local nonrigid movement of the cortical surface in 3D. The results suggest small differences exist in surface-registered locations when compared to positions measured independently with a coregistered stylus and when the full iSV surface was aligned with model-updated MR. The effectiveness and efficiency of the registration method is also minimally disruptive to surgical workflow.

Parametric-based brain Magnetic Resonance Elastography using a Rayleigh damping material model

The three-parameter Rayleigh damping (RD) model applied to time-harmonic Magnetic Resonance Elastography (MRE) has potential to better characterise fluid-saturated tissue systems. However, it is not uniquely identifiable at a single frequency. One solution to this problem involves simultaneous inverse problem solution of multiple input frequencies over a broad range. As data is often limited, an alternative elegant solution is a parametric RD reconstruction, where one of the RD parameters (μI or ρI) is globally constrained allowing accurate identification of the remaining two RD parameters. This research examines this parametric inversion approach as applied to in vivo brain imaging. Overall, success was achieved in reconstruction of the real shear modulus (μR) that showed good correlation with brain anatomical structures. The mean and standard deviation shear stiffness values of the white and gray matter were found to be 3±0.11kPa and 2.2±0.11kPa, respectively, which are in good agreement with values established in the literature or measured by mechanical testing. Parametric results with globally constrained μI indicate that selecting a reasonable value for the μI distribution has a major effect on the reconstructed ρI image and concomitant damping ratio (ξd). More specifically, the reconstructed ρI image using a realistic μI=333Pa value representative of a greater portion of the brain tissue showed more accurate differentiation of the ventricles within the intracranial matter compared to μI=1000Pa, and ξd reconstruction with μI=333Pa accurately captured the higher damping levels expected within the vicinity of the ventricles. Parametric RD reconstruction shows potential for accurate recovery of the stiffness characteristics and overall damping profile of the in vivo living brain despite its underlying limitations. Hence, a parametric approach could be valuable with RD models for diagnostic MRE imaging with single frequency data.

A dynamic mechanical analysis technique for porous media

Dynamic mechanical analysis (DMA) is a common way to measure the mechanical properties of materials as functions of frequency. Traditionally, a viscoelastic mechanical model is applied and current DMA techniques fit an analytical approximation to measured dynamic motion data by neglecting inertial forces and adding empirical correction factors to account for transverse boundary displacements. Here, a finite-element (FE) approach to processing DMA data was developed to estimate poroelastic material properties. Frequency-dependent inertial forces, which are significant in soft media and often neglected in DMA, were included in the FE model. The technique applies a constitutive relation to the DMA measurements and exploits a nonlinear inversion to estimate the material properties in the model that best fit the model response to the DMA data. A viscoelastic version of this approach was developed to validate the approach by comparing complex modulus estimates to the direct DMA results. Both analytical and FE poroelastic models were also developed to explore their behavior in the DMA testing environment. All of the models were applied to tofu as a representative soft poroelastic material that is a common phantom in elastography imaging studies. Five samples of three different stiffnesses were tested from 1-14 Hz with rough platens placed on the top and bottom surfaces of the material specimen under test to restrict transverse displacements and promote fluid-solid interaction. The viscoelastic models were identical in the static case, and nearly the same at frequency with inertial forces accounting for some of the discrepancy. The poroelastic analytical method was not sufficient when the relevant physical boundary constraints were applied, whereas the poroelastic FE approach produced high quality estimates of shear modulus and hydraulic conductivity. These results illustrated appropriate shear modulus contrast between tofu samples and yielded a consistent contrast in hydraulic conductivity as well.

Brain shift during bur hole–based procedures using interventional MRI

Object

Brain shift during minimally invasive, bur hole–based procedures such as deep brain stimulation (DBS) electrode implantation and stereotactic brain biopsy is not well characterized or understood. We examine shift in various regions of the brain during a novel paradigm of DBS electrode implantation using interventional imaging throughout the procedure with high-field interventional MRI.

Methods

Serial MR images were obtained and analyzed using a 1.5-T magnet prior to, during, and after the placement of DBS electrodes via frontal bur holes in 44 procedures. Three-dimensional coordinates in MR space of unique superficial and deep brain structures were recorded, and the magnitude, direction, and rate of shift were calculated. Measurements were recorded to the nearest 0.1 mm.

Results

Shift ranged from 0.0 to 10.1 mm throughout all structures in the brain. The greatest shift was seen in the frontal lobe, followed by the temporal and occipital lobes. Shift was also observed in deep structures such as the anterior and posterior commissures and basal ganglia; shift in the pallidum and subthalamic region ipsilateral to the bur hole averaged 0.6 mm, with 9% of patients having over 2 mm of shift in deep brain structures. Small amounts of shift were observed during all procedures; however, the initial degree of shift and its direction were unpredictable.

Conclusions

Brain shift is continual and unpredictable and can render traditional stereotactic targeting based on preoperative imaging inaccurate even in deep brain structures such as those used for DBS.

Spatially-resolved hydraulic conductivity estimation via poroelastic magnetic resonance elastography

Poroelastic magnetic resonance elastography is an imaging technique that could recover mechanical and hydrodynamical material properties of in vivo tissue. To date, mechanical properties have been estimated while hydrodynamical parameters have been assumed homogeneous with literature-based values. Estimating spatially-varying hydraulic conductivity would likely improve model accuracy and provide new image information related to a tissue's interstitial fluid compartment. A poroelastic model was reformulated to recover hydraulic conductivity with more appropriate fluid-flow boundary conditions. Simulated and physical experiments were conducted to evaluate the accuracy and stability of the inversion algorithm. Simulations were accurate (property errors were < 2%) even in the presence of Gaussian measurement noise up to 3%. The reformulated model significantly decreased variation in the shear modulus estimate (p << 0.001) and eliminated the homogeneity assumption and the need to assign hydraulic conductivity values from literature. Material property contrast was recovered experimentally in three different tofu phantoms and the accuracy was improved through soft-prior regularization. A frequency-dependence in hydraulic conductivity contrast was observed suggesting that fluid-solid interactions may be more prominent at low frequency. In vivo recovery of both structural and hydrodynamical characteristics of tissue could improve detection and diagnosis of neurological disorders such as hydrocephalus and brain tumors.

Near Real-Time Computer Assisted Surgery for Brain Shift Correction Using Biomechanical Models

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.

Model-based correction of tissue compression for tracked ultrasound in soft tissue image-guided surgery

Acquisition of ultrasound data negatively affects image registration accuracy during image-guided therapy because of tissue compression by the probe. We present a novel compression correction method that models sub-surface tissue displacement resulting from application of a tracked probe to the tissue surface. Patient landmarks are first used to register the probe pose to pre-operative imaging. The ultrasound probe geometry is used to provide boundary conditions to a biomechanical model of the tissue. The deformation field solution of the model is inverted to non-rigidly transform the ultrasound images to an estimation of the tissue geometry before compression. Experimental results with gel phantoms indicated that the proposed method reduced the tumor margin modified Hausdorff distance (MHD) from 5.0 ± 1.6 to 1.9 ± 0.6 mm, and reduced tumor centroid alignment error from 7.6 ± 2.6 to 2.0 ± 0.9 mm. The method was applied to a clinical case and reduced the tumor margin MHD error from 5.4 ± 0.1 to 2.6 ± 0.1 mm and the centroid alignment error from 7.2 ± 0.2 to 3.5 ± 0.4 mm.

Neck tumor dissection improved with 3-dimensional ultrasound image guidance: technical case report

BACKGROUND AND IMPORTANCE

Three-dimensional ultrasound navigation has been performed to assist in resection of cranial and spinal tumors, but to the best of our knowledge, no one has described the use of real-time 3-dimensional ultrasound navigation in the resection of neck tumors beyond biopsy.

CLINICAL PRESENTATION

This case report describes the use of 3-dimensional ultrasonic navigation in assisting with resection of a large neck paraganglioma. The 3-dimensional ultrasonic navigation improved real-time visualization of the carotid arteries, the trachea, and other vital structures.

CONCLUSION

The use of 3-dimensional ultrasound navigation should be considered in aiding resection of large neck tumors because it can allow more efficient and safer tumor resection.

Volumetric intraoperative brain deformation compensation: model development and phantom validation

During neurosurgery, nonrigid brain deformation may affect the reliability of tissue localization based on preoperative images. To provide accurate surgical guidance in these cases, preoperative images must be updated to reflect the intraoperative brain. This can be accomplished by warping these preoperative images using a biomechanical model. Due to the possible complexity of this deformation, intraoperative information is often required to guide the model solution. In this paper, a linear elastic model of the brain is developed to infer volumetric brain deformation associated with measured intraoperative cortical surface displacement. The developed model relies on known material properties of brain tissue, and does not require further knowledge about intraoperative conditions. To provide an initial estimation of volumetric model accuracy, as well as determine the model's sensitivity to the specified material parameters and surface displacements, a realistic brain phantom was developed. Phantom results indicate that the linear elastic model significantly reduced localization error due to brain shift, from > 16 mm to under 5 mm, on average. In addition, though in vivo quantitative validation is necessary, preliminary application of this approach to images acquired during neocortical epilepsy cases confirms the feasibility of applying the developed model to in vivo data.

Toward a preoperative planning tool for brain tumor resection therapies

BACKGROUND

Neurosurgical procedures involving tumor resection require surgical planning such that the surgical path to the tumor is determined to minimize the impact on healthy tissue and brain function. This work demonstrates a predictive tool to aid neurosurgeons in planning tumor resection therapies by finding an optimal model-selected patient orientation that minimizes lateral brain shift in the field of view. Such orientations may facilitate tumor access and removal, possibly reduce the need for retraction, and could minimize the impact of brain shift on image-guided procedures.

METHODS

In this study, preoperative magnetic resonance images were utilized in conjunction with pre- and post-resection laser range scans of the craniotomy and cortical surface to produce patient-specific finite element models of intraoperative shift for 6 cases. These cases were used to calibrate a model (i.e., provide general rules for the application of patient positioning parameters) as well as determine the current model-based framework predictive capabilities. Finally, an objective function is proposed that minimizes shift subject to patient position parameters. Patient positioning parameters were then optimized and compared to our neurosurgeon as a preliminary study.

RESULTS

The proposed model-driven brain shift minimization objective function suggests an overall reduction of brain shift by 23 % over experiential methods.

CONCLUSIONS

This work recasts surgical simulation from a trial-and-error process to one where options are presented to the surgeon arising from an optimization of surgical goals. To our knowledge, this is the first realization of an evaluative tool for surgical planning that attempts to optimize surgical approach by means of shift minimization in this manner.

Serial FEM/XFEM-Based Update of Preoperative Brain Images Using Intraoperative MRI

Current neuronavigation systems cannot adapt to changing intraoperative conditions over time. To overcome this limitation, we present an experimental end-to-end system capable of updating 3D preoperative images in the presence of brain shift and successive resections. The heart of our system is a nonrigid registration technique using a biomechanical model, driven by the deformations of key surfaces tracked in successive intraoperative images. The biomechanical model is deformed using FEM or XFEM, depending on the type of deformation under consideration, namely, brain shift or resection. We describe the operation of our system on two patient cases, each comprising five intraoperative MR images, and we demonstrate that our approach significantly improves the alignment of nonrigidly registered images.

Pre- and intraoperative tractographic evaluation of corticospinal tract shift

BACKGROUND

Magnetic resonance with diffusion tensor image (DTI) may be able to estimate trajectories compatible with subcortical tracts close to brain lesions. A limit of DTI is brain shifting (movement of the brain after dural opening and tumor resection).

OBJECTIVE

To calculate the brain shift of trajectories compatible with the corticospinal tract (CST) in patients undergoing glioma resection and predict the shift directions of CST.

METHODS

DTI was acquired in 20 patients and carried out through 12 noncollinear directions. Dedicated software "merged" all sequences acquired with tractographic processing and the whole dataset was sent to the neuronavigation system. Preoperative, after dural opening (in 11) and tumor resection (in all) DTI acquisitions were performed to evaluate CST shifting. The extent of shifting was considered as the maximum distance between the preoperative and intraoperative contours of the trajectories.

RESULTS

An outward shift of CST was observed in 8 patients and an inward shift in 10 patients during surgery. In the remaining 2 patients, no intraoperative displacement was detected. Only peritumoral edema showed a statistically significant correlation with the amount of shift. In those patients in which DTI was acquired after dural opening as well (11 patients), an outward shifting of CST was evident in that phase.

CONCLUSION

The use of intraoperative DTI demonstrated brain shifting of the CST. DTI evaluation of white matter tracts can be used during surgical procedures only if updated with intraoperative acquisitions.

3D XFEM-based modeling of retraction for preoperative image update

Outcomes for neurosurgery patients can be improved by enhancing intraoperative navigation and guidance. Current navigation systems do not accurately account for intraoperative brain deformation. So far, most studies of brain deformation have focused on brain shift, whereas this paper focuses on the brain deformation due to retraction. The heart of our system is a 3D nonrigid registration technique using a biomechanical model driven by the deformations of key surfaces tracked between two intraoperative images. The key surfaces, e.g., the whole-brain region boundary and the lips of the retraction cut, thus deform due to the combination of gravity and retractor deployment. The tissue discontinuity due to retraction is handled via the eXtended Finite Element Method (XFEM), which has the appealing feature of being able to handle arbitrarily shaped discontinuity without any remeshing. Our approach is shown to significantly improve the alignment of intraoperative MRI.

2D XFEM-based modeling of retraction and successive resections for preoperative image update

This paper considers an approach to improving outcomes for neurosurgery patients by enhancing intraoperative navigation and guidance. Currently, intraoperative navigation systems do not accurately account for brain shift or tissue resection. We describe how preoperative images can be incrementally updated to take into account any type of brain tissue deformation that may occur during surgery, and thus to improve the accuracy of image-guided navigation systems. For this purpose, we have developed a non-rigid image registration technique using a biomechanical model, which deforms based on the Finite Element Method (FEM). While the FEM has been used successfully for dealing with deformations such as brain shift, it has difficulty with tissue discontinuities. Here, we describe a novel application of the eXtended Finite Element Method (XFEM) in the field of image-guided surgery in order to model brain deformations that imply tissue discontinuities. In particular, this paper presents a detailed account of the use of XFEM for dealing with retraction and successive resections, and demonstrates the feasibility of the approach by considering 2D examples based on intraoperative MR images. To evaluate our results, we compute the modified Hausdorff distance between Canny edges extracted from images before and after registration. We show that this distance decreases after registration, and thus demonstrate that our approach improves alignment of intraoperative images.

Standard guidelines for publication of deep brain stimulation studies in Parkinson's disease (Guide4DBS-PD)

While the use of deep brain stimulation (DBS) for the treatment of neurological disorders has risen substantially over the last decade, it is often difficult to compare the results from different studies due to the lack of consistent reporting of key study parameters. We present guidelines to standardize the reporting of clinical studies of DBS for Parkinson's disease (PD). These guidelines provide a minimal set of required data elements to facilitate the interpretation and comparison of results across published clinical studies. The guidelines, summarized in the format of a checklist, may also have utility in the planning of clinical studies of DBS for PD as well as other neurological and psychiatric disorders.

Model-updated image-guided liver surgery: preliminary results using surface characterization

The current protocol for image guidance in open abdominal liver tumor removal surgeries involves a rigid registration between the patient's operating room space and the pre-operative diagnostic image-space. Systematic studies have shown that the liver can deform up to 2 cm during surgeries in a non-rigid fashion thereby compromising the accuracy of these surgical navigation systems. Compensating for intra-operative deformations using mathematical models has shown promising results. In this work, we follow up the initial rigid registration with a computational approach that is geared towards minimizing the residual closest point distances between the un-deformed pre-operative surface and the rigidly registered intra-operative surface. We also use a surface Laplacian equation based filter that generates a realistic deformation field. Preliminary validation of the proposed computational framework was performed using phantom experiments and clinical trials. The proposed framework improved the rigid registration errors for the phantom experiments on average by 43%, and 74% using partial and full surface data, respectively. With respect to clinical data, it improved the closest point residual error associated with rigid registration by 54% on average for the clinical cases. These results are highly encouraging and suggest that computational models can be used to increase the accuracy of image-guided open abdominal liver tumor removal surgeries.

Intraoperative brain shift compensation: accounting for dural septa

Biomechanical models that describe soft tissue deformation provide a relatively inexpensive way to correct registration errors in image-guided neurosurgical systems caused by nonrigid brain shift. Quantifying the factors that cause this deformation to sufficient precision is a challenging task. To circumvent this difficulty, atlas-based methods have been developed recently that allow for uncertainty, yet still capture the first-order effects associated with deformation. The inverse solution is driven by sparse intraoperative surface measurements, which could bias the reconstruction and affect the subsurface accuracy of the model prediction. Studies using intraoperative MR have shown that the deformation in the midline, tentorium, and contralateral hemisphere is relatively small. The dural septa act as rigid membranes supporting the brain parenchyma and compartmentalizing the brain. Accounting for these structures in models may be an important key to improving subsurface shift accuracy. A novel method to segment the tentorium cerebelli will be described, along with the procedure for modeling the dural septa. Results in seven clinical cases show a qualitative improvement in subsurface shift accuracy making the predicted deformation more congruous with previous observations in the literature. The results also suggest a considerably more important role for hyperosmotic drug modeling for the intraoperative shift correction environment.

A sparse intraoperative data-driven biomechanical model to compensate for brain shift during neuronavigation

BACKGROUND AND PURPOSE

Intraoperative brain deformation is an important factor compromising the accuracy of image-guided neurosurgery. The purpose of this study was to elucidate the role of a model-updated image in the compensation of intraoperative brain shift.

MATERIALS AND METHODS

An FE linear elastic model was built and evaluated in 11 patients with craniotomies. To build this model, we provided a novel model-guided segmentation algorithm. After craniotomy, the sparse intraoperative data (the deformed cortical surface) were tracked by a 3D LRS. The surface deformation, calculated by an extended RPM algorithm, was applied on the FE model as a boundary condition to estimate the entire brain shift. The compensation accuracy of this model was validated by the real-time image data of brain deformation acquired by intraoperative MR imaging.

RESULTS

The prediction error of this model ranged from 1.29 to 1.91 mm (mean, 1.62 ± 0.22 mm), and the compensation accuracy ranged from 62.8% to 81.4% (mean, 69.2 ± 5.3%). The compensation accuracy on the displacement of subcortical structures was higher than that of deep structures (71.3 ± 6.1%:66.8 ± 5.0%, P < .01). In addition, the compensation accuracy in the group with a horizontal bone window was higher than that in the group with a nonhorizontal bone window (72.0 ± 5.3%:65.7 ± 2.9%, P < .05).

CONCLUSIONS

Combined with our novel model-guided segmentation and extended RPM algorithms, this sparse data-driven biomechanical model is expected to be a reliable, efficient, and convenient approach for compensation of intraoperative brain shift in image-guided surgery.

Contrast detection in fluid-saturated media with magnetic resonance poroelastography

PURPOSE

Recent interest in the poroelastic behavior of tissues has led to the development of magnetic resonance poroelastography (MRPE) as an alternative to single-phase MR elastographic image reconstruction. In addition to the elastic parameters (i.e., Lamé's constants) commonly associated with magnetic resonance elastography (MRE), MRPE enables estimation of the time-harmonic pore-pressure field induced by external mechanical vibration.

METHODS

This study presents numerical simulations that demonstrate the sensitivity of the computed displacement and pore-pressure fields to a priori estimates of the experimentally derived model parameters. In addition, experimental data collected in three poroelastic phantoms are used to assess the quantitative accuracy of MR poroelastographic imaging through comparisons with both quasistatic and dynamic mechanical tests.

RESULTS

The results indicate hydraulic conductivity to be the dominant parameter influencing the deformation behavior of poroelastic media under conditions applied during MRE. MRPE estimation of the matrix shear modulus was bracketed by the values determined from independent quasistatic and dynamic mechanical measurements as expected, whereas the contrast ratios for embedded inclusions were quantitatively similar (10%-15% difference between the reconstructed images and the mechanical tests).

CONCLUSIONS

The findings suggest that the addition of hydraulic conductivity and a viscoelastic solid component as parameters in the reconstruction may be warranted.

A brain-deformation framework based on a linear elastic model and evaluation using clinical data

In image-guided neurosurgery, brain tissue displacement and deformation during neurosurgical procedures are a major source of error. In this paper, we implement and evaluate a linear-elastic-model-based framework for correction of brain shift using clinical data from five brain tumor patients. The framework uses a linear elastic model to simulate brain-shift behavior. The model is driven by cortical surface deformations, which are tracked using a surface-tracking algorithm combined with a laser-range scanner. The framework performance was evaluated using displacements of anatomical landmarks, tumor contours and self-defined evaluation parameters. The results show that tumor deformations predicted by the present framework agreed well with the ones observed intraoperatively, especially in the parts of the larger deformations. On average, a brain shift of 3.9 mm and a tumor margin shift of 4.2 mm were corrected to 1.2 and 1.3 mm, respectively. The entire correction process was performed in less than 5 min. The data from this study suggest that the technique is a suitable candidate for intraoperative brain-deformation correction.

A fast and efficient method to compensate for brain shift for tumor resection therapies measured between preoperative and postoperative tomograms

In this paper, an efficient paradigm is presented to correct for brain shift during tumor resection therapies. For this study, high resolution preoperative (pre-op) and postoperative (post-op) MR images were acquired for eight in vivo patients, and surface/subsurface shift was identified by manual identification of homologous points between the pre-op and immediate post-op tomograms. Cortical surface deformation data were then used to drive an inverse problem framework. The manually identified subsurface deformations served as a comparison toward validation. The proposed framework recaptured 85% of the mean subsurface shift. This translated to a subsurface shift error of 0.4 +/- 0.4 mm for a measured shift of 3.1 +/- 0.6 mm. The patient's pre-op tomograms were also deformed volumetrically using displacements predicted by the model. Results presented allow a preliminary evaluation of correction both quantitatively and visually. While intraoperative (intra-op) MR imaging data would be optimal, the extent of shift measured from pre- to post-op MR was comparable to clinical conditions. This study demonstrates the accuracy of the proposed framework in predicting full-volume displacements from sparse shift measurements. It also shows that the proposed framework can be extended and used to update pre-op images on a time scale that is compatible with surgery.