MATLAB图像处理之图像增强实战

pout = imread(‘pout.tif’);%读取pout灰度图像
tire = imread(‘tire.tif’);%读取tire灰度图像
[X map] = imread(‘shadow.tif’);%读取索引图像
shadow = ind2rgb(X,map); % 转化为真彩色图像
width = 210;%转化为统一宽度，以便进行对比
images = {pout, tire, shadow};
for k = 1:3
dim = size(images{k});
images{k} = imresize(images{k},…
[widthdim(1)/dim(2) width],‘bicubic’);
end
pout = images{1};
tire = images{2};
shadow = images{3};
pout_imadjust = imadjust(pout);%使用imadjust函数对图像进行增强
pout_histeq = histeq(pout); %使用histeq函数对图像进行增强
pout_adapthisteq = adapthisteq(pout); %使用adapthisteq函数进行增强
figure; subplot(221);
imshow(pout);%显示原图像
title(‘Original’); subplot(222);
imshow(pout_imadjust); %显示imadjust函数增强的图像
title(‘Imadjust’); subplot(223);
imshow(pout_histeq); %显示histeq函数增强的图像
title(‘Histeq’); subplot(224);
imshow(pout_adapthisteq); %显示adapthisteq函数增强的图像
title(‘Adapthisteq’);
tire_imadjust = imadjust(tire); %使用imadjust函数对图像进行增强
tire_histeq = histeq(tire); %使用histeq函数对图像进行增强
tire_adapthisteq = adapthisteq(tire); %使用adapthisteq函数进行增强
figure, subplot(221)
imshow(tire);%显示原图像
title(‘Original’); subplot(222);
imshow(tire_imadjust);%显示imadjust函数增强的图像
title(‘Imadjust’); subplot(223);
imshow(tire_histeq);%显示histeq函数增强的图像
title(‘Histeq’); subplot(224),
imshow(tire_adapthisteq);%显示adapthisteq函数增强的图像
title(‘Adapthisteq’); figure; subplot(121);
imhist(pout), %显示pout图像的直方图
title(‘pout.tif’); subplot(122),
imhist(tire), %显示tire图像的直方图
title(‘tire.tif’);
srgb2lab = makecform(‘srgb2lab’);%rgb彩色空间变为Lab空间结构
lab2srgb = makecform(‘lab2srgb’); %Lab变为rgb空间结构
shadow_lab = applycform(shadow, srgb2lab);%图像变换到Lab空间
max_luminosity = 100;%规定最大的光照值
L = shadow_lab(:,:,1)/max_luminosity;%归一化
shadow_imadjust = shadow_lab;
shadow_imadjust(:,:,1) = imadjust(L)…
max_luminosity; %使用imadjust函数进行增强
shadow_imadjust = applycform(shadow_imadjust, …
lab2srgb);%变换到rgb彩色空间
shadow_histeq = shadow_lab;
shadow_histeq(:,:,1) = histeq(L)…
max_luminosity; %使用histeq函数进行增强
shadow_histeq = applycform(shadow_histeq, …
lab2srgb);%变换到rgb彩色空间
shadow_adapthisteq = shadow_lab;
shadow_adapthisteq(:,:,1) = adapthisteq(L)*…
max_luminosity; %使用adapthisteq函数进行增强
shadow_adapthisteq = applycform(shadow_adapthisteq, …
lab2srgb);%变换到rgb彩色空间
figure; subplot(221)
imshow(shadow);%显示原图像shadow
title(‘Original’); subplot(222),
imshow(shadow_imadjust); %显示imadjust函数增强的图像
title(‘Imadjust’); subplot(223),
imshow(shadow_histeq); %显示histeq函数增强的图像
title(‘Histeq’); subplot(224),
imshow(shadow_adapthisteq); %显示adapthisteq函数增强的图像
title(‘Adapthisteq’);
mypout = myhisteq(pout);%使用myhisteq函数对pout图像进行增强
mytire = myhisteq(tire); %使用myhisteq函数对tire图像进行增强
figure,
subplot(121)
imshow(mypout);%显示myhisteq函数对pout增强后的图像
subplot(122)
imshow(mytire);%显示myhisteq函数对tire增强后的图像