Dose reduction and image enhancement in micro-CT using deep learning.

BACKGROUND: In preclinical settings, micro-computed tomography (CT) provides a powerful tool to acquire high resolution anatomical images of rodents and offers the advantage to in vivo non-invasively assess disease progression and therapy efficacy. Much higher resolutions are needed to achieve scale-equivalent discriminatory capabilities in rodents as those in humans. High resolution imaging however comes at the expense of increased scan times and higher doses. Specifically, with preclinical longitudinal imaging, there are concerns that dose accumulation may affect experimental outcomes of animal models. PURPOSE: Dose reduction efforts under the ALARA (as low as reasonably achievable) principles are thus a key point of attention. However, low dose CT acquisitions inherently induce higher noise levels which deteriorate image quality and negatively impact diagnostic performance. Many denoising techniques already exist, and deep learning (DL) has become increasingly popular for image denoising, but research has mostly focused on clinical CT with limited studies conducted on preclinical CT imaging. We investigate the potential of convolutional neural networks (CNN) for restoring high quality micro-CT images from low dose (noisy) images. The novelty of the CNN denoising frameworks presented in this work consists of utilizing image pairs with realistic CT noise present in the input as well as the target image used for the model training; a noisier image acquired with a low dose protocol is matched to a less noisy image acquired with a higher dose scan of the same mouse. METHODS: Low and high dose ex vivo micro-CT scans of 38 mice were acquired. Two CNN models, based on a 2D and 3D four-layer U-Net, were trained with mean absolute error (30 training, 4 validation and 4 test sets). To assess denoising performance, ex vivo mice and phantom data were used. Both CNN approaches were compared to existing methods, like spatial filtering (Gaussian, Median, Wiener) and iterative total variation image reconstruction algorithm. Image quality metrics were derived from the phantom images. A first observer study (n = 23) was set-up to rank overall quality of differently denoised images. A second observer study (n = 18) estimated the dose reduction factor of the investigated 2D CNN method. RESULTS: Visual and quantitative results show that both CNN algorithms exhibit superior performance in terms of noise suppression, structural preservation and contrast enhancement over comparator methods. The quality scoring by 23 medical imaging experts also indicates that the investigated 2D CNN approach is consistently evaluated as the best performing denoising method. Results from the second observer study and quantitative measurements suggest that CNN-based denoising could offer a 2-4× dose reduction, with an estimated dose reduction factor of about 3.2 for the considered 2D network. CONCLUSIONS: Our results demonstrate the potential of DL in micro-CT for higher quality imaging at low dose acquisition settings. In the context of preclinical research, this offers promising future prospects for managing the cumulative severity effects of radiation in longitudinal studies.

Performance evaluation of a micro-CT system for laboratory animal imaging with iterative reconstruction capabilities.

