Fully Automatic Segmentation of AP Pelvis X-rays via Random Forest Regression and Hierarchical Sparse Shape Composition

Knowledge of landmarks and contours in anteroposterior (AP) pelvis X-rays is invaluable for computer aided diagnosis, hip surgery planning and image-guided interventions. This paper presents a fully automatic and robust approach for landmarking and segmentation of both pelvis and femur in a conventional AP X-ray. Our approach is based on random forest regression and hierarchical sparse shape composition. Experiments conducted on 436 clinical AP pelvis x-rays show that our approach achieves an average point-to-curve error around 1.3 mm for femur and 2.2 mm for pelvis, both with success rates around 98%. Compared to existing methods, our approach exhibits better performance in both the robustness and the accuracy.

Three-dimensional, task-specific robot therapy of the arm after stroke: a multicentre, parallel-group randomised trial

BACKGROUND: Arm hemiparesis secondary to stroke is common and disabling. We aimed to assess whether robotic training of an affected arm with ARMin--an exoskeleton robot that allows task-specific training in three dimensions-reduces motor impairment more effectively than does conventional therapy. METHODS: In a prospective, multicentre, parallel-group randomised trial, we enrolled patients who had had motor impairment for more than 6 months and moderate-to-severe arm paresis after a cerebrovascular accident who met our eligibility criteria from four centres in Switzerland. Eligible patients were randomly assigned (1:1) to receive robotic or conventional therapy using a centre-stratified randomisation procedure. For both groups, therapy was given for at least 45 min three times a week for 8 weeks (total 24 sessions). The primary outcome was change in score on the arm (upper extremity) section of the Fugl-Meyer assessment (FMA-UE). Assessors tested patients immediately before therapy, after 4 weeks of therapy, at the end of therapy, and 16 weeks and 34 weeks after start of therapy. Assessors were masked to treatment allocation, but patients, therapists, and data analysts were unmasked. Analyses were by modified intention to treat. This study is registered with ClinicalTrials.gov, number NCT00719433. FINDINGS: Between May 4, 2009, and Sept 3, 2012, 143 individuals were tested for eligibility, of whom 77 were eligible and agreed to participate. 38 patients assigned to robotic therapy and 35 assigned to conventional therapy were included in analyses. Patients assigned to robotic therapy had significantly greater improvements in motor function in the affected arm over the course of the study as measured by FMA-UE than did those assigned to conventional therapy (F=4.1, p=0.041; mean difference in score 0.78 points, 95% CI 0.03-1.53). No serious adverse events related to the study occurred. INTERPRETATION: Neurorehabilitation therapy including task-oriented training with an exoskeleton robot can enhance improvement of motor function in a chronically impaired paretic arm after stroke more effectively than conventional therapy. However, the absolute difference between effects of robotic and conventional therapy in our study was small and of weak significance, which leaves the clinical relevance in question.

Three-dimensional electron beam dose calculations

The MDAH pencil-beam algorithm developed by Hogstrom et al (1981) has been widely used in clinics for electron beam dose calculations for radiotherapy treatment planning. The primary objective of this research was to address several deficiencies of that algorithm and to develop an enhanced version. Two enhancements have been incorporated into the pencil-beam algorithm; one models fluence rather than planar fluence, and the other models the bremsstrahlung dose using measured beam data. Comparisons of the resulting calculated dose distributions with measured dose distributions for several test phantoms have been made. From these results it is concluded (1) that the fluence-based algorithm is more accurate to use for the dose calculation in an inhomogeneous slab phantom, and (2) the fluence-based calculation provides only a limited improvement to the accuracy the calculated dose in the region just downstream of the lateral edge of an inhomogeneity. The source of the latter inaccuracy is believed primarily due to assumptions made in the pencil beam's modeling of the complex phantom or patient geometry.^ A pencil-beam redefinition model was developed for the calculation of electron beam dose distributions in three dimensions. The primary aim of this redefinition model was to solve the dosimetry problem presented by deep inhomogeneities, which was the major deficiency of the enhanced version of the MDAH pencil-beam algorithm. The pencil-beam redefinition model is based on the theory of electron transport by redefining the pencil beams at each layer of the medium. The unique approach of this model is that all the physical parameters of a given pencil beam are characterized for multiple energy bins. Comparisons of the calculated dose distributions with measured dose distributions for a homogeneous water phantom and for phantoms with deep inhomogeneities have been made. From these results it is concluded that the redefinition algorithm is superior to the conventional, fluence-based, pencil-beam algorithm, especially in predicting the dose distribution downstream of a local inhomogeneity. The accuracy of this algorithm appears sufficient for clinical use, and the algorithm is structured for future expansion of the physical model if required for site specific treatment planning problems. ^

