Description of the desired feature

When trying to create rolling windows with a size larger than the minimum dimension of the region, the function raises a very strange error:

For example, if our region is 20e3 x 10e3 and we try to generate windows with a size of 15e3:

import verde as vd

coordinates = vd.grid_coordinates(region=(-10e3, 10e3, -5e3, 5e3), shape=(30, 30

Validation methods such as cross-validation could leverage multiple threads via Transducers.jl. In this case, we should replace pmap calls (see example) by threaded maps.

This is a good issue for first time contributors.

It would be nice if it supports Thin Plate Spline (TPS) interpolation.