Image quality of opportunistic breast examinations in photon-counting computed tomography: A phantom study.

PURPOSE: To compare the breast imaging performance of a clinical whole-body photon-counting CT (PCCT) to that of a dedicated breast CT (BCT) to determine the image quality of opportunistic breast examinations in clinical PCCT. MATERIALS AND METHODS: To quantify image quality for breast cancer applications, acquisitions of a breast phantom including representations of calcifications, fibers, and masses were performed using a clinical PCCT and a dedicated BCT. When imaging with the PCCT, the phantom was also combined with a thorax phantom to simulate realistic patient positioning, while only the breast phantom was imaged in the BCT. Images in BCT were acquired at 7.0 mGy (CTDI16cm) and using 2.6 mGy-25.0 mGy in the PCCT. Spatial resolution between the BCT and PCCT images was matched and data were reconstructed using the default methods of each system. The dose-normalized contrast-to-noise ratio (CNRD) of masses and the structural visibility of fibers and calcifications were evaluated as figures of merit for all reconstructions. RESULTS: CNRD between masses and background was 0.56 mGy-½, on average with BCT and varied between 0.39 mGy-½ to 1.46 mGy-½ with PCCT over all dose levels, phantom configurations, and reconstruction algorithms. Calcifications down to a size of 0.29 mm and fibers down to a size of 0.23 mm could be reliably identified in the images of both systems. CONCLUSIONS: Clinical PCCT provides an image quality superior to that obtained with BCT in terms of CNRD and allows for the identification of calcifications and fibers at comparable dose levels.

Threshold-dependent iodine imaging and spectral separation in a whole-body photon-counting CT system.

OBJECTIVE: To evaluate the dual-energy (DE) performance and spectral separation with respect to iodine imaging in a photon-counting CT (PCCT) and compare it to dual-source CT (DSCT) DE imaging. METHODS: A semi-anthropomorphic phantom extendable with fat rings equipped with iodine vials is measured in an experimental PCCT. The system comprises a PC detector with two energy bins (20 keV, T) and (T, eU) with threshold T and tube voltage U. Measurements using the PCCT are performed at all available tube voltages (80 to 140 kV) and threshold settings (50-90 keV). Further measurements are performed using a conventional energy-integrating DSCT. Spectral separation is quantified as the relative contrast media ratio R between the energy bins and low/high images. Image noise and dose-normalized contrast-to-noise ratio (CNRD) are evaluated in resulting iodine images. All results are validated in a post-mortem angiography study. RESULTS: R of the PC detector varies between 1.2 and 2.6 and increases with higher thresholds and higher tube voltage. Reference R of the EI DSCT is found as 2.20 on average overall phantoms. Maximum CNRD in iodine images is found for T = 60/65/70/70 keV for 80/100/120/140 kV. The highest CNRD of the PCCT is obtained using 140 kV and is decreasing with decreasing tube voltage. All results could be confirmed in the post-mortem angiography study. CONCLUSION: Intrinsically acquired DE data are able to provide iodine images similar to conventional DSCT. However, PCCT thresholds should be chosen with respect to tube voltage to maximize image quality in retrospectively derived image sets. KEY POINTS: • Photon-counting CT allows for the computation of iodine images with similar quality compared to conventional dual-source dual-energy CT. • Thresholds should be chosen as a function of the tube voltage to maximize iodine contrast-to-noise ratio in derived image sets. • Image quality of retrospectively computed image sets can be maximized using optimized threshold settings.

Photon-counting detectors in computed tomography: from quantum physics to clinical practice.

Over the last decade, a fundamentally new type of computed tomography (CT) detectors has proved its superior capabilities in both physical and preclinical evaluations and is now approaching the stage of clinical practice. These detectors are able to discriminate single photons and quantify their energy and are hence called photon-counting detectors. Among the promising benefits of this technology are improved radiation dose efficiency, increased contrast-to-noise ratio, reduced metal artifacts, improved spatial resolution, simultaneous multi-energy acquisitions, and the prospect of multi-phase imaging within a single acquisition using multiple contrast agents. Taking the conventional energy-integrating detectors as a reference, the authors demonstrate the technical principles of this new technology and provide phantom and patient images acquired by a whole-body photon-counting CT. These images serve as a basis for discussing the potential future of clinical CT.

