Search of low-contrast liver lesions in abdominal CT: the importance of scrolling behavior

Purpose: Visual search using volumetric images is becoming the standard in medical imaging. However, we do not fully understand how eye movement strategies mediate diagnostic performance. A recent study on computed tomography (CT) images showed that the search strategies of radiologists could be classified based on saccade amplitudes and cross-quadrant eye movements [eye movement index (EMI)] into two categories: drillers and scanners. Approach: We investigate how the number of times a radiologist scrolls in a given direction during analysis of the images (number of courses) could add a supplementary variable to use to characterize search strategies. We used a set of 15 normal liver CT images in which we inserted 1 to 5 hypodense metastases of two different signal contrast amplitudes. Twenty radiologists were asked to search for the metastases while their eye-gaze was recorded by an eye-tracker device (EyeLink1000, SR Research Ltd., Mississauga, Ontario, Canada). Results: We found that categorizing radiologists based on the number of courses (rather than EMI) could better predict differences in decision times, percentage of image covered, and search error rates. Radiologists with a larger number of courses covered more volume in more time, found more metastases, and made fewer search errors than those with a lower number of courses. Our results suggest that the traditional definition of drillers and scanners could be expanded to include scrolling behavior. Drillers could be defined as scrolling back and forth through the image stack, each time exploring a different area on each image (low EMI and high number of courses). Scanners could be defined as scrolling progressively through the stack of images and focusing on different areas within each image slice (high EMI and low number of courses). Conclusions: Together, our results further enhance the understanding of how radiologists investigate three-dimensional volumes and may improve how to teach effective reading strategies to radiology residents.

A delivery quality assurance tool based on the actual leaf open times in tomotherapy

PURPOSE

To validate a delivery quality assurance (DQA) protocol for tomotherapy based on the measurement of the leaf open times (LOTs). In addition, to show the correlation between the mean relative LOT discrepancy and the dose deviation in the planning target volume (PTV).

MATERIALS AND METHODS

We used a LOT measurement algorithm presented in a previous work on our two tomotherapy treatment units (TOMO1 and TOMO2). We generated TomoPhant plans with intentional random LOT discrepancies following Gaussian distributions of -6%, -4%, -2%, 2%, 4%, and 6%. We irradiated them on the Cheese Phantom with two ion chambers and collected the raw data on both our treatment units. Using the raw data, we measured the actual LOTs and verified that the induced discrepancies were highlightable. Then, we calculated the actual dose using Accuray's standalone dose calculator and verified that the calculated dose agreed with the ion chamber measurement. We randomly chose 60 clinical treatment plans, delivered them in air, and collected the raw detector data. We measured the actual LOTs from the raw data and calculated the corresponding dose distributions using Accuray's standalone dose calculator. We assessed the Pearson coefficient correlation of the deviation between expected and actual dose in the PTV (a) with the mean relative LOT discrepancy and (b) with the γ-index pass rate for different tolerances.

RESULTS

The mean relative discrepancy between actual (measured by the algorithm) and expected LOTs on the modified TomoPhant plans was 1.10 ± 0.05% on TOMO1 and 0.02 ± 0.03% on TOMO2, respectively. The agreement between measured and calculated dose was 0.2 ± 0.3% on TOMO1 and 0.1 ± 0.3% on TOMO2, respectively. On clinical plans, the means of the relative LOT discrepancies ranged from -3.0 % to 1.4%. The dose deviation in the PTVs ranged from -1.6% to 2.4%. The Pearson coefficient correlation between the mean relative LOT discrepancy and the dose deviation in the PTV was 0.76 (P ≈ 10^-15 ) on TOMO1 and 0.65 (P ≈ 10^-10 ) on TOMO2, respectively. There was no correlation between the γ-index pass rate and the dose deviation in the PTV.

CONCLUSION

The method made it possible to measure and to correctly highlight the LOT discrepancies on the TomoPhant plans. The dose subsequently calculated was accurate. On clinical plans, the mean LOT discrepancy correlated with the dose deviation in the PTV. This makes the mean LOT discrepancy a handy indicator of the plan quality.

Measurement of the useful field of view for single slices of different imaging modalities and targets

Purpose: With three-dimensional (3-D) images displayed as stacks of 2-D images, radiologists rely more heavily on vision away from their fixation point to visually process information, guide eye movements, and detect abnormalities. Thus the ability to detect targets away from the fixation point, commonly characterized as the useful field of view (UFOV), becomes critical for these 3-D imaging modalities. We investigate how the UFOV, defined as the eccentricity, in which detection performance degrades to a given probability, varies across imaging modalities and targets. Approach: We measure the detectability of different targets at various distances from gaze locations for single slices of liver computed tomography (CT), 2-D digital mammograms (DM), and single slices of digital breast tomosynthesis (DBT) cases. Observers with varying expertise were instructed to maintain their gaze at a point while a short display of the image was flashed and an eye tracker verified observer's steady fixation. Display times were 200 and 1000 ms for CT images and 500 ms for DM and DBT images. Results: We find variations in the UFOV from 9 to 12 deg for liver CT to as small as 2.5 to 5 deg for calcification clusters in breast images (DM and DBT). We compare our results to those reported in the literature for lung nodules and discuss the differences across methods used to measure the UFOV, their dependence on case selection/task difficulty, viewing conditions, and observer expertise. We propose a complementary measure defined in terms of performance degradation relative to the peak foveal performance (relative-UFOV) to circumvent UFOV's variations with case selection/task difficulty. Conclusion: Our results highlight the variations in the UFOV across imaging modalities, target types, observer expertise, and measurement methods and suggest an additional relative-UFOV measure to more thoroughly characterize the detection performance away from point of fixation.

In air and in vivo measurement of the leaf open time in tomotherapy using the on-board detector pulse-by-pulse data

PURPOSE

