图像处理学习之图像增强

实际操作

例1：读取 dicom 图像，完成以下操作： 利用线性函数增强图像对比度、降低对比度、增加亮度和 降低亮度，并分别对比显示。

I=dicomread('CT000010');
I1=5*I;
I2=0.5*I;
I3=I+10000;
I4=I-4000;
subplot(2,3,1);imshow(I);title('原始图像');
subplot(2,3,2);imshow(I1);title('增强对比度图像');
subplot(2,3,3);imshow(I2);title('减小对比度图像');
subplot(2,3,4);imshow(I3);title('增加亮度图像');
subplot(2,3,5);imshow(I4);title('减小亮度图像');

例2：分别利用对数变换、反对数变换、幂次变化对图像进行处理，对比显示实验结果。

I=dicomread('CT000010');
v=20;
I1=mat2gray(double(I));
Ilog1=log(1+v*I1)/(log(v+1));
Ilog2=log(1+I1);
Ivlog=((v+1).^I1-1)/v;
Ipower1=double(I1).^2;
Ipower2=double(I1).^0.5;
subplot(2,3,1);imshow(I);title('原始图像');
subplot(2,3,2);imshow(Ilog1);title('对数变换1');
subplot(2,3,3);imshow(Ilog2);title('对数变换2');
subplot(2,3,4);imshow(Ivlog);title('反对数变换');
subplot(2,3,5);imshow(Ipower1);title('幂次变换＞1');
subplot(2,3,6);imshow(Ipower2);title('幂次变换＜1');

例3：用以下两种方法开窗显示 dicom 图像：
（1）读取图像的窗 宽窗位，并将其从 CT 值转换成灰度值，代码需包含窗宽、 窗位的读取及转换过程。利用窗宽、窗位显示 CT 图像。
（2） 将图像从灰度值转换为 CT 值。

%开窗显示法的2种方法
I0=dicomread('CT000010');
Iinfo=dicominfo('CT000010');
ww=Iinfo.WindowWidth;
wc=Iinfo.WindowCenter;
wl=wc-ww/2;
wh=wc+ww/2;
slope=Iinfo.RescaleSlope;
inter=Iinfo.RescaleIntercept;
%方法一：CT值转灰度值
cwl=(wl-inter)/slope;
cwh=(wh-inter)/slope;
%方法二：灰度值转CT值
Ict=I0*slope+inter;
subplot(1,3,1);imshow(I0);title('原始图像CT000010');
subplot(1,3,2);imshow(I0,[cwl,cwh]);title('开窗显示法1');
subplot(1,3,3);imshow(Ict,[wl,wh]);title('开窗显示法2');

例3：分别利用 mat2gray、imadjust 对图像进行变换，拉伸其对比 度，并将其另存为 bmp 格式图像。

%利用 mat2gray、imadjust 对图像进行变换
I0=dicomread('CT000010');
Iinfo=dicominfo('CT000010');
subplot(1,3,1);imshow(I0);title('原始图像CT000010');
ww=Iinfo.WindowWidth;
wc=Iinfo.WindowCenter;
wl=wc-ww/2;
wh=wc+ww/2;
slope=Iinfo.RescaleSlope;
inter=Iinfo.RescaleIntercept;
%CT值转灰度值
cwl=(wl-inter)/slope;
cwh=(wh-inter)/slope;
%使用函数mat2gray
Imat2=mat2gray(I0,[cwl,cwh]);
imwrite(Imat2,'Imat2.bmp');%另存为 bmp 格式图像
%使用函数imadjust
I0=double(I0);
I1=double(I0+abs(min(I0(:))));
I2=I1/(max(I1(:)));
l=(cwl+abs(min(I0(:))))/max(I1(:));
h=(cwh+abs(min(I0(:))))/max(I1(:));
Iad=imadjust(I2,[l,h],[0,1]);
imwrite(Iad,'Iad.bmp');%另存为 bmp 格式图像
subplot(1,3,2);imshow(Imat2);title('mat2gray处理后的图像');
subplot(1,3,3);imshow(Iad);title('imadjust处理后的图像');

运行结果：

原来的工作区间：

运行成功后在工作区间能看到保存的bmp格式图片