Potential of ultra-high-resolution photon-counting CT of bone metastases: initial experiences in breast cancer patients.

Conventional CT scanners use energy-integrating detectors (EIDs). Photon-counting detector (PCD) computed tomography (CT) utilizes a CT detector technology based on smaller detector pixels capable of counting single photons and in addition discriminating their energy. Goal of this study was to explore the potential of higher spatial resolution for imaging of bone metastases. Four female patients with histologically confirmed breast cancer and bone metastases were included between July and October 2019. All patients underwent conventional EID CT scans followed by a high resolution non-contrast experimental PCD CT scan. Ultra-high resolution (UHR) reconstruction kernels were used to reconstruct axial slices with voxel sizes of 0.3 mm × 0.3 mm (inplane) × 1 mm (z-direction). Four radiologists blinded for patient identity assessed the images and compared the quality to conventional CT using a qualitative Likert scale. In this case series, we present images of bone metastases in breast cancer patients using an experimental PCD CT scanner and ultra-high-resolution kernels. A tendency to both a smaller inter-reader variability in the structural assessment of lesion sizes and in the readers' opinion to an improved visualization of lesion margins and content was observed. In conclusion, while further studies are warranted, PCD CT has a high potential for therapy monitoring in breast cancer.

A semi-automated quantitative comparison of metal artifact reduction in photon-counting computed tomography by energy-selective thresholding.

An evaluation of energy thresholding and acquisition mode for metal artifact reduction in Photon-counting detector CT (PCD-CT) compared to conventional energy-integrating detector CT (EID-CT) was performed. Images of a hip prosthesis phantom placed in a water bath were acquired on a scanner with PCD-CT and EID-CT (tube potentials: 100, 120 and 140 kVp) and energy thresholds (above 55-75 keV) in Macro and Chess mode. Only high-energy threshold images (HTI) were used. Metal artifacts were quantified by a semi-automated segmentation algorithm, calculating artifact volumes, means and standard deviations of CT numbers. Images of a human cadaver with hip prosthesis were acquired on the PCD-CT in Macro mode as proof-of-concept. Images at 140 kVp showed less metal artifacts than 120 kVp or 100 kVp. HTI (70, 75 keV) had fewer artifacts than low energy thresholds (55, 60, 65 keV). Fewer artifacts were observed in the Macro-HTI (8.9-13.3%) for cortical bone compared to Chess-HTI (9.4-19.1%) and EID-CT (10.7-19.0%) whereas in bone marrow Chess-HTI (19.9-45.1%) showed less artifacts compared to Macro-HTI (21.9-38.3%) and EID-CT (36.4-54.9%). Noise for PCD-CT (56-81 HU) was higher than EID-CT (33-36 HU) irrespective of tube potential. High-energy thresholding could be used for metal artifact reduction in PCD-CT, but further investigation of acquisition modes depending on target structure is required.

Neural networks-based regularization for large-scale medical image reconstruction.

In this paper we present a generalized Deep Learning-based approach for solving ill-posed large-scale inverse problems occuring in medical image reconstruction. Recently, Deep Learning methods using iterative neural networks (NNs) and cascaded NNs have been reported to achieve state-of-the-art results with respect to various quantitative quality measures as PSNR, NRMSE and SSIM across different imaging modalities. However, the fact that these approaches employ the application of the forward and adjoint operators repeatedly in the network architecture requires the network to process the whole images or volumes at once, which for some applications is computationally infeasible. In this work, we follow a different reconstruction strategy by strictly separating the application of the NN, the regularization of the solution and the consistency with the measured data. The regularization is given in the form of an image prior obtained by the output of a previously trained NN which is used in a Tikhonov regularization framework. By doing so, more complex and sophisticated network architectures can be used for the removal of the artefacts or noise than it is usually the case in iterative NNs. Due to the large scale of the considered problems and the resulting computational complexity of the employed networks, the priors are obtained by processing the images or volumes as patches or slices. We evaluated the method for the cases of 3D cone-beam low dose CT and undersampled 2D radial cine MRI and compared it to a total variation-minimization-based reconstruction algorithm as well as to a method with regularization based on learned overcomplete dictionaries. The proposed method outperformed all the reported methods with respect to all chosen quantitative measures and further accelerates the regularization step in the reconstruction by several orders of magnitude.

