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

Purpose: Flat-detector, 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. Methods: The rich sources of prior information in IGRT are incorporated into a hidden Markov random field (MRF) 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 (OAR). The voxel labels are estimated using the iterated conditional modes (ICM) algorithm. Results: The accuracy of the method has been evaluated using 27 CBCT scans of an electron density phantom (CIRS, Inc. model 062). The mean voxel-wise misclassification rate was 6.2%, with Dice similarity coefficient of 0.73 for liver, muscle, breast and adipose tissue. Conclusions: 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.

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.

Cone signals and activity in myopia and emmetropia

Myopia (short-sightedness) is a common ocular disorder of children and young adults. Studies primarily using animal models have shown that the retina controls eye growth and the outer retina is likely to have a key role. One theory is that the proportion of L (long-wavelength-sensitive) and M (medium-wavelength-sensitive) cones is related to myopia development; with a high L/M cone ratio predisposing individuals to myopia. However, not all dichromats (persons with red-green colour vision deficiency) with extreme L/M cone ratios have high refractive errors. We predict that the L/M cone ratio will vary in individuals with normal trichromatic colour vision but not show a systematic difference simply due to refractive error. The aim of this study was to determine if L/M cone ratios in the central 30° are different between myopic and emmetropic young, colour normal adults. Information about L/M cone ratios was determined using the multifocal visual evoked potential (mfVEP). The mfVEP can be used to measure the response of visual cortex to different visual stimuli. The visual stimuli were generated and measurements performed using the Visual Evoked Response Imaging System (VERIS 5.1). The mfVEP was measured when the L and M cone systems were separately stimulated using the method of silent substitution. The method of silent substitution alters the output of three primary lights, each with physically different spectral distributions to control the excitation of one or more photoreceptor classes without changing the excitation of the unmodulated photoreceptor classes. The stimulus was a dartboard array subtending 30° horizontally and 30° vertically on a calibrated LCD screen. The m-sequence of the stimulus was 215-1. The N1-P1 amplitude ratio of the mfVEP was used to estimate the L/M cone ratio. Data were collected for 30 young adults (22 to 33 years of age), consisting of 10 emmetropes (+0.3±0.4 D) and 20 myopes (–3.4±1.7 D). The stimulus and analysis techniques were confirmed using responses of two dichromats. For the entire participant group, the estimated central L/M cone ratios ranged from 0.56 to 1.80 in the central 3°-13° diameter ring and from 0.94 to 1.91 in the more peripheral 13°-30° diameter ring. Within 3°-13°, the mean L/M cone ratio of the emmetropic group was 1.20±0.33 and the mean was similar, 1.20±0.26, for the myopic group. For the 13°-30° ring, the mean L/M cone ratio of the emmetropic group was 1.48±0.27 and it was slightly lower in the myopic group, 1.30±0.27. Independent-samples t-test indicated no significant difference between the L/M cone ratios of the emmetropic and myopic group for either the central 3°-13° ring (p=0.986) or the more peripheral 13°-30° ring (p=0.108). The similar distributions of estimated L/M cone ratios in the sample of emmetropes and myopes indicates that there is likely to be no association between the L/M cone ratio and refractive error in humans.

Evaluation of MegaVoltage Cone Beam CT image quality with an unmodified Elekta Precise Linac and EPID : a feasibility study

In order to increase the accuracy of patient positioning for complex radiotherapy treatments various 3D imaging techniques have been developed. MegaVoltage Cone Beam CT (MVCBCT) can utilise existing hardware to implement a 3D imaging modality to aid patient positioning. MVCBCT has been investigated using an unmodified Elekta Precise linac and 15 iView amorphous silicon electronic portal imaging device (EPID). Two methods of delivery and acquisition have been investigated for imaging an anthropomorphic head phantom and quality assurance phantom. Phantom projections were successfully acquired and CT datasets reconstructed using both acquisition methods. Bone, tissue and air were 20 clearly resolvable in both phantoms even with low dose (22 MU) scans. The feasibility of MegaVoltage Cone beam CT was investigated using a standard linac, amorphous silicon EPID and a combination of a free open source reconstruction toolkit as well as custom in-house software written in Matlab. The resultant image quality has 25 been assessed and presented. Although bone, tissue and air were resolvable 2 in all scans, artifacts are present and scan doses are increased when compared with standard portal imaging. The feasibility of MVCBCT with unmodified Elekta Precise linac and EPID has been considered as well as the identification of possible areas for future development in artifact correction techniques to 30 further improve image quality.

Assessment of cone beam CT registration for prostate radiation therapy : fiducial marker and soft tissue methods

