Updated breast CT dose coefficients (DgN CT ) using patient-derived breast shapes and heterogeneous fibroglandular distributions

Andrew M. Hernandez, Amy E. Becker, John M Boone

Research output: Contribution to journalArticle

1 Citation (Scopus)

Abstract

Purpose: The purpose of this work was to generate uncompressed heterogeneous breast phantom models using size-dependent fibroglandular distributions derived from a large cohort of breast CT (bCT) datasets, and to compare differences in normalized glandular dose coefficients for bCT “DgN CT ” when the realistic heterogeneous model is considered relative to the simple, homogeneous model used in the past. Methods: A cohort of 274 segmented bCT datasets were used to quantify the fibroglandular tissue distribution within the breast parenchyma. Each dataset was interpolated to an isotropic voxel size of 0.25 mm and the breast center-of-mass was aligned for all coronal slices. Each aligned dataset was converted into two binarized volumes representing voxels containing only glandular tissue “G(x,y,z)” and voxels containing glandular or adipose tissue “AG(x,y,z)”. The datasets were classified by volume in accordance with previously reported size-dependent, breast-shaped phantoms. Within the five groups — each containing on average 55 datasets, all G(x,y,z) and AG(x,y,z) volumes were summed separately representing the cumulative distribution of glandular tissue or breast parenchyma (glandular + adipose), respectively. G(x,y,z) was divided by AG(x,y,z) on a voxel-by-voxel basis resulting in a glandular fraction “GF(x,y,z)” distribution for each phantom size. The GF(x,y,z) distributions were used to construct heterogeneous mathematical phantoms for the small, median, and large breast sizes with a 1.5 mm skin thickness — based on previously reported measurements from bCT, and a 5 mm skin thickness for comparison with outdated assumptions of skin thickness. A subset of 15 bCT datasets from the cohort (five for each breast size) were used to construct voxelized patient models for validation of the heterogeneous phantom models. Monte Carlo techniques were used to estimate monoenergetic DgN(E) C T values for photon energies from 9 to 70 keV (in 1 keV intervals) using the mathematical phantoms composed of either heterogeneous or homogeneous breast parenchyma. Polyenergetic (pDgN CT ) coefficients were determined by weighting the DgN(E) C T values by x-ray spectra tuned to the beam characteristic of breast CT. Dose coefficients were compared between the two breast compositions for each volume class, breast density, and skin thickness. Results: For photon energies ≲45 keV, the homogeneous model overestimates DgN(E) values relative to the realistic heterogeneous model sorted into five volume classes. The 5 mm skin thickness underestimates DgN(E) values relative to the realistic 1.5 mm thickness for lower energies and the differences diminish up to 70 keV. Averaged across all phantom sizes the homogeneous model overestimates pDgN CT by 5.7% and 23.3% for the 60 kV W/Cu and 49 kV W/Al spectra, respectively. The heterogeneous model was also found to be in agreement with the voxelized bCT patient models with pDgN CT differences less than 2.3% and 5.2% for the 60 kV W/Cu and 49 kV W/Al spectra, respectively, across all phantom sizes. Conclusion: Anatomically accurate heterogeneous phantom models were developed using bCT image-derived fibroglandular tissue distributions. These new models improve the accuracy of bCT dosimetry, and in conjunction with previous models for mammography, may help in providing a more universally accepted breast dosimetry model.

Original languageEnglish (US)
JournalMedical Physics
DOIs
StatePublished - Jan 1 2019

Fingerprint

Breast
Tissue Distribution
Skin
Photons
Mammography
Datasets
Adipose Tissue

Keywords

  • breast CT
  • breast dosimetry
  • glandular distribution
  • MCNP
  • skin thickness

ASJC Scopus subject areas

  • Biophysics
  • Radiology Nuclear Medicine and imaging

Cite this

Updated breast CT dose coefficients (DgN CT ) using patient-derived breast shapes and heterogeneous fibroglandular distributions . / Hernandez, Andrew M.; Becker, Amy E.; Boone, John M.

In: Medical Physics, 01.01.2019.

Research output: Contribution to journalArticle