BACKGROUND: In recent years, there has been a rapid proliferation in micro-computed tomography (micro-CT) systems becoming more available for routine preclinical research, with applications in many areas, including bone, lung, cancer, and cardiac imaging. Micro-CT provides the means to non-invasively acquire detailed anatomical information, but high-resolution imaging comes at the cost of longer scan times and higher doses, which is not desirable given the potential risks related to x-ray radiation. To achieve dose reduction and higher throughputs without compromising image quality, fewer projections can be acquired. This is where iterative reconstruction methods can have the potential to reduce noise since these algorithms can better handle sparse projection data, compared to filtered backprojection PURPOSE: We evaluate the performance characteristics of a compact benchtop micro-CT scanner that provides iterative reconstruction capabilities with GPU-based acceleration. We thereby investigate the potential benefit of iterative reconstruction for dose reduction. METHODS: Based on a series of phantom experiments, the benchtop micro-CT system was characterized in terms of image uniformity, noise, low contrast detectability, linearity, and spatial resolution. Whole-body images of a plasticized ex vivo mouse phantom were also acquired. Different acquisition protocols (general-purpose versus high-resolution, including low dose scans) and different reconstruction strategies (analytic versus iterative algorithms: FDK, ISRA, ISRA-TV) were compared. RESULTS: Signal uniformity was maintained across the radial and axial field-of-view (no cupping effect) with an average difference in Hounsfield units (HU) between peripheral and central regions below 50. For low contrast detectability, regions with at least ∆HU of 40 to surrounding material could be discriminated (for rods of 2.5 mm diameter). A high linear correlation (R2  = 0.997) was found between measured CT values and iodine concentrations (0-40 mg/ml). Modulation transfer function (MTF) calculations on a wire phantom evaluated a resolution of 10.2 lp/mm at 10% MTF that was consistent with the 8.3% MTF measured on the 50 µm bars (10 lp/mm) of a bar-pattern phantom. Noteworthy changes in signal-to-noise and contrast-to-noise values were found for different acquisition and reconstruction protocols. Our results further showed the potential of iterative reconstruction to deliver images with less noise and artefacts. CONCLUSIONS: In summary, the micro-CT system that was evaluated in the present work was shown to provide a good combination of performance characteristics between image uniformity, low contrast detectability, and resolution in short scan times. With the iterative reconstruction capabilities of this micro-CT system in mind (ISRA and ISRA-TV), the adoption of such algorithms by GPU-based acceleration enables the integration of noise reduction methods which here demonstrated potential for high-quality imaging at reduced doses.

Guide to Plant-PET Imaging Using 11CO2.

Due to its high sensitivity and specificity for tumor detection, positron emission tomography (PET) has become a standard and widely used molecular imaging technique. Given the popularity of PET, both clinically and preclinically, its use has been extended to study plants. However, only a limited number of research groups worldwide report PET-based studies, while we believe that this technique has much more potential and could contribute extensively to plant science. The limited application of PET may be related to the complexity of putting together methodological developments from multiple disciplines, such as radio-pharmacology, physics, mathematics and engineering, which may form an obstacle for some research groups. By means of this manuscript, we want to encourage researchers to study plants using PET. The main goal is to provide a clear description on how to design and execute PET scans, process the resulting data and fully explore its potential by quantification via compartmental modeling. The different steps that need to be taken will be discussed as well as the related challenges. Hereby, the main focus will be on, although not limited to, tracing 11CO2 to study plant carbon dynamics.

Roadmap toward the 10 ps time-of-flight PET challenge.

Since the seventies, positron emission tomography (PET) has become an invaluable medical molecular imaging modality with an unprecedented sensitivity at the picomolar level, especially for cancer diagnosis and the monitoring of its response to therapy. More recently, its combination with x-ray computed tomography (CT) or magnetic resonance (MR) has added high precision anatomic information in fused PET/CT and PET/MR images, thus compensating for the modest intrinsic spatial resolution of PET. Nevertheless, a number of medical challenges call for further improvements in PET sensitivity. These concern in particular new treatment opportunities in the context personalized (also called precision) medicine, such as the need to dynamically track a small number of cells in cancer immunotherapy or stem cells for tissue repair procedures. A better signal-to-noise ratio (SNR) in the image would allow detecting smaller size tumours together with a better staging of the patients, thus increasing the chances of putting cancer in complete remission. Moreover, there is an increasing demand for reducing the radioactive doses injected to the patients without impairing image quality. There are three ways to improve PET scanner sensitivity: improving detector efficiency, increasing geometrical acceptance of the imaging device and pushing the timing performance of the detectors. Currently, some pre-localization of the electron-positron annihilation along a line-of-response (LOR) given by the detection of a pair of annihilation photons is provided by the detection of the time difference between the two photons, also known as the time-of-flight (TOF) difference of the photons, whose accuracy is given by the coincidence time resolution (CTR). A CTR of about 10 picoseconds FWHM will ultimately allow to obtain a direct 3D volume representation of the activity distribution of a positron emitting radiopharmaceutical, at the millimetre level, thus introducing a quantum leap in PET imaging and quantification and fostering more frequent use of 11C radiopharmaceuticals. The present roadmap article toward the advent of 10 ps TOF-PET addresses the status and current/future challenges along the development of TOF-PET with the objective to reach this mythic 10 ps frontier that will open the door to real-time volume imaging virtually without tomographic inversion. The medical impact and prospects to achieve this technological revolution from the detection and image reconstruction point-of-views, together with a few perspectives beyond the TOF-PET application are discussed.

