Biomechanical wall properties of human intracranial aneurysms resected following surgical clipping (IRRAs Project).

BACKGROUND AND PURPOSE: Individual rupture risk assessment of intracranial aneurysms is a major issue in the clinical management of asymptomatic aneurysms. Aneurysm rupture occurs when wall tension exceeds the strength limit of the wall tissue. At present, aneurysmal wall mechanics are poorly understood and thus, risk assessment involving mechanical properties is inexistent. Aneurysm computational hemodynamics studies make the assumption of rigid walls, an arguable simplification. We therefore aim to assess mechanical properties of ruptured and unruptured intracranial aneurysms in order to provide the foundation for future patient-specific aneurysmal risk assessment. This work also challenges some of the currently held hypotheses in computational flow hemodynamics research. METHODS: A specific conservation protocol was applied to aneurysmal tissues following clipping and resection in order to preserve their mechanical properties. Sixteen intracranial aneurysms (11 female, 5 male) underwent mechanical uniaxial stress tests under physiological conditions, temperature, and saline isotonic solution. These represented 11 unruptured and 5 ruptured aneurysms. Stress/strain curves were then obtained for each sample, and a fitting algorithm was applied following a 3-parameter (C(10), C(01), C(11)) Mooney-Rivlin hyperelastic model. Each aneurysm was classified according to its biomechanical properties and (un)rupture status. RESULTS: Tissue testing demonstrated three main tissue classes: Soft, Rigid, and Intermediate. All unruptured aneurysms presented a more Rigid tissue than ruptured or pre-ruptured aneurysms within each gender subgroup. Wall thickness was not correlated to aneurysmal status (ruptured/unruptured). An Intermediate subgroup of unruptured aneurysms with softer tissue characteristic was identified and correlated with multiple documented risk factors of rupture. CONCLUSION: There is a significant modification in biomechanical properties between ruptured aneurysm, presenting a soft tissue and unruptured aneurysms, presenting a rigid material. This finding strongly supports the idea that a biomechanical risk factor based assessment should be utilized in the to improve the therapeutic decision making.

New electrophysiological mapping combined with MRI in parkinsonian's subthalamic region.

The subthalamic nucleus (STN) is the main target for deep brain stimulation in Parkinson's disease. We analysed the relationships between magnetic resonance imaging (MRI) anatomy and spontaneous neuronal activity to confirm the potential of microelectrode recordings to assist in determining the optimal surgical target. Ten bilateral surgeries were performed after 1.5-T (T2-weighted) anatomical MRI identification of the STN, zona incerta (ZI), Forel's field H2 (H2) and substantia nigra (SN). Spontaneous neuronal activity was recorded simultaneously along the distal 10 mm on a central track (optimally covering the STN) and a 2-mm anterior track. We calculated off-line mean firing rate and burst frequency on 248 neurons clustered according to anatomical structure. Subjective visual analysis of signal was also realized on-line, during surgery, to classify patterns of activity. Mean firing rate and burst frequency increased from H2-ZI to SN. The mean firing rate was higher in SN only using paired comparison (SN vs. its neighbours). The burst frequency was lower in H2 than in SN; using comparison with neighbours, it was lower in H2 and ZI. An irregular high activity (type 2C) was more often detected in STN and SN than in H2 and ZI. Anatomical boundaries and unitary recordings appear to be linked, supporting the ability of MRI to provide a detailed anatomy. Electrophysiological mapping combined with MRI is a useful tool for precise targeting in the subthalamic region.

Evaluation of a stochastic approach to 2D-3D intensity-based registration.

this paper assesses a new 2D-3D rigid registration method based on stochastic clustering and an enhanced estimator for mutual information. It combines precision, accuracy and acceptable computation time. Spine datasets (fluoroscopy and computed tomography) with their gold standard transformations are used. Both optimization method and similarity measure are assessed separately using standardized evaluation methodologies. Sub-millimeter accuracy and the high convergence rate obtained within one minute are compared to other quasi-global optimization processes such as particle filtering.

In vitro evaluation of 2D-digital subtraction angiography versus 3D-time-of-flight in assessment of intracranial cerebral aneurysm filling after endovascular therapy.

BACKGROUND AND PURPOSE: The aim of this study was to evaluate 2D-digital subtraction angiographic (DSA) and 3D-time-of-flight (TOF) MR imaging in assessment of aneurysmal residue by using a pulsating silicon aneurysm model. For each imaging system, we studied intra- and interobserver reproducibility and the agreement between interpretations and reference measurements. We also examined how each imaging technique affected the operator's therapeutic decision. METHODS: Two silicon aneurysm models depicting subarachnoidal aneurysms were used, one with a wide neck and one with a narrow neck. Each aneurysm model was placed in series on a pulsed flow circuit and was filled with Guglielmi detachable coils to simulate a clinical case. Each aneurysm was then gradually filled with silicon gel in increments of 10%, up to 100% to simulate different levels of occlusion (residual neck or dog ear, partial, complete) at each filling level. For each level of filling, we performed conventional 2D-DSA and 3D-TOF MR imaging. We submitted the images for examination by 2 senior medical staff with 2 readings per image. A combined reading of the 2 images was submitted to each expert to determine whether the 2 examinations were complementary. RESULTS: The 2D-DSA analysis showed good reproducibility (k = 0.8 and k = 0.57) and agreement (k = 0.71) in describing "complete" treatments. The distinction between a "residual neck" and "partial treatment," however, was not reliable. The 2D-DSA provided a good description of the coil and silicon protrusion into the parent artery. The 3D-TOF analysis of the residual aneurysm, however, was not reproducible, though it was more effective than the 2D-DSA in evaluation of partially wide-necked aneurysms (k = 0.68 MR imaging vs k = 0.041 2D-DSA; P = .018). At the same filling level, the 2D-DSA analysis indicated repeat treatment more often than 3D-TOF analysis (P = .059). CONCLUSION: The 2D-DSA remains the gold standard, but MR imaging is more effective in evaluating a "partial treatment." The 2D-DSA analysis indicated repeat treatment more often than the 3D-TOF for the same occlusion level. The distinction between "partial treatment" and a "residual neck" was not reliable with either method of evaluation.

