图像工程概论实验报告——二值化直方图均衡化

图像工程概论

实验报告

课程报告1：

一、实验目的：

二、实验内容

把这幅图像分成同样大小的10幅人脸图片

然后分别对第一行5幅人脸图像的第3 第4 第5 第二行5幅人脸图像的第1 第5 进行如下处理：

1.进行大津法阈值分割的二值化处理

2.进行直方图均衡化处理

三、实验具体代码以及结果

（一）大津法阈值分割的二值化处理

1.实验代码

%图像分割部分

clear;

clc;

pic0=imread('G: \image\histogram matching.bmp');

figure(1),imshow(pic0);

wdiv=5;

hdiv=2;

[hd,wd,l]=size(pic0);

sub_wd=floor(wd/wdiv);

sub_hd=floor(hd/hdiv);

for row=1:2

for col=1:5

pic1=pic0((row-1)*sub_hd+1:row*sub_hd+1,(col-1)*sub_wd+1:col*sub_ wd+1);

figure(2),imshow(pic1);

imwrite(pic1,['G: \image\'num2str(row) '-'num2str(col) '.jpg'],'jpg'); %存储图像

end

end

%————————图像处理部分

cd 'G:\pic'

x1=imread('1-3.jpg'); %目标读取图像第一行第三幅人脸图像

x2=imread('1-4.jpg'); %第一行第四幅人脸图像

x3=imread('1-5.jpg'); %第一行第五幅人脸图像

x4=imread('2-1.jpg'); %第二行第一幅人脸图像

x5=imread('2-5.jpg'); %第二行第五幅人脸图像

%matlab 自动确定阈值的方法，大津法，类间方差

level1=graythresh(x1);

level2=graythresh(x2);

level3=graythresh(x3);

level4=graythresh(x4);

level5=graythresh(x5);

%用得到的阈值直接对图像进行二值化处理并显示

BW1=im2bw(x1,level1);

BW2=im2bw(x2,level2);

BW3=im2bw(x3,level3);

BW4=im2bw(x4,level4);

BW5=im2bw(x5,level5);

subplot(1,5,1),imshow(BW1);

subplot(1,5,2),imshow(BW2);

subplot(1,5,3),imshow(BW3);

subplot(1,5,4),imshow(BW4);

subplot(1,5,5),imshow(BW5);

2.实验结果显示

分割结果截图：

进行大津法阈值分割的二值化处理后截图：

（二）直方图均衡化处理

1.实验代码

%图像分割部分

clear;

clc;

pic0=imread('G: \image\histogram matching.bmp'); %figure(1),imshow(pic0);

wdiv=5;

hdiv=2;

[hd,wd,l]=size(pic0);

sub_wd=floor(wd/wdiv);

sub_hd=floor(hd/hdiv);

for row=1:2

for col=1:5

pic1=pic0((row-1)*sub_hd+1:row*sub_hd+1,(col-1)*sub_wd+1:col*sub_ wd+1);

figure(2),imshow(pic1);

imwrite(pic1,['G: \image\'num2str(row) '-'num2str(col) '.jpg'],'jpg'); %存储图像

end

end

%图像处理部分

cd 'G:\pic'

x1=imread('1-3.jpg'); %目标读取图像第一行第三幅人脸图像

x2=imread('1-4.jpg'); %第一行第四幅人脸图像

x3=imread('1-5.jpg'); %第一行第五幅人脸图像

x4=imread('2-1.jpg'); %第二行第一幅人脸图像

x5=imread('2-5.jpg'); %第二行第五幅人脸图像

%直方图均衡化处理

%调用直方图均衡化函数 histeq（）

%均衡化处理前后的灰度级直方图分布，采用两行十列显示

subplot(2,5,1),imhist(rgb2gray(x1)); %第一行第三幅人脸图像subplot(2,5,6),imhist(histeq(rgb2gray(x1)));

subplot(2,5,2),imhist(rgb2gray(x2)); %第一行第四幅人脸图像subplot(2,5,7),imhist(histeq(rgb2gray(x2)));

subplot(2,5,3),imhist(rgb2gray(x3)); %第一行第五幅人脸图像subplot(2,5,8),imhist(histeq(rgb2gray(x3)));

subplot(2,5,4),imhist(rgb2gray(x4)); %第二行第一幅人脸图像subplot(2,5,9),imhist(histeq(rgb2gray(x4)));

subplot(2,5,5),imhist(rgb2gray(x5)); %第二行第五幅人脸图像subplot(2,5,10),imhist(histeq(rgb2gray(x5)));

%显示均衡化处理后的图像

figure(11),imshow(histeq(rgb2gray(x1)));

figure(12),imshow(histeq(rgb2gray(x2)));

figure(13),imshow(histeq(rgb2gray(x3)));

figure(14),imshow(histeq(rgb2gray(x4)));

figure(15),imshow(histeq(rgb2gray(x5)));

2.实验结果显示

①直方图