State of the art in total body PET.

The idea of a very sensitive positron emission tomography (PET) system covering a large portion of the body of a patient already dates back to the early 1990s. In the period 2000-2010, only some prototypes with long axial field of view (FOV) have been built, which never resulted in systems used for clinical research. One of the reasons was the limitations in the available detector technology, which did not yet have sufficient energy resolution, timing resolution or countrate capabilities for fully exploiting the benefits of a long axial FOV design. PET was also not yet as widespread as it is today: the growth in oncology, which has become the major application of PET, appeared only after the introduction of PET-CT (early 2000).The detector technology used in most clinical PET systems today has a combination of good energy and timing resolution with higher countrate capabilities and has now been used since more than a decade to build time-of-flight (TOF) PET systems with fully 3D acquisitions. Based on this technology, one can construct total body PET systems and the remaining challenges (data handling, fast image reconstruction, detector cooling) are mostly related to engineering. The direct benefits of long axial FOV systems are mostly related to the higher sensitivity. For single organ imaging, the gain is close to the point source sensitivity which increases linearly with the axial length until it is limited by solid angle and attenuation of the body. The gains for single organ (compared to a fully 3D PET 20-cm axial FOV) are limited to a factor 3-4. But for long objects (like body scans), it increases quadratically with scanner length and factors of 10-40 × higher sensitivity are predicted for the long axial FOV scanner. This application of PET has seen a major increase (mostly in oncology) during the last 2 decades and is now the main type of study in a PET centre. As the technology is available and the full body concept also seems to match with existing applications, the old concept of a total body PET scanner is seeing a clear revival. Several research groups are working on this concept and after showing the potential via extensive simulations; construction of these systems has started about 2 years ago. In the first phase, two PET systems with long axial FOV suitable for large animal imaging were constructed to explore the potential in more experimental settings. Recently, the first completed total body PET systems for human use, a 70-cm-long system, called PennPET Explorer, and a 2-m-long system, called uExplorer, have become reality and first clinical studies have been shown. These results illustrate the large potential of this concept with regard to low-dose imaging, faster scanning, whole-body dynamic imaging and follow-up of tracers over longer periods. This large range of possible technical improvements seems to have the potential to change the current clinical routine and to expand the number of clinical applications of molecular imaging. The J-PET prototype is a prototype system with a long axial FOV built from axially arranged plastic scintillator strips.This paper gives an overview of the recent technical developments with regard to PET scanners with a long axial FOV covering at least the majority of the body (so called total body PET systems). After explaining the benefits and challenges of total body PET systems, the different total body PET system designs proposed for large animal and clinical imaging are described in detail. The axial length is one of the major factors determining the total cost of the system, but there are also other options in detector technology, design and processing for reducing the cost these systems. The limitations and advantages of different designs for research and clinical use are discussed taking into account potential applications and the increased cost of these systems.

18F-FDG micro-PET/CT for intra-operative margin assessment during breast-conserving surgery.

Rationale: Positive surgical margins for invasive breast cancer (BC) treated with breast-conserving surgery (BCS) are defined as ink on tumor. The rate of positive margins is approximately 20%, since a time- and cost-effective method for margin assessment is lacking. In this study, we investigated margin status by intra-operative imaging using high-resolution 18 F-fluoro-deoxyglucose (FDG) positron emission tomography (PET) and X-ray computed tomography (CT).Methods: Twenty patients were enrolled and received 4 MBq/kg of FDG prior to surgery. Intra-operative imaging of the specimens was performed by the MOLECUBES ß-CUBE (PET) and X-CUBE (CT). Margin status was assessed by three surgeons and compared with an algorithm. The sensitivity and specificity were calculated by using histopathological assessment as a gold standard.Results: A region with high FDG uptake was visualized in all specimens. Automated analysis showed a sensitivity of 90%, a specificity of 60%, and an area under the curve (AUC) of 0.86 after ROC analysis. Margin assessment by the surgeons resulted in a mean sensitivity and specificity of 79% and 72%, respectively.Conclusions: This proof-of-concept study demonstrates that high-resolution FDG-PET/CT can facilitate intra-operative margin assessment during BCS. This technique achieves good sensitivity and specificity and may therefore reduce re-operation rates in the future.

