Home Quantitative Measurement of Iodine Concentration in the Liver Using Abdominal C-Arm Computed Tomography
Post
Cancel

Quantitative Measurement of Iodine Concentration in the Liver Using Abdominal C-Arm Computed Tomography

Rationale and Objectives

C-arm computed tomography is an option on a C-arm angiographic system capable of acquiring projections while rotating the C-arm around the patient and reconstructing cross-sectional images with improved contrast resolution of 5 to 10 Hounsfield units. Typical abdominal C-arm computed tomographic (CCT) images, however, exhibit artifacts with spatially varying and drifting pixel values. Considering liver tumor oncologic procedures, the aim of this study was to evaluate the accuracy of liver iodine concentration (IC) estimated from CCT images under such challenging conditions.

Materials and Methods

The proposed method estimates the IC in a region of interest (ROI) using pixel values of CCT images measured at the ROI and a nonenhanced background. Two approaches to measure the background value were tested: one approach, L-BG, measured a corresponding local background value near each ROI, and the other, G-BG, used one global background value for the entire object. The accuracy of estimations using CCT and computed tomographic scanners was evaluated; an elliptical cylinder water phantom with iodine solution inserts and seven patient data sets with transcatheter arterial chemoembolization were used.

Results

With the least “truncation” (the edge of the object being outside the field of view) of 27 mm, the IC was accurately estimated with CCT images ( n = 9; root-mean-square error [RMSE], 1.60–1.63 mg/mL; normalized RMSE, 11.8%; r 2 = 0.97; P < .001), with the true concentration ranging from 2.32 to 31.82 mg/mL. With truncations of up to 100 mm ( n = 88), the estimation by L-BG remained accurate independent of the amount of truncation (RMSE, 1.58 mg/mL; normalized RMSE, 12.5%; r 2 = 0.06; P = .02), whereas the estimation by G-BG reduced the accuracy (RMSE, 4.61 mg/mL; normalized RMSE, 34.3%; r 2 = 0.10; P = .003). Clinical data ( n = 37) showed that the estimation from CCT images using the L-BG method agreed well with that from computed tomographic images (RMSE, 2.86 mg/mL; normalized RMSE, 38.7%; r 2 = 0.76; P < .001).

Conclusion

The liver IC can be accurately estimated with abdominal CCT images.

Many image-guided interventional procedures remain subjective, even with the introduction of digital flat-panel detectors. Many interventional procedures that use x-ray images involve iodine, for example, nonionized iodine contrast agent material ( ), iodine-mixed lipiodol oil for transcatheter arterial chemoembolization (TACE) ( ), and embolic particles in an iodized solution ( ). Measuring the concentration of iodine, which is likely associated with drug concentration or blood perfusion, can lead to objective, evidence-driven procedures.

In organs, iodine is distributed three dimensionally ( ). Thus, two-dimensional projection images such as fluoroscopic or angiographic images, which show the total amount of iodine along each x-ray path, are not sufficient for accurately estimating iodine concentration (IC)—the amount of iodine per unit volume—and its spatial distribution ( ). Therefore, in this study, we used a recently developed C-arm computed tomographic (CCT) system ( ) to quantify the IC. CCT imaging is an option on C-arm angiographic systems that acquires projections while rotating the C-arm around the patient and reconstructs cross-sectional images with improved contrast resolution of 5 to 10 Hounsfield units ( ). There are two major challenges with abdominal CCT imaging: (1) to quantify the IC from a single scan and (2) to overcome the effect of artifacts unique to abdominal applications. Thus, we proposed combining two methods, each of which solves one problem. We focused on measuring the IC in the liver, because our primary clinical application was TACE for liver tumors. In principal, however, the method can be applied to other applications and organs.

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Materials and methods

Get Radiology Tree app to read full this article<

CCT System

Get Radiology Tree app to read full this article<

CT Scanner

Get Radiology Tree app to read full this article<

Phantom

Get Radiology Tree app to read full this article<