We developed an algorithm to measure the leaf open times (LOT) from the on-board detector (OBD) pulse-by-pulse data in tomotherapy. We assessed the feasibility of measuring the LOTs in dynamic jaw mode and validated the algorithm on machine QA and clinical data. Knowledge of the actual LOTs is a basis toward calculating the delivered dose and performing efficient phantom-less delivery quality assurance (DQA) controls of the multileaf collimator (MLC). In tomotherapy, the quality of the delivered dose depends on the correct performance of the MLC, hence on the accuracy of the LOTs.

MATERIALS AND METHODS

In the detector signal, the period of time during which a leaf is open corresponds to a high intensity region. The algorithm described here locally normalizes the detector signal and measures the FWHM of the high intensity regions. The Daily QA module of the TomoTherapy Quality Assurance (TQA) tool measures LOT errors. The Daily QA detector data were collected during 9 days on two tomotherapy units. The errors yielded by the method were compared to these reported by the Daily QA module. In addition, clinical data were acquired on the two units (25 plans in total), in air without attenuation material in the beam path and in vivo during a treatment fraction. The study included plans with fields of all existing sizes (1.05, 2.51, 5.05 cm). The collimator jaws were in dynamic mode (TomoEDGE^TM ). The feasibility of measuring the LOTs was assessed with respect to the jaw aperture.

RESULTS

The mean discrepancy between LOTs measured by the algorithm and those measured by TQA was of 0 ms, with a standard deviation of 0.3 ms. The LOT measured by the method had thus an uncertainty of 1 ms with a confidence level of 99%. In 5.05 cm dynamic jaw procedures, the detector is in the beam umbra at the beginning and at the end of the delivery. In such procedures, the algorithm could not measure the LOTs at jaw apertures between 7 and maximum 12.4 mm. Otherwise, no measurement error due to the jaw movement was observed. No LOT measurement difference between air and in vivo data was observed either.

CONCLUSION

The method we propose is reliable. It can equivalently measure the LOTs from data acquired in air or in vivo. It handles fully the static procedures and the 2.51 cm dynamic procedures. It handles partially the 5.05 cm dynamic procedures. The limitation was evaluated with respect to the jaw aperture.

Objective comparison of high-contrast spatial resolution and low-contrast detectability for various clinical protocols on multiple CT scanners

PURPOSE

We sought to compare objectively computed tomography (CT) scanner performance for three clinically relevant protocols using a task-based image quality assessment method in order to assess the potential for radiation dose reduction.

METHODS

Four CT scanners released between 2003 and 2007 by different manufacturers were compared with four CT scanners released between 2012 and 2014 by the same manufacturers using ideal linear model observers (MO): prewhitening (PW) MO and channelized Hotelling (CHO) MO with Laguerre-Gauss channels for high-contrast spatial resolution and low-contrast detectability (LCD) performance, respectively. High-contrast spatial resolution was assessed using a custom-made phantom that enabled the computation of the target transfer function (TTF) and noise power spectrum (NPS). Low-contrast detectability was assessed using a commercially available anthropomorphic abdominal phantom providing equivalent diameters of 24, 29.6, and 34.6 cm. Three protocols were reviewed: a head (trauma) and an abdominal (urinary stones) protocol were applied to assess high-contrast spatial resolution performance; and another abdominal (focal liver lesions) protocol was applied for LCD. The liver protocol was tested using fixed and modulated tube currents. The PW MO was proposed for assessing high-contrast detectability performance of the various CT scanners.

RESULTS

Compared with older generation CT scanners, three newer systems displayed significant improvements in high-contrast detectability over that of their predecessors. A fourth, newer system had lower performance. The CHO MO was appropriate for assessing LCD performance and revealed that an excellent level of image quality could be obtained with newer scanners at significantly lower dose levels.

CONCLUSIONS

This study shows that MO can objectively benchmark CT scanners using a task-based image quality method, thus helping to estimate the potential for further dose reductions offered by the latest systems. Such an approach may be useful for adequately and quantitatively comparing clinically relevant image quality among various scanners.

Inter-laboratory comparison of channelized hotelling observer computation

PURPOSE

The task-based assessment of image quality using model observers is increasingly used for the assessment of different imaging modalities. However, the performance computation of model observers needs standardization as well as a well-established trust in its implementation methodology and uncertainty estimation. The purpose of this work was to determine the degree of equivalence of the channelized Hotelling observer performance and uncertainty estimation using an intercomparison exercise.

MATERIALS AND METHODS

Image samples to estimate model observer performance for detection tasks were generated from two-dimensional CT image slices of a uniform water phantom. A common set of images was sent to participating laboratories to perform and document the following tasks: (a) estimate the detectability index of a well-defined CHO and its uncertainty in three conditions involving different sized targets all at the same dose, and (b) apply this CHO to an image set where ground truth was unknown to participants (lower image dose). In addition, and on an optional basis, we asked the participating laboratories to (c) estimate the performance of real human observers from a psychophysical experiment of their choice. Each of the 13 participating laboratories was confidentially assigned a participant number and image sets could be downloaded through a secure server. Results were distributed with each participant recognizable by its number and then each laboratory was able to modify their results with justification as model observer calculation are not yet a routine and potentially error prone.

RESULTS

Detectability index increased with signal size for all participants and was very consistent for 6 mm sized target while showing higher variability for 8 and 10 mm sized target. There was one order of magnitude between the lowest and the largest uncertainty estimation.

CONCLUSIONS

This intercomparison helped define the state of the art of model observer performance computation and with thirteen participants, reflects openness and trust within the medical imaging community. The performance of a CHO with explicitly defined channels and a relatively large number of test images was consistently estimated by all participants. In contrast, the paper demonstrates that there is no agreement on estimating the variance of detectability in the training and testing setting.

Dose assessment following an overexposure of a worker at a Swiss nuclear power plant

AIM

The aim of this work was to assess the doses received by a diver exposed to a radiation source during maintenance work in the fuel transfer pool at a Swiss nuclear power plant, and to define whether the statutory limit was breached or not.

METHOD