A dosimetry procedure for organs-at-risk in 177Lu peptide receptor radionuclide therapy of patients with neuroendocrine tumours.

PURPOSE: Peptide receptor radionuclide therapy with 177Lu-DOTATATE has become a standard treatment modality in neuroendocrine tumours (NETs). No consensus has yet been reached however regarding the absorbed dose threshold for lesion response, the absorbed dose limit to organs-at-risk, and the optimal fractionation and activity to be administered. This is partly due to a lack of uniform and comparable dosimetry protocols. The present article details the development of an organ-at-risk dosimetry procedure, which could be implemented and used routinely in a clinical context. METHODS: Forty-seven patients with NETs underwent 177Lu-DOTATATE therapy. Three SPECT/CT images were acquired at 4, 24 and 144-192 h post-injection. Three blood samples were obtained together with the SPECT/CT acquisitions and 2 additional samples were obtained around 30 min and 1 h post-injection. A bi-exponential fit was used to compute the source organ time-integrated activity coefficients. Coefficients were introduced into OLINDA/EXM software to compute organ-at-risk absorbed doses. Median values for all patients were computed for absorbed dose coefficient D/A0 and for late effective half-life T1/2eff for kidneys, spleen and red marrow. RESULTS: Dosimetry resulted in a median[interquartile range] of 0.78[0.35], 1.07[0.58] and 0.028[0.010] Gy/GBq for D/A0 and of 55[9], 71[9] and 52[18] h for T1/2eff for kidneys, spleen and red marrow respectively. CONCLUSIONS: A dosimetry procedure for organs-at-risk in 177Lu-DOTATATE therapy based on serial SPECT/CT images and blood samples can be implemented routinely in a clinical context with limited patient burden. The results obtained were in accordance with those of other centres.

Performance evaluation of the LightPath imaging system for intra-operative Cerenkov luminescence imaging.

The performances of an intra-operative optical imaging system for Cerenkov luminescence imaging of resected tumor specimens were evaluated with phantom studies. The spatial resolution, the linearity of the measured signal with the activity concentration and the minimum detectable activity concentration were considered. A high linearity was observed over a broad range of activity concentration (R2⩾0.99 down to ∼40 kBq/ml of 18F-FDG). For 18F-FDG activity distributions 2 mm deep in biological tissue, the measured detection limit was 8 kBq/ml and a spatial resolution of 2.5 mm was obtained. The detection limit of the imaging system is comparable with clinical activity concentrations in tumor specimens, and the spatial resolution is compatible with clinical requirements.

Quantitative Imaging for Targeted Radionuclide Therapy Dosimetry - Technical Review.

Targeted radionuclide therapy (TRT) is a promising technique for cancer therapy. However, in order to deliver the required dose to the tumor, minimize potential toxicity in normal organs, as well as monitor therapeutic effects, it is important to assess the individualized internal dosimetry based on patient-specific data. Advanced imaging techniques, especially radionuclide imaging, can be used to determine the spatial distribution of administered tracers for calculating the organ-absorbed dose. While planar scintigraphy is still the mainstream imaging method, SPECT, PET and bremsstrahlung imaging have promising properties to improve accuracy in quantification. This article reviews the basic principles of TRT and discusses the latest development in radionuclide imaging techniques for different theranostic agents, with emphasis on their potential to improve personalized TRT dosimetry.

EEG source connectivity to localize the seizure onset zone in patients with drug resistant epilepsy.

