Attenuation Correction Algorithm Research Articles

Current Applications of PET/MR: Part I: Technical Basics and Preclinical/Clinical Applications.

Positron emission tomography/magnetic resonance (PET/MR) imaging has gone through major hardware improvements in recent years, making it a reliable state-of-the-art hybrid modality in clinical practice. At the same time, image reconstruction, attenuation correction, and motion correction algorithms have significantly evolved to provide high-quality images. Part I of the current review discusses technical basics, pre-clinical applications, and clinical applications of PET/MR in radiation oncology and head and neck imaging. PET/MR offers a broad range of advantages in preclinical and clinical imaging. In the preclinic, small and large animal-dedicated devices were developed, making PET/MR capable of delivering new insight into animal models in diseases and facilitating the development of methods that inform clinical PET/MR. Regarding PET/MR's clinical applications in radiation medicine, PET and MR already play crucial roles in the radiotherapy process. Their combination is particularly significant as it can provide molecular and morphological characteristics that are not achievable with other modalities. In addition, the integration of PET/MR information for therapy planning with linear accelerators is expected to provide potentially unique biomarkers for treatment guidance. Furthermore, in clinical applications in the head and neck region, it has been shown that PET/MR can be an accurate modality in head and neck malignancies for staging and resectability assessment. Also, it can play a crucial role in diagnosing residual or recurrent diseases, reliably distinguishing from oedema and fibrosis. PET/MR can furthermore help with tumour characterization and patient prognostication. Lastly, in head and neck carcinoma of unknown origin, PET/MR, with its diagnostic potential, may obviate multiple imaging sessions in the near future.

Evaluating a radiotherapy deep learning synthetic CT algorithm for PET-MR attenuation correction in the pelvis

BackgroundPositron emission tomography–magnetic resonance (PET-MR) attenuation correction is challenging because the MR signal does not represent tissue density and conventional MR sequences cannot image bone. A novel zero echo time (ZTE) MR sequence has been previously developed which generates signal from cortical bone with images acquired in 65 s. This has been combined with a deep learning model to generate a synthetic computed tomography (sCT) for MR-only radiotherapy. This study aimed to evaluate this algorithm for PET-MR attenuation correction in the pelvis.MethodsTen patients being treated with ano-rectal radiotherapy received a 18\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$^{18}$$\\end{document}F-FDG-PET-MR in the radiotherapy position. Attenuation maps were generated from ZTE-based sCT (sCTAC) and the standard vendor-supplied MRAC. The radiotherapy planning CT scan was rigidly registered and cropped to generate a gold standard attenuation map (CTAC). PET images were reconstructed using each attenuation map and compared for standard uptake value (SUV) measurement, automatic thresholded gross tumour volume (GTV) delineation and GTV metabolic parameter measurement. The last was assessed for clinical equivalence to CTAC using two one-sided paired t tests with a significance level corrected for multiple testing of p≤0.05/7=0.007\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$p \\le 0.05/7 = 0.007$$\\end{document}. Equivalence margins of ±3.5%\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\pm 3.5\\%$$\\end{document} were used.ResultsMean whole-image SUV differences were −0.02% (sCTAC) compared to −3.0% (MRAC), with larger differences in the bone regions (−0.5% to −16.3%). There was no difference in thresholded GTVs, with Dice similarity coefficients ≥0.987\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\ge 0.987$$\\end{document}. However, there were larger differences in GTV metabolic parameters. Mean differences to CTAC in SUVmax\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${\\mathrm {SUV}}_{\\max}$$\\end{document} were 1.0±0.8%\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$1.0 \\pm 0.8\\%$$\\end{document} (± standard error, sCTAC) and -4.6±0.9%\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$-4.6 \\pm 0.9\\%$$\\end{document} (MRAC), and 1.0±0.7%\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$1.0 \\pm 0.7\\%$$\\end{document} (sCTAC) and -4.3±0.8%\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$-4.3 \\pm 0.8\\%$$\\end{document} (MRAC) in SUVmean\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${\\mathrm {SUV}}_{\\rm mean}$$\\end{document}. The sCTAC was statistically equivalent to CTAC within a ±3.5%\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\pm 3.5\\%$$\\end{document} equivalence margin for SUVmax\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${\\mathrm {SUV}}_{\\max}$$\\end{document} and SUVmean\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${\\mathrm {SUV}}_{\\rm mean}$$\\end{document} (p=0.007\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$p = 0.007$$\\end{document} and p=0.002\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$p = 0.002$$\\end{document}), whereas the MRAC was not (p=0.88\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$p = 0.88$$\\end{document} and p=0.83\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$p = 0.83$$\\end{document}).ConclusionAttenuation correction using this radiotherapy ZTE-based sCT algorithm was substantially more accurate than current MRAC methods with only a 40 s increase in MR acquisition time. This did not impact tumour delineation but did significantly improve the accuracy of whole-image and tumour SUV measurements, which were clinically equivalent to CTAC. This suggests PET images reconstructed with sCTAC would enable accurate quantitative PET images to be acquired on a PET-MR scanner.

