Abstract
Background/Aim: Multiparametric dual energy comptuted tomography (CT) imaging allows for multidimensional tissue characterization beyond the measurement of Hounsfield units. The purpose of this study was to evaluate multiple imaging parameters for different abdominal organs in dual energy CT (DECT) and analyze the effects of the contrast agent on these different parameters and provide normal values for characterization of parenchymatous organs. Patients and Methods: This retrospective analysis included a total of 484 standardized DECT scans of the abdomen. Hounsfield Units (HU), rho (electron density relative to water), Zeff (effective atomic number) and FF (fat fraction) were evaluated for liver, spleen, kidney, muscle, fat-tissue. Independent generalized estimation equation models were fitted. Results: In DECT imaging there is only little difference in mean HUmixed for parenchymatous abdominal organs. Analysis including Zeff, rho and FF allows for better discrimination while a large overlap remains for liver, spleen and muscle. Including multidimensional analysis and the effects of contrast medium further enhances tissue characterization. Small differences remain for liver and spleen. Conclusion: Organ characterization using multiparametric dual energy CT analysis is possible. An increased number of parameters obtained from DECT improves organ characterization. To our knowledge this is the first attempt to provide normal values for characterization of parenchymatous organs.
The differentiation of tissue and organs on single-energy computed-tomography (CT) is performed according to differences in X-ray attenuation, which in the case of CT is measured as Hounsfield units (HU) relative to the attenuation of pure water (1, 2). The CT number is influenced by the effective atomic number (Zeff) and the electron density of each material (ρe). As the effective atomic number is dependent to spectral properties, a change in X-ray energy leads to a change in the resulting Hounsfield units for a tissue accordingly (3-6).
The basic idea of dual-energy CT (DECT) is to apply two different X-ray energies and hence take advantage of the different spectral properties. This method has become established over recent years and offers the possibility of new and better tissue characterizations (7-15). By acquiring two CT datasets either simultaneously or one immediately after another, DECT enables the concurrent acquisition of two attenuation maps at low- and high-energy spectrums (10, 15-17). Although different tissues may show quite similar attenuation (in terms of CT numbers) at a certain energy level (7, 8), they may show large differences in attenuation at other energy levels because of their individual electron binding energies (9, 10, 12). The two main reasons for this effect are Compton-scattering and photoelectric-absorption (5, 6, 11-13, 15, 18). Compton scatter is nearly independent of the photon-energy, depending primarily on the electron density of the material. It occurs predominately at high energies. However, photoelectric-absorption is highly dependent on the electron binding energy, which in turn is highly dependent on the atomic number of the examined element (3-6, 13, 18).
The probability of photoelectric-absorption increases substantially when the energy of a photon approaches the individual binding energy of a material. On the basis of these effects, different tissues may show large differences in their probabilities for photoelectric interactions, and the individual spectral characteristics may assist in the characterization of different tissues on DECT (9-12, 14-18).
The purpose of this study was to evaluate multiparametric analysis of different tissues on DECT-scans implemented in a standardized clinical routine setting over the course of a year. The tissues evaluated were abdominal organs (liver, spleen, kidney) and muscle, as all these materials have quite similar elemental compositions and may be challenging for automated segmentation methods using single-energy CT images. Additionally, abdominal fat tissue was measured, which should be easier to separate from soft tissue because of its lower density (1). By using dedicated software (syngo.via® Dual Energy, Siemens Healthineers, Forchheim, Germany), CT images acquired at two different energy levels can be decomposed into two new material-specific images by the application of spectral processing techniques. The results can be presented as either Zeff and ρe images, or as one image representing high atomic number material (e.g., iodine) and the other image representing soft tissue, which also permits analysis of the fat fraction. A future idea is to create automated segmentation techniques for abdominal imaging, using dedicated trained software for the automatic recognition of abdominal organs.
We, furthermore, analyzed the effect of a contrast agent on the different DECT parameters of each organ, as this information could contribute to future studies researching pathological changes in organs. It is likely that angiopathic-ischemic lesions with reduced blood circulation or blood clots in the vascular system of an organ will lead to low contrast agent uptake in the respective organ.
Patients and Methods
This retrospective study was performed using routinely acquired clinical DECT-scans of human participants. The requirement for informed patient consent for this retrospective study was waived by the local ethics committee (EKNZ 2018-01641). In addition, written informed consent was obtained for the research use of the clinical data of each participant.
Study population. This study was performed at a public hospital using data acquired from patients between February 2018 and February 2019. The institution has subspecialty-trained radiologists dedicated to the analysis of CT and DECT images. The inclusion criteria for subject selection were patients above 18 years who provided written informed consent and underwent one of two defined standardized DECT-protocols. All patients were referred for DECT due to suspected kidney/ureter stones or work-up of hematuria. Sixteen DECT scans from 16 patients were excluded because the patients received two non-contrast and two contrast-enhanced DECT scans during the study period. In total, imaging data from 331 different patients were included in the study.
Imaging modalities. Dual-energy CT was performed with a SOMATOM Definition Flash second-generation dual-source CT scanner (Siemens Healthineers). Two standardized DECT protocols for the abdomen (non-contrast and contrast-enhanced) were included in the study. The standard-of-care imaging was a non-contrast DECT scan of the abdomen in patients with suspected kidney/ureter stones and a combined non-contrast and venous-phase contrast DECT scan of the abdomen for work-up of hematuria. Axial spiral datasets of the abdomen were acquired during a single breath-hold.
The scanning and reconstruction parameters for the unenhanced images were: kV 100/Sn 140, reference mAs 70, rotation time 0.5 s, pitch 0.7, collimation 32×0.6 mm, and reconstruction 1.5 mm. The scanning and reconstruction parameters for the contrast-enhanced images were: kV 100/Sn 140, reference mAs 94, rotation time 0.5 s, pitch 1.2, collimation 32×0.6 mm, and reconstruction 1 mm.
Image evaluation. All dual-energy images were evaluated using syngo.via® software (syngo Dual Energy, Siemens Healthcare GmbH 2009-2018, Version 05.01.0000.0030, Erlangen, Germany). Standardized segments of each organ were measured by one reader using standardized regions of interest (ROIs). The measured DECT-parameters included HU values of the mixed 100 kV and 140 Sn kV images (HUmixed; representing a 120 kV image), electron density relative to water (ρe), effective atomic number (Zeff), and fat fraction (FF). Image analysis and measurement processing were performed under the supervision of a board-certified radiologist.
The readers were free to choose optimal window settings for the image analysis. After visual assessment, the parameters were measured using ROIs set in the liver, spleen, kidney, muscle, and fat tissue. All parameters were evaluated on DECT-scans with and without contrast agent. For each tissue, at least two ROIs of between 0.5 and 4 cm2 were placed on axial image views. Contact with macroscopic vessels or organ capsules were avoided. In the liver, three ROIs were placed in the periphery of the right hepatic lobe (hepatic segments VIII, VI, V) and one ROI was placed in the left hepatic lobe (III). In the spleen, one ROI was set in the upper half and one in the lower half. In each kidney, three ROIs were set in the cortex between the upper and lower pole, avoiding contact with the renal pelvis. For muscle measurements, the psoas, obturatorius, sartorius, or gluteus muscles were measured on both sides, avoiding muscles showing signs of fatty degeneration or atrophy in the image reporting. Fat tissue was measured in subcutaneous adipose tissue, preferably in the lumbar region. As several ROIs were placed in each organ, the measurements were averaged for each participant and organ.
Figure 1 shows the application profiles used in syngo.via® software. Measurements of attenuation were taken on default view images. Afterwards, the specific DECT-parameters (ρe, Zeff, FF) of each tissue were extracted using the application profiles “Rho/Z” and “fat map” (subcategories of the application profile “Liver VNC”). Both profiles create material decomposition images and display all structures in a color-coded overlay. The application class “Liver VNC” is based on an iodine subtraction using a material decomposition algorithm encoding the iodine of the contrast agent (19, 20). Consequently, every voxel is decomposed into iodine, fat, and soft tissue maps. The Liver VNC application offers the possibility of quantifying the fat fraction of each material in a certain ROI using the sub-algorithm “fat map”. The application class “Rho/Z” provides Zeff directly, where Rho is encoded in HU-like units (HURho), and is converted into electron density values by applying the linear relationship from Saito et al., as shown below (21):
Statistical analysis. Statistical analysis was performed using R software version 4.0.0 (2020-04-24). Summary statistics for each organ were created separately, with and without contrast agent, and included the number of observations on which the calculations were based, minimum, 1st and 3rd quartile, median, mean, maximum, standard deviation, and interquartile range. Potential interactions between organs and contrast agent were taken into account by performing analyses separately with and without contrast agent. Independent generalized estimation equation models (GEEs) were fitted to quantify the differences between the organs with respect to the parameters HUmixed, FF, ρe, and Zeff, and to assess the effects of contrast agent on the parameters. The reported (unadjusted) 95% Wald-type confidence intervals are based on the sandwich variance estimator, thus taking intraindividual correlation into account. As there were at most two measurements, the choice of correlation structure was irrelevant.
Results
Study participants. The full data set consisted of a total of 2,420 measurements from 331 patients. One hundred and seventy-eight patients underwent DECT with contrast agent (n=40) or without contrast agent (n=138), whereas 153 patients underwent DECT scans with and without contrast agent during the same session. In total, we report on 291 non-contrast DECT scans and 193 contrast-enhanced DECT scans. Measurements of the spleen were not possible in some cases because of splenectomy. Figure 2 presents a flow chart of the patient selection process.
Data analysis. All parameters were successfully measured in all regions of interest in each organ or tissue. The measured HURho values of the “Rho/Z” application profile were successfully transformed into the electron density relative to water (ρe) using the formula described by Saito et al. (21). The color-coded fat overlay images of the application profile “Liver VNC ‐ fat map” showed visually clear differences between the different included tissues, with different levels of fat content. Abdominal organs were displayed in blue (no or low fat) and fat tissue in yellow to black (fat).
Statistical analysis. The summary statistics for the DECT scans without contrast agent, representing the number of observations on which the calculations are based, minimum, 1st and 3rd quartiles, median, mean, maximum, standard deviation, and interquartile range, are shown in Table I. Summary statistics for DECT scans with contrast agent are shown in Table II.
The effect of contrast agent on the variables was separately assessed for each organ, and the results are shown in Table III. We report the estimated difference between non-contrast and contrast-enhanced DECT scans for each parameter and organ, as well as the 95% confidence intervals. Taking every variable into account, muscle and fat tissue showed the lowest differences, whereas kidney and spleen showed the highest differences when non-contrast and contrast-enhanced results were compared.
Between-organ differences in the parameters were quantified separately for DECT scans performed with and without contrast agent. Results of the analyses of scans without contrast agent are shown in Table IV, whereas those with contrast agent are shown in Table V. Again, we report the estimated difference as well as the 95% confidence intervals. Comparisons of the different organs within non-contrast DECT acquisitions revealed only low differences in HUmixed between liver and muscle, and in Zeff between liver, muscle, and spleen. On contrast-enhanced DECT scans, low differences were shown in HUmixed between liver and spleen, in Zeff between liver and spleen, and in FF between liver, spleen, and muscle. The value ranges given above showed large overlaps for the different tissues within the same acquisition. Overlaps are visualized in graphs in Figure 3. We also inserted the measured values of HUmixed, Zeff and FF for liver, spleen, kidney and muscle in a 3D graph (Figure 2). Especially the overlaps of liver, spleen and muscle can be regarded here. In Figure 3 and Figure 4, fat tissue was not included in the visualization since all values of fat were highly significant different from all other tissues.
Comparing the mean values in non-contrast DECT scans of HUmixed only showed no significant differences by comparing spleen to muscle (p=0.99). Mean values of ρe showed no significant differences by comparing liver to muscle (p=0.24), Zeff by liver to spleen (p=0.32) and liver to muscle (p=0.16). For FF and all other tissue comparisons results showed statistically relevant differences (p<0.05).
Comparison of the mean HUmixed values in contrast-enhanced scans showed no significant differences between liver to spleen (p=0.96), whereas all other comparisons were highly significant (p<0.001).
Discussion
This retrospective analysis of clinical dual-energy CT scans (DECT) reports a statistical summary of CT values for different abdominal tissues with and without the use of contrast agent in a 12-month cohort. Furthermore, we evaluated the estimated differences between abdominal organs, as well the effects of a contrast agent on different CT parameters. A current focus of interest in radiology is the idea of automated segmentation techniques that can automatically recognize different organs and tissues. Therefore, it is of great importance to gather data from many clinical study populations and determine the generally observed distributions to derive normal values. Furthermore, we were able to document the effects of contrast agent on the DECT values of abdominal organs. We are convinced that angiopathic-ischemic lesions with reduced blood circulation, or blood clots in the vascular system of an organ, lead to low uptake of contrast agent by the respective organ. Thus, estimating the effects of contrast agent on different DECT variables is of major interest.
By applying the DECT method with low- and high-energy spectra, we obtained additional parameters for tissue characterization, such as the CT number for mixed images at different tube voltages (HUmixed), the electron density relative to water (ρe), the effective atomic number (Zeff), and the fat fraction (FF). In contrast, only HU values are generated on single-energy CT, which is a disadvantage, because tissues are generally likely to show more differences if more variables are involved. Thus, we examined DECT scans to obtain more knowledge on organic structures and improve organ characterization.
We could determine the effects of the contrast agent on each DECT parameter in each organ. Whereas muscle and fat tissue showed lowest changes in contrast-enhanced scans for each parameter due to their low intake of contrast agent, the kidneys and spleen showed the highest changes of each included organ. This is most likely a result of both organs having a high blood circulation, but also because the CT protocol used was dedicated for work-up of the urinary track, and contrast agent was therefore administered using a split bolus technique. Consequently, our results support the assumption that pathologically ischemic organic changes are detectable using our DECT protocol. Depending on the units of the measured variables, Zeff and HUmixed were the most affected by contrast agent because both parameters are affected by the atomic number of elements. Thus, the effects of the iodine-based (atomic number of iodine, Z=53) contrast agent were measurable in the collected DECT variables. We generated the 95% confidence intervals for each organ and DECT variable before and after contrast agent application, which allows the values obtained from patients to be checked, and those patients with values lower or higher than the confidence interval should be observed more closely if pathological changes in the vascular system feeding or within the organ exist. We also measured an effect of contrast agent on the FF. Although this parameter is derived from the material decomposition algorithm of the syngo.via® application profile “Liver VNC”, which extracts the iodine and creates virtual non-contrast images, we still noticed differences in the FF in response to contrast agent administration. We assume that these differences are an error due to contrast agent application, with the software not fully extracting the iodine attenuation from the virtual non-contrast images.
Our study also revealed that the variable ρe is not an ideal parameter to work with alone for tissue characterisation. Application of the formula described in Saito et al. (21) results in very small parameter values. Thus, differences can to some extent only be recognized to the third decimal place. Furthermore, because ρe is directly derieved from HU, it contributes no further information about the organic structure of tissue.
Former studies have already shown that DECT facilitates the differentiation of different materials (7, 8, 10, 15-17, 22, 23). However, these earlier studies referred predominantly to the differentiation of materials with different compositions, for example kidney stones from uric acid or calcium (17, 24, 25). The differentiation of soft tissues is an ambitious idea and may not yet be possible (10, 22) because soft tissues like liver, kidney, spleen, muscle, and fat tissue consist primarily of hydrogen, carbon, nitrogen and oxygen, which all have similar low average atomic numbers and consequently show no major differences in attenuation of X-rays of different energy levels (16, 22). However, even when materials show similar attenuation, they can still show different electron densities and elemental compositions, which may lead to their separation based on material qualities (7, 8). Therefore, we tested whether soft tissues can be separated from each other using the aforementioned DECT parameters or the averaged Hounsfield units obtained from mixing of the two different energy spectra. By measuring values in a patient cohort collected over the course of a year, we were able to estimate differences and confidence intervals, within which the true differences between organs lie. After comparing all collected variable values for the spleen, liver, and muscle, we found only low differences on non-contrast images. Nevertheless, the contrast-enhanced DECT scans showed greater differences, with the variable values extracted from muscle showing great differences from the values of liver and spleen because of the low contrast uptake of muscle. Thus, automated organ detection would be possible by applying contrast agent. Nevertheless, when comparing the liver to spleen, we only detected small differences in FF and Zeff in non-contrast images, and also in HUmixed, FF, and Zeff in contrast-enhanced images. The estimated differences in HUmixed between liver and spleen were higher in non-contrast images than in contrast-enhanced images. Nevertheless, we are convinced that there is potential to increase the differences between these tissues in contrast-enhanced images, by developing the CT protocol and the contrast agent application. For example, the spleen and liver should show greater differences in images acquired during arterial phase CT.
A recent study by Hunemohr et al. also measured HU, ρe, and Zeff on DECT scans (16). They performed scans at 120 kV and measured the parameters three times in kidney, liver, muscle, and three different fat tissues. The HU values ranged from 43 to 63 HU for liver, 41 to 46 HU for kidney, 40 to 44 HU for muscle, and ‒98 to ‒55 HU for fat tissue. The values for ρe were 1.05-1.07 for liver, 1.05 for kidney, 1.05 for muscle, and 0.93-0.97 for fat tissue. Zeff ranged from 7.26 to 7.30 for liver, 7.27 to 7.36 for kidney, 7.21 to 7.31 for muscle, and 6.10 to 6.48 for fat tissue. Their values for ρe and Zeff are similar to those in our study because the tube voltage of 120 kV corresponds with our mixed 100 kV and 140 Sn kV images. Nevertheless, the study sample of Hunemohr et al. was smaller than ours and did not aim to generate normal distributions. By taking measurements from a large sample, as in our study, it can be assumed that patients with values around the median are likely to correspond to healthy adults. Additionally, we calculated the first and third quartiles, as pathological organic changes could possibly be present in patients with values below or above these limits.
Our study has several limitations. First even if we fixed tube voltages, the scanner was allowed for tube current modulation since we used clinical protocols for our retrospective analyses. The aim of the algorithm is to maintain a global image noise level, but we could not take into consideration dose variations and its effect on the different parameters evaluated. Second our standardized scan protocol included a split bolus administration of contrast medium, e.g., two separate following contrast peaks that might have influenced organ attenuation due to overlap of early and late contrast phases. Third our normal distribution dataset was generated using vendor specific software tools that may be difficult to transfer to other technical equipment.
To our knowledge our work is the first attempt to generate normal values for parenchymatous abdominal organs using a standardized dual energy CT imaging protocol. Our data was aimed at multiparametric evaluation that enables multidimensional analysis. Such data may be crucial for the development of deep learning powered algorithms that allow for automatic tissue characterization. In future studies, we aim to examine whether measurements of the variables HUmixed, Zeff, and FF are suitable for diagnosing disease. Studies involving the collaboration of different institutes of pathology should be pursued. As a histological examination is the gold standard for diagnosis in many diseases, we should perform correlation analyses between the histological results and radiological examinations. It is conceivable that the microscopic changes a disease can cause may also affect the measured parameters. Diseases leading to increased organ density, such as fibrosis or cirrhosis of the liver, should consequently lead to higher CT numbers and ρe, whereas diseases with lower organ density such as liver steatosis, should lead to lower CT numbers and ρe, but a higher FF.
Conclusion
Organ characterization using multiparametric dual energy CT analysis is possible. The measurable CT parameters showed differences before and after contrast agent application and may, therefore, be useful for detecting pathological vascular changes. The increased number of parameters obtained from DECT in comparison to single-energy CT led to improved organ characterization and shows promise for future automated segmentation techniques. To our knowledge this is the first attempt to provide normal values for characterization of parenchymatous organs.
Footnotes
This article is freely accessible online.
↵* These Authors contributed equally to this study.
Authors’ Contributions
Diana Kreul (first author): Study design, literature research, data analysis, editing and writing of the article. Tilo Niemann (corresponding author): Literature research, editing, writing and proofreading of the article. Rahel Kubik: Data analysis and proofreading of the article. John Froehlich and Michael Thali: Co-editing, proofreading of the article. All Authors made pertinent contributions to the article, and proofread and approved the final article before submission.
Conflicts of Interest
The Authors declare no conflicts of interest.
- Received June 1, 2021.
- Revision received June 19, 2021.
- Accepted July 27, 2021.
- Copyright © 2021 International Institute of Anticancer Research (Dr. George J. Delinasios), All rights reserved