Electrical source imaging (ESI) from interictal scalp EEG is increasingly validated and used as a valuable tool in the presurgical evaluation of epilepsy as a reflection of the irritative zone. ESI of ictal scalp EEG to localize the seizure onset zone (SOZ) remains challenging. We investigated the value of an approach for ictal imaging using ESI and functional connectivity analysis (FC). Ictal scalp EEG from 111 seizures in 27 patients who had Engel class I outcome at least 1 year following resective surgery was analyzed. For every seizure, an artifact-free epoch close to the seizure onset was selected and ESI using LORETA was applied. In addition, the reconstructed sources underwent FC using the spectrum-weighted Adaptive Directed Transfer Function. This resulted in the estimation of the SOZ in two ways: (i) the source with maximal power after ESI, (ii) the source with the strongest outgoing connections after combined ESI and FC. Next, we calculated the distance between the estimated SOZ and the border of the resected zone (RZ) for both approaches and called this the localization error ((i) LEpow and (ii) LEconn respectively). By comparing LEpow and LEconn, we assessed the added value of FC. The source with maximal power after ESI was inside the RZ (LEpow = 0 mm) in 31% of the seizures and estimated within 10 mm from the border of the RZ (LEpow ≤ 10 mm) in 42%. Using ESI and FC, these numbers increased to 72% for LEconn = 0 mm and 94% for LEconn ≤ 10 mm. FC provided a significant added value to ESI alone (p < 0.001). ESI combined with subsequent FC is able to localize the SOZ in a non-invasive way with high accuracy. Therefore it could be a valuable tool in the presurgical evaluation of epilepsy.

Seizure Onset Zone Localization from Ictal High-Density EEG in Refractory Focal Epilepsy.

Epilepsy surgery is the most efficient treatment option for patients with refractory epilepsy. Before surgery, it is of utmost importance to accurately delineate the seizure onset zone (SOZ). Non-invasive EEG is the most used neuroimaging technique to diagnose epilepsy, but it is hard to localize the SOZ from EEG due to its low spatial resolution and because epilepsy is a network disease, with several brain regions becoming active during a seizure. In this work, we propose and validate an approach based on EEG source imaging (ESI) combined with functional connectivity analysis to overcome these problems. We considered both simulations and real data of patients. Ictal epochs of 204-channel EEG and subsets down to 32 channels were analyzed. ESI was done using realistic head models and LORETA was used as inverse technique. The connectivity pattern between the reconstructed sources was calculated, and the source with the highest number of outgoing connections was selected as SOZ. We compared this algorithm with a more straightforward approach, i.e. selecting the source with the highest power after ESI as the SOZ. We found that functional connectivity analysis estimated the SOZ consistently closer to the simulated EZ/RZ than localization based on maximal power. Performance, however, decreased when 128 electrodes or less were used, especially in the realistic data. The results show the added value of functional connectivity analysis for SOZ localization, when the EEG is obtained with a high-density setup. Next to this, the method can potentially be used as objective tool in clinical settings.

Electrical source imaging of interictal spikes using multiple sparse volumetric priors for presurgical epileptogenic focus localization.

Electrical source imaging of interictal spikes observed in EEG recordings of patients with refractory epilepsy provides useful information to localize the epileptogenic focus during the presurgical evaluation. However, the selection of the time points or time epochs of the spikes in order to estimate the origin of the activity remains a challenge. In this study, we consider a Bayesian EEG source imaging technique for distributed sources, i.e. the multiple volumetric sparse priors (MSVP) approach. The approach allows to estimate the time courses of the intensity of the sources corresponding with a specific time epoch of the spike. Based on presurgical averaged interictal spikes in six patients who were successfully treated with surgery, we estimated the time courses of the source intensities for three different time epochs: (i) an epoch starting 50 ms before the spike peak and ending at 50% of the spike peak during the rising phase of the spike, (ii) an epoch starting 50 ms before the spike peak and ending at the spike peak and (iii) an epoch containing the full spike time period starting 50 ms before the spike peak and ending 230 ms after the spike peak. To identify the primary source of the spike activity, the source with the maximum energy from 50 ms before the spike peak till 50% of the spike peak was subsequently selected for each of the time windows. For comparison, the activity at the spike peaks and at 50% of the peaks was localized using the LORETA inversion technique and an ECD approach. Both patient-specific spherical forward models and patient-specific 5-layered finite difference models were considered to evaluate the influence of the forward model. Based on the resected zones in each of the patients, extracted from post-operative MR images, we compared the distances to the resection border of the estimated activity. Using the spherical models, the distances to the resection border for the MSVP approach and each of the different time epochs were in the same range as the LORETA and ECD techniques. We found distances smaller than 23 mm, with robust results for all the patients. For the finite difference models, we found that the distances to the resection border for the MSVP inversions of the full spike time epochs were generally smaller compared to the MSVP inversions of the time epochs before the spike peak. The results also suggest that the inversions using the finite difference models resulted in slightly smaller distances to the resection border compared to the spherical models. The results we obtained are promising because the MSVP approach allows to study the network of the estimated source-intensities and allows to characterize the spatial extent of the underlying sources.