Attenuation correction and truncation completion for breast PET/MR imaging using deep learning.

Objective. Simultaneous PET/MR scanners combine the high sensitivity of MR imaging with the functional imaging of PET. However, attenuation correction of breast PET/MR imaging is technically challenging. The purpose of this study is to establish a robust attenuation correction algorithm for breast PET/MR images that relies on deep learning (DL) to recreate the missing portions of the patient's anatomy (truncation completion), as well as to provide bone information for attenuation correction from only the PET data.Approach. Data acquired from 23 female subjects with invasive breast cancer scanned with18F-fluorodeoxyglucose PET/CT and PET/MR localized to the breast region were used for this study. Three DL models, U-Net with mean absolute error loss (DLMAE) model, U-Net with mean squared error loss (DLMSE) model, and U-Net with perceptual loss (DLPerceptual) model, were trained to predict synthetic CT images (sCT) for PET attenuation correction (AC) given non-attenuation corrected (NAC) PETPET/MRimages as inputs. The DL and Dixon-based sCT reconstructed PET images were compared against those reconstructed from CT images by calculating the percent error of the standardized uptake value (SUV) and conducting Wilcoxon signed rank statistical tests.Main results. sCT images from the DLMAEmodel, the DLMSEmodel, and the DLPerceptualmodel were similar in mean absolute error (MAE), peak-signal-to-noise ratio, and normalized cross-correlation. No significant difference in SUV was found between the PET images reconstructed using the DLMSEand DLPerceptualsCTs compared to the reference CT for AC in all tissue regions. All DL methods performed better than the Dixon-based method according to SUV analysis.Significance. A 3D U-Net with MSE or perceptual loss model can be implemented into a reconstruction workflow, and the derived sCT images allow successful truncation completion and attenuation correction for breast PET/MR images.

List-mode quantitative joint reconstruction of activity and attenuation maps in Time-of-Flight PET

Quantitative PET imaging requires accurate determination of patient-specific attenuation maps, which remains challenging on hybrid PET/MRI systems owing to the lack of a direct relationship between MR image intensity and attenuation coefficients. The aim of the present study is to develop a list-mode based algorithm for accurate and robust attenuation correction of PET data using time-of-flight (TOF) emission information. We analyze and address the challenges of list-mode emission-based maximum-likelihood joint estimation of activity and attenuation (LM-MLAA) in state-of-the-art PET imaging. The proposed method exploits a rapid on-the-fly system matrix calculation algorithm based on elliptic integrals while updating the attenuation map from accumulating list-mode coincidences to achieve accelerated image reconstruction. The scattering compensation is incorporated inside it using an iterative approach in such that the current estimation of attenuation map used on a course grid sampling scattering points to make an estimate of scattering. The performance of the proposed LM-MLAA approach was evaluated on Monte Carlo simulations of a phantom at different time resolutions. The contrast and noise for hot and cold regions on reconstructed images at different time resolutions were analysed. The estimated attenuation map exhibits resilience against noise, effectively eliminates high-frequency cross-talk even in the absence of prior information on attenuation coefficients, and enables discrimination among different anatomical regions in the reconstructed image. The error in the mean estimated attenuation coefficients after 50 iterations was ∼ 2% in water and ∼ -14% in Teflon regions for TOF resolutions corresponding to those of most current commercial PET systems (∼ 500 ps). The proposed LM-MLAA framework can be used for joint reconstruction of activity and attenuation maps from list-mode emission data as standalone or a complementary approach to existing in multimodality imaging such as PET/MRI, where direct measurement attenuation maps is not possible.

Attenuation Correction of the X-Band Dual-Polarization Phased Array Radar Based on Observed Raindrop Size Distribution Characteristics

