Commissioning an Anthropomorphic Spine and Lung Phantom for Remote Dose Verification of Institutions Participating in RTOG 0631

The RPC developed a new phantom to ensure comparable and consistent radiation administration in spinal radiosurgery clinical trials. This study assessed the phantom’s dosimetric and anatomic utility. The ‘spine phantom’ is a water filled thorax with anatomy encountered in spinal radiosurgery: target volume, vertebral column, spinal canal, esophagus, heart, and lungs. The dose to the target volume was measured with axial and sagittal planes of radiochromic film and thermoluminescent dosimeters (TLD). The dose distributions were measured with the radiochromic film calibrated to the absolute dose measured by the TLD. Four irradiations were administered: a four angle box plan, a seven angle conformal plan, a seven angle IMRT plan, and a nine angle IMRT plan (denoted as IMRT plan #1 and plan #2, respectively). In each plan, at least 95% of the defined tumor volume received 8 Gy. For each irradiation the planned and administered dose distributions were registered via pinpricks, and compared using point dose measurements, dose profiles, isodose distributions, and gamma analyses. Based on previous experience at the RPC, a gamma analysis was considering passing if greater than 95% of pixels passed the criteria of 5% dose difference and 3 mm distance-to-agreement. Each irradiation showed acceptable agreement in the qualitative assessments and exceeded the 95% passing rate at the 5% / 3 mm criteria, except IMRT plan #1, which was determined to have been poorly localized during treatment administration. The measured and planned dose distributions demonstrated acceptable agreement at the 5% / 3 mm criteria, and the spine phantom was determined to be a useful tool for the remote assessment of an institution’s treatment planning and dose delivery regimen.

IMRT credentialing for prospective trials using institutional virtual phantoms: results of a joint European Organization for the Research and Treatment of Cancer and Radiological Physics Center project.

BACKGROUND AND PURPOSE Intensity-modulated radiotherapy (IMRT) credentialing for a EORTC study was performed using an anthropomorphic head phantom from the Radiological Physics Center (RPC; RPC(PH)). Institutions were retrospectively requested to irradiate their institutional phantom (INST(PH)) using the same treatment plan in the framework of a Virtual Phantom Project (VPP) for IMRT credentialing. MATERIALS AND METHODS CT data set of the institutional phantom and measured 2D dose matrices were requested from centers and sent to a dedicated secure EORTC uploader. Data from the RPC(PH) and INST(PH) were thereafter centrally analyzed and inter-compared by the QA team using commercially available software (RIT; ver.5.2; Colorado Springs, USA). RESULTS Eighteen institutions participated to the VPP. The measurements of 6 (33%) institutions could not be analyzed centrally. All other centers passed both the VPP and the RPC ±7%/4 mm credentialing criteria. At the 5%/5 mm gamma criteria (90% of pixels passing), 11(92%) as compared to 12 (100%) centers pass the credentialing process with RPC(PH) and INST(PH) (p = 0.29), respectively. The corresponding pass rate for the 3%/3 mm gamma criteria (90% of pixels passing) was 2 (17%) and 9 (75%; p = 0.01), respectively. CONCLUSIONS IMRT dosimetry gamma evaluations in a single plane for a H&N prospective trial using the INST(PH) measurements showed agreement at the gamma index criteria of ±5%/5 mm (90% of pixels passing) for a small number of VPP measurements. Using more stringent, criteria, the RPC(PH) and INST(PH) comparison showed disagreement. More data is warranted and urgently required within the framework of prospective studies.

Coded Aperture Flow

Real cameras have a limited depth of field. The resulting defocus blur is a valuable cue for estimating the depth structure of a scene. Using coded apertures, depth can be estimated from a single frame. For optical flow estimation between frames, however, the depth dependent degradation can introduce errors. These errors are most prominent when objects move relative to the focal plane of the camera. We incorporate coded aperture defocus blur into optical flow estimation and allow for piecewise smooth 3D motion of objects. With coded aperture flow, we can establish dense correspondences between pixels in succeeding coded aperture frames. We compare several approaches to compute accurate correspondences for coded aperture images showing objects with arbitrary 3D motion.

Localization and Segmentation of 3D Intervertebral Discs in MR Images by Data Driven Estimation

This paper addresses the problem of fully-automatic localization and segmentation of 3D intervertebral discs (IVDs) from MR images. Our method contains two steps, where we first localize the center of each IVD, and then segment IVDs by classifying image pixels around each disc center as foreground (disc) or background. The disc localization is done by estimating the image displacements from a set of randomly sampled 3D image patches to the disc center. The image displacements are estimated by jointly optimizing the training and test displacement values in a data-driven way, where we take into consideration both the training data and the geometric constraint on the test image. After the disc centers are localized, we segment the discs by classifying image pixels around disc centers as background or foreground. The classification is done in a similar data-driven approach as we used for localization, but in this segmentation case we are aiming to estimate the foreground/background probability of each pixel instead of the image displacements. In addition, an extra neighborhood smooth constraint is introduced to enforce the local smoothness of the label field. Our method is validated on 3D T2-weighted turbo spin echo MR images of 35 patients from two different studies. Experiments show that compared to state of the art, our method achieves better or comparable results. Specifically, we achieve for localization a mean error of 1.6-2.0 mm, and for segmentation a mean Dice metric of 85%-88% and a mean surface distance of 1.3-1.4 mm.

Shaping the PSF to nearly top-hat profile: CHEOPS laboratory results

Spreading the PSF over a quite large amount of pixels is an increasingly used observing technique in order to reach extremely precise photometry, such as in the case of exoplanets searching and characterization via transits observations. A PSF top-hat profile helps to minimize the errors contribution due to the uncertainty on the knowledge of the detector flat field. This work has been carried out during the recent design study in the framework of the ESA small mission CHEOPS. Because of lack of perfect flat-fielding information, in the CHEOPS optics it is required to spread the light of a source into a well defined angular area, in a manner as uniform as possible. Furthermore this should be accomplished still retaining the features of a true focal plane onto the detector. In this way, for instance, the angular displacement on the focal plane is fully retained and in case of several stars in a field these look as separated as their distance is larger than the spreading size. An obvious way is to apply a defocus, while the presence of an intermediate pupil plane in the Back End Optics makes attractive to introduce here an optical device that is able to spread the light in a well defined manner, still retaining the direction of the chief ray hitting it. This can be accomplished through an holographic diffuser or through a lenslet array. Both techniques implement the concept of segmenting the pupil into several sub-zones where light is spread to a well defined angle. We present experimental results on how to deliver such PSF profile by mean of holographic diffuser and lenslet array. Both the devices are located in an intermediate pupil plane of a properly scaled laboratory setup mimicking the CHEOPS optical design configuration. © (2014) COPYRIGHT Society of Photo-Optical Instrumentation Engineers (SPIE). Downloading of the abstract is permitted for personal use only.

RETICULAR PSEUDODRUSEN ON INFRARED IMAGING ARE TOPOGRAPHICALLY DISTINCT FROM SUBRETINAL DRUSENOID DEPOSITS ON EN FACE OPTICAL COHERENCE TOMOGRAPHY.

PURPOSE: To evaluate the quantitative and topographic relationship between reticular pseudodrusen (RPD) on infrared reflectance (IR) and subretinal drusenoid deposits (SDD) on en face volumetric spectral domain optical coherence tomography. METHODS: Reticular pseudodrusen were marked on IR images by a masked observer. Subretinal drusenoid deposits were visualized on en face sections of spectral domain optical coherence tomography below the external limiting membrane and identified by a semiautomated technique. Control RPD lesions were generated in a random distribution for each IR image. Binary maps of control and experimental RPD and SDD were merged and analyzed in terms of topographic localization and quantitative drusen load comparison. RESULTS: A total of 54 eyes of 41 patients diagnosed with RPD were included in this study. The average number of RPD lesions on IR images was 320 ± 44.62 compared with 127 ± 26.02 SDD lesions on en face (P < 0.001). The majority of RPD lesions did not overlap with SDD lesions and were located >30 μm away (92%). The percentage of total SDD lesions overlapping RPD was 2.91 ± 0.87% compared with 1.73 ± 0.68% overlapping control RPD lesions (P < 0.05). The percentage of total SDD lesions between 1 and 3 pixels of the nearest RPD lesion was 5.08 ± 1.40% compared with 3.33 ± 1.07% between 1 and 3 pixels of the nearest control RPD lesion (P < 0.05). CONCLUSION: This study identified significantly more RPD lesions on IR compared with SDD lesions on en face spectral domain optical coherence tomography and found that a large majority of SDD (>90% of lesions) were >30 μm away from the nearest RPD. Together, our findings indicate that RPD and SDD are two entities that are only occasionally topographically associated, suggesting that at some stage in their development, they may be pathologically related.

Temporally Coherent Disparity Maps Using CRFs with Fast 4D Filtering

State of the art methods for disparity estimation achieve good results for single stereo frames, but temporal coherence in stereo videos is often neglected. In this paper we present a method to compute temporally coherent disparity maps. We define an energy over whole stereo sequences and optimize their Conditional Random Field (CRF) distributions using mean-field approximation. We introduce novel terms for smoothness and consistency between the left and right views, and perform CRF optimization by fast, iterative spatio-temporal filtering with linear complexity in the total number of pixels. Our results rank among the state of the art while having significantly less flickering artifacts in stereo sequences.

Regularizing Image Reconstruction for Gradient-Domain Rendering with Feature Patches

We present a novel algorithm to reconstruct high-quality images from sampled pixels and gradients in gradient-domain rendering. Our approach extends screened Poisson reconstruction by adding additional regularization constraints. Our key idea is to exploit local patches in feature images, which contain per-pixels normals, textures, position, etc., to formulate these constraints. We describe a GPU implementation of our approach that runs on the order of seconds on megapixel images. We demonstrate a significant improvement in image quality over screened Poisson reconstruction under the L1 norm. Because we adapt the regularization constraints to the noise level in the input, our algorithm is consistent and converges to the ground truth.

The Effect of Broadleaf Canopies on Survey-grade Horizontal GPS/GLONASS Measurements

A study was conducted to empirically determine the degradation of survey-grade GPS horizontal position measurements due to the effects of broadleaf forest canopies. The measurements were taken using GPS/GLONASS-capable receivers measuring C/A and P-codes, and carrier phase. Fourteen survey markers were chosen in central Connecticut to serve as reference markers for the study. These markers had varying degrees of sky obstruction due to overhanging tree canopies. Sky obstruction was measured by photographing the sky with a 35mm reflex camera fitted with a hemispherical lens. The negative was scanned and the image mapped using an equal- area projection to remove the distortion caused by the lens. The resulting digital image was thresholded to produce a black-and-white image in which a count of the black pixels is a measure of sky-area obstruction. The locations of the markers were determined independently before the study. During the study, each marker was occupied for four 20-minute sessions over the period of one week in mid-July, 1999. The location of the study markers produced relatively long baselines, as compared with similar studies. We compared the accuracy of GPS-only vs. GPS&GLONASS as a function of sky obstruction. Based on our results, GLONASS observations did not improve or degrade the accuracy of the position measurements. There is a loss of 2mm of accuracy per percent of sky obstruction for both GPS only and GPS&GLONASS.

The design and evaluation of a 2D dose verification system for intensity modulated radiotherapy

The usage of intensity modulated radiotherapy (IMRT) treatments necessitates a significant amount of patient-specific quality assurance (QA). This research has investigated the precision and accuracy of Kodak EDR2 film measurements for IMRT verifications, the use of comparisons between 2D dose calculations and measurements to improve treatment plan beam models, and the dosimetric impact of delivery errors. New measurement techniques and software were developed and used clinically at M. D. Anderson Cancer Center. The software implemented two new dose comparison parameters, the 2D normalized agreement test (NAT) and the scalar NAT index. A single-film calibration technique using multileaf collimator (MLC) delivery was developed. EDR2 film's optical density response was found to be sensitive to several factors: radiation time, length of time between exposure and processing, and phantom material. Precision of EDR2 film measurements was found to be better than 1%. For IMRT verification, EDR2 film measurements agreed with ion chamber results to 2%/2mm accuracy for single-beam fluence map verifications and to 5%/2mm for transverse plane measurements of complete plan dose distributions. The same system was used to quantitatively optimize the radiation field offset and MLC transmission beam modeling parameters for Varian MLCs. While scalar dose comparison metrics can work well for optimization purposes, the influence of external parameters on the dose discrepancies must be minimized. The ability of 2D verifications to detect delivery errors was tested with simulated data. The dosimetric characteristics of delivery errors were compared to patient-specific clinical IMRT verifications. For the clinical verifications, the NAT index and percent of pixels failing the gamma index were exponentially distributed and dependent upon the measurement phantom but not the treatment site. Delivery errors affecting all beams in the treatment plan were flagged by the NAT index, although delivery errors impacting only one beam could not be differentiated from routine clinical verification discrepancies. Clinical use of this system will flag outliers, allow physicists to examine their causes, and perhaps improve the level of agreement between radiation dose distribution measurements and calculations. The principles used to design and evaluate this system are extensible to future multidimensional dose measurements and comparisons. ^

Design, Calibration, and Evaluation of Depth-of-Interaction-Capable PET Detector Modules

High-resolution, small-bore PET systems suffer from a tradeoff between system sensitivity, and image quality degradation. In these systems long crystals allow mispositioning of the line of response due to parallax error and this mispositioning causes resolution blurring, but long crystals are necessary for high system sensitivity. One means to allow long crystals without introducing parallax errors is to determine the depth of interaction (DOI) of the gamma ray interaction within the detector module. While DOI has been investigated previously, newly available solid state photomultipliers (SSPMs) well-suited to PET applications and allow new modules for investigation. Depth of interaction in full modules is a relatively new field, and so even if high performance DOI capable modules were available, the appropriate means to characterize and calibrate the modules are not. This work presents an investigation of DOI capable arrays and techniques for characterizing and calibrating those modules. The methods introduced here accurately and reliably characterize and calibrate energy, timing, and event interaction positioning. Additionally presented is a characterization of the spatial resolution of DOI capable modules and a measurement of DOI effects for different angles between detector modules. These arrays have been built into a prototype PET system that delivers better than 2.0 mm resolution with a single-sided-stopping-power in excess of 95% for 511 keV g's. The noise properties of SSPMs scale with the active area of the detector face, and so the best signal-to-noise ratio is possible with parallel readout of each SSPM photodetector pixel rather than multiplexing signals together. This work additionally investigates several algorithms for improving timing performance using timing information from multiple SSPM pixels when light is distributed among several photodetectors.

(Table 1) Effects of desiccation on photochemical processes and nitrogenase activity in Nostoc commune s.l. colonies

The effects of desiccation on photochemical processes and nitrogenase activity were evaluated in Nostoc commune s.l. colonies in situ from a wet thufur meadow at Petuniabukta, Billefjorden, Central Svalbard, during the 2009 arctic summer. The colonies were collected in the fully hydrated state, and were subjected to slow desiccation at ambient temperatures (5 - 8°C) and low light (30 - 80 µmol/m**2/s). For each colony the weight, area, photochemical performance, and nitrogenase activity were determined at the beginning, as well as on every day during the first four days of the experiment; thereafter, on every second day until desiccation was complete. The photochemical performance was evaluated from variable chlorophyll fluorescence parameters (FV/FM, Phi(PSII) , qP, and NPQ), and the nitrogenase activity was estimated by an acetylene-ethylene reduction assay. A significant decrease in the photochemically active area was recorded from the third day, when the colony had lost approximately 40% of its original weight indicating some changes in the extracellular matrix, and stopped on the 14th to 18th day. No effects of the desiccation on the main photochemical parameters (FV/FM, Phi(PSII), qP) were observed up to the sixth to eighth days of desiccation. Slightly lower values of FV/FM and Phi(PSII) recorded in fully-hydrated colonies could be caused by impaired diffusion of CO2 into cells. The steep reduction of photochemical activity occurred between the eighth and tenth day of the experiment, when the colony had lost approximately 80% of its fully-hydrated weight. The nitrogenase activity was highest on the first day, probably due to improved diffusion of N2 into cells, then declined, but was detectable until the sixth day of the experiment. Since Nostoc commune s.l. colonies were capable of photosynthesis and nitrogen fixation to the level of ca. 60% of its fully-hydrated weight, even partly-hydrated colonies contribute substantially to carbon and nitrogen cycling in the High Arctic wet meadow tundra ecosystem.