Background: We recently developed a gene-expression-based HOT score to identify the hot/cold phenotype of head and neck squamous cell carcinomas (HNSCCs), which is associated with the response to immunotherapy. Our goal was to determine whether radiomic profiling from computed tomography (CT) scans can distinguish hot and cold HNSCC. Method: We included 113 patients from The Cancer Genome Atlas (TCGA) and 20 patients from the Groupe Hospitalier Pitié-Salpêtrière (GHPS) with HNSCC, all with available pre-treatment CT scans. The hot/cold phenotype was computed for all patients using the HOT score. The IBEX software (version 4.11.9, accessed on 30 march 2020) was used to extract radiomic features from the delineated tumor region in both datasets, and the intraclass correlation coefficient (ICC) was computed to select robust features. Machine learning classifier models were trained and tested in the TCGA dataset and validated using the area under the receiver operator characteristic curve (AUC) in the GHPS cohort. Results: A total of 144 radiomic features with an ICC >0.9 was selected. An XGBoost model including these selected features showed the best performance prediction of the hot/cold phenotype with AUC = 0.86 in the GHPS validation dataset. Conclusions and Relevance: We identified a relevant radiomic model to capture the overall hot/cold phenotype of HNSCC. This non-invasive approach could help with the identification of patients with HNSCC who may benefit from immunotherapy.