LAUSR

Radiomics have been extensively investigated as quantitative biomarkers that can enhance the utility of imaging studies and aid the clinical decision making process. A major challenge to the clinical translation of radiomics is their variability as a result of different imaging and reconstruction protocols. In this work, we present a novel radiomics standardization framework capable of modeling and recovering the underlying radiomic feature in images that have been corrupted by the effects of spatial resolution and noise. We focus on two classes of radiomics based on pixel value distributions-i.e. histograms and gray-level co-occurrence matrices (GLCMs). We developed a model that predicts these distributions in the presence of system blur and noise, and used that model to invert these physical effects and recover the underlying distributions. Specifically, the effect of blur on histogram and GLCM is highly image-dependent, while additive noise convolves the histogram/GLCM of the noiseless image with those of the noise. The recovery method therefore consists of two deconvolution operations: the first in the image domain to remove the effect of system blur, the second in the histogram/GLCM domain to remove the effect of noise. The performance of the proposed recovery strategy was investigated using a set of texture phantoms and an emulated computed tomography imaging chain with a range of realistic blur and noise levels. The proposed method was able to obtain histogram and GLCM estimates that closely resemble the ground truth. The method performed well across imaging conditions and significantly lowered the variability associated with different imaging protocols. This improvement also translated to better classification accuracy, where recovered radiomic values result in greater separation of radiomic clusters for two different texture phantoms as compared to values derived from the original blurred and noisy images. In summary, the novel radiomics standardization framework demonstrates high potential for mitigating radiomic variability as a result of the imaging system and can potentially be integrated as a preprocessing step towards more robust and reproducible radiomic models.