Michael Ljungberg

A Monte Carlo program for the simulation of scintillation camera characteristics

There is a need for mathematical modelling for the evaluation of important parameters for photon imaging systems. A Monte Carlo program which simulates medical imaging nuclear detectors has been developed. Different materials can be chosen for the detector, a cover and a phantom. Cylindrical, spherical, rectangular and more complex phantom and source shapes can be simulated. Photoelectric, incoherent, coherent interactions and pair production are simulated. Different detector parameters, e.g. the energy pulse-height distribution and pulse pile-up due to finite decay time of the scintillation light emission, can be calculated. An energy resolution of the system is simulated by convolving the energy imparted with an energy-dependent Gaussian function. An image matrix of the centroid of the events in the detector can be simulated. Simulation of different collimators permits studies of spatial resolution and sensitivity. Comparisons of our results with experimental data and other published results have shown good agreement. The usefulness of the Monte Carlo code for the accurately simulation of important parameters in scintillation camera systems, stationary as well as SPECT (single-photon emission computed tomography) systems, has been demonstrated.

EANM/ESC procedural guidelines for myocardial perfusion imaging in nuclear cardiology

The European procedural guidelines for radionuclide imaging of myocardial perfusion and viability are presented in 13 sections covering patient information, radiopharmaceuticals, injected activities and dosimetry, stress tests, imaging protocols and acquisition, quality control and reconstruction methods, gated studies and attenuation-scatter compensation, data analysis, reports and image display, and positron emission tomography. If the specific recommendations given could not be based on evidence from original, scientific studies, we tried to express this state-of-art. The guidelines are designed to assist in the practice of performing, interpreting and reporting myocardial perfusion SPET. The guidelines do not discuss clinical indications, benefits or drawbacks of radionuclide myocardial imaging compared to non-nuclear techniques, nor do they cover cost benefit or cost effectiveness.

EANM/ESC guidelines for radionuclide imaging of cardiac function

Radionuclide imaging of cardiac function represents a number of well-validated techniques for accurate determination of right (RV) and left ventricular (LV) ejection fraction (EF) and LV volumes. These first European guidelines give recommendations for how and when to use first-pass and equilibrium radionuclide ventriculography, gated myocardial perfusion scintigraphy, gated PET, and studies with non-imaging devices for the evaluation of cardiac function. The items covered are presented in 11 sections: clinical indications, radiopharmaceuticals and dosimetry, study acquisition, RV EF, LV EF, LV volumes, LV regional function, LV diastolic function, reports and image display and reference values from the literature of RVEF, LVEF and LV volumes. If specific recommendations given cannot be based on evidence from original, scientific studies, referral is given to “prevailing or general consensus”. The guidelines are designed to assist in the practice of referral to, performance, interpretation and reporting of nuclear cardiology studies for the evaluation of cardiac performance.

MIRD Pamphlet No. 23: Quantitative SPECT for Patient-Specific 3-Dimensional Dosimetry in Internal Radionuclide Therapy

In internal radionuclide therapy, a growing interest in voxel-level estimates of tissue-absorbed dose has been driven by the desire to report radiobiologic quantities that account for the biologic consequences of both spatial and temporal nonuniformities in these dose estimates. This report presents an overview of 3-dimensional SPECT methods and requirements for internal dosimetry at both regional and voxel levels. Combined SPECT/CT image-based methods are emphasized, because the CT-derived anatomic information allows one to address multiple technical factors that affect SPECT quantification while facilitating the patient-specific voxel-level dosimetry calculation itself. SPECT imaging and reconstruction techniques for quantification in radionuclide therapy are not necessarily the same as those designed to optimize diagnostic imaging quality. The current overview is intended as an introduction to an upcoming series of MIRD pamphlets with detailed radionuclide-specific recommendations intended to provide best-practice SPECT quantification–based guidance for radionuclide dosimetry.

Lu-177-[DOTA0,Tyr3] Octreotate Therapy in Patients With Disseminated Neuroendocrine Tumors: Analysis of Dosimetry With Impact on Future Therapeutic Strategy

177Lu‐(DOTA0,Tyr3) octreotate is a new treatment modality for disseminated neuroendocrine tumors. According to a consensus protocol, the calculated maximally tolerated absorbed dose to the kidney should not exceed 27 Gy. In commonly used dosimetry methods, planar imaging is used for determination of the residence time, whereas the kidney mass is determined from a computed tomography (CT) scan.

Evaluation of quantitative (90)Y SPECT based on experimental phantom studies.

