MATLAB数字图像处理（一）局部直方图均衡化

利用MATLAB做图像的局部直方图均衡化：

clc
clear;
picture=imread('yun.png'); %读取图像
figure(1)
imshow(picture),title('原象'); %显示原始图象
picture=rgb2gray(picture);
%%
%建立模板
n=25;
a(1:n,1:n)=1;  %模板小于图像
image=wextend('2D','sym',picture,n);%%扩展
%%
[M,N]=size(image);%%扩展后的图像大小
W=double(image);
jubu_image=W;
for i=n+1:M-n      
    for j=n+1:N-n 
      m=W(i:i+(n-1),j:j+(n-1)).*a(1:n,1:n);  %取出x1中从(i,j)开始的n行n列元素与模板相乘
      k=histeq(uint8(m));%对局部进行直方图均衡
       lll=double(k);
        jubu_image(i,j)=lll(1,1); %将均衡后中心点的像素值赋给原图对应点的元素，未被赋值的元素取原值
    end
end
Picture=jubu_image(n+1:M-n,n+1:N-n);%%均衡化后，取原图像大小
Picture=uint8(Picture);
figure(2)
imshow(Picture),title('均衡后的图像'); %显示均衡后的图象