求负片&灰度扩展&幂律变换

提示：文章写完后，目录可以自动生成，如何生成可参考右边的帮助文档

 

文章目录

• 前言
• 一、imadjust、stretchlim函数简介
• 二、代码展示
• 三、结果展示
• 总结

 

---

前言

在进行医学图像处理时灰度变换是非常重要的操作手段，本节我将学习imadjust在灰度变换方面的应用，另外还将引入stretchlim函数作为灰度变换区间的辅助。

 

一、imadjust、stretchlim函数简介

imadjust函数是matlab中常用的灰度变换函数，在进行求负片、灰度扩展和幂律变换时非常有用。

imadjust函数用法非常简单，具体的用法格式是：

g1=imadjust(f,[low_in high_in],[low_out high_out],gamma);

在这里g1指的是变换后的图像，具体的操作指的是将[low_in high_in]之间的值映射到[low_out high_out]之间，实现灰度扩展；gamma指的是扩展的方式，一般默认采用gamma=1的线性变换方式，当然在幂律变换过程中gamma可以不为1.

如果gamma>1,输出图像关注的是较暗的部分；gamma<1，输出图像关注的是较亮的部分。

在实际操作过程中其实确定我们想要变换的灰度区间是一个相对来说很难精准确定的事，因此我们可以引入stretchlim函数，这个函数可以帮助我们自动确定要变换的参数区间而不需要我们操作，方便了imadjust函数的使用

二、代码展示

I1=imread('breast.tif');
I2=imadjust(I1,[0,1],[1,0],1);%求负片
I3=imadjust(I1,[0.5,0.75],[0,1],1);%将[0.5_0.75]映射至[0，1]
I4=imadjust(I1,[],[],2);%进行r=2的幂律变换
I5=imadjust(I1,stretchlim(I1),[],1);%自动确定变换区间
I6=imadjust(I1,stretchlim(I1),[1,0],1);%自动确定变换区间并求负片
subplot(231),imshow(I1),title('原图');
subplot(232),imshow(I2),title('图像反转');
subplot(233),imshow(I3),title('[0.5-0.75]扩展');
subplot(234),imshow(I4),title('幂律变换r=2');
subplot(235),imshow(I5),title('自动确定变换灰度值');
subplot(236),imshow(I6),title('自动确定变换灰度值并求负片');

三、结果展示

---

总结

可以看到imadjust函数在灰度扩展等方面具有强大的处理功能，而利用stretchlim函数处理输入区间之后输出的图像比我们自己确定区间对比效果要更好一些。