Three-dimensional printing creates models for surgical planning of aortic valve replacement after previous coronary bypass grafting.

PURPOSE Resternotomy for aortic valve replacement in patients with previous coronary artery bypass grafting and an internal mammary artery graft may be a surgical problem. Thus, we are exploring the effect of using rapid prototyping techniques for surgical planning and intraoperative orientation during aortic valve replacement after previous coronary artery bypass grafting (CABG). DESCRIPTION As a proof of concept, we studied a patient who had undergone CABG 5 years earlier. At that time the patient received a left internal mammary artery graft to the left anterior descending artery and a venous graft to the right coronary artery. Now the patient required aortic valve replacement due to symptomatic aortic valve stenosis. The left internal mammary artery bypass and the right coronary artery bypass were patent and showed good flow in the angiography. The patient was examined by 128-slice computed tomography. The image data were visualized and reconstructed. Afterwards, a replica showing the anatomic structures was fabricated using a rapid prototyping machine. EVALUATION Using data derived from 128-slice computed tomography angiography linked to proprietary software, we were able to create three-dimensional reconstructions of the vascular anatomy after the previous CABG. The models were sterilized and taken to the operating theatre for orientation during the surgical procedure. CONCLUSIONS Stereolithographic replicas are helpful for choosing treatment strategies in surgical planning and for intraoperative orientation during reoperations of patients with previous CABG.

Extension and Partitioning in an Oblique Subduction Zone, New Zealand: Constraints from Three-Dimensional Numerical Modeling

Contraction, strike slip, and extension displacements along the Hikurangi margin northeast of the North Island of New Zealand coincide with large lateral gradients in material properties. We use a finite- difference code utilizing elastic and elastic-plastic rheologies to build large- scale, three-dimensional numerical models which investigate the influence of material properties on velocity partitioning within oblique subduction zones. Rheological variation in the oblique models is constrained by seismic velocity and attenuation information available for the Hikurangi margin. We compare the effect of weakly versus strongly coupled subduction interfaces on the development of extension and the partitioning of velocity components for orthogonal and oblique convergence and include the effect of ponded sediments beneath the Raukumara Peninsula. Extension and velocity partitioning occur if the subduction interface is weak, but neither develops if the subduction interface is strong. The simple mechanical model incorporating rheological variation based on seismic observations produces kinematics that closely match those published from the Hikurangi margin. These include extension within the Taupo Volcanic Zone, uplift over ponded sediments, and dextral contraction to the south.

Three-Dimensional Mechanics of Yakutat Convergence in the Southern Alaskan Plate Corner

Three-dimensional numerical models are used to investigate the mechanical evolution of the southern Alaskan plate corner where the Yakutat and the Pacific plates converge on the North American plate. The evolving model plate boundary consists of Convergent, Lateral, and Subduction subboundaries with flow separation of incoming material into upward or downward trajectories forming dual, nonlinear advective thermal/mechanical anomalies that fix the position of major subaerial mountain belts. The model convergent subboundary evolves into two teleconnected orogens: Inlet and Outlet orogens form at locations that correspond with the St. Elias and the Central Alaska Range, respectively, linked to the East by the Lateral boundary. Basins form parallel to the orogens in response to the downward component of velocity associated with subduction. Strain along the Lateral subboundary varies as a function of orogen rheology and magnitude and distribution of erosion. Strain-dependent shear resistance of the plate boundary associated with the shallow subduction zone controls the position of the Inlet orogen. The linkages among these plate boundaries display maximum shear strain rates in the horizontal and vertical planes where the Lateral subboundary joins the Inlet and Outlet orogens. The location of the strain maxima shifts with time as the separation of the Inlet and Outlet orogens increases. The spatiotemporal predictions of the model are consistent with observed exhumation histories deduced from thermochronology, as well as stratigraphic studies of synorogenic deposits. In addition, the complex structural evolution of the St Elias region is broadly consistent with the predicted strain field evolution. Citation: Koons, P. O., B. P. Hooks, T. Pavlis, P. Upton, and A. D. Barker (2010), Three-dimensional mechanics of Yakutat convergence in the southern Alaskan plate corner, Tectonics, 29, TC4008, doi: 10.1029/2009TC002463.