X-band dual-polarization phased array radar (XPAR-D) possesses high resolution and plays a significant role in detecting meso- and micro-scale convective systems. However, the precipitation attenuation it endures necessitates an effective correction method. This study selected radar data from XPAR-D at the peak of Maofeng Mountain in Guangzhou during 16–17 May 2020 from three precipitation stages after quality control. Attenuation coefficients were calculated for different precipitation types through scattering simulations of raindrop size distribution (RSD) data. Next, an attenuation correction algorithm (MZH-KDP method) was proposed for the radar reflectivity factor (ZH) according to different raindrop types and compared to the ZH-KDP method currently in use. The results indicate that the attenuation amount of XPAR-D echoes depends on the attenuation path and echo intensity. When the attenuation path is shorter and the echo intensity is weaker, the amount of attenuation and correction is smaller. Difficulties arise when there are noticeable deviations, which are challenging to resolve using attenuation correction methods. Longer attenuation paths and stronger echoes highlight the advantages of the MZH-KDP method, while the ZH-KDP method tends to overcorrect the bias. The MZH-KDP method outperforms the ZH-KDP method for different precipitation types. The superior correction capability of the MZH-KDP method provides a significant advantage in improving the performance of XPAR-D for the detection of extreme weather.

Attenuation Correction in Weather Radars for Snow

Weather radars play a prominent role in remote sensing of the atmosphere. Various fields, such as meteorology and hydrology, rely on accurate weather radar data as input for their models. Different hydrometeors present during a weather event influence the amount of attenuation encountered by the radar signal. Attenuation correction for dual-polarization weather radar data is necessary to improve the radar products and get accurate measurements. Most of the existing attenuation correction research is associated with rain hydrometeors. Currently, research that addresses the attenuation correction of snow in weather radars is limited. Although it is known that attenuation of radar signals when it encounters rain is much greater than that for snow, attenuation for all hydrometeors needs to be addressed for accurate radar estimates. In this research work, the attenuation of different hydrometeors is studied using signal simulations. Various factors which influence attenuation, such as the elevation angle and particle size distribution, are considered, and the results are presented. An attenuation correction algorithm that uses the hydrometeor classification and specific differential phase products from the DROPS2.0 algorithm is introduced. Signal simulations are employed to obtain the relationship between specific attenuation and specific differential phase for different hydrometeors used in the proposed algorithm. The attenuation correction method is applied to X-band and Ku-band radar data. Path integrated attenuation of about 8 dB was observed in the snow case discussed from Ku band radar data. The method proposed for attenuation shows promising results at both frequency bands.

Evaluation of KDP Estimation Algorithm Performance in Rain Using a Known-Truth Framework

AbstractAccurate estimation of specific differential phase (KDP) is necessary for rain rate estimation, attenuation correction, and hydrometeor classification algorithms. There are numerous published methods to process polarimetric radar observations of propagation differential phase shift (ΦDP) and estimate KDP, but the corresponding KDP estimate uncertainty is unquantified. This study provides guidance on how commonly used KDP estimation algorithms perform in various environments. We create numerous synthetic (“true”) KDP profiles, integrate over them to obtain “smoothed” ΦDP, and then add noise typical of S-band operational weather radar measurements. Each algorithm is applied to our noisy ΦDP profiles and compared to the true KDP profile such that the errors and uncertainty are quantified. The synthetic KDP profiles are Gaussian in shape, which allows systematic variations in their magnitude and width to determine how each algorithm performs in smooth, slowly changing KDP profiles, as well as steep profiles. Results demonstrate that algorithm performance is dependent on the ΦDP field received. These results are further supported by an error analysis of each algorithm for two more complicated synthetic KDP profiles. Some KDP algorithms allow users to change various tuning parameters; a subset of these tuning parameters is tested to provide guidance on how changing these parameters impacts algorithm performance. We then provide evidence that our known-truth framework provides insight into algorithm performance in observed data through two case studies.

Deep learning-guided estimation of attenuation correction factors from time-of-flight PET emission data.

Attenuation correction (AC) is essential for quantitative PET imaging. In the absence of concurrent CT scanning, for instance on hybrid PET/MRI systems or dedicated brain PET scanners, an accurate approach for synthetic CT generation is highly desired. In this work, a novel framework is proposed wherein attenuation correction factors (ACF) are estimated from time-of-flight (TOF) PET emission data using deep learning. In this approach, referred to as called DL-EM), the different TOF sinogram bins pertinent to the same slice are fed into a multi-input channel deep convolutional network to estimate a single ACF sinogram associated with the same slice. The clinical evaluation of the proposed DL-EM approach consisted of 68 clinical brain TOF PET/CT studies, where CT-based attenuation correction (CTAC) served as reference. A two-tissue class consisting of background-air and soft-tissue segmentation of the TOF PET non-AC images (SEG) as a proxy of the technique used in the clinic was also included in the comparative evaluation. Qualitative and quantitative PET analysis was performed through SUV bias maps quantification in 63 different brain regions. The DL-EM approach resulted in 6.1±9.7% relative mean absolute error (RMAE) in bony structures compared to SEG AC method with RMAE of 16.1±8.2% (p-value <0.001). Considering the entire head region, DL-EM led to a root mean square error (RMSE) of 0.3±0.01 outperforming the SEG method with RMSE of 0.8±0.02 SUV (p-value <0.001). The region-wise analysis of brain PET studies revealed less than 7% absolute SUV bias for the DL-EM approach, whereas the SEG method resulted in more than 14% absolute SUV bias (p-value <0.05). Qualitative assessment and quantitative PET analysis demonstrated the superior performance of the DL-EM approach over the segmentation-based technique with clinically acceptable SUV bias. The results obtained using the DL-EM approach are comparable to state-of-the-art MRI-guided AC methods. Yet, this approach enables the extraction of interesting features about patient-specific attenuation which could be employed not only as a stand-alone AC approach but also as complementary/prior information in other AC algorithms.