The Role of Skull Modeling in EEG Source Imaging for Patients with Refractory Temporal Lobe Epilepsy.

We investigated the influence of different skull modeling approaches on EEG source imaging (ESI), using data of six patients with refractory temporal lobe epilepsy who later underwent successful epilepsy surgery. Four realistic head models with different skull compartments, based on finite difference methods, were constructed for each patient: (i) Three models had skulls with compact and spongy bone compartments as well as air-filled cavities, segmented from either computed tomography (CT), magnetic resonance imaging (MRI) or a CT-template and (ii) one model included a MRI-based skull with a single compact bone compartment. In all patients we performed ESI of single and averaged spikes marked in the clinical 27-channel EEG by the epileptologist. To analyze at which time point the dipole estimations were closer to the resected zone, ESI was performed at two time instants: the half-rising phase and peak of the spike. The estimated sources for each model were validated against the resected area, as indicated by the postoperative MRI. Our results showed that single spike analysis was highly influenced by the signal-to-noise ratio (SNR), yielding estimations with smaller distances to the resected volume at the peak of the spike. Although averaging reduced the SNR effects, it did not always result in dipole estimations lying closer to the resection. The proposed skull modeling approaches did not lead to significant differences in the localization of the irritative zone from clinical EEG data with low spatial sampling density. Furthermore, we showed that a simple skull model (MRI-based) resulted in similar accuracy in dipole estimation compared to more complex head models (based on CT- or CT-template). Therefore, all the considered head models can be used in the presurgical evaluation of patients with temporal lobe epilepsy to localize the irritative zone from low-density clinical EEG recordings.

Use of a ray-based reconstruction algorithm to accurately quantify preclinical microSPECT images.

This work aimed to measure the in vivo quantification errors obtained when ray-based iterative reconstruction is used in micro-single-photon emission computed tomography (SPECT). This was investigated with an extensive phantom-based evaluation and two typical in vivo studies using 99mTc and 111In, measured on a commercially available cadmium zinc telluride (CZT)-based small-animal scanner. Iterative reconstruction was implemented on the GPU using ray tracing, including (1) scatter correction, (2) computed tomography-based attenuation correction, (3) resolution recovery, and (4) edge-preserving smoothing. It was validated using a National Electrical Manufacturers Association (NEMA) phantom. The in vivo quantification error was determined for two radiotracers: [99mTc]DMSA in naive mice (n  =  10 kidneys) and [111In]octreotide in mice (n  =  6) inoculated with a xenograft neuroendocrine tumor (NCI-H727). The measured energy resolution is 5.3% for 140.51 keV (99mTc), 4.8% for 171.30 keV, and 3.3% for 245.39 keV (111In). For 99mTc, an uncorrected quantification error of 28 ± 3% is reduced to 8 ± 3%. For 111In, the error reduces from 26 ± 14% to 6 ± 22%. The in vivo error obtained with 99mTc-dimercaptosuccinic acid ([99mTc]DMSA) is reduced from 16.2 ± 2.8% to -0.3 ± 2.1% and from 16.7 ± 10.1% to 2.2 ± 10.6% with [111In]octreotide. Absolute quantitative in vivo SPECT is possible without explicit system matrix measurements. An absolute in vivo quantification error smaller than 5% was achieved and exemplified for both [99mTc]DMSA and [111In]octreotide.