Figure 1, C-arm computed tomographic (CCT) images and computed tomographic images of the elliptical phantom (a–c) and a patient (d, e) . CCT images of the elliptical phantom obtained by the C-arm system with two different amounts of truncation (the distance from the edge of phantom to the field of view): (a) 27 mm at the least amount of truncation and (b) 125 mm. Increased pixel values near the edge of the field of view (arrows) and decreased values near the isocenter are artifacts. (b) Methods for measuring the background values. For the L-BG method, the corresponding local background value was measured within 2 cm from the region of interest (ROI; as indicated by the circles near S1 and S2). In contrast, for the G-BG method, the global background value was measured near the isocenter (at the square ). (c) A reconstructed CT image of the elliptical phantom. A CCT image (d) and a CT image (e) of a patient at the same slice show the lipiodol distribution. (d) The corresponding background value was measured at a nonenhanced tissue lesion (circles) as close to the ROIs (diamonds) as possible. (e) A background value was measured at the liver parenchyma with no enhancement (circle) .

Get Radiology Tree app to read full this article<

Patients

Get Radiology Tree app to read full this article<

Estimation of IC

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Truncation Correction

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Collection and Processing of the Phantom Data

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Collection and Processing of the Patient Data

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Statistical Analysis

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Results

Linearity

Get Radiology Tree app to read full this article<

Table 1

Summaries of the Estimated Iodine Concentrations with Various Iodine Densities (phantom data set 1), with Various Levels of Truncation (phantom data set 2), and with Clinical Settings (patient data)

Phantom Data Set 2 Phantom Data Set 1 ( n = 9) Truncation < 150 mm ( n = 119) Truncation < 100 mm ( n = 88) Patient Data ( n = 37) Modality and Method RMSE (mg/mL) ⁎ r 2 ∥ RMSE (mg/mL) ⁎ r 2 ∥ RMSE (mg/mL) ⁎ r 2 ∥ RMSE (mg/mL) ⁎ r 2 ‡ CCT imaging L-BG 1.63 (12.5%) 0.97 † 3.12 (23.2%) 0.25 † 1.58 (11.8%) 0.06 § 2.86 (38.7%) 0.76 † G-BG 1.60 (12.5%) 0.97 † 11.05 (82.1%) 0.09 ¶ 4.61 (34.3%) 0.10 ¶ CT imaging 0.51 (3.9%) 1.00 †

CCT, C-arm computed tomographic; CT, computed tomographic; RMSE, root-mean-square error.

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Truncation

Get Radiology Tree app to read full this article<

Figure 2, The result of phantom experiment data set 2: iodine concentration values estimated from C-arm computed tomographic (CCT) images by two methods, L-BG and G-BG, with various amounts of truncation. The dashed line indicates the true concentration (13.46 mg/mL). The discrepancies increase with truncation with the G-BG method. Less than 100 mm of truncation is a more realistic value for clinical conditions.

Get Radiology Tree app to read full this article<

Patient Data

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Figure 3, The result of patient data: iodine concentration values obtained from C-arm computed tomographic (CCT) images compared with those obtained from computed tomographic (CT) images.

Get Radiology Tree app to read full this article<

Discussion

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Acknowledgments

Get Radiology Tree app to read full this article<

Appendix A

Get Radiology Tree app to read full this article<

μx=wIρIμ(ρ)I+wBGρBGμ(ρ)BG, μ

x

=

w

I

ρ

I

μ

(

ρ

)

I

+

w

BG

ρ

BG

μ

(

ρ

)

B

G

,

where w I and w BG are the relative concentration in terms of a fraction by weight, ρ I and ρ BG are the native density (in grams per cubic centimeter), and μ (ρ) I and μ (ρ) BG are the mass attenuation coefficient (in square centimeters per gram) of iodine and background material, respectively. In a case in which there is only one substance (iodine) dissolved in the background material (water or tissue), the enhancement from the background due to the iodine is related to the linear attenuation coefficients of the ROI and background by

μx= −μBG=[wIρIμ(ρ)I+wBGρBGμ(ρ)BG]−[ρBGμ(ρ)BG]=wI×[ρIμ(ρ)I−ρBGμ(ρ)BG], μ

x

= −

μ

BG

=

[

w

I

ρ

I

μ

(

ρ

)

I

+

w

BG

ρ

BG

μ

(

ρ

)

B

G

]

[

ρ

BG

μ

(

ρ

)

B

G

]

=

w

I

×

[

ρ

I

μ

(

ρ

)

I

ρ

BG

μ

(

ρ

)

B

G

]

,

because w I = 1 − w BG . Thus, once μ x is known, w I can be estimated from a single measurement by

wI=(μx−μBG)/[ρIμ(ρ)I−ρBGμ(ρ)BG], w

I

