引用
原理参考该博客
%.基于局部均方差增强图像
% blockMean 3*3邻域的均值
% blockVar 3*3邻域的方差
% PixelVal 3*3邻域中心像素点的灰度值
% GlobalVar 整个图像的方差
% blockMean + GlobalVar/blockVar*(PixelVal - blockMean)
clc;
clearvars;
close all;
img = imread('D:/picture/lena.jpg');
gray = rgb2gray(img);
GlobalMean = mean2(gray);%整幅图像的平均值
GlobalVar = std2(gray).^2;
BlockSize = 3;
len = floor(BlockSize/2);
grayPad = padarray(gray,[len,len], 'symmetric');
[grayPadRow, grayPadCol] = size(grayPad);
enhancePad = grayPad;
for ii = len+1:grayPadRow-len
for jj = len+1:grayPadCol - len
block = grayPad(ii-len:ii+len,jj-len:jj+len);
blockMean = mean2(block);
blockVar = std2(block).^2;
CG = GlobalVar/blockVar;
if CG>4
CG = 4;
end
enhancePad(ii, jj) = uint8( blockMean + CG*(grayPad(ii, jj) - blockMean));
end
end
enhanceGray = enhancePad(len+1:end-len, len+1:end-len);
figure,imshow(gray); title('original image');
figure,imshow(enhanceGray); title('enhanced image');