Deep-JASC: joint attenuation and scatter correction in whole-body 18F-FDG PET using a deep residual network

We demonstrate the feasibility of direct generation of attenuation and scatter-corrected images from uncorrected images (PET-nonASC) using deep residual networks in whole-body 18F-FDG PET imaging. Two- and three-dimensional deep residual networks using 2D successive slices (DL-2DS), 3D slices (DL-3DS) and 3D patches (DL-3DP) as input were constructed to perform joint attenuation and scatter correction on uncorrected whole-body images in an end-to-end fashion. We included 1150 clinical whole-body 18F-FDG PET/CT studies, among which 900, 100 and 150 patients were randomly partitioned into training, validation and independent validation sets, respectively. The images generated by the proposed approach were assessed using various evaluation metrics, including the root-mean-squared-error (RMSE) and absolute relative error (ARE %) using CT-based attenuation and scatter-corrected (CTAC) PET images as reference. PET image quantification variability was also assessed through voxel-wise standardized uptake value (SUV) bias calculation in different regions of the body (head, neck, chest, liver-lung, abdomen and pelvis). Our proposed attenuation and scatter correction (Deep-JASC) algorithm provided good image quality, comparable with those produced by CTAC. Across the 150 patients of the independent external validation set, the voxel-wise REs (%) were - 1.72 ± 4.22%, 3.75 ± 6.91% and - 3.08 ± 5.64 for DL-2DS, DL-3DS and DL-3DP, respectively. Overall, the DL-2DS approach led to superior performance compared with the other two 3D approaches. The brain and neck regions had the highest and lowest RMSE values between Deep-JASC and CTAC images, respectively. However, the largest ARE was observed in the chest (15.16 ± 3.96%) and liver/lung (11.18 ± 3.23%) regions for DL-2DS. DL-3DS and DL-3DP performed slightly better in the chest region, leading to AREs of 11.16 ± 3.42% and 11.69 ± 2.71%, respectively (p value < 0.05). The joint histogram analysis resulted in correlation coefficients of 0.985, 0.980 and 0.981 for DL-2DS, DL-3DS and DL-3DP approaches, respectively. This work demonstrated the feasibility of direct attenuation and scatter correction of whole-body 18F-FDG PET images using emission-only data via a deep residual network. The proposed approach achieved accurate attenuation and scatter correction without the need for anatomical images, such as CT and MRI. The technique is applicable in a clinical setting on standalone PET or PET/MRI systems. Nevertheless, Deep-JASC showing promising quantitative accuracy, vulnerability to noise was observed, leading to pseudo hot/cold spots and/or poor organ boundary definition in the resulting PET images.

Interpretation of myocardial perfusion SPECT with attenuation correction. Part 2