In SPECT imaging of pure beta emitters, such as (90)Y, the acquired spectrum is very complex, which increases the demands on the imaging protocol and the reconstruction. In this work, we have evaluated the quantitative accuracy of bremsstrahlung SPECT with focus on the reconstruction algorithm including model-based attenuation, scatter and collimator-detector response (CDR) compensations. The scatter and CDR compensation methods require pre-calculated point-spread functions, which were generated with the SIMIND MC program. The SIMIND program is dedicated for simulation of scintillation camera imaging and only handles photons. The aim of this work was therefore twofold. The first aim was to implement simulation of bremsstrahlung imaging into the SIMIND code and to validate simulations against experimental measurements. The second was to investigate the quality of bremsstrahlung SPECT imaging and to evaluate the possibility of quantifying the activity in differently shaped sources. In addition, a feasibility test was performed on a patient that underwent treatment with (90)Y-Ibritumomab tiuxetan (Zevalin). The MCNPX MC program was used to generate bremsstrahlung photon spectra which were used as source input in the SIMIND program. The obtained bremsstrahlung spectra were separately validated by experimental measurement using a HPGe detector. Validation of the SIMIND generated images was done by a comparison to gamma camera measurements of a syringe containing (90)Y. Results showed a slight deviation between simulations and measurements in image regions outside the source, but the agreement was sufficient for the purpose of generating scatter and CDR kernels. For the bremsstrahlung SPECT experiment, the RSD torso phantom with (90)Y in the liver insert was measured with and without background activities. Projection data were obtained using a GE VH/Hawkeye system. Image reconstruction was performed by using the OSEM algorithm with and without different combinations of model-based attenuation, scatter and CDR compensations. The reconstructed images were then evaluated in terms of the accuracy of the total activity estimate in the liver insert. It was found that the activity in a large source such as the liver was estimated with a bias of around -70%, when no compensations were included in the reconstruction, whereas when compensations were included the bias obtained was between -10 and 16%. It is concluded that although the (90)Y bremsstrahlung spectrum is continuous with no pronounced peak and the count rate is low, it is possible to achieve reasonably accurate activity estimates from bremsstrahlung SPECT images if proper compensations are applied in the reconstruction. This conclusion was also confirmed by the patient study.

A Monte Carlo investigation of artifacts caused by liver uptake in single-photon emission computed tomography perfusion imaging with technetium 99m-labeled agents

BackgroundSignificant hepatobiliary accumulation of technetium 99m-labeled cardiac perfusion agents has been shown to cause alterations in the apparent localization of the agents in the cardiac walls. A Monte Carlo study was conducted to investigate the hypothesis that the cardiac count changes are due to the inconsistencies in the projection data input to reconstruction, and that correction of the causes of these inconsistencies before reconstruction, or including knowledge of the physics underlying them in the reconstruction algorithm, would virtually eliminate these artifacts.Methods and ResultsThe SIMIND Monte Carlo package was used to simulate 64×64 pixel projection images at 128 angles of the three-dimensional mathematical cardiac-torso (MCAT) phantom. Simulations were made of (1) a point source in the liver, (2) cardiac activity only, and (3) hepatic activity only. The planar projections and reconstructed point spread functions (PSFs) of the point source in the liver were investigated to study the nature of the inconsistencies introduced into the projections by imaging, and how these affect the distribution of counts in the reconstructed slices. Bull’s eye polar maps of the counts at the center of the left ventricular wall of filtered back-projection (FBP) and maximum-likelihood expectation-maximization (MLEM) reconstructions of projections with solely cardiac activity, and with cardiac activity plus hepatic activity scaled to have twice the cardiac concentration, were compared to determine the magnitude and location of apparent changes in cardiac activity when hepatic activity is present. Separate simulations were made to allow the investigation of stationary spatial resolution, distance-dependent spatial resolution, attenuation, and scatter. The point source projections showed significant inconsistencies as a function of projection angle with the largest effect being caused by attenuation. When consistent projections were simulated, no significant impact of hepatic activity on cardiac counts was noted with FBP, or 100 iterations of MLEM. With inconsistent projections, reconstruction of 180 degrees resulted in greater apparent cardiac count losses than did 360 degrees reconstruction for both FBP and MLEM. The incorporation of attenuation correction in MLEM reconstruction reduced the changes in cardiac counts to that seen in simulations in which attenuation was not included, but resulted in increased apparent localization of activity in the posterior wall of the left ventricle when scatter was present in the simulated images.ConclusionsThe apparent alterations in cardiac counts when significant hepatic localization is present is due to the inconsistency of the projections inherent in imaging. Prior correction of these, or accounting for them in the reconstruction algorithm, will virtually eliminate them as causes of artifactual changes in localization. Attenuation correction and scatter correction are both required to overcome the major sources of apparent count changes in the heart associated with hepatic uptake.

EANM procedural guidelines for radionuclide myocardial perfusion imaging with SPECT and SPECT/CT : 2015 revision