Onsite measurements were carried out and different scenarios were simulated using the MicroShield Software and the MCNPX Monte Carlo radiation transport code to estimate the activity of the irradiating object as well as the doses to the limbs and the effective dose delivered to the operator.

RESULTS

The activity of the object was estimated to 1.8 TBq. From the various dose estimations, a conservative value of 7.5 Sv was proposed for the equivalent dose to the skin on the hands and an effective dose of 28 mSv.

CONCLUSION

The use of different experimental and calculation methods allowed us to accurately estimate the activity of the object and the dose delivered to the diver, useful information for making a decision on the most appropriate scheme of follow up for the patient.

CONVERTING SPECIFIC ACTIVITY INTO AMBIENT DOSE EQUIVALENT: UPDATED COEFFICIENTS FOR IN SITU GAMMA SPECTROMETRY

In situ gamma spectrometry is a valuable tool to assess the radionuclides released in the environment and the associated dose. This requires prior establishment of coefficients allowing the conversion of the specific activity into ambient equivalent dose. The aim of this work is to calculate updated conversion factors for monoenergetic photons and for a series of radionuclides of interest. The calculation was performed using the Monte Carlo (MC) method, the GEANT4 MC code, various activity distribution models and up-to-date nuclear decay data. A new set of conversion factors is established in the energy range extending from  <100 keV to 8.5 MeV. The coefficients calculated in this work were compared to the data published in the literature.

High dose-per-pulse electron beam dosimetry: Usability and dose-rate independence of EBT3 Gafchromic films

PURPOSE

The aim of this study was to assess the suitability of Gafchromic EBT3 films for reference dose measurements in the beam of a prototype high dose-per-pulse linear accelerator (linac), capable of delivering electron beams with a mean dose-rate (Ḋ_m ) ranging from 0.07 to 3000 Gy/s and a dose-rate in pulse (Ḋ_p ) of up to 8 × 10⁶ Gy/s. To do this, we evaluated the overall uncertainties in EBT3 film dosimetry as well as the energy and dose-rate dependence of their response.

MATERIAL AND METHODS

Our dosimetric system was composed of EBT3 Gafchromic films in combination with a flatbed scanner and was calibrated against an ionization chamber traceable to primary standard. All sources of uncertainties in EBT3 dosimetry were carefully analyzed using irradiations at a clinical radiotherapy linac. Energy dependence was investigated with the same machine by acquiring and comparing calibration curves for three different beam energies (4, 8 and 12 MeV), for doses between 0.25 and 30 Gy. Ḋ_m dependence was studied at the clinical linac by changing the pulse repetition frequency (f) of the beam in order to vary Ḋ_m between 0.55 and 4.40 Gy/min, while Ḋ_p dependence was probed at the prototype machine for Ḋ_p ranging from 7 × 10³ to 8 × 10⁶ Gy/s. Ḋ_p dependence was first determined by studying the correlation between the dose measured by films and the charge of electrons measured at the exit of the machine by an induction torus. Furthermore, we compared doses from the films to independently calibrated thermo-luminescent dosimeters (TLD) that have been reported as being dose-rate independent up to such high dose-rates.

RESULTS

We report that uncertainty below 4% (k = 2) can be achieved in the dose range between 3 and 17 Gy. Results also demonstrated that EBT3 films did not display any detectable energy dependence for electron beam energies between 4 and 12 MeV. No Ḋ_m dependence was found either. In addition, we obtained excellent consistency between films and TLDs over the entire Ḋ_p range attainable at the prototype linac confirming the absence of any dose-rate dependence within the investigated range (7 × 10³ to 8 × 10⁶ Gy/s). This aspect was further corroborated by the linear relationship between the dose-per-pulse (D_p ) measured by films and the charge per pulse (C_p ) measured at the prototype linac exit.

CONCLUSION

Our study shows that the use of EBT3 Gafchromic films can be extended to reference dosimetry in pulsed electron beams with a very high dose rate. The measurement results are associated with an overall uncertainty below 4% (k = 2) and are dose-rate and energy independent.

BENCHMARKING OF CT FOR PATIENT EXPOSURE OPTIMISATION

Patient dose optimisation in computed tomography (CT) should be done using clinically relevant tasks when dealing with image quality assessments. In the present work, low-contrast detectability for an average patient morphology was assessed on 56 CT units, using a model observer applied on images acquired with two specific protocols of an anthropomorphic phantom containing spheres. Images were assessed using the channelised Hotelling observer (CHO) with dense difference of Gaussian channels. The results were computed by performing receiver operating characteristics analysis (ROC) and using the area under the ROC curve (AUC) as a figure of merit. The results showed a small disparity at a volume computed tomography dose index (CTDIvol) of 15 mGy depending on the CT units for the chosen image quality criterion. For 8-mm targets, AUCs were 0.999 ± 0.018 at 20 Hounsfield units (HU) and 0.927 ± 0.054 at 10 HU. For 5-mm targets, AUCs were 0.947 ± 0.059 and 0.702 ± 0.068 at 20 and 10 HU, respectively. The robustness of the CHO opens the way for CT protocol benchmarking and optimisation processes.

PATIENT EXPOSURE OPTIMISATION THROUGH TASK-BASED ASSESSMENT OF A NEW MODEL-BASED ITERATIVE RECONSTRUCTION TECHNIQUE

The goal of the present work was to report and investigate the performances of a new iterative reconstruction algorithm, using a model observer. For that, a dedicated low-contrast phantom containing different targets was scanned at four volume computed tomography dose index (CTDI_vol) levels on a Siemens SOMATOM Force computed tomography (CT). The acquired images were reconstructed using the ADMIRE algorithm and were then assessed by three human observers who performed alternative forced choice experiments. Next, a channelised hotelling observer model was applied on the same set of images. The comparison between the two was performed using the percentage correct as a figure of merit. The results indicated a strong agreement between human and model observer as well as an improvement in the low-contrast detection when switching from an ADMIRE strength of 1-3. Good results were also observed even in situations where the target was hard to detect, suggesting that patient dose could be further reduced and optimised.