Objective: to investigate the possibilities of a combined single photon emission computed tomography/computed tomography (SPECT/CT) system during myocardial perfusion studies and to examine the features of their interpretation. Subjects and methods. In April 2013 to April 2019, the Department of Radionuclide Diagnosis, National Medical Research Center of Cardiology, performed myocardial perfusion SPECT with attenuation correction (AC) in 3144 patients with various cardiovascular diseases. Based on the experience gained, the authors expanded and adjusted the principles of processing and interpreting the results of this study, which they had first set forth in 2014. Special emphasis is placed on the interpretation of apical and septal perfusion defects and diffuse perfusion irregularity, on the demonstration of the role of current reconstruction algorithms in obtaining highest-quality myocardial images, and on the assessment of self-potentialities of lowdose CT findings at myocardial SPECT/CT, as well as on proposals to standardize quantitative parameters for perfusion assessment at myocardial SPECT/CT. Results. AC affected in large measure the interpretation of myocardial perfusion SPECT in most cases. Apical perfusion defect appeared on AC images was observed in 50% of patients without established coronary heart disease and could be interpreted as a normal variant (apical thinning). The diffuse irregularity of perfusion radiopharmaceutical tracer accumulation in the application of the current reconstruction algorithms stood out in a separate pattern as a sign of microcirculatory disturbances. Low-dose CT data allow visualization of a severe thoracic abnormality at imaging and should be reflected in the conclusion. The perfusion study descriptions based on an analysis of AC and nAC images had a greater interoperator consistency than those based on an analysis of only nAC images (w = 0.915 and 0.809, respectively; p = 0.027). At the same time, the use of arithmetic means of the extent of transient ischemia (Reversibility Extent) is better consistent with visual analysis (= 0.08Ѓ}1.46%, p = 0.49). Conclusion. The use of AC and current iterative algorithms at myocardial perfusion SPECT should be part of a mandatory study protocol, since this affects to a large extent the diagnostic value of the technique. The borderline values of disturbed perfusion parameters at AC should be reconsidered.

Integrated Correction Algorithm for X Band Dual-Polarization Radar Reflectivity Based on CINRAD/SA Radar

The values of ratio a of the linear relationship between specific attenuation and specific differential phase vary significantly in convective storms as a result of resonance scattering. The best-linear-fit ratio a at X band is determined using the modified attenuation correction algorithm based on differential phase and attenuation, as well as the premise that reflectivity is unattenuated in S band radar detection. Meanwhile, the systemic reflectivity bias between the X band radar and S band radar and water layer attenuation (ZW) on the wet antenna cover of the X band radar are also considered. The good performance of the modified correction algorithm is demonstrated in a moderate rainfall event. The data were collected by four X band dual-polarization (X-POL) radar sites, namely, BJXCP, BJXFS, BJXSY, and BJXTZ, and a China’s New Generation Weather Radar (CINRAD/SA radar) site, BJSDX, in Beijing on 20 July 2016. Ratio a is calculated for each volume scan of the X band radar, with a mean value of 0.26 dB deg−1 varying from 0.20 to 0.31 dB deg−1. The average values of systemic reflectivity bias between the X band radar (at BJXCP, BJXFS, BJXSY, and BJXTZ) and S band radar (at BJSDX) are 0, −3, 2, and 0 dB, respectively. The experimentally determined ZW is in substantial agreement with the theoretically calculated ones, and their values are an order of magnitude smaller than rain attenuation. The comparison of the modified attenuation correction algorithm and the empirical-fixed-ratio correction algorithm is further evaluated at the X-POL radar. It is shown that the modified attenuation correction algorithm in the present paper provides higher correction accuracy for rain attenuation than the empirical-fixed-ratio correction algorithm.

Applications of QC and Merged Doppler Spectral Density Data from Ka-Band Cloud Radar to Microphysics Retrieval and Comparison with Airplane in Situ Observation

The new Chinese Ka-band solid-state transmitter cloud radar (CR) uses four operational modes with different pulse widths and coherent integration and non-coherent integration numbers to meet long-term cloud measurement requirements. The CR and an instrument-equipped aircraft were used to observe clouds and precipitation on the east side of Taihang Mountain in Hebei Province in 2018. To resolve the data quality problems caused by attenuation in the precipitation area; we focused on developing an algorithm for attenuation correction based on rain drop size distribution (DSD) retrieved from the merged Doppler spectral density data of the four operational modes following data quality control (QC). After dealiasing Doppler velocity and removal of range sidelobe artifacts; we merged the four types of Doppler spectral density data. Vertical air speed and DSD are retrieved from the merged Doppler spectral density data. Finally, we conducted attenuation correction of Doppler spectral density data and recalculated Doppler moments such as reflectivity; radial velocity; and spectral width. We evaluated the consistencies of reflectivity spectra from the four operational modes and DSD retrieval performance using airborne in situ observation. We drew three conclusions: First, the four operational modes observed similar reflectivity and velocity for clouds and low-velocity solid hydrometeors; however; three times of coherent integration underestimated Doppler reflectivity spectra for velocities greater than 2 m s−1. Reflectivity spectra were also underestimated for low signal-to-noise ratios in the low-sensitivity operational mode. Second, QC successfully dealiased Doppler velocity and removed range sidelobe artifacts; and merging of the reflectivity spectra mitigated the effects of coherent integration and pulse compression on radar data. Lastly, the CR observed similar DSD and liquid water content vertical profiles to airborne in situ observations. Comparing CR and aircraft data yielded uncertainty due to differences in observation space and temporal and spatial resolutions of the data.

