MATLAB图像处理入门

一、读写图像、直方图 

1、读取并显示图像

imread与imshow函数

clc;clear; close all;
I = imread('pout.tif');
imshow(I);

结果:

MATLAB图像处理入门_第1张图片

如果我们想显示图像坐标轴的话只需在imshow后加上下面这行代码:

set(gca,'Visible','on');

MATLAB图像处理入门_第2张图片

 我们还可以用whos命令查看图像的信息:

在命令行窗口输入whos I查看上面那种图片的信息:

MATLAB图像处理入门_第3张图片

 2、直方图

显示图像的直方图用hist函数,将直方图均衡化用histeq函数

我们在上述代码基础上继续输入下面代码:

figure();
imhist(I);

结果:

MATLAB图像处理入门_第4张图片

 可以看到图像的对比度比较低,我们进行直方图均衡,继续输入下面代码:

I2 = histeq(I);
imshow(I2);
figure();
imhist(I2);

结果: 

MATLAB图像处理入门_第5张图片

 3、写入图像

我们使用imwrite函数将上述均衡化的图像I2写入磁盘:

imwrite(I2, 'output.png');

可以看到左侧文件区多了output.png

我们在命令行窗口使用imfinfo函数能看见该图片更为详细的信息:

MATLAB图像处理入门_第6张图片

 拓展:

使用montage函数能一次显示多张图片在一个窗口,并且这些图像不需要是同一尺寸和大小,我们只需要将在当前目录下或已经在matlab路径里的图片的名字放在一个字典里:

运行下面代码:

clc;clear; close all;
files = {'catdog.jpg', 'jiedao.jfif', 'Fig0638(a)(lenna_RGB).tif', 'output.png', 'Fig0335(a)(ckt_board_saltpep_prob_pt05).tif'};
montage(files)

结果:

MATLAB图像处理入门_第7张图片

 我们甚至可以看到不是同一类型的(灰度图和RGB图)都能一起展示!

二、图像类型转换:

在MATLAB的官方文档里有着这么几个图像类型转换的函数:

MATLAB图像处理入门_第8张图片

 我们选择几个常用的讲解:

1、img2gray将 RGB 图像转换为灰度图像

代码:

clc;clear; close all;
RGB = imread('example.tif');
figure;
imshow(RGB);
I = im2gray(RGB);
figure;
imshow(I);

结果:

MATLAB图像处理入门_第9张图片

im2gray 函数与 rgb2gray 基本相同,不同之处是它可以接受灰度图像作为输入并原样返回它们。如果输入图像是灰度图像,则 rgb2gray 函数返回错误。

2、mat2gray将矩阵转换为灰度图像 

 代码:

clc;clear; close all;
mat = [0:255]' * ones(1,256);
mat = mat';
img = mat2gray(mat);
imshow(img);

结果:

MATLAB图像处理入门_第10张图片

 3、imsplit将多通道图像分割成其独立的通道

代码:

clc;clear; close all;
I = imread('peppers.png');
imshow(I);
[r,g,b] = imsplit(I);
montage({r,g,b},'Size',[1 3]);

结果:

MATLAB图像处理入门_第11张图片

4、rgb2hsv将RGB图像转为HSV图像 

代码:

clc;clear; close all;
rgbImage = imread('peppers.png');
imshow(rgbImage);
hsvImage = rgb2hsv(rgbImage);
[h,s,v] = imsplit(hsvImage);
montage({h,s,v},'Size',[1 3])

结果:

MATLAB图像处理入门_第12张图片

还有很多颜色空间之间的转换就不一一列举了:

MATLAB图像处理入门_第13张图片

 5、灰度图像二值化

MATLAB图像处理入门_第14张图片

这里就不演示了

6、数据类型转换

 其实这些都可以直接转换,比如我们想将uint16的矩阵I转换为uint8直接I = uint8(I)即可,其他数据类型转换也是同理!

示例:

MATLAB图像处理入门_第15张图片

三、几何变换和图像配准

常见几何变换

MATLAB图像处理入门_第16张图片

1、imcrop裁剪图像 

Icropped = imcrop(I,rect)裁剪图像

clc;clear; close all;
I = imread('circuit.tif');
I2 = imcrop(I,[75 68 130 112]);
axis image;
subplot(1,2,1);
imshow(I);
title('Original Image');
subplot(1,2,2);
imshow(I2);
title('Cropped Image');
axis image;