OBJECTIVE TASK-BASED ASSESSMENT OF LOW-CONTRAST DETECTABILITY IN ITERATIVE RECONSTRUCTION

Evaluating image quality by using receiver operating characteristic studies is time consuming and difficult to implement. This work assesses a new iterative algorithm using a channelised Hotelling observer (CHO). For this purpose, an anthropomorphic abdomen phantom with spheres of various sizes and contrasts was scanned at 3 volume computed tomography dose index (CTDIvol) levels on a GE Revolution CT. Images were reconstructed using the iterative reconstruction method adaptive statistical iterative reconstruction-V (ASIR-V) at ASIR-V 0, 50 and 70 % and assessed by applying a CHO with dense difference of Gaussian and internal noise. Both CHO and human observers (HO) were compared based on a four-alternative forced-choice experiment, using the percentage correct as a figure of merit. The results showed accordance between CHO and HO. Moreover, an improvement in the low-contrast detection was observed when switching from ASIR-V 0 to 50 %. The results underpin the finding that ASIR-V allows dose reduction.

Eye-tracking of nodule detection in lung CT volumetric data

PURPOSE

Signal detection on 3D medical images depends on many factors, such as foveal and peripheral vision, the type of signal, and background complexity, and the speed at which the frames are displayed. In this paper, the authors focus on the speed with which radiologists and naïve observers search through medical images. Prior to the study, the authors asked the radiologists to estimate the speed at which they scrolled through CT sets. They gave a subjective estimate of 5 frames per second (fps). The aim of this paper is to measure and analyze the speed with which humans scroll through image stacks, showing a method to visually display the behavior of observers as the search is made as well as measuring the accuracy of the decisions. This information will be useful in the development of model observers, mathematical algorithms that can be used to evaluate diagnostic imaging systems.

METHODS

The authors performed a series of 3D 4-alternative forced-choice lung nodule detection tasks on volumetric stacks of chest CT images iteratively reconstructed in lung algorithm. The strategy used by three radiologists and three naïve observers was assessed using an eye-tracker in order to establish where their gaze was fixed during the experiment and to verify that when a decision was made, a correct answer was not due only to chance. In a first set of experiments, the observers were restricted to read the images at three fixed speeds of image scrolling and were allowed to see each alternative once. In the second set of experiments, the subjects were allowed to scroll through the image stacks at will with no time or gaze limits. In both static-speed and free-scrolling conditions, the four image stacks were displayed simultaneously. All trials were shown at two different image contrasts.

RESULTS

The authors were able to determine a histogram of scrolling speeds in frames per second. The scrolling speed of the naïve observers and the radiologists at the moment the signal was detected was measured at 25-30 fps. For the task chosen, the performance of the observers was not affected by the contrast or experience of the observer. However, the naïve observers exhibited a different pattern of scrolling than the radiologists, which included a tendency toward higher number of direction changes and number of slices viewed.

CONCLUSIONS

The authors have determined a distribution of speeds for volumetric detection tasks. The speed at detection was higher than that subjectively estimated by the radiologists before the experiment. The speed information that was measured will be useful in the development of 3D model observers, especially anthropomorphic model observers which try to mimic human behavior.

Anthropomorphic model observer performance in three-dimensional detection task for low-contrast computed tomography

X-ray medical imaging is increasingly becoming three-dimensional (3-D). The dose to the population and its management are of special concern in computed tomography (CT). Task-based methods with model observers to assess the dose-image quality trade-off are promising tools, but they still need to be validated for real volumetric images. The purpose of the present work is to evaluate anthropomorphic model observers in 3-D detection tasks for low-contrast CT images. We scanned a low-contrast phantom containing four types of signals at three dose levels and used two reconstruction algorithms. We implemented a multislice model observer based on the channelized Hotelling observer (msCHO) with anthropomorphic channels and investigated different internal noise methods. We found a good correlation for all tested model observers. These results suggest that the msCHO can be used as a relevant task-based method to evaluate low-contrast detection for CT and optimize scan protocols to lower dose in an efficient way.

Monte Carlo simulation of a whole-body counter using IGOR phantoms

Whole-body counting is a technique of choice for assessing the intake of gamma-emitting radionuclides. An appropriate calibration is necessary, which is done either by experimental measurement or by Monte Carlo (MC) calculation. The aim of this work was to validate a MC model for calibrating whole-body counters (WBCs) by comparing the results of computations with measurements performed on an anthropomorphic phantom and to investigate the effect of a change in phantom's position on the WBC counting sensitivity. GEANT MC code was used for the calculations, and an IGOR phantom loaded with several types of radionuclides was used for the experimental measurements. The results show a reasonable agreement between measurements and MC computation. A 1-cm error in phantom positioning changes the activity estimation by >2%. Considering that a 5-cm deviation of the positioning of the phantom may occur in a realistic counting scenario, this implies that the uncertainty of the activity measured by a WBC is ∼10-20%.

Exposure of the Swiss population to computed tomography

Background

The frequency of CT procedures has registered a significant increase over the last decade, which led at the international level to an increasing concern on the radiological risk associated with the use of CT especially in paediatrics. This work aimed at investigating the use of computed tomography in Switzerland, following the evolution of CT frequency and dose data over a decade and comparing it to data reported in other countries.

Methods

The frequency and dose data related to CT are obtained by means of a nationwide survey. National frequencies were established by projecting the collected data, using the ratio of the number of CT units belonging to the respondents to the total number of CT units in the country. The effective doses per examination were collected during an auditing campaign.

Results

In 2008 about 0.8 Million CT procedures (~ 100 CT examinations / 1000 population) were performed in the country, leading to a collective effective dose of more than 6000 man.Sv (0.8 mSv/caput). In a decade the frequency of CT examinations averaged over the population and the associated average effective dose per caput increased by a factor of 2.2 and 2.9 respectively.

