I am developing in Python and using GDAL from OSGEO to manipulate and interact with rasters and shapefiles.
I want to take a shapefile that has point features and interpolate it into a surface raster. Right now I am using the method 'RasterizeLayer' which burns a value from the point feature into the raster (which is set with all nodata values) but leaves all untouched pixels as a 'nodata' value. I am therefore left with a checkerboard type raster.
What I have after using RasterizeLayer:
What I want for a final product:
I believe the function I am looking for is known as 'Spline_sa()' from the arcgisscripting import.
Does GDAL have a similar function, or is there a different method to get my desired output?