=

(

μ

x

μ

BG

)

/

[

ρ

I

μ

(

ρ

)

I

ρ

BG

μ

(

ρ

)

B

G

]

,

where the values for the mass density of iodine, water, and liver tissue, available in the literature (20,21), are 4.93, 1.00, and 1.06 g/cm 3 , respectively. The reference x-ray photon energies of CCT and CT scanners in our study were found to be 79.0 and 78.6 keV, respectively. The mass attenuation coefficients of iodine, water, and liver tissue for CCT images were 3.713, 0.185, and 0.183 cm 2 /g, respectively. The corresponding values for CT images were 3.795, 0.185, and 0.184 cm 2 /g. To obtain the relative concentration w I , we use Equation A3 , whereby μ x and μ BG are obtained from the values of appropriate image pixels. (See Appendix C for CCT imaging and Appendix B for CT imaging.)

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

IC=1000×ρI×wI. IC

=

1000

×

ρ

I

×

w

I

.

Get Radiology Tree app to read full this article<

Appendix B

Get Radiology Tree app to read full this article<

μx=μw×(HUX/1000+1);μBG=μw×(HUBG/1000+1), μ

x

=

μ

w

×

(

HU

X

/

1000

+

1

)

;

μ

BG

=

μ

w

×

(

HU

BG

/

1000

+

1

)

,

where μ w = ρ w × μ (ρ)w = 0.185 cm −1 in our study. Then, we use Equation A3 to calculate w I (see Appendix A ).

Get Radiology Tree app to read full this article<

Appendix C

Get Radiology Tree app to read full this article<

fx=ax×μx+bx;fBG=aBG×μBG+bBG, f

x

=

a

x

×

μ

x

+

b

x

;

f

BG

=

a

BG

×

μ

BG

+

b

BG

,

where a__x , b__x and a BG , b BG are scaling and offset parameters at ROI and background, respectively. We assume that the additive offsets, b__x and b BG , are much smaller than the first terms of Equation C1 and that the scalings, a__x and a BG , are similar to each other. The ratio of μ x and μ BG can then be approximated reasonably well by

μx/μBG≈fx/fBG. μ

x

/

μ

BG

f

x

/

f

BG

.

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

wI=μBG×(μx/μBG−1)/[ρIμ(ρ)I−ρBGμ(ρ)BG], w

I

=

μ

BG

×

(

μ

x

/

μ

BG

1

)

/

[

ρ

I

μ

(

ρ

)

I

ρ

BG

μ

(

ρ

)

B

G

]

,

where μ BG is the linear attenuation coefficient of the background lesion (water or liver parenchyma), which is obtained by ρ BG μ (ρ) BG , with values available in the literature.

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