Conclusions

Although the contribution of CT to the total medical X-rays is 6% in terms of the frequency, it represents 68% in terms of the collective effective dose. These results are comparable to those reported in a number of countries in Europe and America with similar health level.

Model-based iterative reconstruction in pediatric chest CT: assessment of image quality in a prospective study of children with cystic fibrosis

BACKGROUND

The potential effects of ionizing radiation are of particular concern in children. The model-based iterative reconstruction VEO(TM) is a technique commercialized to improve image quality and reduce noise compared with the filtered back-projection (FBP) method.

OBJECTIVE

To evaluate the potential of VEO(TM) on diagnostic image quality and dose reduction in pediatric chest CT examinations.

MATERIALS AND METHODS

Twenty children (mean 11.4 years) with cystic fibrosis underwent either a standard CT or a moderately reduced-dose CT plus a minimum-dose CT performed at 100 kVp. Reduced-dose CT examinations consisted of two consecutive acquisitions: one moderately reduced-dose CT with increased noise index (NI = 70) and one minimum-dose CT at CTDIvol 0.14 mGy. Standard CTs were reconstructed using the FBP method while low-dose CTs were reconstructed using FBP and VEO. Two senior radiologists evaluated diagnostic image quality independently by scoring anatomical structures using a four-point scale (1 = excellent, 2 = clear, 3 = diminished, 4 = non-diagnostic). Standard deviation (SD) and signal-to-noise ratio (SNR) were also computed.

RESULTS

At moderately reduced doses, VEO images had significantly lower SD (P < 0.001) and higher SNR (P < 0.05) in comparison to filtered back-projection images. Further improvements were obtained at minimum-dose CT. The best diagnostic image quality was obtained with VEO at minimum-dose CT for the small structures (subpleural vessels and lung fissures) (P < 0.001). The potential for dose reduction was dependent on the diagnostic task because of the modification of the image texture produced by this reconstruction.

CONCLUSIONS

At minimum-dose CT, VEO enables important dose reduction depending on the clinical indication and makes visible certain small structures that were not perceptible with filtered back-projection.

Swiss population exposure to radiation by interventional radiology in 2008

The aim of this study was to investigate the radiation exposure of the Swiss population to interventional procedures. A nationwide survey was conducted in Switzerland. The annual effective dose per capita due to interventional procedures was found to be 0.14 mSv, corresponding to 12% of the total dose. Coronary angiography and percutaneous coronary interventions were found to be the most frequent and the most irradiating interventional procedures, accounting for 52% of the total examination frequency and 64% of the dose delivered to the population. Switzerland stands at the same level as other countries in terms of effective dose per capita due to interventional radiology.

Calculation of correction factors for ionization chamber measurements with small fields in low-density media

The quantity of interest for high-energy photon beam therapy recommended by most dosimetric protocols is the absorbed dose to water. Thus, ionization chambers are calibrated in absorbed dose to water, which is the same quantity as what is calculated by most treatment planning systems (TPS). However, when measurements are performed in a low-density medium, the presence of the ionization chamber generates a perturbation at the level of the secondary particle range. Therefore, the measured quantity is close to the absorbed dose to a volume of water equivalent to the chamber volume. This quantity is not equivalent to the dose calculated by a TPS, which is the absorbed dose to an infinitesimally small volume of water. This phenomenon can lead to an overestimation of the absorbed dose measured with an ionization chamber of up to 40% in extreme cases. In this paper, we propose a method to calculate correction factors based on the Monte Carlo simulations. These correction factors are obtained by the ratio of the absorbed dose to water in a low-density medium □D(w,Q,V1)(low) averaged over a scoring volume V₁ for a geometry where V₁ is filled with the low-density medium and the absorbed dose to water □D(w,QV2)(low) averaged over a volume V₂ for a geometry where V₂ is filled with water. In the Monte Carlo simulations, □D(w,QV2)(low) is obtained by replacing the volume of the ionization chamber by an equivalent volume of water, according to the definition of the absorbed dose to water. The method is validated in two different configurations which allowed us to study the behavior of this correction factor as a function of depth in phantom, photon beam energy, phantom density and field size.

Exposure of the Swiss population by medical x-rays: 2008 review

Nationwide surveys on radiation dose to the population from medical radiology are recommended in order to follow the trends in population exposure and ensure radiation protection.The last survey in Switzerland was conducted in 1998,and the annual effective dose from medical radiology was estimated to be 1 mSv y j(-1) per capita. The purpose of this work was to follow the trends in diagnostic radiology between 1998 and 2008 in Switzerland and determine the contribution of different modalities and types of examinations to the collective effective dose from medical x-rays. For this reason, an online database(www.raddose.ch) was developed. All healthcare providers who hold a license to run an x-ray unit in the country were invited to participate in the survey. More than 225 examinations, covering eight radiological modalities, were included in the survey. The average effective dose for each examination was reassessed. Data from about 3,500 users were collected (42% response rate). The survey showed that the annual effective dose was 1.2 mSv/capita in 2008. The most frequent examinations are conventional and dental radiographies (88%). The contribution of computed tomography was only 6% in terms of examination frequency but 68% in terms of effective dose. The comparison with other countries showed that the effective dose per capita in Switzerland was in the same range as in other countries with similar healthcare systems, although the annual number of examinations performed in Switzerland was higher.

Usefulness of specific calibration coefficients for gamma-emitting sources measured by radionuclide calibrators in nuclear medicine

PURPOSE

In nuclear medicine, the activity of a radionuclide is measured with a radionuclide calibrator that often has a calibration coefficient independent of the container type and filling.

METHODS

To determine the effect of the container on the accuracy of measuring the activity injected into a patient, The authors simulated a commercial radionuclide calibrator and 18 container types most typically used in clinical practice. The instrument sensitivity was computed for various container thicknesses and filling levels. Monoenergetic photons and electrons as well as seven common radionuclides were considered.

