Multicenter CT phantoms public dataset for radiomics reproducibility tests

Proton FLASH-arc therapy (PFAT): A feasibility study for meeting FLASH dose-rate requirements in the clinic

BACKGROUND AND PURPOSE

Proton arc therapy and FLASH radiotherapy (FLASH-RT) each offer unique advantages in proton therapy. However, clinical translation of FLASH-RT faces challenges in defining and delivering high dose rates. We propose the use of proton FLASH-arc therapy (PFAT) to leverage the benefits of arc while addressing FLASH delivery concerns by spatially fractionating dose delivery to healthy tissue.

MATERIALS AND METHODS

Treatment plans for an abdominal phantom and a clinical brain case were designed in OpenTPS, using monoenergetic beams within a 360-degree gantry rotation. Beams were optimized to achieve target coverage while maximizing spatial fractionation in non-target regions. The temporal dose delivery to healthy-tissue voxels, or in specified organs-at-risk (OARs), was constrained via selective spot removal in the beamlets matrix. The dose, LET, number of spots per voxel, and voxel-wise average dose rate were calculated for each PFAT plan and compared to a corresponding IMPT scenario.

RESULTS

PFAT plans demonstrated comparable dose conformity to IMPT, with LET hotspots shifted towards the target center. The number of spots influencing healthy-tissue voxels was reduced, leading to regions of substantially higher dose rates in many points outside the target. OAR dose-rate optimization in the brain plan resulted in dose rates exceeding 40 Gy/s in the majority of points in the brainstem.

CONCLUSION

The PFAT technique combines the advantages of FLASH and arc therapy, providing improved LET distributions and enhanced biological effect in the target, while achieving high dose rates in healthy tissue, thus reducing healthy tissue damage. This feasibility study demonstrates the capability of PFAT, setting the foundation for further optimization and application in diverse patient cases and complex geometries.

Emulating the Delivery of Sawtooth Proton Arc Therapy Plans on a Cyclotron-Based Proton Beam Therapy System

Purpose: To evaluate and compare the deliverability of 'sawtooth' proton arc therapy (PAT) plans relative to static intensity modulated proton therapy (IMPT) at a cyclotron-based clinical facility. Methods: The delivery of single and dual arc Sawtooth PAT plans for an abdominal CT phantom and multiple clinical cases of brain, head and neck (H&N) and base of skull (BoS) targets was emulated under the step-and-shoot and continuous PAT delivery regimes and compared to that of a corresponding static IMPT plan. Results: Continuous PAT delivery increased the time associated with beam delivery and gantry movement in single/dual PAT plans by 4.86/7.34 min (brain), 7.51/12.40 min (BoS) and 6.59/10.57 min (H&N) on average relative to static IMPT. Step-and-shoot PAT increased this delivery time further by 4.79 min on average as the delivery was limited by gantry motion. Conclusions: The emulator can approximately model clinical sawtooth PAT delivery but requires experimental validation. No clear benefit was observed regarding beam-on time for sawtooth PAT relative to static IMPT.

A reference framework for standardization and harmonization of CT radiomics features on cadaveric sample

Radiomics features (RFs) serve as quantitative metrics to characterize shape, density/intensity, and texture patterns in radiological images. Despite their promise, RFs exhibit reproducibility challenges across acquisition settings, thus limiting implementation into clinical practice. In this investigation, we evaluate the effects of different CT scanners and CT acquisition protocols (KV, mA, field-of-view, and reconstruction kernel settings) on RFs extracted from lumbar vertebrae of a cadaveric trunk. Employing univariate and multivariate Generalized Linear Models (GLM), we evaluated the impact of each acquisition parameter on RFs. Our findings indicate that variations in mA had negligible effects on RFs, while alterations in kV resulted in exponential changes in several RFs, notably First Order (94.4%), GLCM (87.5%), and NGTDM (100%). Moreover, we demonstrated that a tailored GLM model was superior to the ComBat algorithm in harmonizing CT images. GLM achieved R2 > 0.90 in 21 RFs (19.6%), contrasting ComBat's mean R2 above 0.90 in only 1 RF (0.9%). This pioneering study unveils the effects of CT acquisition parameters on bone RFs in cadaveric specimens, highlighting significant variations across parameters and scanner datasets. The proposed GLM model presents a robust solution for mitigating these differences, potentially advancing harmonization efforts in Radiomics-based studies across diverse CT protocols and vendors.

Independent validation of CT radiomics models in colorectal liver metastases: predicting local tumour progression after ablation

OBJECTIVES

Independent internal and external validation of three previously published CT-based radiomics models to predict local tumor progression (LTP) after thermal ablation of colorectal liver metastases (CRLM).

MATERIALS AND METHODS

Patients with CRLM treated with thermal ablation were collected from two institutions to collect a new independent internal and external validation cohort. Ablation zones (AZ) were delineated on portal venous phase CT 2-8 weeks post-ablation. Radiomics features were extracted from the AZ and a 10 mm peri-ablational rim (PAR) of liver parenchyma around the AZ. Three previously published prediction models (clinical, radiomics, combined) were tested without retraining. LTP was defined as new tumor foci appearing next to the AZ up to 24 months post-ablation.

