Feature request for resampling of rectilinear grids

[mjo22]

Hello, I was wondering what the timeline is on adding interpolation routines for nd rectilinear grids. I am currently a user of jax.scipy.map_coordinates for data on a 3D cartesian grid, and I’m hoping to transition soon to this package for a routine that is similarly fast and supports cubic splines.

[f0uriest]

1, 2, and 3d are currently supported. I'm not sure when I'll have time to look into nd.

For the case you describe I think the existing interp3d method would work, just pass integer ranges for the x, y, z coordinates. IE, if your data is shape(Nx, Ny, Nz, ...) then x = np.arange(Nx) etc.