RESULTS

The quality of the simulation with gamma-emitting sources was validated by an agreement with measurements better than 4% in five selected radionuclides. The results show that the measured activity can vary by more than a factor of 2 depending on the type of container. The filling level and the thickness of the container wall only have a marginal effect for radionuclides of high energy but could induce differences up to 4%.

CONCLUSIONS

The authors conclude that radionuclide calibrators should be tailored to the uncertainty required by clinical applications. For most clinical cases, and at least for the low-energy gamma and x-ray emitters, measurements should be performed with calibration coefficients specific to the container type.

Primary activity measurements with a 4πβ-4πγ coincidence counting system

The radioactive concentrations of (166m)Ho, (134)Cs and (133)Ba solutions have been standardised using a 4πβ-4πγ coincidence counting system we have recently set up. The detection in the beta channel is performed using various geometries of a UPS-89 plastic scintillator optically coupled to a selected low-noise 1in. diameter photomultiplier tube. The light-tight thin capsule that encloses this beta detector is housed within the well of a 5in.×5in. NaI(Tl) monocrystal detector. The beta detection efficiency can be varied either by optical filtering or electronic discrimination when the electrons loose all their energy in the plastic scintillator. This 4πβ-4πγ coincidence system improves on our 4πβ(PC)-γ system in that its sample preparation is less labour intensive, it yields larger beta- and gamma-counting efficiencies thus enabling the standardisation of low activity sources with good statistics in reasonable time, and it makes standardising short-lived radionuclides easier. The resulting radioactive concentrations of (166m)Ho, (134)Cs and (133)Ba are found to agree with those measured with other primary measurement methods thus validating our 4πβ-4πγ coincidence counting system.

Variability of radioiodine measurements in the thyroid

Monte Carlo simulations were carried out to study the response of a thyroid monitor for measuring intake activities of (125)I and (131)I. The aim of the study was 3-fold: to cross-validate the Monte Carlo simulation programs, to study the response of the detector using different phantoms and to study the effects of anatomical variations. Simulations were performed using the Swiss reference phantom and several voxelised phantoms. Determining the position of the thyroid is crucial for an accurate determination of radiological risks. The detector response using the Swiss reference phantom was in fairly good agreement with the response obtained using adult voxelised phantoms for (131)I, but should be revised for a better calibration for (125)I and for any measurements taken on paediatric patients.

Characterisation of the PSI whole body counter by radiographic imaging

A joint project between the Paul Scherrer Institut (PSI) and the Institute of Radiation Physics was initiated to characterise the PSI whole body counter in detail through measurements and Monte Carlo simulation. Accurate knowledge of the detector geometry is essential for reliable simulations of human body phantoms filled with known activity concentrations. Unfortunately, the technical drawings provided by the manufacturer are often not detailed enough and sometimes the specifications do not agree with the actual set-up. Therefore, the exact detector geometry and the position of the detector crystal inside the housing were determined through radiographic images. X-rays were used to analyse the structure of the detector, and (60)Co radiography was employed to measure the core of the germanium crystal. Moreover, the precise axial alignment of the detector within its housing was determined through a series of radiographic images with different incident angles. The hence obtained information enables us to optimise the Monte Carlo geometry model and to perform much more accurate and reliable simulations.

Individual monitoring of internal exposure for nuclear medicine workers in Switzerland

Monitoring of internal exposure for nuclear medicine workers requires frequent measurements due to the short physical half-lives of most radionuclides used in this field. The aim of this study was to develop screening measurements performed at the workplace by local staff using standard laboratory instrumentation, to detect whether potential intake has occurred. Such measurements do not enable to determine the committed effective dose, but are adequate to verify that a given threshold is not exceeded. For radioiodine, i.e. (123)I, (124)I, (125)I and (131)I, a calibrated surface contamination monitor is placed in front of the thyroid to detect whether the activity threshold has been exceeded. For radionuclides with very short physical half-lives (≤ 6 h), such as (99m)Tc and those used in positron emission tomography imaging, i.e. (11)C, (15)O, (18)F and (68)Ga, screening procedures consist in performing daily measurements of the ambient dose rate in front of the abdomen. Other gamma emitters used for imaging, i.e. (67)Ga, (111)In and (201)Tl, are measured with a scintillation detector located in front of the thorax. For pure beta emitters, i.e. (90)Y and (169)Er, as well as beta emitters with low-intensity gamma rays, i.e. (153)Sm, (177)Lu, (186)Re and (188)Re, the procedure consists in measuring hand contamination immediately after use. In Switzerland, screening procedures have been adopted by most nuclear medicine services since such measurements enable an acceptable monitoring while taking into account practical and economic considerations.

Predicting hematologic toxicity in patients undergoing radioimmunotherapy with 90Y-ibritumomab tiuxetan or 131I-tositumomab

This study aimed at identifying clinical factors for predicting hematologic toxicity after radioimmunotherapy with (90)Y-ibritumomab tiuxetan or (131)I-tositumomab in clinical practice.

METHODS

Hematologic data were available from 14 non-Hodgkin lymphoma patients treated with (90)Y-ibritumomab tiuxetan and 18 who received (131)I-tositumomab. The percentage baseline at nadir and 4 wk post nadir and the time to nadir were selected as the toxicity indicators for both platelets and neutrophils. Multiple linear regression analysis was performed to identify significant predictors (P < 0.05) of each indicator.

RESULTS

For both platelets and neutrophils, pooled and separate analyses of (90)Y-ibritumomab tiuxetan and (131)I-tositumomab data yielded the time elapsed since the last chemotherapy as the only significant predictor of the percentage baseline at nadir. The extent of bone marrow involvement was not a significant factor in this study, possibly because of the short time elapsed since the last chemotherapy of the 7 patients with bone marrow involvement. Because both treatments were designed to deliver a comparable bone marrow dose, this factor also was not significant. None of the 14 factors considered was predictive of the time to nadir. The R(2) value for the model predicting percentage baseline at nadir was 0.60 for platelets and 0.40 for neutrophils. This model predicted the platelet and neutrophil toxicity grade to within ±1 for 28 and 30 of the 32 patients, respectively. For the 7 patients predicted with grade I thrombocytopenia, 6 of whom had actual grade I-II, dosing might be increased to improve treatment efficacy.