RESULTS

The internal cohort included 39 patients with 68 CRLM and the external cohort 52 patients with 78 CRLM. 34/146 CRLM developed LTP after a median follow-up of 24 months (range 5-139). The median time to LTP was 8 months (range 2-22). The combined clinical-radiomics model yielded a c-statistic of 0.47 (95%CI 0.30-0.64) in the internal cohort and 0.50 (95%CI 0.38-0.62) in the external cohort, compared to 0.78 (95%CI 0.65-0.87) in the previously published original cohort. The radiomics model yielded c-statistics of 0.46 (95%CI 0.29-0.63) and 0.39 (95%CI 0.28-0.52), and the clinical model 0.51 (95%CI 0.34-0.68) and 0.51 (95%CI 0.39-0.63) in the internal and external cohort, respectively.

CONCLUSION

The previously published results for prediction of LTP after thermal ablation of CRLM using clinical and radiomics models were not reproducible in independent internal and external validation.

CLINICAL RELEVANCE STATEMENT

Local tumour progression after thermal ablation of CRLM cannot yet be predicted with the use of CT radiomics of the ablation zone and peri-ablational rim. These results underline the importance of validation of radiomics results to test for reproducibility in independent cohorts.

KEY POINTS

• Previous research suggests CT radiomics models have the potential to predict local tumour progression after thermal ablation in colorectal liver metastases, but independent validation is lacking. • In internal and external validation, the previously published models were not able to predict local tumour progression after ablation. • Radiomics prediction models should be investigated in independent validation cohorts to check for reproducibility.

Investigation of a potential upstream harmonization based on image appearance matching to improve radiomics features robustness: a phantom study

Objective. Radiomics is a promising valuable analysis tool consisting in extracting quantitative information from medical images. However, the extracted radiomics features are too sensitive to variations in used image acquisition and reconstruction parameters. This limited robustness hinders the generalizable validity of radiomics-assisted models. Our aim is to investigate a possible harmonization strategy based on matching image quality to improve feature robustness.Approach.We acquired CT scans of a phantom with two scanners across different dose levels and percentages of Iterative Reconstruction algorithms. The detectability index was used as a comprehensive task-based image quality metric. A statistical analysis based on the Intraclass Correlation Coefficient was performed to determine if matching image quality/appearance could enhance the robustness of radiomics features extracted from the phantom images. Additionally, an Artificial Neural Network was trained on these features to automatically classify the scanner used for image acquisition.Main results.We found that the ICC of the features across protocols providing a similar detectability index improves with respect to the ICC of the features across protocols providing a different detectability index. This improvement was particularly noticeable in features relevant for distinguishing between scanners.Significance.This preliminary study demonstrates that a harmonization based on image quality/appearance matching could improve radiomics features robustness and heterogeneous protocols can be used to obtain a similar image appearance in terms of the detectability index. Thus protocols with a lower dose level could be selected to reduce the amount of radiation dose delivered to the patient and simultaneously obtain a more robust quantitative analysis.

CT Cadaveric dataset for Radiomics features stability assessment in lumbar vertebrae

Radiomics features (RFs) studies have showed limitations in the reproducibility of RFs in different acquisition settings. To date, reproducibility studies using CT images mainly rely on phantoms, due to the harness of patient exposure to X-rays. The provided CadAIver dataset has the aims of evaluating how CT scanner parameters effect radiomics features on cadaveric donor. The dataset comprises 112 unique CT acquisitions of a cadaveric truck acquired on 3 different CT scanners varying KV, mA, field-of-view, and reconstruction kernel settings. Technical validation of the CadAIver dataset comprises a comprehensive univariate and multivariate GLM approach to assess stability of each RFs extracted from lumbar vertebrae. The complete dataset is publicly available to be applied for future research in the RFs field, and could foster the creation of a collaborative open CT image database to increase the sample size, the range of available scanners, and the available body districts.

Radiogenomics: a key component of precision cancer medicine

Radiogenomics, focusing on the relationship between genomics and imaging phenotypes, has been widely applied to address tumour heterogeneity and predict immune responsiveness and progression. It is an inevitable consequence of current trends in precision medicine, as radiogenomics costs less than traditional genetic sequencing and provides access to whole-tumour information rather than limited biopsy specimens. By providing voxel-by-voxel genetic information, radiogenomics can allow tailored therapy targeting a complete, heterogeneous tumour or set of tumours. In addition to quantifying lesion characteristics, radiogenomics can also be used to distinguish benign from malignant entities, as well as patient characteristics, to better stratify patients according to disease risk, thereby enabling more precise imaging and screening. Here, we have characterised the radiogenomic application in precision medicine using a multi-omic approach. we outline the main applications of radiogenomics in diagnosis, treatment planning and evaluations in the field of oncology with the aim of developing quantitative and personalised medicine. Finally, we discuss the challenges in the field of radiogenomics and the scope and clinical applicability of these methods.

