interpolate_smooth_2d.py

from scipy import interpolate
import numpy as np
# x_resolution = 0.38
# y_resolution = 0.38

# z = matrix_data
# xsize,ysize = z.shape

# x = np.arange(0,xsize+extra,1)*x_resolution
# y = np.arange(0,ysize+extra,1)*y_resolution

def interpolate_2d_data(x,y, z,x_resolution, y_resolution, new_N=None):
    extra = 0

    x_min, x_max = np.min(x),np.max(x)
    y_min, y_max = np.min(y),np.max(y)

    xx, yy = np.meshgrid(x, y)

    f = interpolate.interp2d(x, y, z, kind='cubic')
    # Interpolate the matrix data

    xnew = np.linspace(x_min, x_max, new_N)
    ynew = np.linspace(y_min, y_max, new_N)

    znew = f(xnew, ynew)
    return (xnew,ynew,znew)