Reactive mass transport modelling of a three-dimensional vertical fault zone with a finger-like convective flow regime

For a three-dimensional vertically-oriented fault zone, we consider the coupled effects of fluid flow, heat transfer and reactive mass transport, to investigate the patterns of fluid flow, temperature distribution, mineral alteration and chemically induced porosity changes. We show, analytically and numerically, that finger-like convection patterns can arise in a vertically-oriented fault zone. The onset and patterns of convective fluid flow are controlled by the Rayleigh number which is a function of the thermal properties of the fluid and the rock, the vertical temperature gradient, and the height and the permeability of the fault zone. Vigorous fluid flow causes low temperature gradients over a large region of the fault zone. In such a case, flow across lithological interfaces becomes the most important mechanism for the formation of sharp chemical reaction fronts. The degree of rock buffering, the extent and intensity of alteration, the alteration mineralogy and in some cases the formation of ore deposits are controlled by the magnitude of the flow velocity across these compositional interfaces in the rock. This indicates that alteration patterns along compositional boundaries in the rock may provide some insights into the convection pattern. The advective mass and heat exchanges between the fault zone and the wallrock depend on the permeability contrast between the fault zone and the wallrock. A high permeability contrast promotes focussed convective flow within the fault zone and diffusive exchange of heat and chemical reactants between the fault zone and the wallrock. However, a more gradual permeability change may lead to a regional-scale convective flow system where the flow pattern in the fault affects large-scale fluid flow, mass transport and chemical alteration in the wallrocks

Rivastigmine effects on EEG spectra and three-dimensional LORETA functional imaging in Alzheimer's disease

OBJECTIVE The objective of the study is to investigate the electrocortical and the global cognitive effects of 3 months rivastigmine medication in a group of mild to moderate Alzheimer's disease patients. MATERIALS AND METHODS Multichannel EEG and cognitive performances measured with the Mini Mental State Examination in a group of 16 patients with mild to moderate Alzheimer's Disease were collected before and 3 months after the onset of rivastigmine medication. RESULTS Spectral analysis of the EEG data showed a significant power decrease in the delta and theta frequency bands during rivastigmine medication, i.e., a shift of the power spectrum towards 'normalization'. Three-dimensional low resolution electromagnetic tomography (LORETA) functional imaging localized rivastigmine effects in a network that includes left fronto-parietal regions, posterior cingulate cortex, bilateral parahippocampal regions, and the hippocampus. Moreover, a correlation analysis between differences in the cognitive performances during the two recordings and LORETA-computed intracortical activity showed, in the alpha1 frequency band, better cognitive performance with increased cortical activity in the left insula. CONCLUSION The results point to a 'normalization' of the EEG power spectrum due to medication, and the intracortical localization of these effects showed an increase of cortical activity in frontal, parietal, and temporal regions that are well-known to be affected in Alzheimer's disease. The topographic convergence of the present results with the memory network proposed by Vincent et al. (J. Neurophysiol. 96:3517-3531, 2006) leads to the speculation that in our group of patients, rivastigmine specifically activates brain regions that are involved in memory functions, notably a key symptom in this degenerative disease.

Three-dimensional study of nasopalatine canal morphology: a descriptive retrospective analysis using cone-beam computed tomography