3D-printed iodine-ink CT phantom for radiomics feature extraction - advantages and challenges

BACKGROUND

To test and validate novel CT techniques, such as texture analysis in radiomics, repeat measurements are required. Current anthropomorphic phantoms lack fine texture and true anatomic representation. 3D-printing of iodinated ink on paper is a promising phantom manufacturing technique. Previously acquired or artificially created CT data can be used to generate realistic phantoms.

PURPOSE

To present the design process of an anthropomorphic 3D-printed iodine ink phantom, highlighting the different advantages and pitfalls in its use. To analyze the phantom's X-ray attenuation properties, and the influences of the printing process on the imaging characteristics, by comparing it to the original input dataset.

METHODS

Two patient CT scans and artificially generated test patterns were combined in a single dataset for phantom printing and cropped to a size of 26 × 19 × 30 cm3 . This DICOM dataset was printed on paper using iodinated ink. The phantom was CT-scanned and compared to the original image dataset used for printing the phantom. The water-equivalent diameter of the phantom was compared to that of a patient cohort (N = 104). Iodine concentrations in the phantom were measured using dual-energy CT. 86 radiomics features were extracted from 10 repeat phantom scans and the input dataset. Features were compared using a histogram analysis and a PCA individually and overall, respectively. The frequency content was compared using the normalized spectrum modulus.

RESULTS

Low density structures are depicted incorrectly, while soft tissue structures show excellent visual accordance with the input dataset. Maximum deviations of around 30 HU between the original dataset and phantom HU values were observed. The phantom has X-ray attenuation properties comparable to a lightweight adult patient (∼54 kg, BMI 19 kg/m2 ). Iodine concentrations in the phantom varied between 0 and 50 mg/ml. PCA of radiomics features shows different tissue types separate in similar areas of PCA representation in the phantom scans as in the input dataset. Individual feature analysis revealed systematic shift of first order radiomics features compared to the original dataset, while some higher order radiomics features did not. The normalized frequency modulus |f(ω)| of the phantom data agrees well with the original data. However, all frequencies systematically occur more frequently in the phantom compared to the maximum of the spectrum modulus than in the original data set, especially for mid-frequencies (e.g., for ω = 0.3942 mm-1 , |f(ω)|original  = 0.09 * |fmax |original and |f(ω)|phantom  = 0.12 * |fmax |phantom ).

CONCLUSIONS

3D-iodine-ink-printing technology can be used to print anthropomorphic phantoms with a water-equivalent diameter of a lightweight adult patient. Challenges include small residual air enclosures and the fidelity of HU values. For soft tissue, there is a good agreement between the HU values of the phantom and input data set. Radiomics texture features of the phantom scans are similar to the input data set, but systematic shifts of radiomics features in first order features, due to differences in HU values, need to be considered. The paper substrate influences the spatial frequency distribution of the phantom scans. This phantom type is of very limited use for dual-energy CT analyses.

A customized anthropomorphic 3D-printed phantom to reproducibility assessment in computed tomography: an oncological case study

Introduction

Studies on computed tomography (CT) reproducibility at different acquisition parameters have to take into account radiation dose administered and related ethical issues. 3D-printed phantoms provide the possibility to investigate these features deeply and to foster CT research, also taking advantage by outperforming new generation scanners. The aim of this study is to propose a new anthropomorphic 3D-printed phantom for chest lesions, tailored on a real patient CT scan, to investigate the variability of volume and Hounsfield Unit (HU) measurements at different CT acquisition parameters.

Methods

The chest CT of a 75-year-old patient with a paramediastinal lung lesion was segmented based on an eight-compartment approach related to HU ranges (air lung, lung interstitium, fat, muscle, vascular, skin, bone, and lesion). From each mask produced, the 3D.stl model was exported and linked to a different printing infill value, based on a preliminary test and HU ratios derived from the patient scan. Fused deposition modeling (FDM) technology printing was chosen with filament materials in polylactic acid (PLA). Phantom was acquired at 50 mAs and three different tube voltages of 80, 100, and 120 kVp on two different scanners, namely, Siemens Somatom Force (Siemens Healthineers, Erlangen, Germany; same setting of real patient for 80 kVp acquisition) and GE 750 HD CT (GE Healthcare, Chicago, IL). The same segmentation workflow was then applied on each phantom acquisition after coregistration pipeline, and Dice Similarity Coefficient (DSC) and HU averages were extracted and compared for each compartment.

Results

DSC comparison among real patient versus phantom scans at different kVp, and on both CT scanners, demonstrated a good overlap of different compartments and lesion vascularization with a higher similarity for lung and lesion masks for each setting (about 0.9 and 0.8, respectively). Although mean HU was not comparable with real data, due to the PLA material, the proportion of intensity values for each compartment remains respected.

Discussion

The proposed approach demonstrated the reliability of 3D-printed technology for personalized approaches in CT research, opening to the application of the same workflow to other oncological fields.

Bridging the simulation-to-real gap for AI-based needle and target detection in robot-assisted ultrasound-guided interventions

BACKGROUND

