[MATLAB]数字图像处理—直方图均衡、直方图匹配、局部直方图增强、直方图分割

数字图像处理课程的第三次作业，包括图像灰度直方图的绘制、基于直方图的各种空域处理。参考课本是冈萨雷斯的《数字图像处理》

clear;
clc;
file=dir（‘.\*.bmp‘）;
figure（‘Numbertitle‘‘off‘‘Name‘‘图像直方图‘）;
for i=1:length（file）
    [img{i}colormap{i}]=imread（file（i）.name）;  %读出每幅图的索引和调色板
    img{i}=ind2gray（img{i}colormap{i}）;    %将索引图转为灰度图
    %img{i}=im2double（img{i}）;
    subplot（44i）;
    imhist（img{i}）;
    title（file（i）.name）;
end
figure（‘Numbertitle‘‘off‘‘Name‘‘原始图像‘）;
for i=1:length（file）
    subplot（44i）;
    imshow（img{i}）;
    title（file（i）.name）;
end
figure（‘Numbertitle‘‘off‘‘Name‘‘直方图均衡化后的图像‘）;
for i=1:length（file）
    img_histeq{i}=histeq（img{i}）;
    subplot（44i）;
    imshow（img_histeq{i}）;
    title（file（i）.name）;
end
figure（‘Numbertitle‘‘off‘‘Name‘‘均衡化后的直方图‘）;
for i=1:length（file）
    subplot（44i）;
    imhist（img_histeq{i}）;
    title（file（i）.name）;
end
%%
city_elain=histmatch（img{1}imhist（img{4}））;
figure（‘Numbertitle‘‘off‘‘Name‘‘直方图匹配1‘）
subplot（231）;imshow（img{1}）;title（‘citywall‘）;
subplot（232）;imshow（img{4}）;title（‘elain‘）;
subplot（233）;imshow（city_elain）;title（‘city\_elain‘）;
subplot（234）;imhist（img{1}）;
subplot（235）;imhist（img{4}）;
subplot（236）;imhist（city_elain）;

lena_woman=histmatch（img{8}imhist（img{12}））;
figure（‘Numbertitle‘‘off‘‘Name‘‘直方图匹配2‘）
subplot（231）;imshow（img{8}）;title（‘lena‘）;
subplot（232）;imshow（img{12}）;title（‘woman‘）;
subplot（233）;imshow（lena_woman）;title（‘lena\_woman‘）;
subplot（234）;imhist（img{8}）;
subplot（235）;imhist（img{12}）;
subplot（236）;imhist（lena_woman）;
%%
function img_out=histmatch（img_inhist）
[MN]=size（img_in）;
hist_img=imhist（img_in）;
%求T1：
a=zeros（2561）;
T1=zeros（2561‘uint8‘）;
a（1）=hist_img（1）;
for i=2:256
    a（i）=a（i-1）+hist_img（i）;
    T1（i）=uint8（round（a（i）*255/（M*N）））;
end
%求T2：
a=zeros（2561）;
T2=zeros（2561‘uint8‘）;
a（1）=hist（1）;
for i=2:256
    a（i）=a（i-1）+hist（i）;
    T2（i）=uint8（round（a（i）*255/sum（hist）））;
end
%求T3：
T3=zeros（2561）;
a=zeros（2561）;
for i=1:256
    for j=1:256
        a（j）=abs（T1（i）-T2（j））;
    end
    [~p]=min（a）;
    T3（i）=uint8（p-1）;
end
%outimage:
for i=1:M
    for j=1:N
        img_out（ij）=T3（uint32（img_in（ij））+1）;
    end
end
img_out=uint8（img_out）;
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件         519  2019-03-21 19:10  README.txt
     文件        2257  2019-03-15 01:02  task3_1.m
     文件        1705  2019-03-16 21:30  task3_2.m
     文件        1005  2019-03-15 17:53  task3_3.m