数字图像处理 - Matlab - 实验笔记二

第三章 灰度变换与空间滤波

1. 仿真实现一种灰度变换增强图像

程序代码：

% 灰度增强
clear
X=imread('pic.jpg');
% 把原图转为灰度图
X0=double(rgb2gray(X));
% 获取图像大小
[width,height]=size(X0);
% 显示原灰度图
subplot(2,3,1)
imshow(uint8(X0));
title('原灰度图')
% 预分配内存
X1=zeros(width,height);
X2=zeros(width,height);
X3=zeros(width,height);
X4=zeros(width,height);
% 线性变换
a=0.5;
b=0;
for i=1:width
    for j=1:height
        X1(i,j)=X0(i,j)*a+b;
    end
end
subplot(2,3,2)
imshow(uint8(X1));
title('线性变换')
% 对数形式变换
c1=35;
for i=1:width
    for j=1:height
        X2(i,j)=c1*log(X0(i,j)+1);
    end
end
subplot(2,3,3)
imshow(uint8(X2));
title('对数形式变换')
% 指数形式变换
c2=0.1;
r=1.5;
for i=1:width
    for j=1:height
        X3(i,j)=c2*(X0(i,j)).^r;
    end
end
subplot(2,3,5)
imshow(uint8(X3));
title('指数形式变换')
% 分段线性变换
k1=0.5;
b1=0;
k2=167.5/80;
b2=-191.25;
k3=0.5;
b3=127.5;
for i=1:width
    for j=1:height
        if X0(i,j)<=120
            X4(i,j)=k1*X0(i,j)+b1;
        elseif X0(i,j)<=200
            X4(i,j)=k2*X0(i,j)+b2;
        else
            X4(i,j)=k3*X0(i,j)+b3;
        end
    end
end
subplot(2,3,6)
imshow(uint8(X4));
title('分段线性变换')

输出结果：

2. 仿真实现一幅图像的直方图

程序代码：

function hist = gethistogram(X)
%ge