Artificial intelligence (AI)-powered, robot-assisted, and ultrasound (US)-guided interventional radiology has the potential to increase the efficacy and cost-efficiency of interventional procedures while improving postsurgical outcomes and reducing the burden for medical personnel.

METHODS

To overcome the lack of available clinical data needed to train state-of-the-art AI models, we propose a novel approach for generating synthetic ultrasound data from real, clinical preoperative three-dimensional (3D) data of different imaging modalities. With the synthetic data, we trained a deep learning-based detection algorithm for the localization of needle tip and target anatomy in US images. We validated our models on real, in vitro US data.

RESULTS

The resulting models generalize well to unseen synthetic data and experimental in vitro data making the proposed approach a promising method to create AI-based models for applications of needle and target detection in minimally invasive US-guided procedures. Moreover, we show that by one-time calibration of the US and robot coordinate frames, our tracking algorithm can be used to accurately fine-position the robot in reach of the target based on 2D US images alone.

CONCLUSIONS

The proposed data generation approach is sufficient to bridge the simulation-to-real gap and has the potential to overcome data paucity challenges in interventional radiology. The proposed AI-based detection algorithm shows very promising results in terms of accuracy and frame rate.

RELEVANCE STATEMENT

This approach can facilitate the development of next-generation AI algorithms for patient anatomy detection and needle tracking in US and their application to robotics.

KEY POINTS

• AI-based methods show promise for needle and target detection in US-guided interventions. • Publicly available, annotated datasets for training AI models are limited. • Synthetic, clinical-like US data can be generated from magnetic resonance or computed tomography data. • Models trained with synthetic US data generalize well to real in vitro US data. • Target detection with an AI model can be used for fine positioning of the robot.

Pilot Study for the Assessment of the Best Radiomic Features for Bosniak Cyst Classification Using Phantom and Radiologist Inter-Observer Selection

Since the Bosniak cysts classification is highly reader-dependent, automated tools based on radiomics could help in the diagnosis of the lesion. This study is an initial step in the search for radiomic features that may be good classifiers of benign-malignant Bosniak cysts in machine learning models. A CCR phantom was used through five CT scanners. Registration was performed with ARIA software, while Quibim Precision was used for feature extraction. R software was used for the statistical analysis. Robust radiomic features based on repeatability and reproducibility criteria were chosen. Excellent correlation criteria between different radiologists during lesion segmentation were imposed. With the selected features, their classification ability in benignity-malignity terms was assessed. From the phantom study, 25.3% of the features were robust. For the study of inter-observer correlation (ICC) in the segmentation of cystic masses, 82 subjects were prospectively selected, finding 48.4% of the features as excellent regarding concordance. Comparing both datasets, 12 features were established as repeatable, reproducible, and useful for the classification of Bosniak cysts and could serve as initial candidates for the elaboration of a classification model. With those features, the Linear Discriminant Analysis model classified the Bosniak cysts in terms of benignity or malignancy with 88.2% accuracy.

An Innovative Low-dose CT Inpainting Algorithm based on Limited-angle Imaging Inpainting Model

BACKGROUND

With the popularity of computed tomography (CT) technique, an increasing number of patients are receiving CT scans. Simultaneously, the public's attention to CT radiation dose is also increasing. How to obtain CT images suitable for clinical diagnosis while reducing the radiation dose has become the focus of researchers.

OBJECTIVE

To demonstrate that limited-angle CT imaging technique can be used to acquire lower dose CT images, we propose a generative adversarial network-based image inpainting model-Low-dose imaging and Limited-angle imaging inpainting Model (LDLAIM), this method can effectively restore low-dose CT images with limited-angle imaging, which verifies that limited-angle CT imaging technique can be used to acquire low-dose CT images.

METHODS

In this work, we used three datasets, including chest and abdomen dataset, head dataset and phantom dataset. They are used to synthesize low-dose and limited-angle CT images for network training. During training stage, we divide each dataset into training set, validation set and testing set according to the ratio of 8:1:1, and use the validation set to validate after finishing an epoch training, and use the testing set to test after finishing all the training. The proposed method is based on generative adversarial networks(GANs), which consists of a generator and a discriminator. The generator consists of residual blocks and encoder-decoder, and uses skip connection.

RESULTS

We use SSIM, PSNR and RMSE to evaluate the performance of the proposed method. In the chest and abdomen dataset, the mean SSIM, PSNR and RMSE of the testing set are 0.984, 35.385 and 0.017, respectively. In the head dataset, the mean SSIM, PSNR and RMSE of the testing set are 0.981, 38.664 and 0.011, respectively. In the phantom dataset, the mean SSIM, PSNR and RMSE of the testing set are 0.977, 33.468 and 0.022, respectively. By comparing the experimental results of other algorithms in these three datasets, it can be found that the proposed method is superior to other algorithms in these indicators. Meanwhile, the proposed method also achieved the highest score in the subjective quality score.

CONCLUSIONS

Experimental results show that the proposed method can effectively restore CT images when both low-dose CT imaging techniques and limited-angle CT imaging techniques are used simultaneously. This work proves that the limited-angle CT imaging technique can be used to reduce the CT radiation dose, and also provides a new idea for the research of low-dose CT imaging.

