bug: resample_to_grid function doesn't return correct values when method='nearest'

[wceyhan]

Describe the bug
I am using resample_to_grid raster utility to map elevation data from a raster to the model grid with method='nearest'. When I manually check the resampled data, it is outside the range of the raster cells covering the model element. It seems to me that 'linear', or 'mean' works. I just can't use them due to memory limit.

To Reproduce
Steps to reproduce the behavior:
I can provide the datasets I am using but it returned a similar result when I tested with other raster data.

Expected behavior
The resampled data using method='nearest' will not match the underlying raster cells.

Additional context
I found these two links below. It sounds like the issue originates from the griddata function of scipy. The first link explains what may be going on in griddata and the second link suggests setting rescale=True in griddata. When I modified the resample_to_grid call of the griddata with rescale=True, it works.

• https://stackoverflow.com/questions/78074923/why-does-scipys-interpolate-griddata-return-different-results-for-different-siz
• https://stackoverflow.com/questions/65735210/python-scipy-interpolate-griddata-only-performs-well-in-method-nearest

Thank you!