PURPOSE For dental implant treatment planning and placement, a precise anatomic description of the nasopalatine canal (NC) is necessary. This descriptive retrospective study evaluated dimensions of the NC and buccal bone plate (BBP) and the tridimensional association of the anatomic variants of NC, using cone-beam computed tomography (CBCT). METHODS This study included 230 CBCTs. Sagittal slices were used for measurements of the NC and BBP and to evaluate shape and direction-course of the NC. Coronal slices were used to assess NC shape and axial slices to assess number of incisive foramina and foramina of Stenson. RESULTS Mean NC length was 12.34 ± 2.79 mm, statistically significant differences were detected between genders (p < 0.001). Mean BBP length was 20.87 ± 3.68 mm, statistically significant differences were found for the dental status (p < 0.001) and mean BBP width was 6.83 ± 1.28 mm, significant differences were detected between genders (p < 0.001). Mean nasopalatine angle was 73.33° ± 8.11°, significant differences were found in sagittal and coronal classifications. The most prevalent canal was: cylindrical sagittal shape (48.2 %); slanted-straight direction-course (57.6 %); Ya-type coronal shape (42.4 %); and one foramen incisive with two Stenson's foramina (1-2) (50.9 %). Sagittal shape was associated with sagittal direction-course (p < 0.001). Coronal shape was associated with axial classification (p < 0.001). CONCLUSIONS The NC anatomy is highly variable. Gender is related to the NC length and BBP width, while dental status is related to BBP length. There was an association between the different sagittal classifications of the NC and between the coronal shape and axial classification.

Three-dimensional printing for perioperative planning of complex aortic arch surgery

Purpose In this study, we show the use of three-dimensional printing models for preoperative planning of surgery for patients with complex aortic arch anomalies. Description A 70-year-old man with an extensively arteriosclerotic aneurysm reaching from the ascending aorta to the descending aorta was referred to our center for complete aortic arch replacement. We visualized and reconstructed computed tomography data of the patient and fabricated a flexible three-dimensional model of the aortic arch including the aneurysm. Evaluation This model was very helpful for the preoperative decision making and planning of the frozen elephant trunk procedure owing to the exact and lifelike illustration of the native aortic arch. Conclusions Three-dimensional models are helpful in preoperative planning and postoperative evaluation of frozen elephant trunk procedures in patients with complex aortic anatomy.

Fully automatic segmentation of AP pelvis X-rays via random forest regression with efficient feature selection and hierarchical sparse shape composition

In clinical practice, traditional X-ray radiography is widely used, and knowledge of landmarks and contours in anteroposterior (AP) pelvis X-rays is invaluable for computer aided diagnosis, hip surgery planning and image-guided interventions. This paper presents a fully automatic approach for landmark detection and shape segmentation of both pelvis and femur in conventional AP X-ray images. Our approach is based on the framework of landmark detection via Random Forest (RF) regression and shape regularization via hierarchical sparse shape composition. We propose a visual feature FL-HoG (Flexible- Level Histogram of Oriented Gradients) and a feature selection algorithm based on trace radio optimization to improve the robustness and the efficacy of RF-based landmark detection. The landmark detection result is then used in a hierarchical sparse shape composition framework for shape regularization. Finally, the extracted shape contour is fine-tuned by a post-processing step based on low level image features. The experimental results demonstrate that our feature selection algorithm reduces the feature dimension in a factor of 40 and improves both training and test efficiency. Further experiments conducted on 436 clinical AP pelvis X-rays show that our approach achieves an average point-to-curve error around 1.2 mm for femur and 1.9 mm for pelvis.

The Utilization of Classical Spin Monte Carlo Methods to Simulate the Magnetic Behavior of Extended Three-Dimensional Cubic Networks Incorporating M(II) Ions with an S = 5/2 Ground State Spin

The numerical simulations of the magnetic properties of extended three-dimensional networks containing M(II) ions with an S = 5/2 ground-state spin have been carried out within the framework of the isotropic Heisenberg model. Analytical expressions fitting the numerical simulations for the primitive cubic, diamond, together with (10−3) cubic networks have all been derived. With these empirical formulas in hands, we can now extract the interaction between the magnetic ions from the experimental data for these networks. In the case of the primitive cubic network, these expressions are directly compared with those from the high-temperature expansions of the partition function. A fit of the experimental data for three complexes, namely [(N(CH3)4][Mn(N3)] 1, [Mn(CN4)]n 2, and [FeII(bipy)3][MnII2(ox)3] 3, has been carried out. The best fits were those obtained using the following parameters, J = −3.5 cm-1, g = 2.01 (1); J = −8.3 cm-1, g = 1.95 (2); and J = −2.0 cm-1, g = 1.95 (3).