Radiomics in clinical trials: perspectives on standardization

The term biomarker is used to describe a biological measure of the disease behavior. The existing imaging biomarkers are associated with the known tissue biological characteristics and follow a well-established roadmap to be implemented in routine clinical practice. Recently, a new quantitative imaging analysis approach named radiomics has emerged. It refers to the extraction of a large number of advanced imaging features with high-throughput computing. Extensive research has demonstrated its value in predicting disease behavior, progression, and response to therapeutic options. However, there are numerous challenges to establishing it as a clinically viable solution, including lack of reproducibility and transparency. The data-driven nature also does not offer insights into the underpinning biology of the observed relationships. As such, additional effort is needed to establish it as a qualified biomarker to inform clinical decisions. Here we review the technical difficulties encountered in the clinical applications of radiomics and current effort in addressing some of these challenges in clinical trial designs. By addressing these challenges, the true potential of radiomics can be unleashed.

A Review of Radiomics and Artificial Intelligence and Their Application in Veterinary Diagnostic Imaging

Simple Summary

The goal of this paper is to provide an overview of current radiomic and AI applications in veterinary diagnostic imaging. We discuss the essential elements of AI for veterinary practitioners with the aim of helping them make informed decisions in applying AI technologies to their practices and that veterinarians will play an integral role in ensuring the appropriate uses and suitable curation of data. The expertise of veterinary professionals will be vital to ensuring suitable data and, subsequently, AI that meets the needs of the profession.

Abstract

Great advances have been made in human health care in the application of radiomics and artificial intelligence (AI) in a variety of areas, ranging from hospital management and virtual assistants to remote patient monitoring and medical diagnostics and imaging. To improve accuracy and reproducibility, there has been a recent move to integrate radiomics and AI as tools to assist clinical decision making and to incorporate it into routine clinical workflows and diagnosis. Although lagging behind human medicine, the use of radiomics and AI in veterinary diagnostic imaging is becoming more frequent with an increasing number of reported applications. The goal of this paper is to provide an overview of current radiomic and AI applications in veterinary diagnostic imaging.

Radiomics: A Primer on Processing Workflow and Analysis

Quantitative analysis of medical images can provide objective tools for diagnosis, prognostication, and disease monitoring. Radiomics refers to automated extraction of a large number of quantitative features from medical images for characterization of underlying pathologies. In this review, we will discuss the principles of radiomics, image preprocessing, feature extraction workflow, and statistical analysis. We will also address the limitations and future directions of radiomics.

XNAT-PIC: Extending XNAT to Preclinical Imaging Centers

Molecular imaging generates large volumes of heterogeneous biomedical imagery with an impelling need of guidelines for handling image data. Although several successful solutions have been implemented for human epidemiologic studies, few and limited approaches have been proposed for animal population studies. Preclinical imaging research deals with a variety of machinery yielding tons of raw data but the current practices to store and distribute image data are inadequate. Therefore, standard tools for the analysis of large image datasets need to be established. In this paper, we present an extension of XNAT for Preclinical Imaging Centers (XNAT-PIC). XNAT is a worldwide used, open-source platform for securely hosting, sharing, and processing of clinical imaging studies. Despite its success, neither tools for importing large, multimodal preclinical image datasets nor pipelines for processing whole imaging studies are yet available in XNAT. In order to overcome these limitations, we have developed several tools to expand the XNAT core functionalities for supporting preclinical imaging facilities. Our aim is to streamline the management and exchange of image data within the preclinical imaging community, thereby enhancing the reproducibility of the results of image processing and promoting open science practices.

The impact of phantom design and material-dependence on repeatability and reproducibility of CT-based radiomics features

PURPOSE

To understand the design of radiomics phantom and material-dependence on repeatability and reproducibility of CT radiomics features METHODS: : A radiomics phantom consisting of various materials with uniformity, textural and biological components, was constructed. The phantom was scanned with different manufacturer CT scanners and the scans were repeated multiple times on the same scanner with different acquisition settings as kVp, mAs, orientation, field of view (FOV), slice thickness, pitch, reconstruction kernels and acquisition mode. A total of 72 phantom scans were included. For each scan, 18 different regions of interest (ROI) were contoured and 708 radiomics features were extracted from each ROI via an open source radiomics tool, IBEX. To relate the phantom data to patient data, the radiomics features from different phantom materials were compared with those extracted from 50 patients' images of five disease sites as brain, head-and-neck, breast, liver and lung cases using box-plots comparison and principal component analysis (PCA). The temporal stability of imaging features was then evaluated with respect to a controlled scenario (test-retest) via the intra-class correlation coefficient (ICC). The reproducibility of radiomics features with respect to different scanners or acquisition settings were further evaluated with concordance correlation coefficients (CCC).

RESULTS