CONCLUSION

The elapsed time since the last chemotherapy can be used to predict hematologic toxicity and customize the current dosing method in radioimmunotherapy.

Criteria for establishing shielding of multi-detector computed tomography (MDCT) rooms

The aim of this work is to compare two methods used for determining the proper shielding of computed tomography (CT) rooms while considering recent technological advances in CT scanners. The approaches of the German Institute for Standardisation and the US National Council on Radiation Protection and Measurements were compared and a series of radiation measurements were performed in several CT rooms at the Lausanne University Hospital. The following three-step procedure is proposed for assuring sufficient shielding of rooms hosting new CT units with spiral mode acquisition and various X-ray beam collimation widths: (1) calculate the ambient equivalent dose for a representative average weekly dose length product at the position where shielding is required; (2) from the maximum permissible weekly dose at the location of interest, calculate the transmission factor F that must be taken to ensure proper shielding and (3) convert the transmission factor into a thickness of lead shielding. A similar approach could be adopted to use when designing shielding for fluoroscopy rooms, where the basic quantity would be the dose area product instead of the load of current (milliampere-minute).

Using a NPWE model observer to assess suitable image quality for a digital mammography quality assurance programme

A method of objectively determining imaging performance for a mammography quality assurance programme for digital systems was developed. The method is based on the assessment of the visibility of a spherical microcalcification of 0.2 mm using a quasi-ideal observer model. It requires the assessment of the spatial resolution (modulation transfer function) and the noise power spectra of the systems. The contrast is measured using a 0.2-mm thick Al sheet and Polymethylmethacrylate (PMMA) blocks. The minimal image quality was defined as that giving a target contrast-to-noise ratio (CNR) of 5.4. Several evaluations of this objective method for evaluating image quality in mammography quality assurance programmes have been considered on computed radiography (CR) and digital radiography (DR) mammography systems. The measurement gives a threshold CNR necessary to reach the minimum standard image quality required with regards to the visibility of a 0.2-mm microcalcification. This method may replace the CDMAM image evaluation and simplify the threshold contrast visibility test used in mammography quality.

Monte Carlo simulation of a clearance box monitor used for nuclear power plant decommissioning

When decommissioning a nuclear facility it is important to be able to estimate activity levels of potentially radioactive samples and compare with clearance values defined by regulatory authorities. This paper presents a method of calibrating a clearance box monitor based on practical experimental measurements and Monte Carlo simulations. Adjusting the simulation for experimental data obtained using a simple point source permits the computation of absolute calibration factors for more complex geometries with an accuracy of a bit more than 20%. The uncertainty of the calibration factor can be improved to about 10% when the simulation is used relatively, in direct comparison with a measurement performed in the same geometry but with another nuclide. The simulation can also be used to validate the experimental calibration procedure when the sample is supposed to be homogeneous but the calibration factor is derived from a plate phantom. For more realistic geometries, like a small gravel dumpster, Monte Carlo simulation shows that the calibration factor obtained with a larger homogeneous phantom is correct within about 20%, if sample density is taken as the influencing parameter. Finally, simulation can be used to estimate the effect of a contamination hotspot. The research supporting this paper shows that activity could be largely underestimated in the event of a centrally-located hotspot and overestimated for a peripherally-located hotspot if the sample is assumed to be homogeneously contaminated. This demonstrates the usefulness of being able to complement experimental methods with Monte Carlo simulations in order to estimate calibration factors that cannot be directly measured because of a lack of available material or specific geometries.

Mass detection on mammograms: influence of signal shape uncertainty on human and model observers

We studied the influence of signal variability on human and model observers for detection tasks with realistic simulated masses superimposed on real patient mammographic backgrounds and synthesized mammographic backgrounds (clustered lumpy backgrounds, CLB). Results under the signal-known-exactly (SKE) paradigm were compared with signal-known-statistically (SKS) tasks for which the observers did not have prior knowledge of the shape or size of the signal. Human observers' performance did not vary significantly when benign masses were superimposed on real images or on CLB. Uncertainty and variability in signal shape did not degrade human performance significantly compared with the SKE task, while variability in signal size did. Implementation of appropriate internal noise components allowed the fit of model observers to human performance.

Mammographic texture synthesis: second-generation clustered lumpy backgrounds using a genetic algorithm

Synthetic yet realistic images are valuable for many applications in visual sciences and medical imaging. Typically, investigators develop algorithms and adjust their parameters to generate images that are visually similar to real images. In this study, we used a genetic algorithm and an objective, statistical similarity measure to optimize a particular texture generation algorithm, the clustered lumpy backgrounds (CLB) technique, and synthesize images mimicking real mammograms textures. We combined this approach with psychophysical experiments involving the judgment of radiologists, who were asked to qualify the visual realism of the images. Both objective and psychophysical approaches show that the optimized versions are significantly more realistic than the previous CLB model. Anatomical structures are well reproduced, and arbitrary large databases of mammographic texture with visual and statistical realism can be generated. Potential applications include detection experiments, where large amounts of statistically traceable yet realistic images are needed.

Human linear template with mammographic backgrounds estimated with a genetic algorithm