结果:

MATLAB图像处理入门_第17张图片

 注意matlab的图像坐标系为:

在这里插入图片描述

使用空间参照矩形使裁剪图像居中

I = imread('parkavenue.jpg');
imshow(I)

将目标窗口大小指定为 [宽度高度] 形式的二元素向量。

targetSize = [300 600];

 创建指定裁剪窗口的空间范围的 Rectangle 对象。

r = centerCropWindow2d(size(I),targetSize);

在空间范围的界限处裁剪图像。显示裁剪的区域。 

J = imcrop(I,r);
imshow(J)

结果:

MATLAB图像处理入门_第18张图片

 2、imresize调整图像大小

J = imresize(I,scale) 返回图像 J,它是将 I 的长宽大小缩放 scale 倍之后的图像。输入图像 I 可以是灰度图像、RGB 图像、二值图像或分类图像。

代码:

clc;clear; close all;
I = imread('rice.png');
J = imresize(I, 0.5);
figure;
imshow(I);
title('Original Image');
figure;
imshow(J);
title('Resized Image');

结果:

MATLAB图像处理入门_第19张图片

我们都知道图像缩放需要进行插值,我们可以在imresize函数中自行选择插值方法:

MATLAB图像处理入门_第20张图片

我们在下面的代码中指定最近邻插值(其他代码不变):

J = imresize(I, 0.5,'nearest');

MATLAB图像处理入门_第21张图片

J = imresize(I,[numrows numcols]) 返回图像 J,其行数和列数由向量 [numrows numcols] 指定。

代码:

clc;clear; close all;
RGB = imread('peppers.png');
RGB2 = imresize(RGB, [128 128]);
figure;
imshow(RGB);
title('Original Image');
figure;
imshow(RGB2);
title('Resized Image');

结果:

MATLAB图像处理入门_第22张图片

 3、imroate旋转图像

J = imrotate(I,angle) 将图像 I 围绕其中心点逆时针方向旋转 angle 度。要顺时针旋转图像,请为 angle 指定负值。imrotate 使输出图像 J 足够大,可以包含整个旋转图像。默认情况下,imrotate 使用最近邻点插值,将 J 中位于旋转后的图像外的像素的值设置为 0

 代码:
 

clc;clear; close all;
I = imread('peppers.png');
J = imrotate(I, 45, 'bilinear','loose');
figure(1);
imshow(I);
figure(2);
imshow(J);

结果:

MATLAB图像处理入门_第23张图片

 我们将bbox参数改为crop观察图片大小是否一致:

J = imrotate(I, 45, 'bilinear','crop');

 MATLAB图像处理入门_第24张图片

 可以看到改为'crop'后图片大小一致!

 4、图像金字塔impyramid

B = impyramid(A,direction) computes a Gaussian pyramid reduction or expansion of A by one level. direction determines whether impyramid performs a reduction or an expansion.

MATLAB图像处理入门_第25张图片代码:

clc;clear; close all;
I = imread('cameraman.tif');
I1 = impyramid(I, 'reduce');
I2 = impyramid(I1, 'reduce');
I3 = impyramid(I2, 'reduce');
figure, imshow(I);
figure, imshow(I1);
figure, imshow(I2);
figure, imshow(I3);

结果:

MATLAB图像处理入门_第26张图片

四、图像滤波和增强 

1、图像滤波

1.1、fspecial创建滤波器核、imfilter根据选定的滤波核进行滤波

MATLAB图像处理入门_第27张图片

MATLAB图像处理入门_第28张图片

注:图像滤波也可以用filter2函数,与imfilter的区别请自行help命令查看 

MATLAB图像处理入门_第29张图片 注意上面说的imfilter默认使用相关性而不是卷积,若要使用卷积进行滤波则需要指定'conv'参数!

示例:

读取图像并显示:

I = imread('cameraman.tif');
imshow(I);

MATLAB图像处理入门_第30张图片

创建运动滤波器,并使用它来对图像进行模糊处理。显示模糊处理后的图像。 

H = fspecial('motion',20,45);
MotionBlur = imfilter(I,H,'replicate');
imshow(MotionBlur);

MATLAB图像处理入门_第31张图片

 创建圆形滤波器,并使用它来对图像进行模糊处理。显示模糊处理后的图像。

