qq.m

非线性多分辨率信号的分解

clear;
I2=imread('alumgrns.tif');%读入图像
imshow(I2)%显示图像
I=double(I2);%将数据转换成double型，unit8型不能进行运算
[m,n]=size(I2);%原图的大小
p=ceil(m*0.5);
q=ceil(n*0.5);%新图的大小
for x=1:p-2
   for y=1:q-2
      I1(x,y)=0.25*(I(2*x,2*y)+I(2*x+1,2*y)+I(2*x+1,2*y+1)+I(2*x,2*y+1))+0.25*(I(2*x-1,2*y)+I(2*x-1,2*y+1)+I(2*x,2*y-1)+I(2*x+1,2*y-1)+I(2*x+2,2*y)+I(2*x+2,2*y+1)+I(2*x,2*y+2)+I(2*x+1,2*y+2))-0.5*(I(2*x-1,2*y-1)+I(2*x+2,2*y-1)+I(2*x+2,2*y+2)+I(2*x-1,2*y+2));
   end
end%实现算法
I3=uint8(I1);%将数据转换回unit8型
figure,imshow(I3)%显示新图