Ictal-onset localization through connectivity analysis of intracranial EEG signals in patients with refractory epilepsy.

PURPOSE: Fifteen percent to 25% of patients with refractory epilepsy require invasive video-electroencephalography (EEG) monitoring (IVEM) to precisely delineate the ictal-onset zone. This delineation based on the recorded intracranial EEG (iEEG) signals occurs visually by the epileptologist and is therefore prone to human mistakes. The purpose of this study is to investigate whether effective connectivity analysis of intracranially recorded EEG during seizures provides an objective method to localize the ictal-onset zone. METHODS: In this study data were analyzed from eight patients who underwent IVEM at Ghent University Hospital in Belgium. All patients had a focal ictal onset and were seizure-free following resective surgery. The effective connectivity pattern was calculated during the first 20 s of ictal rhythmic iEEG activity. The out-degree, which is reflective of the number of outgoing connections, was calculated for each electrode contact for every single seizure during these 20 s. The seizure specific out-degrees were summed per patient to obtain the total out-degree. The electrode contact with the highest total out-degree was considered indicative of localization of the ictal-onset zone. This result was compared to the conclusion of the visual analysis of the epileptologist and the resected brain region segmented from postoperative magnetic resonance imaging (MRI). KEY FINDINGS: In all eight patients the electrode contact with the highest total out-degree was among the contacts identified by the epileptologist as the ictal onset. This contact, that we named "the driver," always laid within the resected brain region. Furthermore, the patient-specific connectivity patterns were consistent over the majority of seizures. SIGNIFICANCE: In this study we demonstrated the feasibility of correctly localizing the ictal-onset zone from iEEG recordings by using effective connectivity analysis during the first 20 s of ictal rhythmic iEEG activity.

Development and evaluation of an improved quantitative (90)Y 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.

Comparison of yttrium-90 quantitative imaging by TOF and non-TOF PET in a phantom of liver selective internal radiotherapy.

The aim of this study is to determine the feasibility of achieving quantitative measurement in (90)Y-microspheres liver selective internal radiotherapy (SIRT) by imaging (90)Y with a conventional non-time of flight (TOF) PET device. Instead of the bremsstrahlung x-rays of the ß-decay, the low branch of e(-)- e(+) pair production in the (90)Y-decay was used. The activity distribution in a phantom-simulated liver SIRT was obtained by direct (90)Y-PET imaging. We tested a LYSO TOF PET and two GSO and BGO non-TOF PET scanners using a 3.6-l cylindrical phantom filled with the (90)Y solution containing two sets of hot and cold spheres. The best hot contrast was obtained with the LYSO TOF. It was close to the expected value and remained constant, even for short acquisition times. The LYSO non-TOF was about 10% lower. The GSO performed similarly but degraded for shorter times whilst the BGO was the worst with 40% loss. For the cold spheres, the LYSO TOF and the GSO provided the best results, while the LYSO non-TOF and the BGO were the worst. (90)Y PET imaging in liver SIRT is achievable with LYSO TOF. Conventional LYSO and GSO show a loss of contrast and require longer acquisition times. BGO imaging is not feasible for dosimetry calculation.

Comparing planar image quality of rotating slat and parallel hole collimation: influence of system modeling.