@article{d22786fe2e774a9bbad0c98b708de2aa,
title = "Updated breast CT dose coefficients (DgN CT ) using patient-derived breast shapes and heterogeneous fibroglandular distributions",
abstract = "Purpose: The purpose of this work was to generate uncompressed heterogeneous breast phantom models using size-dependent fibroglandular distributions derived from a large cohort of breast CT (bCT) datasets, and to compare differences in normalized glandular dose coefficients for bCT “DgN CT ” when the realistic heterogeneous model is considered relative to the simple, homogeneous model used in the past. Methods: A cohort of 274 segmented bCT datasets were used to quantify the fibroglandular tissue distribution within the breast parenchyma. Each dataset was interpolated to an isotropic voxel size of 0.25 mm and the breast center-of-mass was aligned for all coronal slices. Each aligned dataset was converted into two binarized volumes representing voxels containing only glandular tissue “G(x,y,z)” and voxels containing glandular or adipose tissue “AG(x,y,z)”. The datasets were classified by volume in accordance with previously reported size-dependent, breast-shaped phantoms. Within the five groups — each containing on average 55 datasets, all G(x,y,z) and AG(x,y,z) volumes were summed separately representing the cumulative distribution of glandular tissue or breast parenchyma (glandular + adipose), respectively. G(x,y,z) was divided by AG(x,y,z) on a voxel-by-voxel basis resulting in a glandular fraction “GF(x,y,z)” distribution for each phantom size. The GF(x,y,z) distributions were used to construct heterogeneous mathematical phantoms for the small, median, and large breast sizes with a 1.5 mm skin thickness — based on previously reported measurements from bCT, and a 5 mm skin thickness for comparison with outdated assumptions of skin thickness. A subset of 15 bCT datasets from the cohort (five for each breast size) were used to construct voxelized patient models for validation of the heterogeneous phantom models. Monte Carlo techniques were used to estimate monoenergetic DgN(E) C T values for photon energies from 9 to 70 keV (in 1 keV intervals) using the mathematical phantoms composed of either heterogeneous or homogeneous breast parenchyma. Polyenergetic (pDgN CT ) coefficients were determined by weighting the DgN(E) C T values by x-ray spectra tuned to the beam characteristic of breast CT. Dose coefficients were compared between the two breast compositions for each volume class, breast density, and skin thickness. Results: For photon energies ≲45 keV, the homogeneous model overestimates DgN(E) values relative to the realistic heterogeneous model sorted into five volume classes. The 5 mm skin thickness underestimates DgN(E) values relative to the realistic 1.5 mm thickness for lower energies and the differences diminish up to 70 keV. Averaged across all phantom sizes the homogeneous model overestimates pDgN CT by 5.7{\%} and 23.3{\%} for the 60 kV W/Cu and 49 kV W/Al spectra, respectively. The heterogeneous model was also found to be in agreement with the voxelized bCT patient models with pDgN CT differences less than 2.3{\%} and 5.2{\%} for the 60 kV W/Cu and 49 kV W/Al spectra, respectively, across all phantom sizes. Conclusion: Anatomically accurate heterogeneous phantom models were developed using bCT image-derived fibroglandular tissue distributions. These new models improve the accuracy of bCT dosimetry, and in conjunction with previous models for mammography, may help in providing a more universally accepted breast dosimetry model.",
keywords = "breast CT, breast dosimetry, glandular distribution, MCNP, skin thickness",
author = "Hernandez, {Andrew M.} and Becker, {Amy E.} and Boone, {John M}",
year = "2019",
month = "1",
day = "1",
doi = "10.1002/mp.13391",
language = "English (US)",
journal = "Medical Physics",
issn = "0094-2405",
publisher = "AAPM - American Association of Physicists in Medicine",

}

TY - JOUR

T1 - Updated breast CT dose coefficients (DgN CT ) using patient-derived breast shapes and heterogeneous fibroglandular distributions

AU - Hernandez, Andrew M.

AU - Becker, Amy E.

AU - Boone, John M

PY - 2019/1/1

Y1 - 2019/1/1

N2 - Purpose: The purpose of this work was to generate uncompressed heterogeneous breast phantom models using size-dependent fibroglandular distributions derived from a large cohort of breast CT (bCT) datasets, and to compare differences in normalized glandular dose coefficients for bCT “DgN CT ” when the realistic heterogeneous model is considered relative to the simple, homogeneous model used in the past. Methods: A cohort of 274 segmented bCT datasets were used to quantify the fibroglandular tissue distribution within the breast parenchyma. Each dataset was interpolated to an isotropic voxel size of 0.25 mm and the breast center-of-mass was aligned for all coronal slices. Each aligned dataset was converted into two binarized volumes representing voxels containing only glandular tissue “G(x,y,z)” and voxels containing glandular or adipose tissue “AG(x,y,z)”. The datasets were classified by volume in accordance with previously reported size-dependent, breast-shaped phantoms. Within the five groups — each containing on average 55 datasets, all G(x,y,z) and AG(x,y,z) volumes were summed separately representing the cumulative distribution of glandular tissue or breast parenchyma (glandular + adipose), respectively. G(x,y,z) was divided by AG(x,y,z) on a voxel-by-voxel basis resulting in a glandular fraction “GF(x,y,z)” distribution for each phantom size. The GF(x,y,z) distributions were used to construct heterogeneous mathematical phantoms for the small, median, and large breast sizes with a 1.5 mm skin thickness — based on previously reported measurements from bCT, and a 5 mm skin thickness for comparison with outdated assumptions of skin thickness. A subset of 15 bCT datasets from the cohort (five for each breast size) were used to construct voxelized patient models for validation of the heterogeneous phantom models. Monte Carlo techniques were used to estimate monoenergetic DgN(E) C T values for photon energies from 9 to 70 keV (in 1 keV intervals) using the mathematical phantoms composed of either heterogeneous or homogeneous breast parenchyma. Polyenergetic (pDgN CT ) coefficients were determined by weighting the DgN(E) C T values by x-ray spectra tuned to the beam characteristic of breast CT. Dose coefficients were compared between the two breast compositions for each volume class, breast density, and skin thickness. Results: For photon energies ≲45 keV, the homogeneous model overestimates DgN(E) values relative to the realistic heterogeneous model sorted into five volume classes. The 5 mm skin thickness underestimates DgN(E) values relative to the realistic 1.5 mm thickness for lower energies and the differences diminish up to 70 keV. Averaged across all phantom sizes the homogeneous model overestimates pDgN CT by 5.7% and 23.3% for the 60 kV W/Cu and 49 kV W/Al spectra, respectively. The heterogeneous model was also found to be in agreement with the voxelized bCT patient models with pDgN CT differences less than 2.3% and 5.2% for the 60 kV W/Cu and 49 kV W/Al spectra, respectively, across all phantom sizes. Conclusion: Anatomically accurate heterogeneous phantom models were developed using bCT image-derived fibroglandular tissue distributions. These new models improve the accuracy of bCT dosimetry, and in conjunction with previous models for mammography, may help in providing a more universally accepted breast dosimetry model.