Among all phantom materials, the biological component had feature values closest to human tissues, especially for tumors in brain and liver. The textural component showed similar ranges of variation to lung lesions, particularly for cartridges of rice, cereal, and the 3D-printed textural phantom with fine and rough-grid. It also showed that certain materials, such as polystyrene foam, plaster and peanuts, did not have comparable values to human tissue and could be excluded for future phantom design. High repeatability was observed in the test-retest study as indicated by an ICC value of 0.998 ± 0.020. All materials were used for feature stability analysis. For the inter-scanner study, shape-related features were the most-reliable category with 94% of features having CCC ≥0.9, while GOH were the least-reliable with only 14.6% meeting the criteria. For the intra-scanner study, the reproducibility of CT-based radiomics features showed material-dependence. In general, the instability of radiomics features introduced by kVp, mAs, pitch, acquisition mode and orientation were relatively mild. However, the homogeneous materials were more vulnerable to those changes compared to materials with textural patterns. Regardless of material compositions, resolution parameters like FOV and slice thickness, could have large impact on feature stability. Switching between standard and bone reconstruction kernels could also result significant changes to feature reproducibility.

CONCLUSION

We have built a radiomics phantom using materials that cover a wide span of tumor textures seen in oncological CT images. The designed phantom presents a preliminary opportunity for investigating reproducibility of radiomics features and the reproducibility can be material dependent. Thus, in the radiomics quality assurance design, it is important to choose appropriate materials that can provide a close range of radiomics features to patients with specific disease sites dependency taken into consideration. This article is protected by copyright. All rights reserved.

Cardiac SPECT radiomic features repeatability and reproducibility: A multi-scanner phantom study

BACKGROUND

The aim of this work was to assess the robustness of cardiac SPECT radiomic features against changes in imaging settings, including acquisition, and reconstruction parameters.

METHODS

Four commercial SPECT and SPECT/CT cameras were used to acquire images of a static cardiac phantom mimicking typical myorcardial perfusion imaging using 185 MBq of 99mTc. The effects of different image acquisition and reconstruction parameters, including number of views, view matrix size, attenuation correction, as well as image reconstruction related parameters (algorithm, number of iterations, number of subsets, type of post-reconstruction filter, and its associated parameters, including filter order and cut-off frequency) were studied. In total, 5,063 transverse views were reconstructed by varying the aforementioned factors. Eighty-seven radiomic features including first-, second-, and high-order textures were extracted from these images. To assess reproducibility and repeatability, the coefficient of variation (COV), as a widely adopted metric, was measured for each of the radiomic features over the different imaging settings.

RESULTS

The Inverse Difference Moment Normalized (IDMN) and Inverse Difference Normalized (IDN) features from the Gray Level Co-occurrence Matrix (GLCM), Run Percentage (RP) from the Gray Level Co-occurrence Matrix (GLRLM), Zone Entropy (ZE) from the Gray Level Size Zone Matrix (GLSZM), and Dependence Entropy (DE) from the Gray Level Dependence Matrix (GLDM) feature sets were the only features that exhibited high reproducibility (COV ≤ 5%) against changes in all imaging settings. In addition, Large Area Low Gray Level Emphasis (LALGLE), Small Area Low Gray Level Emphasis (SALGLE) and Low Gray Level Zone Emphasis (LGLZE) from GLSZM, and Small Dependence Low Gray Level Emphasis (SDLGLE) from GLDM feature sets turned out to be less reproducible (COV > 20%) against changes in imaging settings. The GLRLM (31.88%) and GLDM feature set (54.2%) had the highest (COV < 5%) and lowest (COV > 20%) number of the reproducible features, respectively. Matrix size had the largest impact on feature variability as most of the features were not repeatable when matrix size was modified with 82.8% of them having a COV > 20%.

CONCLUSION

The repeatability and reproducibility of SPECT/CT cardiac radiomic features under different imaging settings is feature-dependent. Different image acquisition and reconstruction protocols have variable effects on radiomic features. The radiomic features exhibiting low COV are potential candidates for future clinical studies.

A Multicentre Evaluation of Dosiomics Features Reproducibility, Stability and Sensitivity

Simple Summary

Dosiomics is born directly as an extension of radiomics: it entails extracting features from the patients’ three-dimensional (3D) radiotherapy dose distribution rather than from conventional medical images to obtain specific spatial and statistical information. Dosiomic studies, in a multicentre setting, require assessing the features’ stability to dose calculation settings and the features’ capability in distinguishing different dose distributions. This study provides the first multicentre evaluation of the dosiomic features in terms of reproducibility, stability and sensitivity across various dose distributions obtained from multiple technologies and techniques and considering different dose calculation algorithms of TPS and two different resolutions of the dose grid. Harmonisation strategies to account for a possible variation in the dose distribution due to these confounding factors should be adopted when investigating a correlation between dosiomic features and clinical outcomes in multicentre studies.

Abstract

