matlab双立方插值法_双三次插值(bicubic interpolation)原理及MATLAB源码实现

本文介绍了双三次插值的原理,并通过MATLAB代码实现了图像的双三次插值放大。实验中,首先将图像转换为灰度图像,然后使用双线性和双三次插值方法进行放大。结果显示,双三次插值在保持图像细节方面表现出更好的效果。同时,通过傅立叶变换比较了原图像与双三次插值图像的傅立叶幅度谱。
摘要由CSDN通过智能技术生成

%

双三次插值具体实现

clc,clear;

fff=imread('E:\Documents\BUPT\DIP\

图片

\lena.bmp');

ff = rgb2gray(fff);%

转化为灰度图像

[mm,nn]=size(ff);

%

将图像隔行隔列抽取元素,得到缩小的图像

f

m=mm/2;n=nn/2;

f = zeros(m,n);

for i=1:m

for j=1:n

f(i,j)=ff(2*i,2*j);

end

end

k=5;

%

设置放大倍数

bijiao1 = imresize(f,k,'bilinear');%

双线性插值结果比较

bijiao = uint8(bijiao1);

a=f(1,:);c=f(m,:);

%

将待插值图像矩阵前后各扩展两行两列

,

共扩展四行四列

b=[f(1,1),f(1,1),f(:,1)',f(m,1),f(m,1)];d=[f(1,n),f(1,n),f(:,n)',f(m,n),f(m,n)];

a1=[a;a;f;c;c];

b1=[b;b;a1';d;d];

ffff=b1';f1=double(ffff);

g1 = zeros(k*m,k*n);

for i=1:k*m

%

利用双三次插值公式对新图象所有像素赋值

u=rem(i,k)/k; i1=floor(i/k)+2;

A=[sw(1+u) sw(u) sw(1-u) sw(2-u)];

for j=1:k*n

v=rem(j,k)/k;j1=floor(j/k)+2;

C=[sw(1+v);sw(v);sw(1-v);sw(2-v)];

B=[f1(i1-1,j1-1) f1(i1-1,j1) f1(i1-1,j1+1) f1(i1-1,j1+2)

f1(i1,j1-1)

f1(i1,j1)

f1(i1,j1+1)

f1(i1,j1+2)

f1(i1+1,j1-1)

f1(i1+1,j1) f1(i1+1,j1+1) f1(i1+1,j1+2)

f1(i1+2,j1-1) f1(i1+2,j1) f1(i1+2,j1+1) f1(i1+2,j1+2)];

g1(i,j)=(A*B*C);

end

end

g=uint8(g1);

imshow(uint8(f)); title('

缩小的图像

');

%

显示缩小的图像

figure,imshow(ff);title('

原图

');

%

显示原图像

figure,imshow(g);title('

双三次插值放大的图像

');

%

显示插值后的图像

figure,imshow(bijiao);title('

双线性插值放大结果

');

%

显示插值后的图像

mse=0;ff=double(ff);g=double(g);

ff2=fftshift(fft2(ff));

%

计算原图像和插值图像的傅立叶幅度谱

g2=fftshift(fft2(g));

figure,subplot(1,2,1),imshow(log(abs(ff2)),[8,10]);title('

原图像的傅立叶幅度谱

');

subplot(1,2,2),imshow(log(abs(g2)),[8,10]);title('

双三次插值图像的傅立叶幅度谱

');

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值