matlab图像处理基础知识0（双线性插值matlab实现--调整水平和垂直放大倍数）

需求说明：图像处理FPGA实现前期算法准备

当前算法：已通过matlab验证成功

function scaler_bilinear()

%-----------------------------0:配置输入输出--------------------------------

m=2; %放大或缩小的高度

n=2; %放大或缩小的宽度

I=imread('F:\book\Digital image processing and machine vision\640480\plane.jpg');

[height,width,l]=size(I); %输入分辨率 640*480

O=zeros(height*m,width*n,l); %输出分辨率 1920*1080

%-----------------------------1:开始输出计算--------------------------------

rot=[m,0,0;0,n,0;0,0,1]; %变换矩阵

for i=1:height*m

for j=1:width*n

pix=[i,j,1]/rot;

float_Y=pix(1)-floor(pix(1)); %水平宽度小数部分 floor 向下取整

float_X=pix(2)-floor(pix(2)); %垂直高度小数部分 floor 向下取整

%-----------------------------2:边界处理--------------------------------

if pix(1)<1 %边界处理

pix(1)=1;

end

if pix(1)>height %边界处理

pix(2)=height;

end

if pix(2)<1 %边界处理

pix(2)=1;

end

if pix(2)>width %边界处理

pix(2)=width;

end

%-----------------------------3:原坐标四个相邻点----------------------------

%四个相邻的点

pix_up_left = [floor(pix(1)) , floor(pix(2))];

pix_up_right = [floor(pix(1)) , ceil(pix(2))]; %floor 向下取整 ceil 向上取整

pix_down_left = [ceil(pix(1)) , floor(pix(2))];

pix_down_right= [ceil(pix(1)) , ceil(pix(2))];

%-----------------------------4:权重计算-----------------------------------

%计算临近四个点的权重

value_up_left = (1-float_X) * (1-float_Y);

value_up_right = float_X * (1-float_Y);

value_down_left = (1-float_X) * float_Y;

value_down_right= float_X * float_Y;

%按权重进行双线性插值

% O(i,j,:) = value_up_left * I(pix_up_left(1) ,pix_up_left(2),:) +...

% value_up_right * I(pix_up_right(1) ,pix_up_right(2),:) +...

% value_down_left * I(pix_down_left(1) ,pix_down_left(2),:) +...

% value_down_right* I(pix_down_right(1),pix_down_right(2),:);

%-----------------------------5:计算输出图像数值----------------------------

O(i,j,1) = value_up_left * I(pix_up_left(1), pix_up_left(2) ,1) +...

value_up_right * I(pix_up_right(1), pix_up_right(2) ,1) +...

value_down_left * I(pix_down_left(1), pix_down_left(2) ,1) +...

value_down_right* I(pix_down_right(1),pix_down_right(2),1);

O(i,j,2) = value_up_left * I(pix_up_left(1), pix_up_left(2) ,2) +...

value_up_right * I(pix_up_right(1), pix_up_right(2) ,2) +...

value_down_left * I(pix_down_left(1), pix_down_left(2) ,2) +...

value_down_right* I(pix_down_right(1),pix_down_right(2),2);

O(i,j,3) = value_up_left * I(pix_up_left(1), pix_up_left(2) ,3) +...

value_up_right * I(pix_up_right(1), pix_up_right(2) ,3) +...

value_down_left * I(pix_down_left(1), pix_down_left(2) ,3) +...

value_down_right* I(pix_down_right(1),pix_down_right(2),3);

end

end

% figure,imshow(uOt8(O))

%-----------------------------6:输出显示------------------------------------

figure,imshow(I); title('输入图像640*480');

figure,imshow(uint8(O));title('输出图像1920*1080');

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

设计来自：时间的诗