In an effort to generate a more representative icing verification field, an autologistic model with covariates is applied to a variety of case studies. This model is used in the context of binary responses considering spatial correlation and covariates, in an attempt to provide inference over a large, sparsely sampled area of interest. It is motivated by the scarcity of icing observations (i.e., pilot reports of icing conditions) and the existence of easily accessible icing predictors. In this case, the covariate is the Integrated Icing Diagnostic Algorithm (IIDA), developed at the National Center for Atmospheric Research. The set-up is Bayesian, and a Gibbs sampler is implemented to estimate the model parameters. The model is evaluated based on its ability to resolve the relationship between the covariate and observations.