H = fspecial('disk',10);
blurred = imfilter(I,H,'replicate'); 
imshow(blurred);

MATLAB图像处理入门_第32张图片

在执行imfilter时要注意的是将图像转换为uint8型以避免滤波核有负数值导致输出有负数:

MATLAB图像处理入门_第33张图片

 在此要与我们图像之间的运算做区分!我们想把两个图像相加再除2即取均值时我们往往先将两个图像转换为double来防止uint8的截断!而这里则是要刻意使用uint8来截断,注意细节!

1.2、imgaussfilt图像二维高斯滤波

MATLAB图像处理入门_第34张图片

 注意,这里的滤波核大小未给出,经查看源码后,滤波核大小默认为3*3,sigma默认为0.5,如果要自己指定滤波核大小,以第三种键值对的方式给出:

示例: B = imgaussfilt(A,'FilterSize',3);

代码演示:

clc;clear; close all;
I = imread('cameraman.tif');
Iblur = imgaussfilt(I,2, 'FilterSize',7); % sigma = 2, kernel_size = 5*5
montage({I,Iblur});
title('Original Image (Left) Vs. Gaussian Filtered Image (Right)');

MATLAB图像处理入门_第35张图片

 2、对比度调整

MATLAB图像处理入门_第36张图片

同样的,我们选择常用函数进行演示!

2.1、imadjust调整图像强度值

MATLAB图像处理入门_第37张图片

MATLAB图像处理入门_第38张图片

调整灰度图像的对比度

clc;clear; close all;
I = imread('pout.tif');
J = imadjust(I);
montage({I,J});

MATLAB图像处理入门_第39张图片

调整灰度图像的对比度,指定对比度限制 

clc;clear; close all;
I = imread('pout.tif');
K = imadjust(I,[0.3 0.7],[]); %low_out和high_out未给出时默认等于low_in,high_in
montage({I,K});

MATLAB图像处理入门_第40张图片

 调整彩色图像的对比度

clc;clear; close all;
RGB = imread('football.jpg');
RGB2 = imadjust(RGB,[.2 .3 0; .6 .7 1],[]);%对应彩色图像三通道的low_in,high_in
montage({RGB,RGB2});

MATLAB图像处理入门_第41张图片

 2.2、imsharpen图像锐化

MATLAB图像处理入门_第42张图片

 图像锐化

clc;clear; close all;
a = imread('hestain.png');
imshow(a)
title('Original Image');
b = imsharpen(a);
figure, imshow(b)
title('Sharpened Image');

MATLAB图像处理入门_第43张图片

 在边界上控制锐化的数量

clc;clear; close all;
a = imread('rice.png');
imshow(a), title('Original Image');
b = imsharpen(a,'Radius',2,'Amount',1);
figure, imshow(b)
title('Sharpened Image');

MATLAB图像处理入门_第44张图片

2.3、histeq直方图均衡化这个上文提过,这里就不多说了

2.4、imnoise向图像添加噪声

MATLAB图像处理入门_第45张图片

clc;clear; close all;
I = imread('eight.tif');
J = imnoise(I,'salt & pepper',0.02);
montage({I,J});
title('Original Image(left) and Noise Image(right)');

MATLAB图像处理入门_第46张图片

3、形态学运算

MATLAB图像处理入门_第47张图片

3.1、imerode腐蚀 

其中的SE:

MATLAB图像处理入门_第48张图片

original_img = imread('Fig0905(a)(wirebond-mask).tif');
subplot(1,4,1);
imshow(original_img);
%三次不同大小的腐蚀操作比较
SE = strel('rectangle',[11 11]);
out = imerode(original_img, SE);
subplot(1,4,2);
imshow(out);

SE = strel('rectangle',[15 15]);
out = imerode(original_img, SE);
subplot(1,4,3);
imshow(out);

SE = strel('rectangle',[45 45]);
out = imerode(original_img, SE);
subplot(1,4,4);
imshow(out);

注意,每一次腐蚀都会导致中间的白圈变小,这里不大明显! 

MATLAB图像处理入门_第49张图片

 下面我们用nhood再演示一下腐蚀操作:

original_img = imread('Fig0905(a)(wirebond-mask).tif');
subplot(1,2,1);
imshow(original_img);

nhood = ones(45);
out = imerode(original_img, nhood);
subplot(1,2,2);
imshow(out);

