医学图像处理实验报告要Word文档下载推荐.docx

1、1.1.2 显示图像一般使用函数imshow显示，其基本语法为：imshow （f, G）其中，f是一个图像数组，G是显示该图像的灰度级数。若G省略，默认为256。语法：imshow f, low high会将所有小于或等于low的值都显示为黑，所有大于或等于high的值都显示为白色。介于low和high之间的值将以默认的级数显示为中等亮度值。当用imshow显示另外一幅图像时，MATLAB会在屏幕上用新图像替换旧图像，为保持第一幅图像并同时显示第二副图像，可以使用如下figure函数：figure, imshow（g）要了解图像文件的其他详细信息，可以使用imfinfo函数，其语法结构为：i

2、mfinfo filename其中，filename是存储在磁盘中的图像全名。1.1.3 保存图像使用函数imwrite可以将图像写在磁盘上，该函数语法为：imwrite （f, filename）下面的命令可以把图像f写为TIFF格式，在默认情况下使用语句imwrite （f,patient10_run1, tif或imwrite （f, patient10_run1.tif会将文件保存到当前的工作目录中。另外，如果是JPEG图像，其语法可以为：filename.jpgquality, q）其中q是一个在0到100之间的整数（q越小，图像退化越严重）。函数imwrite另一种常用但只用于ti

3、f图像的语法为：imwrite （g, filename.tif, compression, parameter, resolution, colres rowres）其中，parameter有三种情况：为none表示无压缩，packbits表示比特压缩包压缩（非二值图像默认参数），ccitt表示ccitt压缩（二值图像默认参数）。colres为列分辨率。1.1.4 数据类型见下表：表格 01 数据类型。前8项数值类，第9项字符类，第10项逻辑类。名称描述double双精度浮点数，范围【-10308,10308】（8比特每像素）unit8无符号8比特整数，范围【0,255】（1比特每像素）un

4、it16无符号16比特整数，范围【0,65535】（2比特每像素）unit32无符号32比特整数，范围【0,4294967295】（4比特每像素）int8有符号8比特整数，范围【-128,127】int16有符号16比特整数，范围【-32768,32767】int32有符号32比特整数，范围【-2147483648,2147483647】single单精度浮点数，范围【-1038，1038】char字符（2比特每像素）logical值为0或者1（1比特每像素）1.1.5 图像类型主要分为4种：亮度图像（Intensity image）二值图像（Binary image）索引图像（Indexed

5、 image）RGB图像（RGB image）注意：一副图像的特征是由数据类和图像类这两者来表征的。例如，“unit8亮度图像”表示一副像素都是unit8数据类的亮度图像。1.1.6 数据类与图像类型之间的转换（1）数据类间的转换通用语法为：B=data_class_name（A）其中，data_class_name可以为表1第一列中的任意一项。（2）图像类和数据类型间的转换表格 02 IPT中用于进行图像类和数据类转换的函数将输入转换为有效的图像数据类im2unint8logical,uint8,uint16,doubleim2uint16mat2graydouble,范围【0,1】im2d

6、oubleim2bwuint8,uint16,double要把一个double类的任意数组转换成取值范围为0,1的归一化double类数组，可以通过函数mat2gray完成，其基本语法为：g=mat2gray（A, Amin, Amax）其中，图像g的取值范围在0（黑色）到1（白色）之间。转换时，小于Amin的值在g中变为0，大于Amax的值变为1。1.1.7 数组索引在MATLAB中，v=1 3 5 7 9代表一维行向量，使用转置运算符（.）可以把行向量转为列向量。要存取元素数据块，可以使用冒号（:）。例如，下列语句代表存取第三个到最后一个元素：v（3:end）ans= 5 7 9符号1:2

7、：end表示索引从1开始计数，步长为2，直到结束。步长可以为负。例如，求一个矩阵所有元素之和：s=sum（A（:）1.1.8 运算符MATLAB运算符一般分为以下三种主要类别：执行数值计算的算术运算符在数量上比较操作数的关系运算符执行函数AND、OR和NOT的逻辑运算符常用的运算符和算术函数如下面几个表所示：表格 03 数组和矩阵运算符运算符MATLAB函数注释和示例+数组和矩阵加plus（A,B）a+b,A+B或a+A-数组和矩阵减minus（A,B）a-b,A-B或a-A.*数组乘times（A,B）C=A.*B,C（I,J）=A（I,J）*B（I*J）*矩阵乘mtimes（A,B）标准矩

8、阵乘./数组右除rdivide（A,B）C=A./B.数组左除ldivide（A,B）C=A.B/矩阵右除mrdivide（A,B）矩阵左除mldivide（A,B）.数组求幂power（A,B）C=A.B,有C（I,J）=A（I,J）B（I,J）矩阵求幂mpower（A,B）.|向量和矩阵转置transpose（A）标准的向量和矩阵转置|向量和矩阵复共轭转置ctanspose（A）标准的向量和矩阵共轭转置一元加uplus（A）+A=0+A一元减uminus（A）-A=0-A:冒号表格 04 IPT支持的图像算术函数函数imadd两幅图像相加或把常数加到图像imsubtract两幅图像相减或从

9、图像中减去常数immultiply两幅图像相乘，在相应的像素对间进行；或图像乘以一个常数imdivide两幅图像相除，在相应的像素对间进行；或图像除以一个常数imadsdiff计算两幅图像间的绝对差imcomplement对图像求补imlicomb计算两幅或多幅的线性组合表格 05 关系运算符小于=小于等于大于大于等于=等于=不等于表格 06 逻辑运算符&AND（与）OR（或）NOT（非）表格 07 逻辑函数注释xor若两个操作数逻辑上不同，则函数xor返回1；否则，返回0all若一个向量中的所有元素都非零，返回1；any若一个向量中的任何元素都非零，返回1；否则，返回0。该函数在矩阵中按列操

10、作1.2 实验内容1.2.1 在“素材”文件夹下，读入“Fig2.19（a）.jpg”并且提取该图像的基本信息，写出实验步骤和实验结果。图 1 图像的读取与显示 f=imshow（Fig2.19（a）.jpg whos f Name Size Bytes Class Attributes f 1x1 8 double 1.2.2 读取图像“Fig1.14（d）.jpg”，使用函数imwrite和subplot函数分别显示质量参数q=50,25,5,0的结果如图所示。写出实验步骤。 f=imread（Fig1.14（d）.jpg subplot（3,2,1）,imshow（f）; subplot

11、（3,2,2）,imshow（bubbles50.jpg imwrite（f,bubbles25.jpg,25） subplot（3,2,3）,imshow（bubbles15.jpg,15） subplot（3,2,4）,imshow（bubbles5.jpg,5） subplot（3,2,5）,imshow（bubbles0.jpg,0） subplot（3,2,6）,imshow（观察哪副图像压缩最大后没有“伪轮廓”。此时键入命令： imfinfo bubbles25.jpg会得到如下信息：ans = Filename: FileModDate:15-十月-2012 15:49:44 FileSize: 13875 Format:jpg FormatVersion: Width: 714 Height: 682 BitDepth: 8 ColorType:grayscale FormatSignature: NumberOfSamples: 1 CodingMethod:Huffman CodingProcess:Sequential Comment: 其中，使用width乘以height乘以BitDepth再除以8，得到：714*682*8/8=486948字节。把这个结果除以FileSize就可以的得到压缩比：486948/138