Iodine contrast-to-noise ratio improvement at unit dose and contrast media volume reduction in whole-body photon-counting CT.

PURPOSE: To assess the dose-normalized iodine contrast-to-noise-ratio (CNRD) improvement and contrast media reduction potential obtained with photon-counting (PC) CT compared to conventional energy-integrating (EI) CT as a function of patient size and tube voltage. METHOD: Images of a semi-anthropomorphic phantom of different sizes (small, medium, large) equipped with vials containing different iodine concentrations were acquired at the SOMATOM CounT prototype CT system using tube voltages of 80 kV-140 kV. CNRD is evaluated in reconstructions obtained using the EI detector, the PC detector using a single bin, and in reconstructions obtained by statistically optimally weighting acquisitions with two bins. Iodine CNRD improvements, potential dose reduction and the potential contrast media volume reduction are reported. RESULTS: In general, iodine CNRD improvement increases with increasing tube voltage for all patient sizes. In particular, if only one energy bin is used, the CNRD improvement is up to 30 % (small: 10 %, medium: 18 %, large: 30 %) and up to 37 % if an optimal weighting of two bins is performed (small: 13 %, medium: 25 %, large: 37 %) which is equivalent to the potential contrast media volume reduction. The improved iodine CNRD of PC compared to EI may allow for a potential radiation dose reduction of up to 46 %. CONCLUSIONS: All patients' iodine contrast at given x-ray dose, and particularly medium and large sized patients acquired at higher tube voltages, may benefit from photon-counting CT. The iodine contrast improvement can be used to reduce patient dose or to reduce the amount of contrast agent that is administered.

Cosmic-ray antinuclei as messengers of new physics: status and outlook for the new decade.

The precise measurement of cosmic-ray antinuclei serves as an important means for identifying the nature of dark matter and other new astrophysical phenomena, and could be used with other cosmic-ray species to understand cosmic-ray production and propagation in the Galaxy. For instance, low-energy antideuterons would provide a "smoking gun" signature of dark matter annihilation or decay, essentially free of astrophysical background. Studies in recent years have emphasized that models for cosmic-ray antideuterons must be considered together with the abundant cosmic antiprotons and any potential observation of antihelium. Therefore, a second dedicated Antideuteron Workshop was organized at UCLA in March 2019, bringing together a community of theorists and experimentalists to review the status of current observations of cosmic-ray antinuclei, the theoretical work towards understanding these signatures, and the potential of upcoming measurements to illuminate ongoing controversies. This review aims to synthesize this recent work and present implications for the upcoming decade of antinuclei observations and searches. This includes discussion of a possible dark matter signature in the AMS-02 antiproton spectrum, the most recent limits from BESS Polar-II on the cosmic antideuteron flux, and reports of candidate antihelium events by AMS-02; recent collider and cosmic-ray measurements relevant for antinuclei production models; the state of cosmic-ray transport models in light of AMS-02 and Voyager data; and the prospects for upcoming experiments, such as GAPS. This provides a roadmap for progress on cosmic antinuclei signatures of dark matter in the coming years.

AAfrag: Interpolation routines for Monte Carlo results on secondary production in proton-proton, proton-nucleus and nucleus-nucleus interactions.

We provide a compilation of predictions of the QGSJET-II-04m model for the production of secondary species (photons, neutrinos, electrons, positrons, and antinucleons) that are covering a wide range of energies of the beam particles in proton-proton, proton-nucleus, nucleus-proton, and nucleus-nucleus reactions. The current version of QGSJET-II-04m has an improved treatment of the production of secondary particles at low energies: the parameters of the hadronization procedure have been fine-tuned, based on a number of recent benchmark experimental data, notably, from the LHCf, LHCb, and NA61 experiments. Our results for the production spectra are made publicly accessible through the interpolation routines AAfrag which are described below. Besides, we comment on the impact of Feynman scaling violation and isospin symmetry effects on antinucleon production.

