The Image Processing Toolbox offers a great function called imresize, which conveniently rescales gridded datasets and automatically applies antialiasing. However, because imresize interpolates data, imresize output tends to lose data near NaN boundaries. This function uses John D'Errico's inpaint_nans function before applying imresize, then re-masks NaN output.