Matlab——图像缩放（插值法）

实验内容：

    用双线性内插法实现位深度为8的灰度图像的缩放。

 

思路：

     输入原图像以及缩放后图像的像素要求（宽度*高度），处理后输出新图像。

    我是用matlab来实现scale(input_img,scale_size)函数的，输入图像路径以及要求实现的宽度scale_w和高度值scale_h即可。通过matlab把原图像转化为img矩阵，这样矩阵img(a, b)元素的值就是原图像高度为a，宽度为b处的像素值。

   设处理前后图像的宽度和高度分别为w、h，scale_w、scale_h。

    如上图所示，row对应的是高度，设为x轴；col对应的是宽度，设为y轴；value对应的是灰度值，设为z轴。

 

    通过缩放比例，可以求得新图像B在(i,j)处对应原图像的点为(x,y)（其中x= i*h/scale_h，y= j*w/scale_w）。

但是这两个值可能为浮点数，而像素中的位置是整数，所以这个点在原图像中对应的可能是“虚”点。所以，我们需

要根据找该“虚”点周围的四个点来进行双线性插值得到新图像的灰度值。我们应该取的是(x,y)邻近的四个像素

(x0,y0) 、(x0,y1) 、(x1,y0) 、(x1,y1)，它们对应的灰度值为f(x0,y0) 、f(x0,y1)、 f(x1,y0)、 f(x1,y1)。

  

    先对x方向进行插值，再对y方向进行插值，计算如下：

    

    

    

 

    反过来，先对y方向：

    

    

    

 

    显然，邻近像素的x1-x0 = y1-y0 = 1
    整理式子得到：

 

    这里u、v分别是x、y的小数部分，即u = x - x0; v = y - y0

       要注意的是：上述公式img(a,b)中的a和b不能小于1，所以当x和y小于1时需要将其重新赋值为1。

 

代码：

function output_img = scale(input_img, scale_size)
%Input - input_img is a two-dimensional matrices storing image
%      - scale_size is a tuple of [width, height] defining the spatial resolution of output
%Output - output_img is the same as input_img
 
img = imread(input_img); %读取输入图片的数据
[h,w] = size(img); %获取行和列，即原图的高度和宽度
 
scale_w = scale_size(1); %根据输入获得缩放后的新宽度
scale_h = scale_size(2); %根据输入获得缩放后的新高度
output_img = zeros(scale_h, scale_w); %初始化
 
for i = 1 : scale_h         %缩放后的图像的（i,j）位置对应原图的（x,y）
    for j = 1 : scale_w
        x = i * h / scale_h;
        y = j * w / scale_w;
        u = x - floor(x);
        v = y - floor(y); %取小数部分
        
        if x < 1           %边界处理
            x = 1;
        end
        
        if y < 1
            y = 1;
        end
        
 
        %用原图的四个真实像素点来双线性插值获得“虚”像素的像素值
        output_img(i, j) = img(floor(x), floor(y)) * (1-u) * (1-v) + ...
                               img(floor(x), ceil(y)) * (1-u) * v + ...
                               img(ceil(x), floor(y)) * u * (1-v) + ...
                               img(ceil(x), ceil(y)) * u * v;
    end
end
 
imwrite(uint8(output_img), '../output_img.png'); %保存处理后的图像
imshow(input_img); %显示原图
figure,imshow(uint8(output_img)) %显示处理后的图像

运行时若图片在代码上一级目录，则输入scale('../a.png', [450,300]);

 

效果图：

原图为384*256：

缩小成192*128：  

放大成450*300：

缩放成500*200：

 

//

 转载：https://blog.csdn.net/Stella_Chan/article/details/78025020?locationNum=5&fps=1