[mtalab图像处理] 灰度的线性变换

 clc;
 clear all;
I= imread('cameraman.tif');        % 读入原图像 
 
I = im2double(I);           % 转换数据类型为double 
[M,N] = size(I);            % 计算图像面积 
 
figure(1);              % 打开新窗口 
imshow(I);              % 显示原图像 
title('原图像'); 
 
figure(2);                  % 打开新窗口 
[H,x] = imhist(I, 64);      % 计算64个小区间的灰度直方图 
stem(x, (H/M/N), '.');      % 显示原图像的直方图  直方图做了归一化（H/M/N）
title('直方图原图像'); 
 
% 增加对比度 
Fa = 2; Fb = -55; 
O = Fa .* I + Fb/255;   %Fa>1时，输出图像的对比度将增大；Fa<1时，输出图像对比度将减小。
                       %Fa=1且Fb非零时，所有像素的灰度值上移或下移，使整个图像更暗或
                       %更亮。Fa<0，暗区变亮，亮区变暗。 实际处理图像是Fa，Fb是通过计算出来的
 
figure(3); 
subplot(2,2,1); 
imshow(O); 
title('Fa = 2 Fb = -55 增加对比度图像'); 
 
figure(4); 
subplot(2,2,1); 
[H,x] = imhist(O, 64); 
stem(x, (H/M/N), '.'); 
title('Fa = 2 Fb = -55 增加对比度直方图'); 
 
% 减小对比度 
Fa = 0.5; Fb = -55; 
O = Fa .* I + Fb/255; 
 
figure(3); 
subplot(2,2,2); 
imshow(O); 
title('Fa = 0.5 Fb = -55 减小对比度图像'); 
 
figure(4); 
subplot(2,2,2); 
[H,x] = imhist(O, 64); 
stem(x, (H/M/N), '.'); 
title('Fa = 0.5 Fb = -55 减小对比度直方图'); 
 
% 线性增加亮度 
Fa = 1; Fb = 55; 
O = Fa .* I + Fb/255; 
 
figure(3); 
subplot(2,2,3); 
imshow(O); 
title('Fa = 1 Fb = 55 线性平移增加亮度图像'); 
 
figure(4); 
subplot(2,2,3); 
[H,x] = imhist(O, 64); 
stem(x, (H/M/N), '.'); 
title('Fa = 1 Fb = 55 线性平移增加亮度直方图'); 
 
% 反相显示   黑变白 白变黑
Fa = -1; Fb = 255; 
O = Fa .* I + Fb/255; 
 
figure(3); 
subplot(2,2,4); 
imshow(O); 
title('Fa = -1 Fb = 255 反相显示图像'); 
 
figure(4); 
subplot(2,2,4); 
[H,x] = imhist(O, 64); 
stem(x, (H/M/N), '.'); 
title('Fa = -1 Fb = 255 反相显示直方图');