Multivariate interpolation to incorporate thematic surface data using inverse distance weighting (IDW)