Intercomparison of attenuation correction algorithms for single-polarized X-band radars

Attenuation due to liquid water is one of the largest uncertainties in radar observations. The effects of attenuation are generally inversely proportional to the wavelength, i.e. observations from X-band radars are more affected by attenuation than those from C- or S-band systems. On the other hand, X-band radars can measure precipitation fields in higher temporal and spatial resolution and are more mobile and easier to install due to smaller antennas.A first algorithm for attenuation correction in single-polarized systems was proposed by Hitschfeld and Bordan (1954) (HB), but it gets unstable in case of small errors (e.g. in the radar calibration) and strong attenuation. Therefore, methods have been developed that restrict attenuation correction to keep the algorithm stable, using e.g. surface echoes (for space-borne radars) and mountain returns (for ground radars) as a final value (FV), or adjustment of the radar constant (C) or the coefficient α. In the absence of mountain returns, measurements from C- or S-band radars can be used to constrain the correction. All these methods are based on the statistical relation between reflectivity and specific attenuation. Another way to correct for attenuation in X-band radar observations is to use additional information from less attenuated radar systems, e.g. the ratio between X-band and C- or S-band radar measurements. Lengfeld et al. (2016) proposed such a method based isotonic regression of the ratio between X- and C-band radar observations along the radar beam.This study presents a comparison of the original HB algorithm and three algorithms based on the statistical relation between reflectivity and specific attenuation as well as two methods implementing additional information of C-band radar measurements. Their performance in two precipitation events (one mainly convective and the other one stratiform) shows that a restriction of the HB is necessary to avoid instabilities. A comparison with vertically pointing micro rain radars (MRR) reveals good performance of two of the methods based in the statistical k-Z-relation: FV and α. The C algorithm seems to be more sensitive to differences in calibration of the two systems and requires additional information from C- or S-band radars.Furthermore, a study of five months of radar observations examines the long-term performance of each algorithm. From this study conclusions can be drawn that using additional information from less attenuated radar systems lead to best results. The two algorithms that use this additional information eliminate the bias caused by attenuation and preserve the agreement with MRR observations.

Implementation of absolute quantification in small-animal SPECT imaging: Phantom and animal studies.

PurposePresence of photon attenuation severely challenges quantitative accuracy in single‐photon emission computed tomography (SPECT) imaging. Subsequently, various attenuation correction methods have been developed to compensate for this degradation. The present study aims to implement an attenuation correction method and then to evaluate quantification accuracy of attenuation correction in small‐animal SPECT imaging.MethodsImages were reconstructed using an iterative reconstruction method based on the maximum‐likelihood expectation maximization (MLEM) algorithm including resolution recovery. This was implemented in our designed dedicated small‐animal SPECT (HiReSPECT) system. For accurate quantification, the voxel values were converted to activity concentration via a calculated calibration factor. An attenuation correction algorithm was developed based on the first‐order Chang's method. Both phantom study and experimental measurements with four rats were used in order to validate the proposed method.ResultsThe phantom experiments showed that the error of −15.5% in the estimation of activity concentration in a uniform region was reduced to +5.1% when attenuation correction was applied. For in vivo studies, the average quantitative error of −22.8 ± 6.3% (ranging from −31.2% to −14.8%) in the uncorrected images was reduced to +3.5 ± 6.7% (ranging from −6.7 to +9.8%) after applying attenuation correction.ConclusionThe results indicate that the proposed attenuation correction algorithm based on the first‐order Chang's method, as implemented in our dedicated small‐animal SPECT system, significantly improves accuracy of the quantitative analysis as well as the absolute quantification.

Benchmarking PET for geoscientific applications: 3D quantitative diffusion coefficient determination in clay rock

The 3D diagonal anisotropic effective diffusion coefficient of Na+, Deff=(Dxx, Dyy, Dzz), was quantified in a clay material in one single experiment/simulation. That is possible due to the combination of the non-invasive observation of Na+ diffusion in Opalinus clay by means of GeoPET method (PET: positron emission tomography) followed by quantitative 3D+t data evaluation by means of the finite element numerical modelling (FEM). The extracted anisotropic effective diffusion coefficient parallel (||) and normal (⊥) to the bedding of the clay rock, Deff=(D||, D⊥, D||) are comparable to those obtained on earlier experimental studies in the same clay material but with different methods. We consider this study as benchmark for the long-standing development of our GeoPET method, that explicitly includes a resolute and physics based attenuation and Compton scatter correction algorithm (Kulenkampff, J., M. Gründig, A. Zakhnini and J. Lippmann-Pipke (2016). "Geoscientific process monitoring with positron emission tomography (GeoPET)." Solid Earth 7: 1217–1231). We suggest GeoPET based fluid flow transport visualization combined with computer based process simulation henceforth as a qualified way for the quantification of three-dimensional, effective transport parameters in geosciences.

