中南民族大学多媒体实验二.doc

年级

班号

计科2003班

?学号

?

专业

?姓名

实验名称

?图像处理

实

验

目

的

或

要

求

实验内容：用代码实现对于图像的灰度处理，椒盐噪声处理，中值滤波处理，sobel算子边缘检测，图片的腐蚀与膨胀，人脸识别及马赛克处理

实验平台：matlab2021a

实

验

过

程

clear;

clc;

closeall;

lib_face(E:\as2.jpeg);

function[]=lib_face(file_name)

RGB_data=imread(file_name);

figure;

subplot(1,2,1);

imshow(RGB_data);

title(原始图像);

%将源图像三个通道值取出

R_data=RGB_data(:,:,1);

G_data=RGB_data(:,:,2);

B_data=RGB_data(:,:,3);

[ROW,COL]=size(R_data);

Y_data=zeros(ROW,COL);

Gy_data=RGB_data;

%RGB图像转换成Ycbcr图像，使用Y分量去显示图像，来实现彩色图像转灰度图

forr=1:ROW

forc=1:COL

Y_data(r,c)=0.299*R_data(r,c)+0.587*G_data(r,c)+0.114*B_data(r,c);

end

end

Gy_data(:,:,1)=Y_data;

Gy_data(:,:,2)=Y_data;

Gy_data(:,:,3)=Y_data;

%灰度

subplot(1,2,2)

imshow(Gy_data);

title(灰度图像);

%0.02为噪声密度，噪声值的图像区域的百分比

imgn1=imnoise(Gy_data,saltpepper,0.02);

%加入椒盐噪声后的图像

figure;

subplot(1,2,1)

imshow(imgn1);%加入椒盐噪声的图像

title(椒盐噪声);

%中值滤波：

%(1)对窗内的每行像素按降序排序，得到最大值、中间值和最小值。

%(2)把三行的最小值即第三列相比较，取其中的最大值。

%(3)把三行的最大值即第一列相比较，取其中的最小值。

%(4)把三行的中间值即第二列相比较，再取一次中间值。

%(5)把前面的到的三个值再做一次排序，获得的中值即该窗口的中值。

Mid_Img=Gy_data;%转存一次

forr=2:ROW-1

forc=2:COL-1

median3x3=[imgn1(r-1,c-1)imgn1(r-1,c)imgn1(r-1,c+1)

imgn1(r,c-1)imgn1(r,c)imgn1(r,c+1)

imgn1(r+1,c-1)imgn1(r+1,c)imgn1(r+1,c+1)];

sort1=sort(median3x3,2,descend);

sort2=sort([sort1(1),sort1(4),sort1(7)],descend);

sort3=sort([sort1(2),sort1(5),sort1(8)],descend);

sort4=sort([sort1(3),sort1(6),sort1(9)],descend);

mid_num=sort([sort2(3),sort3(2),sort4(1)],descend);

Mid_Img(r,c)=mid_num(2);

end

end

subplot(1,2,2)

imshow(Mid_Img);

title(中值滤波);

%soble算子

%(1)计算Gx与Gy与模板每行的乘积。

%(2)两个3x3矩阵的卷积即将每一行每一列对应相乘然后相加。

%(3)求得3*3模板运算后的Gx、Gy。

%(4)求Gx^2+Gy^2的平方根或者直接对Gx和Gy取绝对值后求和。

%(5)设置一个阈值，运算后的像素值大于该阈值输出为255，小于该阈值输出为0。

%在对图像数据进行运算时最好将图像数据转化成double类型的，避免精度损失。