python双三次插值算法代码

以下是 Python 双三次插值算法的代码：

import numpy as np
from scipy import interpolate

def bicubic_interpolation(image, scale):
    height, width, channels = image.shape
    new_height = int(height * scale)
    new_width = int(width * scale)
    new_image = np.zeros((new_height, new_width, channels))

    for c in range(channels):
        for i in range(new_height):
            for j in range(new_width):
                x = i / scale
                y = j / scale
                if x < 1 or x > height - 2 or y < 1 or y > width - 2:
                    new_image[i, j, c] = 0
                else:
                    x1 = int(x)
                    y1 = int(y)
                    dx = x - x1
                    dy = y - y1
                    points = np.zeros((4, 4))
                    for m in range(-1, 3):
                        for n in range(-1, 3):
                            points[m+1, n+1] = image[x1+m, y1+n, c]
                    f = interpolate.interp2d(np.arange(-1, 3), np.arange(-1, 3), points, kind='cubic')
                    new_image[i, j, c] = f(dx, dy)

    return np.clip(new_image, 0, 255).astype(np.uint8)

希望对你有帮助！