Segmentation of cone-beam CT using a hidden Markov random field with informative priors

Cone-beam computed tomography (CBCT) has enormous potential to improve the accuracy of treatment delivery in image-guided radiotherapy (IGRT). To assist radiotherapists in interpreting these images, we use a Bayesian statistical model to label each voxel according to its tissue type. The rich sources of prior information in IGRT are incorporated into a hidden Markov random field model of the 3D image lattice. Tissue densities in the reference CT scan are estimated using inverse regression and then rescaled to approximate the corresponding CBCT intensity values. The treatment planning contours are combined with published studies of physiological variability to produce a spatial prior distribution for changes in the size, shape and position of the tumour volume and organs at risk. The voxel labels are estimated using iterated conditional modes. The accuracy of the method has been evaluated using 27 CBCT scans of an electron density phantom. The mean voxel-wise misclassification rate was 6.2\%, with Dice similarity coefficient of 0.73 for liver, muscle, breast and adipose tissue. By incorporating prior information, we are able to successfully segment CBCT images. This could be a viable approach for automated, online image analysis in radiotherapy.

Constructing a clinical decision-making framework for image-guided radiotherapy using a Bayesian Network

A decision-making framework for image-guided radiotherapy (IGRT) is being developed using a Bayesian Network (BN) to graphically describe, and probabilistically quantify, the many interacting factors that are involved in this complex clinical process. Outputs of the BN will provide decision-support for radiation therapists to assist them to make correct inferences relating to the likelihood of treatment delivery accuracy for a given image-guided set-up correction. The framework is being developed as a dynamic object-oriented BN, allowing for complex modelling with specific sub-regions, as well as representation of the sequential decision-making and belief updating associated with IGRT. A prototype graphic structure for the BN was developed by analysing IGRT practices at a local radiotherapy department and incorporating results obtained from a literature review. Clinical stakeholders reviewed the BN to validate its structure. The BN consists of a sub-network for evaluating the accuracy of IGRT practices and technology. The directed acyclic graph (DAG) contains nodes and directional arcs representing the causal relationship between the many interacting factors such as tumour site and its associated critical organs, technology and technique, and inter-user variability. The BN was extended to support on-line and off-line decision-making with respect to treatment plan compliance. Following conceptualisation of the framework, the BN will be quantified. It is anticipated that the finalised decision-making framework will provide a foundation to develop better decision-support strategies and automated correction algorithms for IGRT.

Predicting the likelihood of QA failure using treatment plan accuracy metrics

This study used automated data processing techniques to calculate a set of novel treatment plan accuracy metrics, and investigate their usefulness as predictors of quality assurance (QA) success and failure. 151 beams from 23 prostate and cranial IMRT treatment plans were used in this study. These plans had been evaluated before treatment using measurements with a diode array system. The TADA software suite was adapted to allow automatic batch calculation of several proposed plan accuracy metrics, including mean field area, small-aperture, off-axis and closed-leaf factors. All of these results were compared the gamma pass rates from the QA measurements and correlations were investigated. The mean field area factor provided a threshold field size (5 cm2, equivalent to a 2.2 x 2.2 cm2 square field), below which all beams failed the QA tests. The small aperture score provided a useful predictor of plan failure, when averaged over all beams, despite being weakly correlated with gamma pass rates for individual beams. By contrast, the closed leaf and off-axis factors provided information about the geometric arrangement of the beam segments but were not useful for distinguishing between plans that passed and failed QA. This study has provided some simple tests for plan accuracy, which may help minimise time spent on QA assessments of treatments that are unlikely to pass.

Using narrow beam profiles to quantify focal spot size, for accurate Monte Carlo simulations of SRS/SRT systems