Dosiomics is a texture analysis method to produce dose features that encode the spatial 3D distribution of radiotherapy dose. Dosiomic studies, in a multicentre setting, require assessing the features’ stability to dose calculation settings and the features’ capability in distinguishing different dose distributions. Dose distributions were generated by eight Italian centres on a shared image dataset acquired on a dedicated phantom. Treatment planning protocols, in terms of planning target volume coverage and dose–volume constraints to the organs at risk, were shared among the centres to produce comparable dose distributions for measuring reproducibility/stability and sensitivity of dosiomic features. In addition, coefficient of variation (CV) was employed to evaluate the dosiomic features’ variation. We extracted 38,160 features from 30 different dose distributions from six regions of interest, grouped by four features’ families. A selected group of features (CV < 3 for the reproducibility/stability studies, CV > 1 for the sensitivity studies) were identified to support future multicentre studies, assuring both stable features when dose distributions variation is minimal and sensitive features when dose distribution variations need to be clearly identified. Dosiomic is a promising tool that could support multicentre studies, especially for predictive models, and encode the spatial and statistical characteristics of the 3D dose distribution.

Quality control and whole-gland, zonal and lesion annotations for the PROSTATEx challenge public dataset

PURPOSE

Radiomic features are promising quantitative parameters that can be extracted from medical images and employed to build machine learning predictive models. However, generalizability is a key concern, encouraging the use of public image datasets. We performed a quality assessment of the PROSTATEx training dataset and provide publicly available lesion, whole-gland, and zonal anatomy segmentation masks.

METHOD

Two radiology residents and two experienced board-certified radiologists reviewed the 204 prostate MRI scans (330 lesions) included in the training dataset. The quality of provided lesion coordinate was scored using the following scale: 0 = perfectly centered, 1 = within lesion, 2 = within the prostate without lesion, 3 = outside the prostate. All clearly detectable lesions were segmented individually slice-by-slice on T2-weighted and apparent diffusion coefficient images. With the same methodology, volumes of interest including the whole gland, transition, and peripheral zones were annotated.

RESULTS

Of the 330 available lesion identifiers, 3 were duplicates (1%). From the remaining, 218 received score = 0, 74 score = 1, 31 score = 2 and 4 score = 3. Overall, 299 lesions were verified and segmented. Independently of lesion coordinate score and other issues (e.g., lesion coordinates falling outside DICOM images, artifacts etc.), the whole prostate gland and zonal anatomy were also manually annotated for all cases.

CONCLUSION

While several issues were encountered evaluating the original PROSTATEx dataset, the improved quality and availability of lesion, whole-gland and zonal segmentations will increase its potential utility as a common benchmark in prostate MRI radiomics.

Radiogenomics in Colorectal Cancer

Simple Summary

Colorectal carcinoma is characterized by intratumoral heterogeneity that can be assessed by radiogenomics. Radiomics, high-throughput quantitative data extracted from medical imaging, combined with molecular analysis, through genomic and transcriptomic data, is expected to lead to significant advances in personalized medicine. However, a radiogenomics approach in colorectal cancer is still in its early stages and many problems remain to be solved. Here we review the progress and challenges in this field at its current stage, as well as future developments.

Abstract

The steady improvement of high-throughput technologies greatly facilitates the implementation of personalized precision medicine. Characterization of tumor heterogeneity through image-derived features—radiomics and genetic profile modifications—genomics, is a rapidly evolving field known as radiogenomics. Various radiogenomics studies have been dedicated to colorectal cancer so far, highlighting the potential of these approaches to enhance clinical decision-making. In this review, a general outline of colorectal radiogenomics literature is provided, discussing the current limitations and suggested further developments.

Harmonization of radiomic feature distributions: impact on classification of hepatic tissue in CT imaging

OBJECTIVES

Following the craze for radiomic features (RF), their lack of reliability raised the question of the generalizability of classification models. Inter-site harmonization of images therefore becomes a central issue. We compared RF harmonization processing designed to detect liver diseases in CT images.

METHODS

We retrospectively analyzed 76 multi-center portal CT series of non-diseased (NDL) and diseased liver (DL) patients. In each series, we positioned volumes of interest in spleen and liver, then extracted 9 RF (histogram and texture). We evaluated two RF harmonization approaches. First, in each series, we computed the Z-score of liver measurements based on those computed in the spleen. Second, we evaluated the ComBat method according to each imaging center; parameters were computed in the spleen and applied to the liver. We compared RF distributions and classification performances before/after harmonization. We classified NDL versus spleen and versus DL tissues.

RESULTS

The RF distributions were all different between liver and spleen (p < 0.05). The Z-score harmonization outperformed for the detection of liver versus spleen: AUC = 93.1% (p < 0.001). For the detection of DL versus NDL, in a case/control setting, we found no differences between the harmonizations: mean AUC = 73.6% (p = 0.49). Using the whole datasets, the performances were improved using ComBat (p = 0.05) AUC = 82.4% and degraded with Z-score AUC = 67.4% (p = 0.008).

CONCLUSIONS

Data harmonization requires to first focus on data structuring to not degrade the performances of subsequent classifications. Liver tissue classification after harmonization of spleen-based RF is a promising strategy for improving the detection of DL tissue.

KEY POINTS

• Variability of acquisition parameter makes radiomics of CT features non-reproducible. • Data harmonization can help circumvent the inter-site variability of acquisition protocols. • Inter-site harmonization must be carefully implemented and requires designing consistent data sets.