Solution of Heliospheric Propagation: Unveiling the Local Interstellar Spectra of Cosmic-ray Species.

Local interstellar spectra (LIS) for protons, helium, and antiprotons are built using the most recent experimental results combined with state-of-the-art models for propagation in the Galaxy and heliosphere. Two propagation packages, GALPROP and HelMod, are combined to provide a single framework that is run to reproduce direct measurements of cosmic-ray (CR) species at different modulation levels and at both polarities of the solar magnetic field. To do so in a self-consistent way, an iterative procedure was developed, where the GALPROP LIS output is fed into HelMod, providing modulated spectra for specific time periods of selected experiments to compare with the data; the HelMod parameter optimization is performed at this stage and looped back to adjust the LIS using the new GALPROP run. The parameters were tuned with the maximum likelihood procedure using an extensive data set of proton spectra from 1997 to 2015. The proposed LIS accommodate both the low-energy interstellar CR spectra measured by Voyager 1 and the high-energy observations by BESS, Pamela, AMS-01, and AMS-02 made from the balloons and near-Earth payloads; it also accounts for Ulysses counting rate features measured out of the ecliptic plane. The found solution is in a good agreement with proton, helium, and antiproton data by AMS-02, BESS, and PAMELA in the whole energy range.

A Supernova at 50 pc: Effects on the Earth's Atmosphere and Biota.

Recent 60Fe results have suggested that the estimated distances of supernovae in the last few million years should be reduced from ∼100 to ∼50 pc. Two events or series of events are suggested, one about 2.7 million years to 1.7 million years ago, and another about 6.5-8.7 million years ago. We ask what effects such supernovae are expected to have on the terrestrial atmosphere and biota. Assuming that the Local Bubble was formed before the event being considered, and that the supernova and the Earth were both inside a weak, disordered magnetic field at that time, TeV-PeV cosmic rays (CRs) at Earth will increase by a factor of a few hundred. Tropospheric ionization will increase proportionately, and the overall muon radiation load on terrestrial organisms will increase by a factor of ∼150. All return to pre-burst levels within 10 kyr. In the case of an ordered magnetic field, effects depend strongly on the field orientation. The upper bound in this case is with a largely coherent field aligned along the line of sight to the supernova, in which case, TeV-PeV CR flux increases are ∼104; in the case of a transverse field they are below current levels. We suggest a substantial increase in the extended effects of supernovae on Earth and in the "lethal distance" estimate; though more work is needed. This paper is an explicit follow-up to Thomas et al. We also provide more detail on the computational procedures used in both works.

TERRESTRIAL EFFECTS OF NEARBY SUPERNOVAE IN THE EARLY PLEISTOCENE.

Recent results have strongly confirmed that multiple supernovae happened at distances of ∼100 pc, consisting of two main events: one at 1.7-3.2 million years ago, and the other at 6.5-8.7 million years ago. These events are said to be responsible for excavating the Local Bubble in the interstellar medium and depositing 60Fe on Earth and the Moon. Other events are indicated by effects in the local cosmic ray (CR) spectrum. Given this updated and refined picture, we ask whether such supernovae are expected to have had substantial effects on the terrestrial atmosphere and biota. In a first look at the most probable cases, combining photon and CR effects, we find that a supernova at 100 pc can have only a small effect on terrestrial organisms from visible light and that chemical changes such as ozone depletion are weak. However, tropospheric ionization right down to the ground, due to the penetration of ⩾TeV CRs, will increase by nearly an order of magnitude for thousands of years, and irradiation by muons on the ground and in the upper ocean will increase twentyfold, which will approximately triple the overall radiation load on terrestrial organisms. Such irradiation has been linked to possible changes in climate and increased cancer and mutation rates. This may be related to a minor mass extinction around the Pliocene-Pleistocene boundary, and further research on the effects is needed.

Signatures of a Two Million Year Old Supernova in the Spectra of Cosmic Ray Protons, Antiprotons, and Positrons.