AB - Purpose: The purpose of this work was to generate uncompressed heterogeneous breast phantom models using size-dependent fibroglandular distributions derived from a large cohort of breast CT (bCT) datasets, and to compare differences in normalized glandular dose coefficients for bCT “DgN CT ” when the realistic heterogeneous model is considered relative to the simple, homogeneous model used in the past. Methods: A cohort of 274 segmented bCT datasets were used to quantify the fibroglandular tissue distribution within the breast parenchyma. Each dataset was interpolated to an isotropic voxel size of 0.25 mm and the breast center-of-mass was aligned for all coronal slices. Each aligned dataset was converted into two binarized volumes representing voxels containing only glandular tissue “G(x,y,z)” and voxels containing glandular or adipose tissue “AG(x,y,z)”. The datasets were classified by volume in accordance with previously reported size-dependent, breast-shaped phantoms. Within the five groups — each containing on average 55 datasets, all G(x,y,z) and AG(x,y,z) volumes were summed separately representing the cumulative distribution of glandular tissue or breast parenchyma (glandular + adipose), respectively. G(x,y,z) was divided by AG(x,y,z) on a voxel-by-voxel basis resulting in a glandular fraction “GF(x,y,z)” distribution for each phantom size. The GF(x,y,z) distributions were used to construct heterogeneous mathematical phantoms for the small, median, and large breast sizes with a 1.5 mm skin thickness — based on previously reported measurements from bCT, and a 5 mm skin thickness for comparison with outdated assumptions of skin thickness. A subset of 15 bCT datasets from the cohort (five for each breast size) were used to construct voxelized patient models for validation of the heterogeneous phantom models. Monte Carlo techniques were used to estimate monoenergetic DgN(E) C T values for photon energies from 9 to 70 keV (in 1 keV intervals) using the mathematical phantoms composed of either heterogeneous or homogeneous breast parenchyma. Polyenergetic (pDgN CT ) coefficients were determined by weighting the DgN(E) C T values by x-ray spectra tuned to the beam characteristic of breast CT. Dose coefficients were compared between the two breast compositions for each volume class, breast density, and skin thickness. Results: For photon energies ≲45 keV, the homogeneous model overestimates DgN(E) values relative to the realistic heterogeneous model sorted into five volume classes. The 5 mm skin thickness underestimates DgN(E) values relative to the realistic 1.5 mm thickness for lower energies and the differences diminish up to 70 keV. Averaged across all phantom sizes the homogeneous model overestimates pDgN CT by 5.7% and 23.3% for the 60 kV W/Cu and 49 kV W/Al spectra, respectively. The heterogeneous model was also found to be in agreement with the voxelized bCT patient models with pDgN CT differences less than 2.3% and 5.2% for the 60 kV W/Cu and 49 kV W/Al spectra, respectively, across all phantom sizes. Conclusion: Anatomically accurate heterogeneous phantom models were developed using bCT image-derived fibroglandular tissue distributions. These new models improve the accuracy of bCT dosimetry, and in conjunction with previous models for mammography, may help in providing a more universally accepted breast dosimetry model.

KW - breast CT

KW - breast dosimetry

KW - glandular distribution

KW - MCNP

KW - skin thickness

UR - http://www.scopus.com/inward/record.url?scp=85061277349&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85061277349&partnerID=8YFLogxK

U2 - 10.1002/mp.13391

DO - 10.1002/mp.13391

M3 - Article

C2 - 30661250

AN - SCOPUS:85061277349

JO - Medical Physics

JF - Medical Physics

SN - 0094-2405

ER -