I am looking for some hints here. There is a shape file with point data of rainfall sampled at a regular grid, and there is a shape file with point data of soil attributes sampled at a regular grip but with smaller pixel size (finer grid). For each point on the finer grid I need to obtain an interpolated value of rainfall data. Basically I was thinking of converting rainfall data into a raster using gdal_grid and interpolating it using idw method.
Are there any Gdal tools which will enable to convert the raster back to a shape file with point data matching locations of the soil properties (given in shape file)?