【数字图像处理】对一个图像进行非线性动态调整和直方图均衡

选取“舅妈真是的”图片如下：

利用imread函数导入matlab，利用rgb2gray转灰度矩阵，利用histeq函数进行处理，并利用subplot与imshow函数展示。非线性动态拉伸参考公式也很容易得到（代码里处理的是不同的两种图片）

源代码如下：

%%读照片
A=imread('灰度图片1-直方图均衡化.jpg');
B=imread('灰度图片2-非线性动态拉伸.jpg');
%%转灰度
gray_A = rgb2gray(A);
gray_B = rgb2gray(B);
%%处理照片 进行动态范围拉伸处理 注意类型转换
[m,n]=size(gray_B);
B_boost=zeros(m,n);
gray_B_double = double(gray_B);
for i=1:m
    for j=1:n
       B_boost(i,j) = log10(1 + gray_B_double(i,j)) ;
    end
end
%%确定c值
B_max = max(gray_B(:));
B_bmax = max(B_boost(:));
c=B_max/B_bmax;
%%完成动态处理，进行类型转换【转换为 8 位无符号整数（取值范围为 0 到 255）】
B_boost=double(c)*B_boost;
B_boost=uint8(B_boost);
A_boost = histeq(gray_A);
%%画图展示
figure;
subplot(2,2,1),imshow(gray_A);
title('原图');
subplot(2,2,2),imshow(A_boost);
title('对比图（直方图均衡化）'); 
subplot(2,2,3),imshow(gray_B);
title('原图');
subplot(2,2,4),imshow(B_boost);
title('对比图（非线性动态范围拉伸）'); 
figure;
subplot(1,2,1),imhist(gray_A);
title('灰度图片1未均衡化前灰度直方图')
xlabel('灰度级别');
ylabel('像素个数');
subplot(1,2,2),imhist(A_boost);
title('灰度图片1均衡化后灰度直方图')
xlabel('灰度级别');
ylabel('像素个数');