This study investigates the variation of photon field penumbra shape with initial electron beam diameter, for very narrow beams. A Varian Millenium MLC (Varian Medical Systems, Palo Alto, USA) and a Brainlab m3 microMLC (Brainlab AB. Feldkirchen, Germany) were used, with one Varian iX linear accelerator, to produce fields that were (nominally) 0.20 cm across. Dose profiles for these fields were measured using radiochromic film and compared with the results of simulations completed using BEAMnrc and DOSXYZnrc, where the initial electron beam was set to FWHM = 0.02, 0.10, 0.12, 0.15, 0.20 and 0.50 cm. Increasing the electron-beam FWHM produced increasing occlusion of the photon source by the closely spaced collimator leaves and resulted in blurring of the simulated profile widths from 0.26 to 0.64 cm, for the MLC, from 0.12 to 0.43 cm, for the microMLC. Comparison with measurement data suggested that the electron spot size in the clinical linear accelerator was between FWHM = 0.10 and 0.15 cm, encompassing the result of our previous output-factor based work, which identified a FWHM of 0.12. Investigation of narrow-beam penumbra variation has been found to be a useful procedure, with results varying noticeably with linear accelerator spot size and allowing FWHM estimates obtained using other methods to be verified.

The influence of Monte Carlo source parameters on detector design and dose perturbation in small field dosimetry

To obtain accurate Monte Carlo simulations of small radiation fields, it is important model the initial source parameters (electron energy and spot size) accurately. However recent studies have shown that small field dosimetry correction factors are insensitive to these parameters. The aim of this work is to extend this concept to test if these parameters affect dose perturbations in general, which is important for detector design and calculating perturbation correction factors. The EGSnrc C++ user code cavity was used for all simulations. Varying amounts of air between 0 and 2 mm were deliberately introduced upstream to a diode and the dose perturbation caused by the air was quantified. These simulations were then repeated using a range of initial electron energies (5.5 to 7.0 MeV) and electron spot sizes (0.7 to 2.2 FWHM). The resultant dose perturbations were large. For example 2 mm of air caused a dose reduction of up to 31% when simulated with a 6 mm field size. However these values did not vary by more than 2 % when simulated across the full range of source parameters tested. If a detector is modified by the introduction of air, one can be confident that the response of the detector will be the same across all similar linear accelerators and the Monte Carlo modelling of each machine is not required.

Analysis on the spin reorientation of GdFeCo/DyFeCo exchange-coupled double-layer films

In this paper, GdFeCo/DyFeCo exchange-coupled double-layer films used for center aperture type magnetically induced super resolution were investigated through experiments and theoretical calculation. The samples were prepared by magnetron sputtering method. The polar Kerr effect was measured to prove the spin reorientation of the readout layer. Theoretical study of magnetization profiles was performed on the basis of the mean-field theory and the continuum model. The theoretical results showed that the magnetization orientation of the readout layer changed gradually from in-plane to out-of-plane with the rise of the temperature. Theoretical analysis explained the experimental results successfully. (c) 2005 Elsevier B.V. All rights reserved.

水生植物叶片的分形特征研究

应用分形几何理论对保安湖部分水生植物叶的分形特征进行研究 ,以寻找其中的自相似性特征 ,计算其分形维数 ,并力图了解各个分维的生态学含义。计算得到各物种分维平均值 ,马来眼子菜为 1 441(1 3 2 2 -1 52 0 ,sd =0 0 547) ,菹草为 1 472 (1 3 81-1 565,sd =0 0 489) ,大茨藻为 1 3 2 5(1 2 0 9-1 40 6,sd =0 0 470 ) ,轮叶黑藻为 1 10 3(1 0 0 1-1 2 17,sd =0 0 456)。按分维值排序

A new genus and species of freshwater monostiliferous hoplonemertean (Nemertea, Enopla) from the People's Republic of China

A new genus and species of freshwater monostiliferous hoplonemertean, Limnemertes poyangensis gen. et sp. nov., from Poyang Lake, People's Republic of China, is described and illustrated. The taxon is compared and contrasted with previously described freshwater hoplonemerteans. This is the fourth species of freshwater nemertean to be described from China and the first recorded from Poyang Lake.

Micro structures, tensile properties and corrosion behavior of die-cast Mg-4Al-based alloys containing La and/or Ce

