RESUMEN
Objective.Contrast-enhanced computed tomography (CECT) is commonly used in the pre-treatment evaluation of liver Y-90 radioembolization feasibility. CECT provides detailed imaging of the liver and surrounding structures, allowing healthcare providers to assess the size, location, and characteristics of liver tumors prior to the treatment. Here we propose a method for translating CECT images to an expected dose distribution for tumor(s) and normal liver tissue.Approach.A pre-procedure CECT is used to obtain an iodine arterial-phase distribution by subtracting the non-contrast CT from the late arterial phase. The liver segments surrounding the targeted tumor are selected using Couinaud's method. The resolution of the resulting images is then degraded to match the resolution of the positron emission tomography (PET) images, which can image the Y-90 activity distribution post-treatment. The resulting images are then used in the same way as PET images to compute doses using the local deposition method. CECT images from three patients were used to test this method retrospectively and were compared with Y-90 PET-based dose distributions through dose volume histograms.Main results.Results show a concordance between predicted and delivered Y-90 dose distributions with less than 10% difference in terms of mean dose, for doses greater than 10% of the 98th percentile (D2%).Significance.CECT-derived predictions of Y-90 radioembolization dose distributions seem promising as a supplementary tool for physicians when assessing treatment feasibility. This dosimetry prediction method could provide a more comprehensive pre-treatment evaluation-offering greater insights than a basic assessment of tumor opacification on CT images.
Asunto(s)
Embolización Terapéutica , Neoplasias Hepáticas , Tomografía Computarizada por Rayos X , Humanos , Neoplasias Hepáticas/radioterapia , Neoplasias Hepáticas/diagnóstico por imagen , Radioisótopos de Itrio/uso terapéutico , Dosis de Radiación , Medios de Contraste , Hígado/diagnóstico por imagen , Hígado/efectos de la radiación , Dosificación Radioterapéutica , Tomografía de Emisión de Positrones , Procesamiento de Imagen Asistido por Computador/métodosRESUMEN
BACKGROUND AND PURPOSE: Magnetic resonance imaging (MRI) scanners are a major contributor to greenhouse gas emissions from the healthcare sector, and efforts to improve energy efficiency and reduce energy consumption rely on quantification of the characteristics of energy consumption. The purpose of this work was to develop a semi-automatic analytical methodology for the characterization of the energy consumption of MRI systems using only the load duration curve (LDC). LDCs are a fundamental tool used across various fields to analyze and understand the behavior of loads over time. METHODS: An electric current transformer sensor and data logger were installed on two 3T MRI scanners from two vendors, termed M1 (outpatient scanner) and M2 (inpatient/emergency scanner). Data was collected for 1 month (7/11/2023 to 8/11/2023). Active power was calculated, assuming a balanced three-phase system, using the average current measured across all three phases, a 480 V reference voltage for both machines, and vendor-provided power factors. An LDC was constructed for each system by sorting the active power values in descending order and computing the cumulative time (in units of percentage) for each data point. The first derivative of the LDC was then computed (LDC'), smoothed by convolution with a window function (sLDC'), and used to detect transitions between different system modes including (in descending power levels): scan, prepared-to-scan, idle, low-power, and off. The final, segmented LDC was used to measure time (% total time), total energy (kWh), and mean power (kW) for each system mode on both scanners. The method was validated by comparing mean power values, computed using the segmented 1-month LDC, for each nonproductive system mode (i.e., prepared-to-scan, idle, lower-power, and off) against power levels measured after a deliberate system shutdown was performed for each scanner (1 day worth of data). RESULTS: The validation revealed differences in mean power values <1.4% for all nonproductive modes and both scanners. In the scan system mode, the mean power values ranged from 29.8 to 37.2 kW and the total energy consumed for 1 month ranged from 11 106 to 14 466 kWh depending on the scanner. Over the course of 1 month, the portion of time the scanners were in nonproductive modes ranged from 76% to 80% across scanners and the nonproductive energy consumption ranged from 8010 to 6722 kWh depending on the scanner. The M1 (outpatient) scanner consumed 99.9 and 183.9 kWh/day in idle mode for weekdays and weekends, respectively, because the scanner spent 23% more time proportionally in idle mode on the weekends. CONCLUSIONS: A semi-automatic method for quantifying energy consumption characteristics of MRI scanners was introduced and validated. This method is relatively simple to implement as it requires only power data from the scanners and avoids the technical challenges associated with extracting and processing scanner log files. The methodology enables quantitative evaluation of the power, time, and energy characteristics of MRI scanners in scan and nonproductive system modes, providing baseline data and the capability of identifying potential opportunities for enhancing the energy efficiency of MRI scanners.
RESUMEN
BACKGROUND: In total knee arthroplasty, unrestricted kinematic alignment aims to restore pre-arthritic lower limb alignment and joint lines. Joint line orientations of the contralateral healthy proximal tibia might be used to evaluate accuracy of tibial component alignment post-operatively if asymmetry is minimal. Our objective was to evaluate left-to-right asymmetry of the proximal tibial epiphysis in posterior tibial slope and varus-valgus orientation as related to unrestricted kinematic alignment principles. METHODS: High resolution CT images (0.5 mm slice thickness) were acquired from bilateral lower limbs of 11 skeletally mature subjects with no skeletal abnormalities. Images were segmented to generate 3D tibia models. Asymmetry was quantified by differences in orientations required to shape-match the proximal epiphysis of the mirror 3D tibia model to the proximal epiphysis of the contralateral 3D tibia model. FINDINGS: Systematic and random differences (i.e. mean ± standard deviation) in tibial slope and varus-valgus orientation were - 0.8° ± 1.2° and - 0.2° ± 0.8°, respectively. Ninety five percent confidence intervals on the means included 0° indicating that systematic differences were minimal. INTERPRETATION: Since random differences due to asymmetry are substantial in relation to random surgical deviations from pre-arthritic joint lines previously reported, post-operative computer tomograms of the contralateral healthy tibia should not be used to directly assess accuracy of tibial component alignment on a group level without correcting for differences in tibial slope and varus-valgus orientation due to asymmetry.
Asunto(s)
Artroplastia de Reemplazo de Rodilla , Osteoartritis de la Rodilla , Humanos , Tibia/diagnóstico por imagen , Tibia/cirugía , Fenómenos Biomecánicos , Articulación de la Rodilla/diagnóstico por imagen , Articulación de la Rodilla/cirugía , Artroplastia de Reemplazo de Rodilla/métodos , Epífisis/diagnóstico por imagen , Epífisis/cirugía , Osteoartritis de la Rodilla/cirugíaRESUMEN
OBJECTIVE: To demonstrate the potential of low-dose ultra-high-resolution CT (UHRCT) images to generate high-quality radiographic images on extremity phantoms and to estimate the radiation dose required for this. MATERIALS AND METHODS: A hand and knee phantom containing real human bones was imaged on an UHRCT scanner at full-dose, half-dose, and quarter-dose levels using a high-resolution extremity protocol. The raw data was reconstructed using both filtered back projection (FBP) and an iterative reconstruction algorithm (AIDR3D). Using custom designed software, each CT volume data set was converted to attenuation coefficients, and then a synthesized radiograph (synDX) was generated by forward projecting the volume data sets from a point source onto a 2D synthetic detector. The signal-to-noise ratio (SNR) was measured in the synDXs across all dose levels and the root-mean-squared error (RMSE) was computed with the FD synDXs as the reference. RESULTS: The proposed workflow generates high-quality synDXs at any arbitrary angle. For FBP, the SNR largely tracked with the radiation dose levels for both the knee and hand phantoms. For the knee phantom, iterative reconstruction provided a 6.1% higher SNR when compared to FBP. The RMSE was overall higher for the lowest dose levels and monotonically decreased with increasing dose. No substantial differences were observed qualitatively in the visualization of skeletal detail of the phantoms. CONCLUSION: The fine detail provided by UHRCT acquisitions of extremities facilitates the ability to generate quality radiographs, potentially eliminating the need for additional scanning on a conventional digital radiography system.
RESUMEN
Currently, there are multiple breast dosimetry estimation methods for mammography and its variants in use throughout the world. This fact alone introduces uncertainty, since it is often impossible to distinguish which model is internally used by a specific imaging system. In addition, all current models are hampered by various limitations, in terms of overly simplified models of the breast and its composition, as well as simplistic models of the imaging system. Many of these simplifications were necessary, for the most part, due to the need to limit the computational cost of obtaining the required dose conversion coefficients decades ago, when these models were first implemented. With the advancements in computational power, and to address most of the known limitations of previous breast dosimetry methods, a new breast dosimetry method, based on new breast models, has been developed, implemented, and tested. This model, developed jointly by the American Association of Physicists in Medicine and the European Federation for Organizations of Medical Physics, is applicable to standard mammography, digital breast tomosynthesis, and their contrast-enhanced variants. In addition, it includes models of the breast in both the cranio-caudal and the medio-lateral oblique views. Special emphasis was placed on the breast and system models used being based on evidence, either by analysis of large sets of patient data or by performing measurements on imaging devices from a range of manufacturers. Due to the vast number of dose conversion coefficients resulting from the developed model, and the relative complexity of the calculations needed to apply it, a software program has been made available for download or online use, free of charge, to apply the developed breast dosimetry method. The program is available for download or it can be used directly online. A separate User's Guide is provided with the software.
Asunto(s)
Neoplasias de la Mama , Mama , Humanos , Femenino , Mama/diagnóstico por imagen , Mamografía/métodos , Radiometría/métodos , Método de Montecarlo , Neoplasias de la Mama/diagnóstico por imagenRESUMEN
BACKGROUND: Breast computed tomography (CT) is an emerging breast imaging modality, and ongoing developments aim to improve breast CT's ability to detect microcalcifications. To understand the effects of different parameters on microcalcification detectability, a virtual clinical trial study was conducted using hybrid images and convolutional neural network (CNN)-based model observers. Mathematically generated microcalcifications were embedded into breast CT data sets acquired at our institution, and parameters related to calcification size, calcification contrast, cluster diameter, cluster density, and image display method (i.e., single slices, slice averaging, and maximum-intensity projections) were evaluated for their influence on microcalcification detectability. PURPOSE: To investigate the individual effects and the interplay of parameters affecting microcalcification detectability in breast CT. METHODS: Spherical microcalcifications of varying diameters (0.04, 0.10, 0.15, 0.20, 0.25, 0.30, 0.35, 0.40 mm) and native intensities were computer simulated to portray the partial volume effects of the imaging system. Calcifications were mathematically embedded into 109 patient breast CT volume data sets as individual calcifications or as clusters of calcifications. Six numbers of calcifications (1, 3, 5, 7, 10, 15) distributed within six cluster diameters (1, 3, 5, 6, 8, 10 mm) were simulated to study the effect of cluster density. To study the role of image display method, 2D regions of interest (ROIs) and 3D volumes of interest (VOIs) were generated using single slice extraction, slice averaging, and maximum-intensity projection (MIP). 2D and 3D CNNs were trained on the ROIs and VOIs, and receiver operating characteristic (ROC) curve analysis was used to evaluate detection performance. The area under the ROC curve (AUC) was used as the primary performance metric. RESULTS: Detection performance decreased with increasing section thickness, and peak detection performance occurred using the native section thickness (0.2 mm) and MIP display. The MIP display method, despite using a single slice, yielded comparable performance to the native section thickness, which employed 50 slices. Reduction in slices did not sacrifice detection accuracy and provided significant computational advantages over multi-slice image volumes. Larger cluster diameters resulted in reduced overall detectability, while smaller cluster diameters led to increased detectability. Additionally, we observed that the presence of more calcifications within a cluster improved the overall detectability, while fewer calcifications decreased it. CONCLUSIONS: As breast CT is still a relatively new breast imaging modality, there is an ongoing need to identify optimal imaging protocols. This work demonstrated the utility of MIP presentation for displaying image volumes containing microcalcification clusters. It is likely that human observers may also benefit from viewing MIPs compared to individual slices. The results of this investigation begin to elucidate how model observers interact with microcalcification clusters in a 3D volume, and will be useful for future studies investigating a broader set of parameters related to breast CT.
Asunto(s)
Enfermedades de la Mama , Calcinosis , Humanos , Mamografía/métodos , Tomografía Computarizada por Rayos X/métodos , Calcinosis/diagnóstico por imagen , Redes Neurales de la ComputaciónRESUMEN
BACKGROUND: Mathematical model observers have been shown to reasonably predict human observer performance and are useful when human observer studies are infeasible. Recently, convolutional neural networks (CNNs) have also been used as substitutes for human observers, and studies have shown their utility as an optimal observer. In this study, a CNN model observer is compared to the pre-whitened matched filter (PWMF) model observer in detecting simulated mass lesions inserted into 253 acquired breast computed tomography (bCT) images from patients imaged at our institution. PURPOSE: To compare CNN and PWMF model observers for detecting signal-known-exactly (SKE) location-known-exactly (LKE) simulated lesions in bCT images with real anatomical backgrounds, and to use these model observers collectively to optimize parameters and understand trends in performance with breast CT. METHODS: Spherical lesions with different diameters (1, 3, 5, 9 mm) were mathematically inserted into reconstructed patient bCT image data sets to mimic 3D mass lesions in the breast. 2D images were generated by extracting the center slice along the axial dimension or by slice averaging across adjacent slices to model thicker sections (0.4, 1.2, 2.0, 6.0, 12.4, 20.4 mm). The role of breast density was retrospectively studied using the range of breast densities intrinsic to the patient bCT data sets. In addition, mass lesions were mathematically inserted into Gaussian images matched to the mean and noise power spectrum of the bCT images to better understand the performance of the CNN in the context of a known ideal observer (the PWMF). The simulated Gaussian and bCT images were divided into training and testing data sets. Each training data set consisted of 91 600 images, and each testing data set consisted of 96 000 images. A CNN and PWMF was trained on the Gaussian training images, and a different CNN and PWMF was trained on the bCT training images. The trained model observers were tested, and receiver operating characteristic (ROC) curve analysis was used to evaluate detection performance. The area under the ROC curve (AUC) was the primary performance metric used to compare the model observers. RESULTS: In the Gaussian background, the CNN performed essentially identically to the PWMF across lesion sizes and section thicknesses. In the bCT background, the CNN outperformed the PWMF across lesion size, breast density, and most section thicknesses. These findings suggest that there are higher-order features in bCT images that are harnessed by the CNN observer but are inaccessible to the PWMF. CONCLUSIONS: The CNN performed equivalently to the ideal observer in Gaussian textures. In bCT background, the CNN captures more diagnostic information than the PWMF and may be a more pertinent observer when conducting optimal performance studies in breast CT images.
Asunto(s)
Procesamiento de Imagen Asistido por Computador , Tomografía Computarizada por Rayos X , Humanos , Estudios Retrospectivos , Procesamiento de Imagen Asistido por Computador/métodos , Tomografía Computarizada por Rayos X/métodos , Redes Neurales de la Computación , Mama/diagnóstico por imagenRESUMEN
BACKGROUND: The use of iodine-based contrast agent for better delineation of tumors in breast CT (bCT) has been shown to be compelling, similar to the tumor enhancement in contrast-enhanced breast MRI. Contrast-enhanced bCT (CE-bCT) is a relatively new tool, and a structured evaluation of different imaging parameters at play has yet to be conducted. In this investigation, data sets of acquired bCT images from 253 patients imaged at our institution were used in concert with simulated mathematically inserted spherical contrast-enhanced lesions to study the role of contrast enhancement on detectability. PURPOSE: To quantitatively evaluate the improvement in lesion detectability due to contrast enhancement across lesion diameter, section thickness, view plane, and breast density using a pre-whitened matched filter (PWMF) model observer. METHODS: The relationship between iodine concentration and Hounsfield units (HU) was measured using spectral modeling. The lesion enhancement from clinical CE-bCT images in 22 patients was evaluated, and the average contrast enhancement (ΔHU) was determined. Mathematically generated spherical mass lesions of varying diameters (1, 3, 5, 9, 11, 15 mm) and contrast enhancement levels (0, 0.25, 0.50, 0.75, 1) were inserted at random locations in 253 actual patient bCT datasets. Images with varying thicknesses (0.4-19.8 mm) were generated by slice averaging, and the role of view plane (coronal and axial planes) was studied. A PWMF was used to generate receiver operating characteristic (ROC) curves across parameters of lesion diameter, contrast enhancement, section thickness, view plane, and breast density. The area under the ROC curve (AUC) was used as the primary performance metric, generated from over 90,000 simulated lesions. RESULTS: An average 20% improvement (ΔAUC = 0.1) in lesion detectability due to contrast enhancement was observed across lesion diameter, section thickness, breast density, and view plane. A larger improvement was observed when stratifying patients based on breast density. For patients with VGF ≤ 40%, detection performance improved up to 20% (until AUC â1), and for patients with denser breasts (VGF > 40%), detection performance improved more drastically, ranging from 20% to 80% for 1- and 5-mm lesions. For the 1 mm lesion, detection performance raised slightly at the 1.2 mm section thickness before falling off as thickness increased. For larger lesions, detection performance was generally unaffected as section thickness increased up until it reached 5.8 mm, where performance began to decline. Detection performance was higher in the axial plane compared to the coronal plane for smaller lesions and thicker sections. CONCLUSIONS: For emerging diagnostic tools like CE-bCT, it is important to optimize imaging protocols for lesion detection. In this study, we found that intravenous contrast can be used to detect small lesions in dense breasts. Optimal section thickness for detectability has dependencies on breast density and lesion size, therefore, display thickness should be adjusted in real-time using display software. These findings may be useful for the development of CE-bCT as well as other x-ray-based breast imaging modalities.
Asunto(s)
Yodo , Tomografía Computarizada por Rayos X , Humanos , Tomografía Computarizada por Rayos X/métodos , Mama/diagnóstico por imagen , Mama/patología , Imagenología Tridimensional/métodos , Mamografía/métodos , Fantasmas de ImagenRESUMEN
BACKGROUND: One method for assessing the accuracy of manual, patient-specific, navigational, and robotic-assisted total knee arthroplasty (TKA) instrumentation is to use a post-operative computer tomogram and determine the deviation of the femoral component alignment relative to the planned alignment in the native (i.e. healthy) contralateral distal femoral epiphysis. However, side-to-side asymmetry might introduce errors which inflate alignment deviations. This study quantified asymmetry in the distal femoral epiphysis. METHODS: High resolution CT images (0.5â¯mm slice thickness) were acquired from bilateral lower limb specimens of 13 skeletally mature subjects with no skeletal abnormalities. Images were segmented to generate 3D femur models. Asymmetry was quantified by differences in positions and orientations required to shape-match the distal epiphysis of the mirror 3D femur model to the distal epiphysis of the contralateral 3D femur model. RESULTS: Asymmetry was due to random rather than systematic differences. Random differences (i.e. standard deviations) in proximal-distal (P-D) and anterior-posterior (A-P) positions were 1.1â¯mm and in varus-valgus (V-V) and internal-external (I-E) orientations were 0.9° and 1.3°, respectively. These represented substantial relative errors of up to 50â¯% in previously reported overall alignment deviations. CONCLUSIONS: Although small in an absolute sense, asymmetry of the distal femur epiphysis introduced substantial relative errors when assessing accuracy of femoral component alignment in TKA. When post-operative computer tomograms are used to assess the accuracy of manual, patient specific, navigational, and robotic-assisted TKA instrumentation, the overall deviation should be corrected for asymmetry to better indicate the accuracy of the surgical technique.
Asunto(s)
Artroplastia de Reemplazo de Rodilla , Prótesis de la Rodilla , Humanos , Articulación de la Rodilla/diagnóstico por imagen , Articulación de la Rodilla/cirugía , Fenómenos Biomecánicos , Artroplastia de Reemplazo de Rodilla/métodos , Fémur/diagnóstico por imagen , Fémur/cirugíaRESUMEN
BACKGROUND: Accurate detection and grading of atheromatous stenotic lesions within the cardiac, renal, and intracranial vasculature is imperative for early recognition of disease and guiding treatment strategies. PURPOSE: In this work, a stenotic lesion phantom was used to compare high resolution and normal resolution modes on the same CT scanner in terms of detection and size discrimination performance. MATERIALS AND METHODS: The phantom is comprised of three acrylic cylinders (each 15.0 cm in diameter and 1.3 cm thick) with a matching array of holes in each module. The outer two modules contain holes that are slightly larger than the corresponding hole in the central module to simulate stenotic narrowing in vasculature. The stack of modules was submerged in an iodine solution simulating contrast-enhanced stenotic lesions with a range of lumen diameters (1.32-10.08 mm) and stenosis severity (0%, 50%, 60%, 70%, and 80%). The phantom was imaged on the Canon Aquilion Precision high-resolution CT scanner in high-resolution (HR) mode (0.25 mm × 0.50 mm detector element size) and normal-resolution (NR) mode (0.50 mm × 0.50 mm) using 120 kV and two dose levels (14 and 21 mGy SSDE) with 30 repeat scans acquired for each combination. Filtered back-projection (FBP) and a hybrid-iterative reconstruction (AIDR) were used with the FC18 kernel, as well as a deep learning algorithm (AiCE) which is only available for HR. A non-prewhitening model observer with an eye filter was implemented to quantify performance for detection and size discrimination tasks in the axial plane. RESULTS: Detection performance improved with increasing diameter, dose, and for AIDR in comparison to FBP for a fixed resolution mode. Performance in the HR mode was generally higher than NR for the smaller lumen diameters (1-5 mm) with decreasing differences as the diameter increased. Performance in NR mode surpassed HR mode for lumen diameters greater than â¼4 mm and â¼5 mm for 14 mGy and 21 mGy, respectively. AiCE provided consistently higher detection performance compared with AIDR-FC18 (48% higher for a 6 mm lumen diameter). Discrimination performance increased with increasing nominal diameter, dose, and for larger differences in stenosis severity. When comparing discrimination performance in HR to NR modes, the largest relative differences occur at the smallest nominal diameters and smallest differences in stenosis severity. The AiCE reconstruction algorithm produced the highest overall discrimination performance values, and these were significantly higher than AIDR-FC18 for nominal diameters of 7.14 and 10.08 mm. CONCLUSIONS: HR mode outperforms NR for detection up to a specific diameter and the results improve with AiCE and for higher dose levels. For the task of size discrimination, HR mode consistently outperforms NR if AIDR-FC18 is used for dose levels of at least 21 mGy, and the results improve with AiCE and for the smallest differences in stenosis severity investigated (50% vs. 60%). High-resolution CT appears to be beneficial for detecting smaller simulated lumen diameters (<5 mm) and is generally advantageous for discrimination tasks related to stenotic lesions, which inherently contain information at higher frequencies, given the right reconstruction algorithm and dose level.
Asunto(s)
Algoritmos , Tomografía Computarizada por Rayos X , Humanos , Constricción Patológica/diagnóstico por imagen , Tomografía Computarizada por Rayos X/métodos , Tomógrafos Computarizados por Rayos X , Fantasmas de Imagen , Dosis de RadiaciónRESUMEN
The healthcare sector generates approximately 10% of the total carbon emissions in the United States. Radiology is thought to be a top contributor to the healthcare carbon footprint due to high energy-consuming devices and waste from interventional procedures. In this article, we provide a background on Radiology's environmental impact, describe why hospitals should add sustainability as a quality measure, and give a framework for radiologists to reduce the carbon footprint through quality improvement and collaboration.
Asunto(s)
Radiología , Humanos , Estados Unidos , Huella de Carbono , Atención a la SaludRESUMEN
PURPOSE: To evaluate the bias to the mean glandular dose (MGD) estimates introduced by the homogeneous breast models in digital breast tomosynthesis (DBT) and to have an insight into the glandular dose distributions in 2D (digital mammography, DM) and 3D (DBT and breast dedicated CT, BCT) x-ray breast imaging by employing breast models with realistic glandular tissue distribution and organ silhouette. METHODS: A Monte Carlo software for DM, DBT and BCT simulations was adopted for the evaluation of glandular dose distribution in 60 computational anthropomorphic phantoms. These computational phantoms were derived from 3D breast images acquired via a clinical BCT scanner. RESULTS: g·c·s·T conversion coefficients based on homogeneous breast model led to a MGD overestimate of 18% in DBT when compared to MGD estimated via anthropomorphic phantoms; this overestimate increased up to 21% for recently computed DgNDBT conversion coefficients. The standard deviation of the glandular dose distribution in BCT resulted 60% lower than in DM and 55% lower than in DBT. The glandular dose peak - evaluated as the average value over the 5% of the gland receiving the highest dose - is 2.8 times the MGD in DM, this factor reducing to 2.6 and 1.6 in DBT and BCT, respectively. CONCLUSIONS: Conventional conversion coefficients for MGD estimates based on homogeneous breast models overestimate MGD by 18%, when compared to MGD estimated via anthropomorphic phantoms. The ratio between the peak glandular dose and the MGD is 2.8 in DM. This ratio is 8% and 75% higher than in DBT and BCT, respectively.
Asunto(s)
Neoplasias de la Mama , Mamografía , Mama/diagnóstico por imagen , Femenino , Humanos , Mamografía/métodos , Método de Montecarlo , Fantasmas de Imagen , Tomografía Computarizada por Rayos X/métodosRESUMEN
PURPOSE: Small airways with inner diameters less than 2 mm are sites of major airflow limitations in patients with chronic obstructive pulmonary disease (COPD) and asthma. The purpose of this study is to investigate the limitations for accurate assessment of small airway dimensions using both high-resolution CT (HRCT) and conventional normal-resolution CT at low dose levels. METHODS: To model the normal human airways from the 3rd to 20th generations, a cylindrical polyurethane phantom with 14 airway tubes of inner diameters (ID) ranging from 0.3 to 3.4 mm and wall thicknesses (WT) ranging from 0.15 to 1.6 mm was placed within an Anthropomorphic QRM-Thorax phantom. The Aquilion Precision (Canon Medical Systems Corporation) HRCT scanner was used to acquire images at 80, 100, and 120 kV using high resolution mode (HR, 0.25 mm × 160 detector configuration) and normal-resolution (NR) mode (0.5 mm × 80 detector configuration). The HR data were reconstructed using a 1024 × 1024 matrix (0.22 × 0.22 × 0.25 mm voxel size) and the NR data were reconstructed using a 512 × 512 matrix (0.43 × 0.43 × 0.50 mm). Two reconstruction algorithms (filtered back projection; FBP and an adaptive iterative dose reduction 3D algorithm; AIDR 3D) and three reconstruction kernels (FC30, FC52, and FC56) were investigated. The C T D I vol dose values ranged from 0.2 to 6.2 mGy. A refined automated full-width half-maximum (FWHM) method was used for the measurement of airway dimensions, where the density profiles were computed by radial oversampling using a polar coordinate system. Both ID and WT were compared to the known dimensions using a regression model, and the root-mean-square error (RMSE) and average error were computed across all 14 airway tubes. RESULTS: The results indicate that the ID can be measured within a 15% error down to approximately 0.8 and 2.0 mm using the HR and NR modes, respectively. The overall RMSE (and average error) of ID measurements for HR and NR were 0.10 mm (-0.70%) and 0.31 mm (-2.63%), respectively. The RMSE (and average error) of WT measurements using HR and NR were 0.10 mm (23.27%) and 0.27 mm (53.56%), respectively. The WT measurement using HR yielded a factor of two improvement in accuracy as compared to NR. CONCLUSIONS: High-resolution CT can provide more accurate measurements of airway dimensions as compared with NR CT, potentially improving quantitative assessment of pathologies such as COPD and asthma. The HR mode acquired and reconstructed with AIDR3D and the FC52 kernel provides most accurate measurement of airway dimensions. Low-dose HR measurements at dose level above 0.9 mGy can provide improved accuracy on both inner diameters and wall thicknesses compared to full dose NR airway phantom measurements.
Asunto(s)
Enfermedad Pulmonar Obstructiva Crónica , Tomografía Computarizada por Rayos X , Algoritmos , Humanos , Pulmón/diagnóstico por imagen , Fantasmas de Imagen , Enfermedad Pulmonar Obstructiva Crónica/diagnóstico por imagen , Dosis de RadiaciónRESUMEN
Purpose: To demonstrate the utility of high-resolution micro-computed tomography ( µ CT ) for determining ground-truth size and shape properties of calcium grains for evaluation of detection performance in breast CT (bCT). Approach: Calcium carbonate grains ( â¼ 200 µ m ) were suspended in 1% agar solution to emulate microcalcifications ( µ Calcs ) within a fibroglandular tissue background. Ground-truth imaging was performed on a commercial µ CT scanner and was used for assessing calcium-grain size and shape, and for generating µ Calc signal profiles. Calcium grains were placed within a realistic breast-shaped phantom and imaged on a prototype bCT system at 3- and 6-mGy mean glandular dose (MGD) levels, and the non-prewhitening detectability was assessed. Additionally, the µ CT -derived signal profiles were used in conjunction with the bCT system characterization (MTF and NPS) to obtain predictions of bCT detectability. Results: Estimated detectability of the calcium grains on the bCT system ranged from 2.5 to 10.6 for 3 mGy and from 3.8 to 15.3 for 6 mGy with large fractions of the grains meeting the Rose criterion for visibility. Segmentation of µ CT images based on morphological operations produced accurate results in terms of segmentation boundaries and segmented region size. A regression model linking bCT detectability to µ Calc parameters indicated significant effects of µ Calc size and vertical position within the breast phantom. Detectability using µ CT -derived detection templates and bCT statistical properties (MTF and NPS) were in good correspondence with those measured directly from bCT ( R 2 > 0.88 ). Conclusions: Parameters derived from µ CT ground-truth data were shown to produce useful characterizations of detectability when compared to estimates derived directly from bCT. Signal profiles derived from µ CT imaging can be used in conjunction with measured or hypothesized statistical properties to evaluate the performance of a system, or system component, that may not currently be available.
RESUMEN
Purpose: A computer-aided diagnosis (CADx) system for breast masses is proposed, which incorporates both handcrafted and convolutional radiomic features embedded into a single deep learning model. Approach: The model combines handcrafted and convolutional radiomic signatures into a multi-view architecture, which retrieves three-dimensional (3D) image information by simultaneously processing multiple two-dimensional mass patches extracted along different planes through the 3D mass volume. Each patch is processed by a stream composed of two concatenated parallel branches: a multi-layer perceptron fed with automatically extracted handcrafted radiomic features, and a convolutional neural network, for which discriminant features are learned from the input patches. All streams are then concatenated together into a final architecture, where all network weights are shared and the learning occurs simultaneously for each stream and branch. The CADx system was developed and tested for diagnosis of breast masses ( N = 284 ) using image datasets acquired with independent dedicated breast computed tomography systems from two different institutions. The diagnostic classification performance of the CADx system was compared against other machine and deep learning architectures adopting handcrafted and convolutional approaches, and three board-certified breast radiologists. Results: On a test set of 82 masses (45 benign, 37 malignant), the proposed CADx system performed better than all other model architectures evaluated, with an increase in the area under the receiver operating characteristics curve (AUC) of 0.05 ± 0.02 , and achieving a final AUC of 0.947, outperforming the three radiologists ( AUC = 0.814 - 0.902 ). Conclusions: In conclusion, the system demonstrated its potential usefulness in breast cancer diagnosis by improving mass malignancy assessment.
RESUMEN
PURPOSE: To present a dataset of computational digital breast phantoms derived from high-resolution three-dimensional (3D) clinical breast images for the use in virtual clinical trials in two-dimensional (2D) and 3D x-ray breast imaging. ACQUISITION AND VALIDATION METHODS: Uncompressed computational breast phantoms for investigations in dedicated breast CT (BCT) were derived from 150 clinical 3D breast images acquired via a BCT scanner at UC Davis (California, USA). Each image voxel was classified in one out of the four main materials presented in the field of view: fibroglandular tissue, adipose tissue, skin tissue, and air. For the image classification, a semi-automatic software was developed. The semi-automatic classification was compared via manual glandular classification performed by two researchers. A total of 60 compressed computational phantoms for virtual clinical trials in digital mammography (DM) and digital breast tomosynthesis (DBT) were obtained from the corresponding uncompressed phantoms via a software algorithm simulating the compression and the elastic deformation of the breast, using the tissue's elastic coefficient. This process was evaluated in terms of glandular fraction modification introduced by the compression procedure. The generated cohort of 150 uncompressed computational breast phantoms presented a mean value of the glandular fraction by mass of 12.3%; the average diameter of the breast evaluated at the center of mass was 105 mm. Despite the slight differences between the two manual segmentations, the resulting glandular tissue segmentation did not consistently differ from that obtained via the semi-automatic classification. The difference between the glandular fraction by mass before and after the compression was 2.1% on average. The 60 compressed phantoms presented an average glandular fraction by mass of 12.1% and an average compressed thickness of 61 mm. DATA FORMAT AND ACCESS: The generated digital breast phantoms are stored in DICOM files. Image voxels can present one out of four values representing the different classified materials: 0 for the air, 1 for the adipose tissue, 2 for the glandular tissue, and 3 for the skin tissue. The generated computational phantoms datasets were stored in the Zenodo public repository for research purposes (http://doi.org/10.5281/zenodo.4529852, http://doi.org/10.5281/zenodo.4515360). POTENTIAL APPLICATIONS: The dataset developed within the INFN AGATA project will be used for developing a platform for virtual clinical trials in x-ray breast imaging and dosimetry. In addition, they will represent a valid support for introducing new breast models for dose estimates in 2D and 3D x-ray breast imaging and as models for manufacturing anthropomorphic physical phantoms.
Asunto(s)
Mama , Mamografía , Mama/diagnóstico por imagen , Simulación por Computador , Humanos , Método de Montecarlo , Fantasmas de Imagen , Tomografía Computarizada por Rayos XRESUMEN
PURPOSE: The modulation transfer function (MTF) is widely used as an objective metric of spatial resolution of medical imaging systems. Despite advances in capability for three-dimensional (3D) isotropic spatial resolution in computed tomography (CT) and cone-beam CT (CBCT), MTF evaluation for such systems is typically reported only in the axial plane, and practical methodology for assessment of fully 3D spatial resolution characteristics is lacking. This work reviews fundamental theoretical relationships of two-dimensional (2D) and 3D spread functions and reports practical methods and test tools for analysis of 3D MTF in CBCT. METHODS: Fundamental aspects of 2D and 3D MTF measurement are reviewed within a common notational framework, and three MTF test tools with analysis code are reported and made available online (https://istar.jhu.edu/downloads/): (a) a multi-wire tool for measurement of the axial plane MTF [denoted as M T F ( f r ; φ = 0 ∘ ) , where φ is the measurement angle out of the axial plane] as a function of position in the axial plane; (b) a wedge tool for measurement of the MTF in any direction in the 3D Fourier domain [e.g., φ = 45°, denoted as M T F ( f r ; φ = 45 ∘ ) ]; and (c) a sphere tool for measurement of the MTF in any or all directions in the 3D Fourier domain. Experiments were performed on a mobile C-arm with CBCT capability, showing that M T F ( f r ; φ = 45 ∘ ) yields an informative one-dimensional (1D) representation of the overall 3D spatial resolution characteristics, capturing important characteristics of the 3D MTF that might be missed in conventional analysis. The effects of anisotropic filters and detector readout mode were investigated, and the extent to which a system can be said to provide "isotropic" resolution was evaluated by quantitative comparison of MTF at various φ . RESULTS: All three test tools provided consistent measurement of M T F ( f r ; φ = 0 ∘ ) , and the wedge and sphere tools demonstrated how measurement of the MTF in directions outside the axial plane ( φ > 0 ∘ ) can reveal spatial resolution characteristics to which conventional axial MTF measurement is blind. The wedge tool was shown to reduce statistical measurement error compared to the sphere tool due to improved sampling, and the sphere tool was shown to provide a basis for measurement of the MTF in any or all directions (outside the null cone) from a single scan. The C-arm system exhibited non-isotropic spatial resolution with conventional non-isotropic 1D apodization filters (i.e., frequency cutoff filters) - which is common in CBCT - and implementation of isotropic 2D apodization yielded quantifiably isotropic MTF. Asymmetric pixel binning modes were similarly shown to impart non-isotropic effects on the 3D MTF, and the overall 3D MTF characteristics were evident in each case with a single, 1D measurement of the 1D M T F ( f r ; φ = 45 ∘ ). CONCLUSION: Three test tools and corresponding MTF analysis methods were presented within a consistent framework for analysis of 3D spatial resolution characteristics in a manner amenable to routine, practical measurements. Experiments on a CBCT C-arm validated many intuitive aspects of 3D spatial resolution and quantified the extent to which a CBCT system may be considered to have isotropic resolution. Measurement of M T F ( f r ; φ = 45 ∘ ) provided a practical 1D measure of the underlying 3D MTF characteristics and is extensible to other CT or CBCT systems offering high, isotropic spatial resolution.
Asunto(s)
Tomografía Computarizada de Haz Cónico , Tomografía Computarizada por Rayos X , Fantasmas de ImagenRESUMEN
PURPOSE: The purpose of this study was to quantify location and direction-dependent variations in the 3D modulation transfer function (MTF) of a high-resolution CT scanner with selectable focal spot sizes and resolution modes. METHODS: The Aquilion Precision CT scanner (Canon Medical Systems) has selectable 0.25 mm or 0.5 mm detectors (by binning) in both the axial (x-y) and detector array width (z) directions. For the x-y and z orientations, detectors are configured (x-y) = 0.5 mm/(z) = 0.5 mm for normal resolution (NR), 0.25/0.5 mm for high resolution (HR), and 0.25/0.25 mm for super high resolution (SHR). Six focal spots (FS1-FS6) range in size from 0.4 (x-y) × 0.5 mm (z) for FS1 to 1.6 × 1.4 mm for FS6. Phantoms fabricated from spherical objects were positioned at radial distances of 0, 4.0, 7.5, 11.0, 14.5, and 18.5 cm. Axial and helical acquisitions were utilized and reconstructed using filtered back projection with the FC18 "Body," FC30 "Bone," and FC81 "Bone Sharp" kernels. The reconstructions were used to measure a 1D slice of the 3D MTF by oversampling the 3D ESF in the axial plane [MTF(fr ); φ = 0°)], 45° out of the axial plane [MTF(fr ); φ = 45°)], in the longitudinal direction [MTF(fr ); φ = 80°)], and along the radial and azimuthal directions within the axial plane. RESULTS: The MTF(fr ); φ = 45°) drops to 10% (f10 ) at 1.20, 1.45, and 2.06 mm-1 for NR, HR, and SHR, respectively, for a helical acquisition with FS1, FC30, and r = 4 cm from the isocenter. The MTF(fr ); φ = 45°) includes contributions of both the axial-plane MTF (f10 = 1.10, 2.04, and 2.01 mm-1 ) and the longitudinal MTF (f10 = 1.17, 1.18, and 1.82 mm-1 ) for the NR, HR, and SHR modes, respectively. For SHR, the axial scan mode showed a 15-25% improvement over helical mode in the longitudinal resolution. Helical pitch, ranging from 0.569 to 1.381, did not appreciably affect the 3D resolution (<2%). The radial MTFs across the axial field of view (FOV) showed dependencies on the focal spot length in z; for example, for SHR with FS2 (0.6 × 0.6 mm), f10 at r = 11 cm was within 17% of the value at r = 4 cm, but for SHR with FS3 (0.6 × 1.3), the reduction in f10 was 46% from 4 to 11 cm from the isocenter. The azimuthal MTF also decreased as r increased but less so for longer gantry rotation times and smaller focal spot dimensions in the axial plane. The longitudinal MTF was minimally affected (<11%) by position in the FOV and was principally affected by the focal spot length in the z-dimension. CONCLUSIONS: The 3D MTF was measured throughout the FOV of a high-resolution CT scanner, quantifying the advantages of different resolution modes and focal spot sizes on the axial-plane and longitudinal MTF. Reconstruction kernels were shown to impact axial-plane resolution, imparting non-isotropic 3D resolution characteristics. Focal spot size (both in x-y and in z) and gantry rotation time play important roles in preserving the high-resolution characteristics throughout the field of view for this new high-resolution CT scanner technology.
Asunto(s)
Huesos , Tomografía Computarizada por Rayos X , Fantasmas de Imagen , Tomógrafos Computarizados por Rayos XRESUMEN
OBJECTIVE. The purpose of this study was to investigate whether systematic bias in attenuation measurements occurs among CT scanners made by four major manufacturers and the relevance of this bias regarding opportunistic screening for osteoporosis. MATERIALS AND METHODS. Data on attenuation measurement accuracy were acquired using the American College of Radiology (ACR) accreditation phantom and were evaluated in a blinded fashion for four CT manufacturers (8500 accreditation submissions for manufacturer A; 18,575 for manufacturer B; 8278 for manufacturer C; and 32,039 for manufacturer D). The attenuation value for water, acrylic (surrogate for trabecular bone), and Teflon (surrogate for cortical bone; Chemours) materials for an adult abdominal CT technique (120 kV, 240 mA, standard reconstruction algorithm) was used in the analysis. Differences in attenuation value across all manufacturers were assessed using the Kruskal-Wallis test followed by a post hoc test for pairwise comparisons. RESULTS. The mean attenuation value for water ranged from -0.3 to 2.7 HU, with highly significant differences among all manufacturers (p < 0.001). For the trabecular bone surrogate, differences in attenuation values across all manufacturers were also highly significant (p < 0.001), with mean values of 120.9 (SD, 3.5), 124.6 (3.3), 126.9 (4.4), and 123.9 (3.4) HU for manufacturers A, B, C, and D, respectively. For the cortical bone surrogate, differences in attenuation values across all manufacturers were also highly significant (p < 0.001), with mean values of 939.0 (14.2), 874.3 (13.3), 897.6 (11.3), and 912.7 (13.4) HU for manufacturers A, B, C, and D, respectively. CONCLUSION. CT scanners made by different manufacturers show systematic offsets in attenuation measurement when compared with each other. Knowledge of these off-sets is useful for optimizing the accuracy of opportunistic diagnosis of osteoporosis.
Asunto(s)
Osteoporosis/diagnóstico por imagen , Fantasmas de Imagen , Tomografía Computarizada por Rayos X/instrumentación , Acreditación , Sesgo , Evaluación Educacional , Humanos , Reproducibilidad de los ResultadosRESUMEN
PURPOSE: To develop and evaluate the diagnostic performance of an algorithm for multi-marker radiomic-based classification of breast masses in dedicated breast computed tomography (bCT) images. METHODS: Over 1000 radiomic descriptors aimed at quantifying mass and border heterogeneity, morphology, and margin sharpness were developed and implemented. These included well-established texture and shape feature descriptors, which were supplemented with additional approaches for contour irregularity quantification, spicule and lobe detection, characterization of degree of infiltration, and differences in peritumoral compartments. All descriptors were extracted from a training set of 202 bCT masses (133 benign and 69 malignant), and their individual diagnostic performance was investigated in terms of area under the receiver operating characteristics (ROC) curve (AUC) of single-feature-based linear discriminant analysis (LDA) classifiers. Subsequently, the most relevant descriptors were selected through a multiple-step feature selection process (including stability analysis, statistical significance, evaluation of feature interaction, and dimensionality reduction), and used to develop a final LDA radiomic model for classification of benign and malignant masses, which was then tested on an independent test set of 82 cases (45 benign and 37 malignant). RESULTS: The majority of the individual radiomic descriptors showed, on the training set, an AUC value deriving from a linear decision boundary higher than 0.65, with the lower limit of the associated 95% confidence interval (C.I.) not overlapping with random chance (AUC = 0.5). The final LDA radiomic model resulted in a test set AUC of 0.90 (95% C.I. 0.80-0.96). CONCLUSIONS: The proposed multi-marker radiomic approach achieved high diagnostic accuracy in bCT mass classification, using a radiomic signature based on different feature types. While future studies with larger datasets are needed to further validate these results, quantitative radiomics applied to bCT shows potential to improve the breast cancer diagnosis pipeline.