Introduction This investigation aimed to assess the consistency and accuracy of radiation therapists (RTs) performing cone beam computed tomography (CBCT) alignment to fiducial markers (FMs) (CBCTFM) and the soft tissue prostate (CBCTST). Methods Six patients receiving prostate radiation therapy underwent daily CBCTs. Manual alignment of CBCTFM and CBCTST was performed by three RTs. Inter-observer agreement was assessed using a modified Bland–Altman analysis for each alignment method. Clinically acceptable 95% limits of agreement with the mean (LoAmean) were defined as ±2.0 mm for CBCTFM and ±3.0 mm for CBCTST. Differences between CBCTST alignment and the observer-averaged CBCTFM (AvCBCTFM) alignment were analysed. Clinically acceptable 95% LoA were defined as ±3.0 mm for the comparison of CBCTST and AvCBCTFM. Results CBCTFM and CBCTST alignments were performed for 185 images. The CBCTFM 95% LoAmean were within ±2.0 mm in all planes. CBCTST 95% LoAmean were within ±3.0 mm in all planes. Comparison of CBCTST with AvCBCTFM resulted in 95% LoA of −4.9 to 2.6, −1.6 to 2.5 and −4.7 to 1.9 mm in the superior–inferior, left–right and anterior–posterior planes, respectively. Conclusions Significant differences were found between soft tissue alignment and the predicted FM position. FMs are useful in reducing inter-observer variability compared with soft tissue alignment. Consideration needs to be given to margin design when using soft tissue matching due to increased inter-observer variability. This study highlights some of the complexities of soft tissue guidance for prostate radiation therapy.

Dosimetry of cone-defined stereotactic radiosurgery fields with a commercial synthetic diamond detector

Purpose Small field x-ray beam dosimetry is difficult due to a lack of lateral electronic equilibrium, source occlusion, high dose gradients and detector volume averaging. Currently there is no single definitive detector recommended for small field dosimetry. The objective of this work was to evaluate the performance of a new commercial synthetic diamond detector, namely the PTW 60019 microDiamond, for the dosimetry of small x-ray fields as used in stereotactic radiosurgery (SRS). Methods Small field sizes were defined by BrainLAB circular cones (4 – 30 mm diameter) on a Novalis Trilogy linear accelerator and using the 6 MV SRS x-ray beam mode for all measurements. Percentage depth doses were measured and compared to an IBA SFD and a PTW 60012 E diode. Cross profiles were measured and compared to an IBA SFD diode. Field factors, Ω_(Q_clin,Q_msr)^(f_clin,f_msr ), were calculated by Monte Carlo methods using BEAMnrc and correction factors, k_(Q_clin,Q_msr)^(f_clin,f_msr ), were derived for the PTW 60019 microDiamond detector. Results For the small fields of 4 to 30 mm diameter, there were dose differences in the PDDs of up to 1.5% when compared to an IBA SFD and PTW 60012 E diode detector. For the cross profile measurements the penumbra values varied, depending upon the orientation of the detector. The field factors, Ω_(Q_clin,Q_msr)^(f_clin,f_msr ), were calculated for these field diameters at a depth of 1.4 cm in water and they were within 2.7% of published values for a similar linear accelerator. The corrections factors, k_(Q_clin,Q_msr)^(f_clin,f_msr ), were derived for the PTW 60019 microDiamond detector. Conclusions We conclude that the new PTW 60019 microDiamond detector is generally suitable for relative dosimetry in small 6 MV SRS beams for a Novalis Trilogy linear equipped with circular cones.

Effect of rod-cone interactions on mesopic visual performance mediated by chromatic and luminance pathways

We studied the effect of rod–cone interactions on mesopic visual reaction time (RT). Rod and cone photoreceptor excitations were independently controlled using a four-primary photostimulator. It was observed that (1) lateral rod–cone interactions increase the cone-mediated RTs; (2) the rod–cone interactions are strongest when rod sensitivity is maximal in a dark surround, but weaker with increased rod activity in a light surround; and (3) the presence of a dark surround nonselectively increased the mean and variability of chromatic (+L-M, S-cone) and luminance (L+M+S) RTs independent of the level of rod activity. The results demonstrate that lateral rod–cone interactions must be considered when deriving mesopic luminous efficiency using RT.

Assessing rod, cone and melanopsin contributions to human pupil flicker responses