We estimated human observer linear templates underlying the detection of a realistic, spherical mass signal with mammographic backgrounds. Five trained naïve observers participated in two-alternative forced-choice (2-AFC) detection experiments with the signal superimposed on synthetic, clustered lumpy backgrounds (CLBs) in one condition and on nonstationary real mammographic backgrounds in another. Human observer linear templates were estimated using a genetic algorithm. A variety of common model observer templates were computed, and their shapes and associated performances were compared with those of the human observer. The estimated linear templates are not significantly different for stationary CLBs and real mammographic backgrounds. The estimated performance of the linear template compared with that of the human observers is within 5% in terms of percent correct (Pc) for the 2-AFC task. Channelized Hotelling models can fit human performance, but the templates differ considerably from the human linear template. Due to different local statistics, detection efficiency is significantly higher on nonstationary real backgrounds than on globally stationary synthetic CLBs. This finding emphasizes that nonstationary backgrounds need to be described by their local statistics.

Semiautomatic mammographic parenchymal patterns classification using multiple statistical features

RATIONALE AND OBJECTIVES

Our project was to investigate a complete methodology for the semiautomatic assessment of digital mammograms according to their density, an indicator known to be correlated to breast cancer risk. The BI-RADS four-grade density scale is usually employed by radiologists for reporting breast density, but it allows for a certain degree of subjective input, and an objective qualification of density has therefore often been reported hard to assess. The goal of this study was to design an objective technique for determining breast BI-RADS density.

MATERIALS AND METHODS

The proposed semiautomatic method makes use of complementary pattern recognition techniques to describe manually selected regions of interest (ROIs) in the breast with 36 statistical features. Three different classifiers based on a linear discriminant analysis or Bayesian theories were designed and tested on a database consisting of 1408 ROIs from 88 patients, using a leave-one-ROI-out technique. Classifications in optimal feature subspaces with lower dimensionality and reduction to a two-class problem were studied as well.

RESULTS

Comparison with a reference established by the classifications of three radiologists shows excellent performance of the classifiers, even though extremely dense breasts continue to remain more difficult to classify accurately. For the two best classifiers, the exact agreement percentages are 76% and above, and weighted kappa values are 0.78 and 0.83. Furthermore, classification in lower dimensional spaces and two-class problems give excellent results.

CONCLUSION

The proposed semiautomatic classifiers method provides an objective and reproducible method for characterizing breast density, especially for the two-class case. It represents a simple and valuable tool that could be used in screening programs, training, education, or for optimizing image processing in diagnostic tasks.

Shielding requirements in helical tomotherapy

Helical tomotherapy is a relatively new intensity-modulated radiation therapy (IMRT) treatment for which room shielding has to be reassessed for the following reasons. The beam-on-time needed to deliver a given target dose is increased and leads to a weekly workload of typically one order of magnitude higher than that for conventional radiation therapy. The special configuration of tomotherapy units does not allow the use of standard shielding calculation methods. A conventional linear accelerator must be shielded for primary, leakage and scatter photon radiations. For tomotherapy, primary radiation is no longer the main shielding issue since a beam stop is mounted on the gantry directly opposite the source. On the other hand, due to the longer irradiation time, the accelerator head leakage becomes a major concern. An analytical model based on geometric considerations has been developed to determine leakage radiation levels throughout the room for continuous gantry rotation. Compared to leakage radiation, scatter radiation is a minor contribution. Since tomotherapy units operate at a nominal energy of 6 MV, neutron production is negligible. This work proposes a synthetic and conservative model for calculating shielding requirements for the Hi-Art II TomoTherapy unit. Finally, the required concrete shielding thickness is given for different positions of interest.

Primary activity measurements with 4πγ NaI(Tl) counting and Monte Carlo calculated efficiencies

The radioactive concentrations of (18)F, (88)Y and (152)Eu solutions and the activity of (222)Rn gas ampoules are measured using a primary method involving 4pigamma NaI(Tl) integral counting with a well-type NaI(Tl) detector and efficiencies computed by Monte Carlo simulations. The simulations use the GEANT code coupled with a routine (sch2for), which generates randomly the decay paths and emissions depending on the decay scheme parameters. The resulting radioactive concentrations of (88)Y, (152)Eu and (18)F are found to agree with those measured with other primary measurement methods, such as 4pi (beta, e, X)-gamma coincidence counting or liquid scintillation counting. Results of the determination of the activity of (222)Rn gas ampoules by this method also match the results of an absolute standardisation technique in which radon is condensed onto a cold surface and its alpha-emissions are detected through an accurately specified solid angle.

Neutron measurements around storage casks containing spent fuel and vitrified high-level radioactive waste at ZWILAG

Spectrometric and dosimetric measurements were made around a cask containing spent fuel and a cask containing high-level radioactive waste at the Swiss intermediate waste and spent fuel storage facility. A Bonner sphere spectrometer, an LB 6411 neutron monitor and an Automess Szintomat 6134A were used to characterise the n-gamma fields at several locations around the two casks. The results of these measurements show that the neutron fluence spectra around the cask containing radioactive waste are harder and higher in intensity than those measured in the vicinity of the spent fuel cask. The ambient dose equivalents measured with the LB 6411 neutron monitor are in good agreement with those obtained using the Bonner spheres, except for locations with soft neutron spectra where the monitor overestimates the neutron ambient dose equivalent by almost 50%.

Number of X-ray examinations performed on paediatric and geriatric patients compared with adult patients

The age of the patient is of prime importance when assessing the radiological risk to patients due to medical X-ray exposures and the total detriment to the population due to radiodiagnostics. In order to take into account the age-specific radiosensitivity, three age groups are considered: children, adults and the elderly. In this work, the relative number of examinations carried out on paediatric and geriatric patients is established, compared with adult patients, for radiodiagnostics as a whole, for dental and medical radiology, for 8 radiological modalities as well as for 40 types of X-ray examinations. The relative numbers of X-ray examinations are determined based on the corresponding age distributions of patients and that of the general population. Two broad groups of X-ray examinations may be defined. Group A comprises conventional radiography, fluoroscopy and computed tomography; for this group a paediatric patient undergoes half the number of examinations as that of an adult, and a geriatric patient undergoes 2.5 times more. Group B comprises angiography and interventional procedures; for this group a paediatric patient undergoes a one-fourth of the number of examinations carried out on an adult, and a geriatric patient undergoes five times more.