I am trying to write a python script with help of GDAL to a red point shapefile and considering a field value as z vales where I need interpolation to get raster with interpolated values. I am using scipy grid data for interpolation.
My interpolation code is:
zi = scipy.griddata((x, y), z, (xi, yi), method='cubic')
It shows interpolation inside points. But, I want to interpolate my points to cover the whole rectangular space.
Does anyone have idea how it works?