Regional Attenuation Correction of Weather Radar Using a Distributed Microwave-Links Network

The complex temporal-spatial variation of raindrop size distribution will affect the precision of precipitation quantitative estimates (QPE) produced from radar data, making it difficult to correct echo attenuation. Given the fact that microwave links can obtain the total path attenuation accurately, we introduce the concept of regional attenuation correction using a multiple-microwave-links network based on the tomographic reconstruction of attenuation coefficients. Derived from the radar-based equation, the effect of rainfall distribution on the propagation of radar and microwave link signals was analyzed. This article focuses on modeling of the tomographic reconstruction of attenuation coefficients and regional attenuation correction algorithms. Finally, a numerical simulation of regional attenuation correction was performed to verify the algorithms employed here. The results demonstrate that the correction coefficient (0.9175) falls between the corrected and initial field of radar reflectivity factor (root mean square error, 2.3476 dBz; average deviation, 0.0113 dBz). Compared with uncorrected data, the accuracy of the corrected radar reflectivity factor was improved by 26.12%, and the corrected rainfall intensity distribution was improved by 51.85% validating the region attenuation correction algorithm. This method can correct the regional attenuation of weather radar echo effectively and efficiently; it can be widely used for the radar attenuation correction and the promotion of quantitative precipitation estimation by weather radar.

A Simple Method for Attenuation Correction in Local X-Band Radar Measurements Using C-Band Radar Data

AbstractThis paper presents a novel, simple method to correct reflectivity measurements of weather radars that operate in attenuation-influenced frequency bands using observations from less attenuated radar systems. In recent years radar systems operating in the X-band frequency range have been developed to provide precipitation fields for areas of special interest in high temporal (≤1 min) and spatial (≤250 m) resolution in complement to nationwide radar networks. However, X-band radars are highly influenced by attenuation. C- and S-band radars typically have coarser resolution (250 m–1 km and 5 min) but are less affected by attenuation.Correcting for attenuation effects in simple (non-Doppler) single-polarized X-band radars remains challenging and is often dependent on restriction parameters, for example, those derived from mountain returns. Therefore, these algorithms are applicable only in limited areas. The method proposed here uses measurements from C-band radars and hence can be applied in all regions covered by nationwide C- (or S-) band radar networks. First, a single scan of X-band radar measurements is used exemplary to identify advantages and disadvantages of the novel algorithm compared to a standard single radar algorithm. The performance of the correction algorithms in different types of precipitation is examined in nine case studies. The proposed method provides very promising results for each type of precipitation. Additionally, it is evaluated in a 5-month comparison with Micro Rain Radar (MRR) observations. The bias between uncorrected X-band radar and MRR data is nearly eliminated by the attenuation correction algorithm, and the RMSE is reduced by 20% while the correlation of ~0.9 between both systems remains nearly constant.

Quantitative analysis of MRI-guided attenuation correction techniques in time-of-flight brain PET/MRI

In quantitative PET/MR imaging, attenuation correction (AC) of PET data is markedly challenged by the need of deriving accurate attenuation maps from MR images. A number of strategies have been developed for MRI-guided attenuation correction with different degrees of success. In this work, we compare the quantitative performance of three generic AC methods, including standard 3-class MR segmentation-based, advanced atlas-registration-based and emission-based approaches in the context of brain time-of-flight (TOF) PET/MRI. Fourteen patients referred for diagnostic MRI and (18)F-FDG PET/CT brain scans were included in this comparative study. For each study, PET images were reconstructed using four different attenuation maps derived from CT-based AC (CTAC) serving as reference, standard 3-class MR-segmentation, atlas-registration and emission-based AC methods. To generate 3-class attenuation maps, T1-weighted MRI images were segmented into background air, fat and soft-tissue classes followed by assignment of constant linear attenuation coefficients of 0, 0.0864 and 0.0975 cm(-1) to each class, respectively. A robust atlas-registration based AC method was developed for pseudo-CT generation using local weighted fusion of atlases based on their morphological similarity to target MR images. Our recently proposed MRI-guided maximum likelihood reconstruction of activity and attenuation (MLAA) algorithm was employed to estimate the attenuation map from TOF emission data. The performance of the different AC algorithms in terms of prediction of bones and quantification of PET tracer uptake was objectively evaluated with respect to reference CTAC maps and CTAC-PET images. Qualitative evaluation showed that the MLAA-AC method could sparsely estimate bones and accurately differentiate them from air cavities. It was found that the atlas-AC method can accurately predict bones with variable errors in defining air cavities. Quantitative assessment of bone extraction accuracy based on Dice similarity coefficient (DSC) showed that MLAA-AC and atlas-AC resulted in DSC mean values of 0.79 and 0.92, respectively, in all patients. The MLAA-AC and atlas-AC methods predicted mean linear attenuation coefficients of 0.107 and 0.134 cm(-1), respectively, for the skull compared to reference CTAC mean value of 0.138cm(-1). The evaluation of the relative change in tracer uptake within 32 distinct regions of the brain with respect to CTAC PET images showed that the 3-class MRAC, MLAA-AC and atlas-AC methods resulted in quantification errors of -16.2 ± 3.6%, -13.3 ± 3.3% and 1.0 ± 3.4%, respectively. Linear regression and Bland-Altman concordance plots showed that both 3-class MRAC and MLAA-AC methods result in a significant systematic bias in PET tracer uptake, while the atlas-AC method results in a negligible bias. The standard 3-class MRAC method significantly underestimated cerebral PET tracer uptake. While current state-of-the-art MLAA-AC methods look promising, they were unable to noticeably reduce quantification errors in the context of brain imaging. Conversely, the proposed atlas-AC method provided the most accurate attenuation maps, and thus the lowest quantification bias.