Purpose: To determine the relative contributions of rods, cones and melanopsin to pupil responses in humans using temporal sinusoidal stimulation for light levels spanning the low mesopic to photopic range. Methods: A four-primary Ganzfeld photostimulator controlled flicker stimulations at seven light levels (-2.7 to 2 log cd/m2) and five frequencies (0.5 to 8Hz). Pupil diameter was measured using a high-resolution eyetracker. Three kinds of sinusoidal photoreceptor modulations were generated using silent substitution: 1) rod modulation, 2) cone modulation, and 3) combined rod and cone modulation in phase (Experiment 1) or phase shifted (Experiment 2) from a fixed rod phase. The melanopsin excitation was computed for each condition. A vector sum model was used to estimate the relative contribution of rods, cones and melanopsin to the pupil response. Results: From Experiment 1, the pupil frequency response peaked at 1Hz at two mesopic light levels for the three modulation conditions. Analyzing the rod-cone phase difference for the combined modulations (Experiment 2) identified a V-shaped response amplitude with a minimum between 135° and 180°. The pupil response phases increased as cone modulation phase increased. The pupil amplitude increased with increasing light level for cone and combined in-phase rod and cone modulation, but not for the rod modulation. Conclusions: These results demonstrate that cone- and rod-pathway contributions are more predominant than melanopsin contribution to the phasic pupil response. The combined rod, cone and melanopsin inputs to the phasic state of the pupil light reflex follow linear summation.

Cone Ratios in Myopia and Emmetropia: A Pilot Study

Purpose There is a suggestion that the long wavelength-sensitive (LWS)-to-middle wavelength-sensitive (MWS) cone ratio in the retina is associated with myopia. The aim was to measure the LWS/MWS amplitude modulation ratio, an estimate of the LWS/MWS cone ratio, in young adult emmetropes and myopes. Methods Multifocal visual evoked potentials were measured when the LWS and MWS cone systems were excited separately using the method of silent substitution. The 30 young adult participants (22 to 33 years) included 10 emmetropes (mean [±SD] refraction, +0.3 [±0.4] diopters [D]) and 20 myopes (mean [±SD] refraction, -3.4 [±1.7] D). Results The LWS/MWS amplitude modulation ratios ranged from 0.56 to 1.80 in the central 3- to 13-degree diameter ring and from 0.94 to 1.91 in the peripheral 13- to 30-degree diameter ring. Within the central ring, the mean (±SD) ratios were 1.20 (±0.26) and 1.20 (±0.33) for the emmetropic and the myopic groups, respectively. For the peripheral ring, the mean (±SD) ratios were 1.48 (±0.27) and 1.30 (±0.27), respectively. There were no significant differences in the ratios between the emmetropic and myopic groups for either the central (p = 0.99) or peripheral (p = 0.08) rings. For the latter, more myopic refractive error was associated with lower LWS/MWS amplitude modulation ratio; the refraction explained 16% (p = 0.02) of variation in ratio. Conclusions The relationship between the LWS/MWS amplitude modulation ratios and refraction at 13 to 30 degrees indicates that a large longitudinal study of changes in refraction in persons with known cone ratio is required to determine if a low LWS/MWS cone ratio is associated with myopia development.

An external field prior for the hidden Potts model with application to cone-beam computed tomography

In images with low contrast-to-noise ratio (CNR), the information gain from the observed pixel values can be insufficient to distinguish foreground objects. A Bayesian approach to this problem is to incorporate prior information about the objects into a statistical model. A method for representing spatial prior information as an external field in a hidden Potts model is introduced. This prior distribution over the latent pixel labels is a mixture of Gaussian fields, centred on the positions of the objects at a previous point in time. It is particularly applicable in longitudinal imaging studies, where the manual segmentation of one image can be used as a prior for automatic segmentation of subsequent images. The method is demonstrated by application to cone-beam computed tomography (CT), an imaging modality that exhibits distortions in pixel values due to X-ray scatter. The external field prior results in a substantial improvement in segmentation accuracy, reducing the mean pixel misclassification rate for an electron density phantom from 87% to 6%. The method is also applied to radiotherapy patient data, demonstrating how to derive the external field prior in a clinical context.

Cone location and correction of keratoconus with rigid gas-permeable contact lenses

Purpose To evaluate the influence of cone location and corneal cylinder on RGP corrected visual acuities and residual astigmatism in patients with keratoconus. Methods In this prospective study, 156 eyes from 134 patients were enrolled. Complete ophthalmologic examination including manifest refraction, Best spectacle visual acuity (BSCVA), slit-lamp biomicroscopy was performed and corneal topography analysis was done. According to the cone location on the topographic map, the patients were divided into central and paracentral cone groups. Trial RGP lenses were selected based on the flat Sim K readings and a ‘three-point touch’ fitting approach was used. Over contact lens refraction was performed, residual astigmatism (RA) was measured and best-corrected RGP visual acuities (RGPVA) were recorded. Results The mean age (±SD) was 22.1 ± 5.3 years. 76 eyes (48.6%) had central and 80 eyes (51.4%) had paracentral cone. Prior to RGP lenses fitting mean (±SD) subjective refraction spherical equivalent (SRSE), subjective refraction astigmatism (SRAST) and BSCVA (logMAR) were −5.04 ± 2.27 D, −3.51 ± 1.68 D and 0.34 ± 0.14, respectively. There were statistically significant differences between central and paracentral cone groups in mean values of SRSE, SRAST, flat meridian (Sim K1), steep meridian (Sim K2), mean K and corneal cylinder (p-values < 0.05). Comparison of BSCVA to RGPVA shows that vision has improved 0.3 logMAR by RGP lenses (p < 0.0001). Mean (±SD) RA was −0.72 ± 0.39 D. There were no statistically significant differences between RGPVAs and RAs of central and paracentral cone groups (p = 0.22) and (p = 0.42), respectively. Pearson's correlation analysis shows that there is a statistically significant relationship between corneal cylinder and BSCVA and RGPVA, However, the relationship between corneal cylinder and residual astigmatism was not significant. Conclusions Cone location has no effect on the RGP corrected visual acuities and residual astigmatism in patients with keratoconus. Corneal cylinder and Sim K values influence RGP-corrected visual acuities but do not influence residual astigmatism.

