用双线性插值放大图片img matlab代码怎么写

以下是一段用双线性插值放大图片的 MATLAB 代码：

function [new_img] = bilinear_interpolation(img, scale_factor)
% img为输入的图片，scale_factor为图片的放大倍数
% 将输入的图片img用双线性插值方法放大
[h, w, n_channels] = size(img);
new_h = round(h * scale_factor);
new_w = round(w * scale_factor);

% 初始化输出图片
new_img = zeros(new_h, new_w, n_channels);

% 缩放因子
sh = h/new_h;
sw = w/new_w;

% 计算在输出图片中每个位置对应的在输入图片中的位置
[xx, yy] = meshgrid(1:new_w, 1:new_h);
v = [xx(:), yy(:)];

% 映射到输入图片中的位置
uv = v.*[sw, sh];

% 找到每个输出像素在输入图像中最近的四个像素
x1 = floor(uv(:,1));
x2 = x1 + 1;
y1 = floor(uv(:,2));
y2 = y1 + 1;

% 边界处理
x1(x1 < 1) = 1;
y1(y1 < 1) = 1;
x2(x2 > w) = w;
y2(y2 > h) = h;

% 计算每个输出像素的值
for c = 1:n_channels
    channel = double(img(:,:,c));
    q11 = channel(sub2ind(size(channel), y1, x1));
    q12 = channel(sub2ind(size(channel), y1, x2));
    q21 = channel(sub2ind(size(channel), y2, x1));
    q22 = channel(sub2ind(size(channel), y2, x2));
    x = uv(:,1) - x1;
    y = uv(:,2) - y1;
    f = q11.*(1-x).*(1-y) + q21.*x.*(1-y) + q12.*(1-x).*y + q22.*x.*y;
    new_img(:,:,c) = reshape(f, [new_h, new_w]);
end

% 转换图片类型以保持一致性
new_img = uint8(new_img);

end

请注意，此代码仅供参考，可能存在错误和不完整性，实际使用时需要进行调整和测试。