数字图像四倍三倍二倍图像处理

A=imread('00102.bmp'); %读取原图片

imshow(A);

fn='00102.bmp';

p=imread(fn);

g=rgb2gray(p); %转化为灰度

subplot(211);imshow(p);

subplot(212);imshow(g);

t=graythresh(a);%取阈值

b=im2bw(a,t);%转化为二值图

figure;imshow(b);title('二值图');

[line,row]=size(a);%读取图像像素

L=1;

R=1;

%4倍减采样

for i=1:2:line;

for j=1:2:row;

b1(L,R)=a(i,j);

R=R+1;%取原图像i列下一行的元素赋给新图像的对应位置

end

L=L+1;%换列

R=1;%从换列后的列里的第一个元素开始取元素

end

figure;imshow(b1);title('4倍采样图');

%16倍减采样

K=1;

M=1;

for i=1:4:line;

for j=1:4:row;

b2(K,M)=a(i,j);

M=M+1;%取原图像i列下一行的元素赋给新图像的对应位置

end

K=K+1;%换列

M=1;%从换列后的列里的第一个元素开始取元素

end

figure;imshow(b2);title('16倍采样图');

%量化成128级，64级，32级

c=(0.5)*a;c1=floor(c);

d=(0.25)*a;d1=floor(c);

e=(0.125)*a;e1=floor(e);

figure

subplot(1,3,1),imshow(c1,[0,127]);title('128级灰度图像');

subplot(1,3,2),imshow(d1,[0,63]);title('64级灰度图像');

subplot(1,3,3),imshow(e1,[0,31]);title('32级灰度图像');