Monitoring of polyethylene wear in nonmetal-backed acetubular cups by digitized anteroposterior pelvic radiography.

The aim of this study was to assess polyethylene wear in a total hip prosthesis by digitized radiography of the whole pelvis in the anteroposterior (AP) plane. The three-dimensional (3-D) pose of the nonmetal-backed acetubular cup, materialized by its metal ring and the femoral head made of metal or ceramic, was estimated using iterative algebraic algorithms with inner bias correction and bootstrapping for variance reduction. Points of interest were obtained by maximizing the correlation between sampled density profiles and 3-D geometric models degraded by the modulation transfer function (MTF) of the radiographic system and the film scanner. The error in the maximal correlation estimate were inferred from noise power spectra (NPS) and allowed the calculation of the point covariance matrix. Both NPS and MTF were modeled for each stage and estimated using least-square fitting of the overall NPS model to the autospectral density function calculated in stationary regions. Comparison of the radiographic time series was made possible by the high accuracy level and 3-D matching from the cup orientation. The feasibility of the full 3-D measurement, the assumption of negligible lateral wear and its influence on AP wear are discussed on simulated and real radiographic data.

Blood flow velocity estimation from x-ray densitometric data: an efficient numerical scheme for the inverse advection problem.

In previously published studies, blood flow velocity from x-ray biplane angiography was measured by solving an inverse advection problem, relating velocity to bolus densities summed across sections. Both spatial and temporal velocity variations were recovered through a computationally expensive parameter estimation algorithm. Here we prove the existence and uniqueness of the solution on three sub-domains of the plane defined by the axial position along the vessel and the time of the angiographic sequence. A fast direct scheme was designed in conjunction with a regularization step stemming from the volume flow conservation law applied on consecutive segments. Its accuracy and immunity towards noise were tested on both simulated and real densitometric data. The relative error between the estimated and expected velocities was less than 5% for more than 90% of the points of the spatiotemporal plane with simulated densities normalized to 1.0 and a Gaussian additive noise of standard deviation 0.01. For densities reconstructed from a biplane angiographic sequence, increase in velocity is used as a functional index for the stenosis ratio and to characterize the sharing of flow at bifurcation.

Three-dimensional tracking of coronary arteries from biplane angiographic sequences using parametrically deformable models.

A new method for coronary artery tracking in biplane digital subtraction is presented. The dynamic tracking of nonrigid objects from two views is achieved using a generalization of parametrically deformable models. Three-dimensional (3-D) Fourier descriptors used for shape representation are obtained from the two-dimensional (2-D) descriptors of the projections. A new constraint inferred from epipolar geometry is applied to the contour model. Direct 3-D tracking is compared with the classical approach in two steps: independent 2-D tracking in each of the two projection planes; 3-D reconstruction using the epipolar constraint. Convergence quality and accuracy of the 3-D reconstruction are analyzed for several sequences showing different displacement amplitudes, deformation rates and image contrasts.

Assessment of stenosis severity using a novel method to estimate spatial and temporal variations of blood flow velocity in biplane coronarography.

The authors present a novel method to estimate absolute blood flow velocity in coronary arteries from biplane angiograms. Spatial and temporal velocity variations are derived giving simultaneously a direct geometric and an indirect functional index of stenosis severity, stenosis ratio and coronary flow reserve. No prior assumption concerning stenosis geometry is made. Deformable models are used to track a coronary artery segment dynamically in three dimensions. A densitometric map is obtained by summing densities across sections at every position along the previously calculated path and at every time of the cardiac cycle. An advection relationship between density and velocity is observed. The spatiotemporal velocity map is a solution of a nonlinear least-squares scheme. A simulation protocol based on simple geometric conformations and blood flow properties is used to assess numerical stability and immunity towards noise. Predicted results for temporal velocity variations are compared with the intracoronary Doppler recordings to test the model assumptions for basal state and hyperaemia examinations of the same patient. The stenosis ratio was accurate to within 3% for a simulated additive Gaussian noise with a standard deviation of 0.14. The limits of agreement between angiographic and Doppler velocities were -11.4 and 11.8 cm s-1 for a peak value of 23 cm s-1 (basal state) and -16.8 and 13.5 cm s-1 for a peak value of 52 cm s-1 (hyperaemia), corresponding to 18 and 3.5% errors on the average peak values and a 16% error on the coronary flow reserve. To summarize, the advection model derivation and its solution are presented. Simulated and experimental results corroborate the validity of the numerical schemes and support clinical applicability.