MATLAB图像处理入门_第50张图片

 3.2、imdilate膨胀

 和腐蚀的参数差不多!

original_img = imread('Fig0907(a)(text_gaps_1_and_2_pixels).tif');
imshow(original_img);

nhood = [0 1 0;1 1 1;0 1 0];
out = imdilate(original_img, nhood);
figure();
imshow(out);

MATLAB图像处理入门_第51张图片

 放大后观察效果:

MATLAB图像处理入门_第52张图片

 3.3、imopen开运算

与腐蚀和膨胀的语法差不多哈

3.4、imclose闭运算 

下面这个示例演示先开运算在闭运算结果

original_img = imread('Fig0911(a)(noisy_fingerprint).tif');
subplot(1,3,1);
imshow(original_img);

nhood = ones(3);
out = imopen(original_img, nhood);
subplot(1,3,2);
imshow(out);

out = imclose(out, nhood);
subplot(1,3,3);
imshow(out);

MATLAB图像处理入门_第53张图片

4、基于关注区域的处理(待更新)

5、图像算术

MATLAB图像处理入门_第54张图片

 5.1、imabsdiff

clc;clear; close all;
I = imread('cameraman.tif');
J = uint8(filter2(fspecial('gaussian'), I));
K = imabsdiff(I,J);
figure;
imshow(K,[]);%将K中像素最小值设为黑色,最大值设为白色

MATLAB图像处理入门_第55张图片

 5.2、imadd

将两个图像相加。将输出指定为 uint16 类型以避免截断结果。

clc;clear; close all;
I = imread('rice.png');
J = imread('cameraman.tif');
K = imadd(I,J,'uint16');
imshow(K,[]);

MATLAB图像处理入门_第56张图片

 5.3、imsubtract

 从一个图像中减去另一个图像或从图像中减去常量

MATLAB图像处理入门_第57张图片

减去图像背景 

clc;clear; close all;
I = imread('rice.png');
%估计背景
background = imopen(I,strel('disk',15));
%从图像中减去背景
J = imsubtract(I,background);
montage({I,J});

MATLAB图像处理入门_第58张图片

五、傅里叶变换 

1、fft2实现图像的快速傅里叶变换

MATLAB图像处理入门_第59张图片

 代码演示:

clc;clear;close all;
im = imread('street.jfif');
figure(1);
imshow(im);
im = rgb2gray(im);
out = fft2(im);
figure(2);
imshow(abs(out),[]); %取实部

MATLAB图像处理入门_第60张图片

可能会有疑问,为什么傅里叶变换黑的啥都没有呢?事实上如果我们将Figure2的左上角放大,其实是这个样子:

MATLAB图像处理入门_第61张图片 所以事实是只有左上角值比较大看起来是白色,而这就是我们的直流分量。事实上这种现象的原因是直接用fft2的话频谱的原点不在图像中心 而在左上角。

如果我们想将原点移到中心,我们需要用fftshift函数:

2、fftshift将频谱原点移到中心

MATLAB图像处理入门_第62张图片

 让我们在上面的代码加上:

figure(3);
out = fftshift(out);
imshow(abs(out),[]);

 这样就对了!

3、ifft2进行逆傅里叶变换

下面我们进行频域的低通滤波再逆傅里叶变换查看结果

clc;clear;close all;
im = imread('Fig0638(a)(lenna_RGB).tif');
figure(1);

im = rgb2gray(im);%转灰度
imshow(im);
title('原图');
[M N] = size(im);
im_f = fft2(im);

lp = fspecial('gaussian',25,2);%高斯低通滤波器(空间域)(补零到与原图一样大小)
lp_f = fft2(lp,M,N);%高斯低通滤波器(频域)
figure(2);
surf(fftshift(abs(lp_f)),'edgecolor','none');%surf函数画曲面图
title('滤波核频域曲面图');

out_f = im_f.*lp_f; %频域相乘
out = ifft2(out_f);%空间域结果
figure(3);
imshow(abs(out),[]);
title('结果');


MATLAB图像处理入门_第63张图片

 这里我们的滤波核是用的自带的,当然我们也可以自定义一个频域低通滤波器:

%D(u,v)<=D0时有值
u = 0:Q-1
v = 0:P-1
%calculate D(u,v)
H = double(D<=D0)

参考资料:

MATLAB官方文档的Image Processing Toolbox

你可能感兴趣的:(matlab,图像处理,开发语言)