The main remaining challenge for a gamma camera is to overcome the existing trade-off between collimator spatial resolution and system sensitivity. This problem, strongly limiting the performance of parallel hole collimated gamma cameras, can be overcome by applying new collimator designs such as rotating slat (RS) collimators which have a much higher photon collection efficiency. The drawback of a RS collimated gamma camera is that, even for obtaining planar images, image reconstruction is needed, resulting in noise accumulation. However, nowadays iterative reconstruction techniques with accurate system modeling can provide better image quality. Because the impact of this modeling on image quality differs from one system to another, an objective assessment of the image quality obtained with a RS collimator is needed in comparison to classical projection images obtained using a parallel hole (PH) collimator. In this paper, a comparative study of image quality, achieved with system modeling, is presented. RS data are reconstructed to planar images using maximum likelihood expectation maximization (MLEM) with an accurate Monte Carlo derived system matrix while PH projections are deconvolved using a Monte Carlo derived point-spread function. Contrast-to-noise characteristics are used to show image quality for cold and hot spots of varying size. Influence of the object size and contrast is investigated using the optimal contrast-to-noise ratio (CNR(o)). For a typical phantom setup, results show that cold spot imaging is slightly better for a PH collimator. For hot spot imaging, the CNR(o) of the RS images is found to increase with increasing lesion diameter and lesion contrast while it decreases when background dimensions become larger. Only for very large background dimensions in combination with low contrast lesions, the use of a PH collimator could be beneficial for hot spot imaging. In all other cases, the RS collimator scores better. Finally, the simulation of a planar bone scan on a RS collimator revealed a hot spot contrast improvement up to 54% compared to a classical PH bone scan.

Comparison of image quality of different iodine isotopes (I-123, I-124, and I-131).

I-131 is a frequently used isotope for radionuclide therapy. This technique for cancer treatment requires a pre-therapeutic dosimetric study. The latter is usually performed (for this radionuclide) by directly imaging the uptake of the therapeutic radionuclide in the body or by replacing it by one of its isotopes, which are more suitable for imaging. This study aimed to compare the image quality that can be achieved by three iodine isotopes: I-131 and I-123 for single-photon emission computed tomography imaging, and I-124 for positron emission tomography imaging. The imaging characteristics of each isotope were investigated by simulated data. Their spectrums, point-spread functions, and contrast-recovery curves were drawn and compared. I-131 was imaged with a high-energy all-purpose (HEAP) collimator, whereas two collimators were compared for I-123: low-energy high-resolution (LEHR) and medium energy (ME). No mechanical collimation was used for I-124. The influence of small high-energy peaks (>0.1%) on the main energy window contamination were evaluated. Furthermore, the effect of a scattering medium was investigated and the triple energy window (TEW) correction was used for spectral-based scatter correction. Results showed that I-123 gave the best results with a LEHR collimator when the scatter correction was applied. Without correction, the ME collimator reduced the effects of high-energy contamination. I-131 offered the worst results. This can be explained by the large amount of septal penetration from the photopeak and by the collimator, which gave a low spatial resolution. I-124 gave the best imaging properties owing to its electronic collimation (high sensitivity) and a short coincidence time window.

18 F-Fluorodeoxyglucose triple-headed coincidence detection imaging in oncology: preliminary results and comparison with dedicated PET.

BACKGROUND: Dual head coincidence (DHC) imaging has been proposed as a valuable and cheaper alternative to oncological PET. The increased sensitivity of the recently-developed triple-headed gamma camera (THC) optimized for coincidence detection has not been either validated or compared to the sensitivity of a dedicated PET system in a clinical study. MATERIAL/METHODS: Differences in contrast-to-noise-ratio (CNR) between dedicated FDG PET, DHC and THC imaging were assessed using a hot sphere phantom. Twenty-two oncological patients were prospectively studied with consecutive whole body FDG PET and FDG THC imaging. The images were visually read by 2 experienced nuclear medicine specialists. The diagnostic sensitivity of FDG THC imaging was assessed using FDG PET as the imaging gold standard. Lesion size was determined using computed tomography or magnetic resonance imaging. RESULTS: The mean gain in CNR for THC as compared to DHC imaging was 35% (10-56%). Of 63 lesions, 58 (92%) on FDG PET were also detected with THC imaging. Considering only lesions smaller than 15 mm, the relative sensitivity was 86% (24/29). All lesions above 15 mm (n=34) were detected using THC imaging. CONCLUSIONS: The first clinical results for THC imaging are promising. The overall relative sensitivity in this limited series was 92%. Only a few lesions smaller than 15 mm were missed.