The locally observed cosmic ray spectrum has several puzzling features, such as the excess of positrons and antiprotons above ~20 GeV and the discrepancy in the slopes of the spectra of cosmic ray protons and heavier nuclei in the TeV-PeV energy range. We show that these features are consistently explained by a nearby source which was active approximately two million years ago and has injected (2-3)×10^{50} erg in cosmic rays. The transient nature of the source and its overall energy budget point to the supernova origin of this local cosmic ray source. The age of the supernova suggests that the local cosmic ray injection was produced by the same supernova that has deposited ^{60}Fe isotopes in the deep ocean crust.

The retrobulbar sinus is superior to the lateral tail vein for the injection of contrast media in small animal cardiac imaging.

Cardiac perfusion studies using computed tomography are a common tool in clinical practice. Recent technical advances and the availability of dedicated small animal scanners allow the transfer of these techniques to the preclinical sector in general and to mouse models of cardiac diseases in particular. This necessitates new requirements for contrast injection techniques as a rapid transport of contrast media from the intravenous access to the animal heart. Clinical contrast agents containing high iodine concentrations are used within small animal studies although they exhibit a high viscosity which might limit their transport within the vasculature. The authors provide a comparison of the transport of contrast media following an injection into the lateral tail vein and an injection into the retrobulbar sinus and discuss the anatomy involved. The temporal evolution of a contrast bolus and its in vivo distribution is visualized. It is demonstrated that injecting contrast agents into the lateral tail vein of mice results in a retrograde blood flow to the liver veins and therefore does not deliver a detectable contrast bolus to the heart, and thus it cannot be used for cardiac perfusion studies. By contrast, boli injected into the retrobulbar sinus are rapidly transported to the heart and provide ventricular contrast enabling perfusion studies similar to those in human patients. The results demonstrate that an injection into the retrobulbar sinus is superior to an injection into the lateral tail vein for the delivery of contrast boli to the animal heart, while all drawbacks of an injection into the lateral tail vein are overcome.

Constrained reconstructions for 4D intervention guidance.

Image-guided interventions are an increasingly important part of clinical minimally invasive procedures. However, up to now they cannot be performed under 4D (3D + time) guidance due to the exceedingly high x-ray dose. In this work we investigate the applicability of compressed sensing reconstructions for highly undersampled CT datasets combined with the incorporation of prior images in order to yield low dose 4D intervention guidance. We present a new reconstruction scheme prior image dynamic interventional CT (PrIDICT) that accounts for specific image features in intervention guidance and compare it to PICCS and ASD-POCS. The optimal parameters for the dose per projection and the numbers of projections per reconstruction are determined in phantom simulations and measurements. In vivo experiments in six pigs are performed in a cone-beam CT; measured doses are compared to current gold-standard intervention guidance represented by a clinical fluoroscopy system. Phantom studies show maximum image quality for identical overall doses in the range of 14 to 21 projections per reconstruction. In vivo studies reveal that interventional materials can be followed in 4D visualization and that PrIDICT, compared to PICCS and ASD-POCS, shows superior reconstruction results and fewer artifacts in the periphery with dose in the order of biplane fluoroscopy. These results suggest that 4D intervention guidance can be realized with today's flat detector and gantry systems using the herein presented reconstruction scheme.

Frequency split metal artefact reduction in pelvic computed tomography.

OBJECTIVES: Artefacts from total hip replacement affect image quality and the visualisation of pelvic lesions on computed tomography (CT). We propose a frequency split (FS) approach in addition to the normalised metal artefact reduction (NMAR) algorithm that aims to suppress artefacts and improves image quality in patients with orthopaedic hardware. METHODS: Data from ten consecutive patients with metal artefacts from uni- (n = 5) and bilateral (n = 4) total hip replacement or osteosynthesis (n = 1) were reconstructed with filtered back projection (FBP), linear interpolation MAR (LIMAR), NMAR, FSLIMAR and FSNMAR and analysed for image quality and severity of artefacts. RESULTS: NMAR and FSNMAR significantly improved the assessment of the pelvic organs, lymph nodes and vessels compared with FBP, LIMAR or FSLIMAR (P < 0.05). Assessment of the metal hardware, joint and capsule was improved with the addition of FS (FSLIMAR, FSNMAR). No algorithm-related artefacts were detected in regions that did not contain metal. CONCLUSIONS: NMAR, FSLIMAR and FSNMAR have the potential to improve image quality in patients with artefacts from metal hardware and to improve the diagnostic accuracy of CT of the organs of the pelvis. Although introducing some algorithm-related artefacts, FSNMAR most accurately displayed adjacent bone and tissue next to metal implants. KEY POINTS: • Orthopaedic metallic hardware often creates serious artefacts in computed tomography, hindering diagnosis. • The normalised metal artefact reduction (NMAR) algorithm was developed to suppress such artefacts. • NMAR improves CT assessment of pelvic organs in patients with orthopaedic hardware. • Addition of the frequency split technique (FSNMAR) helps assess tissue near metal hardware. • NMAR and FSNMAR are robust and computationally effective sinogram interpolation algorithms.