Preclinical models of atherosclerosis. The future of Hybrid PET/MR technology for the early detection of vulnerable plaque.

Cardiovascular diseases are the leading cause of death in developed countries. The aetiology is currently multifactorial, thus making them very difficult to prevent. Preclinical models of atherothrombotic diseases, including vulnerable plaque-associated complications, are now providing significant insights into pathologies like atherosclerosis, and in combination with the most recent advances in new non-invasive imaging technologies, they have become essential tools to evaluate new therapeutic strategies, with which can forecast and prevent plaque rupture. Positron emission tomography (PET)/computed tomography imaging is currently used for plaque visualisation in clinical and pre-clinical cardiovascular research, albeit with significant limitations. However, the combination of PET and magnetic resonance imaging (MRI) technologies is still the best option available today, as combined PET/MRI scans provide simultaneous data acquisition together with high quality anatomical information, sensitivity and lower radiation exposure for the patient. The coming years may represent a new era for the implementation of PET/MRI in clinical practice, but first, clinically efficient attenuation correction algorithms and research towards multimodal reagents and safety issues should be validated at the preclinical level.

Characterization of Mediterranean hail-bearing storms using an operational polarimetric X-band radar

Abstract. This work documents the effective use of X-band radar observations for monitoring severe storms in an operational framework. Two severe hail-bearing Mediterranean storms that occurred in 2013 in southern Italy, flooding two important Sicilian cities, are described in terms of their polarimetric radar signatures and retrieved rainfall fields. The X-band dual-polarization radar operating inside the Catania airport (Sicily, Italy), managed by the Italian Department of Civil Protection, is considered here. A suitable processing is applied to X-band radar measurements. The crucial procedural step relies on the differential phase processing, being preparatory for attenuation correction and rainfall estimation. It is based on an iterative approach that uses a very short-length (1 km) moving window, allowing proper capture of the observed high radial gradients of the differential phase. The parameterization of the attenuation correction algorithm, which uses the reconstructed differential phase shift, is derived from electromagnetic simulations based on 3 years of drop size distribution (DSD) observations collected in Rome (Italy). A fuzzy logic hydrometeor classification algorithm was also adopted to support the analysis of the storm characteristics. The precipitation field amounts were reconstructed using a combined polarimetric rainfall algorithm based on reflectivity and specific differential phase. The first storm was observed on 21 February when a winter convective system that originated in the Tyrrhenian Sea, marginally hit the central-eastern coastline of Sicily, causing a flash flood in Catania. Due to an optimal location (the system is located a few kilometers from the city center), it was possible to retrieve the storm characteristics fairly well, including the amount of rainfall field at the ground. Extemporaneous signal extinction, caused by close-range hail core causing significant differential phase shift in a very short-range path, is documented. The second storm, on 21 August 2013, was a summer mesoscale convective system that originated from a Mediterranean low pressure system lasting a few hours that eventually flooded the city of Syracuse. The undergoing physical process, including the storm dynamics, is inferred by analyzing the vertical sections of the polarimetric radar measurements. The high registered amount of precipitation was fairly well reconstructed, although with a trend toward underestimation at increasing distances. Several episodes of signal extinction were clearly manifested during the mature stage of the observed supercells.