Optimal Mass Transport for Robust Texture Analysis

The emerging field of radiomics, which consists of transforming standard-of-care images to quantifiable scalar statistics, endeavors to reveal the information hidden in these macroscopic images. This field of research has found different applications ranging from phenotyping and tumor classification to outcome prediction and treatment planning. Texture analysis, which often consists of reducing spatial texture matrices to summary scalar features, has been shown to be important in many of the latter applications. However, as pointed out in many studies, some of the derived texture statistics are strongly correlated and tend to contribute redundant information; and are also sensitive to the parameters used in their computation, e.g., the number of gray intensity levels. In the present study, we propose first to consider texture matrices, with an emphasis on gray-level co-occurrence matrix (GLCM), as a non-parametric multivariate objects. The proposed modeling approach avoids evaluating redundant and strongly correlated features and also prevents the feature processing steps. Then, via the Wasserstein distance from optimal mass transport theory, we propose to compare these spatial objects to identify computerized tomography slices with dental artifacts in head and neck cancer. We demonstrate the robustness of the proposed classification approach with respect to the GLCM extraction parameters and the size of the training set. Comparisons with the random forest classifier, which is constructed on scalar texture features, demonstrates the efficiency and robustness of the proposed algorithm.