Filamentary diffusion of cosmic rays on small scales.

We investigate the diffusion of cosmic rays (CRs) close to their sources. Propagating individual CRs in purely isotropic turbulent magnetic fields with maximal scale of spatial variations l(max), we find that CRs diffuse anisotropically at distances r ≤/~ l(max) from their sources. As a result, the CR densities around the sources are strongly irregular and show filamentary structures. We determine the transition time t(*) to standard diffusion as t(*) ~ 10(4) yr(l(max)/150 pc)(ß)(E/PeV)(-γ)(B(rms)/4 µG)(γ), with ß =/~ 2 and γ=0.25-0.5 for a turbulent field with a Kolmogorov power spectrum. We calculate the photon emission due to CR interactions with gas and the resulting irregular source images.

The influence of the heel effect in cone-beam computed tomography: artifacts in standard and novel geometries and their correction.

For decades, the heel effect has been known to cause an angular dependence of the emitted spectrum of an x-ray tube. In radiography, artifacts were observed and attributed to the heel effect. However, no problems due to the heel effect were discerned in multi-slice computed tomography (MSCT) so far. With flat-detector CT (FDCT), involving larger cone angles and different system geometries, the heel effect might cause new artifacts. These artifacts were analyzed in this paper for system geometries different from the ones widely used nowadays. Simulations and measurements were performed. Simulations included symmetric as well as asymmetric detector layouts and different x-ray tube orientations with respect to the detector plane. The measurements were performed on a micro-CT system in an asymmetric detector layout. Furthermore, an analytical correction scheme is proposed to overcome heel effect artifacts. It was shown that the type of artifact greatly depends on the orientation of the x-ray tube and also on the type of detector alignment (i.e. symmetric or different types of asymmetric alignment). Certain combinations exhibited almost no significant artifact while others greatly influenced the quality of the reconstructed images. The proposed correction scheme showed good results that were further improved when also applying a scatter correction. When designing CT systems, care should be taken when placing the tube and the detector. Orientation of the x-ray tube like in most MSCT systems seems advisable in asymmetric detector layouts. However, a different type of tube orientation can be overcome with suitable correction schemes.

Measuring the 13 neutrino mixing angle and the CP phase with neutrino telescopes.

The observed excess of high-energy cosmic rays from the Galactic plane in the energy range around 10(18) eV may be explained by neutron primaries generated in the photodissociation of heavy nuclei. In this scenario, lower-energy neutrons decay before reaching the Earth and produce a detectable flux in a 1 km(3) neutrino telescope. The initial flavor composition of the neutrino flux, phi(nu(e)):phi(nu(mu)):phi(nu(tau))=1:0:0, permits a combined nu(mu)/nu(tau) appearance and nu(e) disappearance experiment. The observable flux ratio phi(nu(mu))/phi(nu(e)+nu(tau) at Earth depends on the 13 mixing angle theta(13) and the leptonic CP phase delta(CP), thus opening a new way to measure these two quantities.

Reconstruction from truncated projections in CT using adaptive detruncation.

If the object exceeds the field of measurement (FOM) of a given CT scanner, severe artifacts may result. In this work, we propose an adaptive detruncation (ADT) method to reconstruct images from medical CT projections which are truncated in the transaxial direction. The truncated projections are extrapolated by estimating the convex hull of the patient. The ADT method allows us not only to achieve artifact-free images in the FOM but also to extend the images beyond the FOM, and can therefore be very attractive, for example, in PET/CT scanners for attenuation correction.