References

  • 1. Gailloud P., Oishi S., Carpenter J., Murphy K.J.: Three-dimensional digital angiography: new tool for simultaneous three-dimensional rendering of vascular and osseous information during rotational angiography. AJNR Am J Neuroradiol 2004; 25: pp. 571-573.

  • 2. Lardo A.C., Cordeiro M.A.S., Silva C., et. al.: Contrast-enhanced multidetector computed tomography viability imaging after myocardial infarction: characterization of myocyte death, microvascular obstruction, and chronic scar. Circulation 2006; 113: pp. 394-404.

  • 3. Geschwind J.F.: Chemoembolization for hepatocellular carcinoma: where does the truth lie?. J Vasc Interv Radiol 2002; 13: pp. 991-994.

  • 4. Ramsey D.E., Kernagis L.Y., Soulen M.C., Geschwind J.-F.H.: Chemoembolization of hepatocellular carcinoma. J Vasc Interv Radiol 2002; 13: pp. S211-S221.

  • 5. Hong K., Khwaja A., Liapi E., Torbenson M.S., Georgiades C.S., Geschwind J.-F.H.: New intra-arterial drug delivery system for the treatment of liver cancer: preclinical assessment in a rabbit model of liver cancer. Clin Cancer Res 2006; 12: pp. 2563-2567.

  • 6. Oto A., Tamm E.P., Szklaruk J.: Multidetector row CT of the liver. Radiol Clin North Am 2005; 43: pp. 827-848.

  • 7. Takayasu K., Arii S., Matsuo N., et. al.: Comparison of CT findings with resected specimens after chemoembolization with iodized oil for hepatocellular carcinoma. AJR Am J Roentgenol 2000; 175: pp. 699-704.

  • 8. Sze D.Y., Strobel N., Fahrig R., Moore T., Busque S., Frisoli J.K.: Transjugular intrahepatic portosystemic shunt creation in a polycystic liver facilitated by hybrid cross-sectional/angiographic imaging. J Vasc Interv Radiol 2006; 17: pp. 711-715.

  • 9. Fahrig R., Dixon R., Payne T., Morin R.L., Ganguly A., Strobel N.: Dose and image quality for a cone-beam C-arm CT system. Med Phys 2006; 33: pp. 4541-4550.

  • 10. Hirota S., Nakao N., Yamamoto S., et. al.: Cone-beam CT with flat-panel-detector digital angiography system: early experience in abdominal interventional procedures. Cardiovasc Interv Radiol 2006; 29: pp. 1034-1038.

  • 11. Wallace M.J., Murthy R., Kamat P.P., et. al.: Impact of C-arm CT on hepatic arterial interventions for hepatic malignancies. J Vasc Interv Radiol 2007; 18: pp. 1500-1507.

  • 12. Alvarez R.E., Macovski A.: Energy-selective reconstructions in x-ray computerised tomography. Phys Med Biol 1976; 21: pp. 733-744.

  • 13. Kelcz F., Joseph P.M., Hilal S.K.: Noise considerations in dual energy CT scanning. Med Phys 1979; 6: pp. 418-425.

  • 14. Gailloud P., Oishi S., Murphy K.: Three-dimensional fusion digital subtraction angiography: new reconstruction algorithm for simultaneous three-dimensional rendering of osseous and vascular information obtained during rotational angiography. AJNR Am J Neuroradiol 2005; 26: pp. 908-911.

  • 15. van de Kraats E.B., Carelsen B., Fokkens W.J., et. al.: Direct navigation on 3D rotational x-ray data acquired with a mobile propeller C-arm: accuracy and application in functional endoscopic sinus surgery. Phys Med Biol 2005; 50: pp. 5769-5781.

  • 16. Noo F., Clackdoyle R., Pack J.D.: A two-step Hilbert transform method for 2D image reconstruction. Phys Med Biol 2004; 49: pp. 3903-3923.

  • 17. Feldkamp L.A., Davis L.C., Kress J.W.: Practical cone-beam algorithm. J Opt Soc Am A 1984; 1: pp. 612-619.

  • 18. Parker D.L.: Optimal short scan convolution reconstruction for fanbeam CT. Med Phys 1982; 9: pp. 254-257.

  • 19. Silver M.D.: A method for including redundant data in computed tomography. Med Phys 2000; 27: pp. 773-774.

  • 20. Tissue substitutes in radiation dosimetry and measurement (ICRU Report 44).1989.International Commission on Radiation Units and MeasurementsBethesda, MD

  • 21. Hubbell J.H., Seltzer S.M.: NIST standard reference database 126.1996.National Institute of Standards and TechnologyGaithersburg, MD

  • 22. Barclay A.B., Eisner R.L., DiBella E.V.R.: PET thorax model database. http://www.emory.edu/CRL/abb/thoraxmodel/Contents.html Accessed September 5 2008

  • 23. Barclay AB, Eisner RL, DiBella EVR. Construction of a thorax model database from clinical PET attenuation scans. Presented at: Annual Meeting of the Society of Nuclear Medicine, Minneapolis, MN, June 12–15, 1995.

  • 24. Schafer D., Borgert J., Rasche V., Grass M.: Motion-compensated and gated cone beam filtered back-projection for 3-D rotational x-ray angiography. IEEE Trans Med Imaging 2006; 25: pp. 898-906.

  • 25. Taguchi K, Kudo H. Motion compensated fan-beam reconstruction for computed tomography using derivative backprojection filtering approach. Presented at: 9th International Conference on Fully Three-Dimensional Reconstruction in Radiology and Nuclear Medicine; Lindau, Germany; July 9–13, 2007.

  • 26. Taguchi K., Kudo H.: A simple motion tracking backprojection for a class of affine transformation. Proc SPIE 2008; 6913: pp. 66.

  • 27. Taguchi K., Kudo H.: Motion compensated fan-beam reconstruction for nonrigid transformation. IEEE Trans Med Imaging 2008; 27: pp. 907-917.

This post is licensed under CC BY 4.0 by the author.