Regridding data from 2D grid to another

Question

I'm working on a project involving netcdf data analysis. To correctly process the data I need to regrid it to a coarser level. The matrix is in fact multidimensional, but I only need spatial regridding, so basically I have to work with a lat/lon matrix.

The process I had in mind would consist, giving a source 2d grid (the one retrieved from the netcdf) and a destination grid (the one that I need for my calculations) the second of which is coarser(have less points),in calculating a sort of weight matrix that contains the coefficients needed to compute each point of destination grid.

This will be helpful because I have a bunch of variables laying on the same grid, so if I calculate a weight matrix, the problem is reduced in a sequence of matrices multiplications.

basically what I need is a function f(lat_source,lon_source,lat_dest, lot_dest) that given the coords of the source/destination grids return a matrix of weights that maps the dest lat/lon onto the source one.

is there something already implemented? Or some reads that could point me in the right direction?

I know there are different interpolation methods that could be used, but, for the sake of the argument, let's suppose a linear interpolation.