Prospects and Challenges of Radiomics by Using Nononcologic Routine Chest CT

Chest CT scans are one of the most common medical imaging procedures. The automatic extraction and quantification of imaging features may help in diagnosis, prognosis of, or treatment decision in cardiovascular, pulmonary, and metabolic diseases. However, an adequate sample size as a statistical necessity for radiomics studies is often difficult to achieve in prospective trials. By exploiting imaging data from clinical routine, a much larger amount of data could be used than in clinical trials. Still, there is only little literature on the implementation of radiomics in clinical routine chest CT scans. Reasons are heterogeneous CT scanning protocols and the resulting technical variability (eg, different slice thicknesses, reconstruction kernels or timings after contrast material administration) in routine CT imaging data. This review summarizes the recent state of the art of studies aiming to develop quantifiable imaging biomarkers at chest CT, such as for osteoporosis, chronic obstructive pulmonary disease, interstitial lung disease, and coronary artery disease. This review explains solutions to overcome heterogeneity in routine data such as the use of imaging repositories, the standardization of radiomic features, algorithmic approaches to improve feature stability, test-retest studies, and the evolution of deep learning for modeling radiomics features. Supplemental material is available for this article. © RSNA, 2020 See also the commentary by Kay in this issue.

Radiomics in medical imaging-"how-to" guide and critical reflection

Radiomics is a quantitative approach to medical imaging, which aims at enhancing the existing data available to clinicians by means of advanced mathematical analysis. Through mathematical extraction of the spatial distribution of signal intensities and pixel interrelationships, radiomics quantifies textural information by using analysis methods from the field of artificial intelligence. Various studies from different fields in imaging have been published so far, highlighting the potential of radiomics to enhance clinical decision-making. However, the field faces several important challenges, which are mainly caused by the various technical factors influencing the extracted radiomic features.The aim of the present review is twofold: first, we present the typical workflow of a radiomics analysis and deliver a practical "how-to" guide for a typical radiomics analysis. Second, we discuss the current limitations of radiomics, suggest potential improvements, and summarize relevant literature on the subject.

A review of original articles published in the emerging field of radiomics

PURPOSE

To determine the characteristics of and trends in research in the emerging field of radiomics through bibliometric and hotspot analyses of relevant original articles published between 2013 and 2018.

METHODS

We evaluated 553 original articles concerning radiomics, published in a total of 61 peer-reviewed journals between 2013 and 2018. The following information was retrieved for each article: radiological subspecialty, imaging technique(s), machine learning technique(s), sample size, study setting and design, statistical result(s), study purpose, software used for feature calculation, funding declarations, author number, first author's affiliation, study origin, and journal name. Qualitative and quantitative analyses were performed for the manually extracted data for identification and visualization of the trends in radiomics research.

RESULTS

The annual growth rate in the number of published papers was 177.82% (p < 0.001). The characteristics and trends of research hotspots in the field of radiomics were clarified and visualized in this study. It was found that the field of radiomics is at a more mature stage for lung, breast, and prostate cancers than for other sites. Radiomics studies primarily focused on radiological characterization (215) and monitoring (182). Logistic regression and LASSO were the two most commonly used techniques for feature selection. Non-clinical researchers without a medical background dominated radiomics studies (70.52%), the vast majority of which only highlighted positive results (97.80%) while downplaying negative findings.

CONCLUSIONS

The reporting of quantifiable knowledge about the characteristics and trajectories of radiomics can inform researchers about the gaps in the field of radiomics and guide its future direction.

Technical Note: Virtual phantom analyses for preprocessing evaluation and detection of a robust feature set for MRI-radiomics of the brain

Purpose

The purpose of the paper was to use a virtual phantom to identify a set of radiomic features from T1‐weighted and T2‐weighted magnetic resonance imaging (MRI) of the brain which is stable to variations in image acquisition parameters and to evaluate the effect of image preprocessing on radiomic features stability.

Methods

Stability to different sources of variability (time of repetition and echo, voxel size, random noise and intensity non‐uniformity) was evaluated for both T1‐weighted and T2‐weighted MRI images. A set of 107 radiomic features, accounting for shape and size, first order statistics, and textural features was used. Feature stability was quantified using intraclass correlation coefficient (ICC). For each source of variability, stability was evaluated before and after preprocessing (Z‐score normalization, resampling, gaussian filtering and bias field correction). Features that have ICC > 0.75 in all the analysis of variability are selected as stable features. Last, the robust feature sets were tested on images acquired with random simulation parameters to assess their generalizability to unseen conditions.

Results

Preprocessing significantly increased the robustness of radiomic features to the different sources of variability. When preprocessing is applied, a set of 67 and 61 features resulted as stable for T1‐weighted and T2‐wieghted images respectively, over 80% of which were confirmed by the analysis on the images acquired with random simulation parameters.

Conclusion

A set of MRI‐radiomic features, robust to changes in TR/TE/PS/ST, was identified. This set of features may be used in radiomic analyses based on T1‐weighted and T2‐weighted MRI images.