中国史之【褒姒得宠】：
褒姒(baosi），西周君主周幽王第二任王后，太子伯服生母，幽王甚为宠爱，有“烽火戏诸侯”之说。前771年，申侯联合犬戎攻打周幽王，周幽王被杀于骊山之下，褒姒被犬戎掳走，从此下落不明，西周灭亡。
——来源：全历史APP
【路漫漫其修远兮，吾将上下而求索】
今天介绍图像的灰度变换实践，包括对数变换和分段线性变换。本内容参考自《实用MATLAB图像和视频处理》第8章。
- 对于matlab2021的安装，可以参考公众号【CV之道】的文章：【matlab安装】手把手图文并茂安装matlab2021（win10版）
1、对数变换
对数变换和反对数变换都是非线性变换，分别用于压缩或扩展一幅图像中像素值的动态范围。对数变换的数学描述如下：
其中，r是原始像素灰度，s是结果像素值，c是一个常数，控制范围。使用（1+r）是为了确保输入像素值大于1，log输出为正数。
书本中，举了一个对数的例子。使用对数函数来改善和显示傅里叶变换的结果。如下，左图是大米图像的傅里叶频谱（仅振幅），其矩阵范围是【0, 2.8591x10的四次方】，当将他显示在一个线性放缩的8比特系统上，仅看到中心的亮点而很难看到其他东西。
我们用matlab来实践一下，代码如下：
%%本代码用于学习如何对图像进行对数变换操作
% 按照公式生成对数函数
x = 0:255;
c = 255 / log(256);
y = c * log(x+1);
figure, subplot(2,2,1), plot(y), title('对数变换函数'), axis tight, axis square
%I = imread('tiger.jpg');
I = imread('radio.tif');
I_log = uint8(y(I+1));
subplot(2,2,2), imshow(I), title('原图');
subplot(2,2,3), imshow(I_log), title('对数变换调整');
% 增强图像亮度
I_br = imadd(I, 100);
subplot(2,2,4), imshow(I_br), title('图像增强亮度');
再测试一张图片看看效果：
2、分段线性变换
分段线性变换可以使用若干个线性方程来描述，每个用于输入图像中灰度值的一个区间。它的主要优点是它们可有任意的复杂度；而缺点是需要更多的用户输入。
下图为书中所给的例子，用来改善输入图像的对比度。
%%本代码用于学习如何对图像进行分段线性变换操作
% 按照公式生成对数函数
y(1:175) = 0:174;
y(176:200) = 255;
y(201:256) = 200:255;
figure, subplot(1,3,1), plot(y), title('分段线性变换函数'), axis tight, axis square
%I = imread('tiger.jpg');
I = imread('micro.tif');
subplot(1,3,2), imshow(I), title('原图');
I2 = uint8(y(I+1));
subplot(1,3,3), imshow(I2), title('分段线性变换调整');
从上图可以看出，调整后的图片在特定范围像素更亮，但仍然有其他干扰。我们试一次另一个分段函数。
%%本代码用于学习如何对图像进行分段线性变换操作
% 按照公式生成对数函数
y(1:175) = 0:174;
y(176:200) = 255;
y(201:256) = 200:255;
figure, subplot(1,3,1), plot(y), title('分段线性变换函数'), axis tight, axis square
%I = imread('tiger.jpg');
I = imread('micro.tif');
subplot(1,3,2), imshow(I), title('原图');
I2 = uint8(y(I+1));
subplot(1,3,3), imshow(I2), title('分段线性变换调整');
z(1:175) = 50;
z(176:200) = 250;
z(201:256) = 50;
I3 = uint8(z(I+1));
figure, subplot(1,2,1), plot(z), xlim([0 255]), ylim([0 255]), axis square;
subplot(1,2,2), imshow(I3)
这样效果就好很多了。
以上是本期内容，建议自己动手来实践一下。如果需要图片和代码，可以关注公众号，回复关键字【matlab代码】即可获得。我会上传到云盘。
【声明】：学习笔记基于互联网上各种学习资源的个人整理。
以上是本期内容，下期介绍图像处理的直方图处理。
我叫小保，一名计算机视觉爱好者、学习者、追随者，欢迎关注我【CV之道】一起学习。