Selectivity of human retinotopic visual cortex to S-cone-opponent, L/M-cone-opponent and achromatic stimulation

Our aim was to make a quantitative comparison of the response of the different visual cortical areas to selective stimulation of the two different cone-opponent pathways [long- and medium-wavelength (L/M)- and short-wavelength (S)-cone-opponent] and the achromatic pathway under equivalent conditions. The appropriate stimulus-contrast metric for the comparison of colour and achromatic sensitivity is unknown, however, and so a secondary aim was to investigate whether equivalent fMRI responses of each cortical area are predicted by stimulus contrast matched in multiples of detection threshold that approximately equates for visibility, or direct (cone) contrast matches in which psychophysical sensitivity is uncorrected. We found that the fMRI response across the two colour and achromatic pathways is not well predicted by threshold-scaled stimuli (perceptual visibility) but is better predicted by cone contrast, particularly for area V1. Our results show that the early visual areas (V1, V2, V3, VP and hV4) all have robust responses to colour. No area showed an overall colour preference, however, until anterior to V4 where we found a ventral occipital region that has a significant preference for chromatic stimuli, indicating a functional distinction from earlier areas. We found that all of these areas have a surprisingly strong response to S-cone stimuli, at least as great as the L/M response, suggesting a relative enhancement of the S-cone cortical signal. We also identified two areas (V3A and hMT+) with a significant preference for achromatic over chromatic stimuli, indicating a functional grouping into a dorsal pathway with a strong magnocellular input.

MegaVoltage Cone Beam Computed Tomography with a standard medical linear accelerator

Accurate patient positioning is vital for improved clinical outcomes for cancer treatments using radiotherapy. This project has developed Mega Voltage Cone Beam CT using a standard medical linear accelerator to allow 3D imaging of the patient position at treatment time with no additional hardware required. Providing 3D imaging functionality at no further cost allows enhanced patient position verification on older linear accelerators and in developing countries where access to new technology is limited.

Reconstruction from incomplete data in cone-beam tomography

A method for reconstruction of an object f(x) x=(x,y,z) from a limited set of cone-beam projection data has been developed. This method uses a modified form of convolution back-projection and projection onto convex sets (POCS) for handling the limited (or incomplete) data problem. In cone-beam tomography, one needs to have a complete geometry to completely reconstruct the original three-dimensional object. While complete geometries do exist, they are of little use in practical implementations. The most common trajectory used in practical scanners is circular, which is incomplete. It is, however, possible to recover some of the information of the original signal f(x) based on a priori knowledge of the nature of f(x). If this knowledge can be posed in a convex set framework, then POCS can be utilized. In this report, we utilize this a priori knowledge as convex set constraints to reconstruct f(x) using POCS. While we demonstrate the effectiveness of our algorithm for circular trajectories, it is essentially geometry independent and will be useful in any limited-view cone-beam reconstruction.

Miniature cone tip resistance of sand with fly ash using triaxial setup

In a continuation of the authors' recent work, the ultimate tip resistance of a miniature cone using triaxial equipment was determined for samples of dry sand mixed with dry fly ash. The effect of (i) the proportion of fly ash, (ii) the relative density of samples, and (iii) the vertical overburden pressure was examined. It was noted that an addition of fly ash in sand for the same range of relative density leads to a significant reduction in the ultimate tip resistance of the cone (q(cu)). This occurs due to a decrease in the friction angle (phi) of the sample with an increase in the fly ash content for a given relative density. For phi greater than about 30 degrees, two widely used correlation curves from published literature, providing the relationships between q(cu) and phi for cohesionless soils, were found to provide satisfactory predictions, even for sand - fly ash mixtures. As was expected, the values of qcu increase continuously with an increase in the relative density of the soil mass and the vertical effective ( overburden) stress on the sample.