图形图像处理识别技术与医学信息分析MATLAB图像处理与分析.pptVIP

MATLAB的图像处理与分析 1. MATLAB软件 软件的工作环境如图9-13所示，右侧的Command Window称为命令窗口，提示符为“”，在命令窗口中可以直接进行简单的算术运算和函数调用。 定义程序文件 在主菜单中选择“File”|“New”|“M-file”打开MATLAB编辑窗口，建立程序文件，如图9-14所示，MATLAB文件的扩展名为“.m”。 运行程序文件(两种方式) 在MATLAB命令窗口中输入文件名 选择“File”|“Open”打开m文件，在编辑窗口下选择“Debug”|“Run”或者按F5键运行该程序，如图9-15所示。 一、显示一张图片 例题1 ：建立程序文件，读入图像“MRI_brain.jpg”并显示。 操作提示： （1） 单击开始 |所有程序| MATLAB | MATLAB ，启动软件； （2） 单击File |New|M-file，打开编辑窗口； （3） 在编辑窗口输入以下程序 X=imread(MRI_brain.jpg, jpg); imshow(X); title(显示一幅MRI图像); （4） 单击File | save将文件保存为 example1.m ； （5） 单击“Debug”|“Run”或者按F5键运行程序， 结果如图9-16所示。 * * 图9-13 MATLAB编程环境 图9-16 1. imread(): 函数imread(‘filename’,’format’)是读入图像，filename是要读入的图像的文件名，包括路径和扩展名，当图像文件和matlab文件在同一目录下时，可省略路径；fotmat是图像的格式，可以支持jpg，bmp，gif，png等众多格式；filename和fotmat都要用单引号； 2. imshow(); Imshow(i):显示一幅图像 3. title(‘text); 显示文本text，文本text用单引号引起来 函数说明: 二、显示灰度图像和直方图 左侧工作区双击文件m02.m打开，则会打开一个Editor，里面就是matlab的源代码，m02.m如下： %%显示灰度图像和直方图 clear all %清空内存数据 G = imread(‘X_chest.bmp,bmp);%读入原始图像，为胸部的X片图 subplot(2,1,1);imshow(G,256);title(原始图像);%显示原始图像 subplot(2,1,2);imhist(G);title(原始图像的直方图);%显示原始图像的直方图 按F5运行得如下结果：见图9-17 图9-17 subplot(2,1,1);imshow(G,256);title(‘原始图像’);%显示原始图像 subplot(2,1,2);imhist(G);title(原始图像的直方图);% 函数说明: 1. 符号%表示其之后的内容为注释 2.函数subplot(m,n,p) 功能是划分绘图窗口，创建m×n排列的子图阵列，p是子图的序号。比如subplot(2,1,1);是创建两行一列的子图阵列，且是第一副图。 3. 函数imshow(I,N) 功能是以N级灰度显示灰度图象I.。imshow(G,256);就是显示256级灰度图象G。 4. 函数imhist(I) 是显示灰度图象I的直方图。 三、图像的旋转 打开文件m03.m: %图像的旋转 clear all G=imread(X_chest.bmp,bmp); %读入原始图像，为胸部的X片图 R=imrotate(G,30,bilinear); %逆时针旋转30度 subplot(1,2,1);imshow(G);title(原图); %显示原图 subplot(1,2,2);imshow(R);title(逆时针旋转30度后的图像);%显示逆时针旋转30度后的图像 按F5运行得如下结果：见图9-18 图9-18 1. 函数imrotate(G,30,’bilinear’); 是将图形G逆时针方向旋转30度，使用’bilinear’双线性法差值，若要顺时针方向旋转30度，则需要在度数30的前面加个负号 函数说明: R=imrotate(G,30,bilinear); %逆时针旋转30度 四、减影处理 打开文件m04.m: %%减影处理 clear all G1=double(imread(D21.png,png)); %读入第一副图 G2=double(imread(D22.png,png)); %读入第二幅图 [m,n]=size(G1);