In order to study the properties of Mg-Al-RE (AE) series alloys, the Mg-4Al-4RE-0.4Mn (RE= La, Ce/La mischmetal or Ce) alloys were developed. Their microstructures, tensile properties and corrosion behavior have been investigated. The results show that the phase compositions of Mg-4Al-4La-0.4Mn alloy consist of alpha-Mg and Al11La3 phases. While two binary Al-RE (RE = Ce/La) phases, Al11RE3 and Al2RE, are formed in Mg-4Al-4Ce/La-0.4Mn alloy, and Al11Ce3 and Al2Ce are formed in Mg-4Al-4Ce-0.4Mn alloy.

Epidemiology and predictors of spinal injury in adult major trauma patients: European cohort study

This is a European cohort study on predictors of spinal injury in adult (≥16 years) major trauma patients, using prospectively collected data of the Trauma Audit and Research Network from 1988 to 2009. Predictors for spinal fractures/dislocations or spinal cord injury were determined using univariate and multivariate logistic regression analysis. 250,584 patients were analysed. 24,000 patients (9.6%) sustained spinal fractures/dislocations alone and 4,489 (1.8%) sustained spinal cord injury with or without fractures/dislocations. Spinal injury patients had a median age of 44.5 years (IQR = 28.8-64.0) and Injury Severity Score of 9 (IQR = 4-17). 64.9% were male. 45% of patients suffered associated injuries to other body regions. Age <45 years (≥45 years OR 0.83-0.94), Glasgow Coma Score (GCS) 3-8 (OR 1.10, 95% CI 1.02-1.19), falls >2 m (OR 4.17, 95% CI 3.98-4.37), sports injuries (OR 2.79, 95% CI 2.41-3.23) and road traffic collisions (RTCs) (OR 1.91, 95% CI 1.83-2.00) were predictors for spinal fractures/dislocations. Age <45 years (≥45 years OR 0.78-0.90), male gender (female OR 0.78, 95% CI 0.72-0.85), GCS <15 (OR 1.36-1.93), associated chest injury (OR 1.10, 95% CI 1.01-1.20), sports injuries (OR 3.98, 95% CI 3.04-5.21), falls >2 m (OR 3.60, 95% CI 3.21-4.04), RTCs (OR 2.20, 95% CI 1.96-2.46) and shooting (OR 1.91, 95% CI 1.21-3.00) were predictors for spinal cord injury. Multilevel injury was found in 10.4% of fractures/dislocations and in 1.3% of cord injury patients. As spinal trauma occurred in >10% of major trauma patients, aggressive evaluation of the spine is warranted, especially, in males, patients <45 years, with a GCS <15, concomitant chest injury and/or dangerous injury mechanisms (falls >2 m, sports injuries, RTCs and shooting). Diagnostic imaging of the whole spine and a diligent search for associated injuries are substantial.

Stable oxygen isotope ratios and grain size analyses of sediment cores from the Propeller Mound, Northeast Atlantic

The first detailed stratigraphic record from a deep-water carbonate mound in the Northeast Atlantic based on absolute datings (U/Th and AMS 14C) and stable oxygen isotope records reveals that its top sediment sequences are condensed by numerous hiatuses. According to stable isotope data, mainly sediments with an intermediate signal are preserved on the mound, while almost all fully glacial and interglacial sediments have either not been deposited or have been eroded later. The resulting hiatuses reduce the Late Pleistocene sediment accumulation at Propeller Mound to amounts smaller than the background sedimentation. The hiatuses most likely result due to the sweeping of the mound in turn with the re-establishment of vigour interglacial circulation patterns after sluggish current regimes during glacials. Thus, within the discussion if internal, fluid-driven or external environmentally driven processes control the evolution of such carbonate mounds, our findings for Propeller Mound clearly point to environmental forcing as the dominant mechanism shaping deep-water carbonate mounds in the NE Atlantic during the Late Pleistocene and Holocene.