MATLAB 图像边缘提取 图像高斯低通滤波 DCT变换 程序

实验二边缘提取

作业要求:

1.参考“空间域图像增强”课的内容, 用Sobel边缘算子提取图像边缘实验程序如下：

RGB=imread('E:\图像信号处理\bmp格式图片\美女.bmp');

isrgb(RGB);

figure;

imshow(RGB);

I=rgb2gray(RGB);

imshow(I);

colorbar('horiz');

isgray(I);

ED=edge(I,'sobel',0.08);

imshow(ED)

运行结果如下：

实验四高斯低通滤波器

作业要求:

参考“频率域图像增强”课的内容, 采用高斯低通滤波器对图片进行滤波操作, 取半径为5, 15, 30, 80和230, 分别输出结果图片

实验程序如下：（d0=5）

I1=imread('E:\图像信号处理\bmp格式图片\荷花.bmp');

I=rgb2gray(I1);

s=fftshift(fft2(I));

[M,N]=size(s);

n=2;

d0=5;

n1=floor(M/2);

n2=floor(N/2);

for i=1:M

for j=1:N

d=sqrt((i-n1)^2+(j-n2)^2);

h=1*exp(-1/2*(d^2/d0^2));

s(i,j)=h*s(i,j);

end

end

s=ifftshift(s);

s=uint8(real(ifft2(s)));

figure;

imshow(s);

title('GLPF滤波（d0=5)');

实验六DCT变换

作业要求:

1.参考“图像压缩”课的内容, 并参考有关资料, 对输入图像进行离散余弦变换

2.通过离散余弦逆变换, 还原出原图片

实验程序如下：

RGB=imread('E:\图像信号处理\bmp格式图片\美女.bmp');

figure(1);

imshow(RGB);

I=rgb2gray(RGB); figure(2);

imshow(I);

J=dct2(I);

figure(3);

imshow(log(abs(J)),[]); colormap(jet(64)); colorbar;

J(abs(J)<11)=0;

K=idct2(J)/255; figure(4);

imshow(K)