Since the publication of the European Association of Nuclear Medicine (EANM) procedural guidelines for radionuclide myocardial perfusion imaging (MPI) in 2005, many small and some larger steps of progress have been made, improving MPI procedures. In this paper, the major changes from the updated 2015 procedural guidelines are highlighted, focusing on the important changes related to new instrumentation with improved image information and the possibility to reduce radiation exposure, which is further discussed in relation to the recent developments of new International Commission on Radiological Protection (ICRP) models. Introduction of the selective coronary vasodilator regadenoson and the use of coronary CT-contrast agents for hybrid imaging with SPECT/CT angiography are other important areas for nuclear cardiology that were not included in the previous guidelines. A large number of minor changes have been described in more detail in the fully revised version available at the EANM home page: http://eanm.org/publications/guidelines/2015_07_EANM_FINAL_myocardial_perfusion_guideline.pdf.

RADIATION DOSIMETRY IN NUCLEAR MEDICINE

Radionuclides are used in nuclear medicine in a variety of diagnostic and therapeutic procedures. A knowledge of the radiation dose received by different organs in the body is essential to an evaluation of the risks and benefits of any procedure. In this paper, current methods for internal dosimetry are reviewed, as they are applied in nuclear medicine. Particularly, the Medical Internal Radiation Dose (MIRD) system for dosimetry is explained, and many of its published resources discussed. Available models representing individuals of different age and gender, including those representing the pregnant woman are described; current trends in establishing models for individual patients are also evaluated. The proper design of kinetic studies for establishing radiation doses for radiopharmaceuticals is discussed. An overview of how to use information obtained in a dosimetry study, including that of the effective dose equivalent (ICRP 30) and effective dose (ICRP 60), is given. Current trends and issues in internal dosimetry, including the calculation of patient-specific doses and in the use of small scale and microdosimetry techniques, are also reviewed.

Development and evaluation of an improved quantitative 90Y bremsstrahlung SPECT method

PURPOSE Yttrium-90 ((90)Y) is one of the most commonly used radionuclides in targeted radionuclide therapy (TRT). Since it decays with essentially no gamma photon emissions, surrogate radionuclides (e.g., (111)In) or imaging agents (e.g., (99m)Tc MAA) are typically used for treatment planning. It would, however, be useful to image (90)Y directly in order to confirm that the distributions measured with these other radionuclides or agents are the same as for the (90)Y labeled agents. As a result, there has been a great deal of interest in quantitative imaging of (90)Y bremsstrahlung photons using single photon emission computed tomography (SPECT) imaging. The continuous and broad energy distribution of bremsstrahlung photons, however, imposes substantial challenges on accurate quantification of the activity distribution. The aim of this work was to develop and evaluate an improved quantitative (90)Y bremsstrahlung SPECT reconstruction method appropriate for these imaging applications. METHODS Accurate modeling of image degrading factors such as object attenuation and scatter and the collimator-detector response is essential to obtain quantitatively accurate images. All of the image degrading factors are energy dependent. Thus, the authors separated the modeling of the bremsstrahlung photons into multiple categories and energy ranges. To improve the accuracy, the authors used a bremsstrahlung energy spectrum previously estimated from experimental measurements and incorporated a model of the distance between (90)Y decay location and bremsstrahlung emission location into the SIMIND code used to generate the response functions and kernels used in the model. This improved Monte Carlo bremsstrahlung simulation was validated by comparison to experimentally measured projection data of a (90)Y line source. The authors validated the accuracy of the forward projection model for photons in the various categories and energy ranges using the validated Monte Carlo (MC) simulation method. The forward projection model was incorporated into an iterative ordered subsets-expectation maximization (OS-EM) reconstruction code to allow for quantitative SPECT reconstruction. The resulting code was validated using both a physical phantom experiment with spherical objects in a warm background and a realistic anatomical phantom simulation. In the physical phantom study, the authors evaluated the method in terms of quantitative accuracy of activity estimates in the spheres; in the simulation study, the authors evaluated the accuracy and precision of activity estimates from various organs and compared them to results from a previously proposed method. RESULTS The authors demonstrated excellent agreement between the experimental measurement and Monte Carlo simulation. In the XCAT phantom simulation, the proposed method achieved much better accuracy in the modeling (error in photon counts was -1.1 %) compared to a previously proposed method (errors were more than 20  %); the quantitative accuracy of activity estimates was excellent for all organs (errors were from -1.6 % to 11.9 %) and comparable to previously published results for (131)I using the same collimator. CONCLUSIONS The proposed (90)Y bremsstrahlung SPECT reconstruction method provided very accurate estimates of organ activities, with accuracies approaching those previously observed for (131)I. The method may be useful in verifying organ doses for targeted radionuclide therapy using (90)Y.