医学图像处理绪论-new.pptVIP

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * Matlab环境 * * 图像的输入输出 读取图像 读入Tiff,bmp,jpeg,gif,png图像 X=imread(‘filename’); 显示读入的图像 imshow(X) 同时显示多幅图像 figure,imshow(g) 保存图像 imwrite(f,’filename’) 读入dicom图像 [X , map] = dicomread(filename) 显示读入的图像 montage(X, map); 保存dicom图像dicomwrite(X, filename) * * 图像类型 亮度图像 像素值表示亮度 uint8 uint16 double 二值图像 索引图像 RGB图像 * * Matlab中的数据类 名称 描述 double Uint8 Uint16 Uint32 Int8 Int16 Int32 Single Char logical 双精度浮点数，范围：-10308~10308（8字节/像素） 无符号整数 [0 255] 1字节/像素 [0 65535] 2字节/像素 [0 4294967295] 4字节/像素 有符号整数 [-128 127] 1字节/像素 [-32768，32767] 2字节/像素 [-2147483648，2147483647] 4字节/像素 单精度浮点数，范围：-1038~1038（4字节/像素） 字符（ 2字节/像素） 值为0或1 1字节/像素 * * 数据类和图像类型之间的转换 数据类间的转换 B=data_class_name(A) 图像类和类型间的转换 * * 数组索引 向量（行向量）创建 向量的元素用方括号括起，并由空格或逗号分开。例如： v=[1 3 5 7 9] W=[1，2，3，4,5] 单个像素的访问: 向量名(下标) 如：v(3) 访问连续元素 ：例如 要存取w的中间三个元素 w(2:4) 存取w中的最后三个元素w(3:end) 访问不连续数据 w(1:2:end) 从1开始计数步长为2(步长可以为负数) end结束 一个向量可以做另一个向量的索引 W([1 4 5]) 从行向量产生列向量 向量转置（.’） x=w.’ w(:) 注意与w(1:end)的区别 * * 矩阵索引 矩阵的创建 用一列被方括号括起来并用分号隔开的行向量表示 例如：A=[1 2 3；4 5 6；7 8 9] 矩阵元素的访问 访问单个元素A(2,3) 访问二维元素块 A(:,3) A(2,:) A(1:2,2:3) A(2:end,end:-2:1) 使用向量作为矩阵的索引 E=A([1 3],[2 3]) 使用逻辑数组 D=logical([1 0 0;0 0 1;0 0 0]);A(D) A(:) 逐列的形式将元素组成列向量 sum(A(:)) sum(A) 区别 * * 一些重要的标准数组 zeros(M,N)生成一个大小为M*N的double类矩阵，其元素均为0 ones(M,N) ture(M,N) false(M,N) rand(M,N) * * M函数编程简介 M文件可以是简单的执行一系列Matlab语句的脚本，也可以是接受变量并产生一个或多个输出的函数 M文件扩展名.m M文件组成： 函数定义行 H1行 帮助文本 函数体 命令 函数定义行的形式为： function[outputs]=name(input) 例：function [s,p]=sumprod(f,g) 调用 〉〉[s,p]=sumprod(f,g) * * 运算符 算术运算符 关系运算符 逻辑运算符 * * 算术运算符 +-两种运算是一致的，其余用圆点（.）来区分两种运算 例如： C=A*B表示传统意义上的矩阵成积 D=A.*B表示数组乘法，要求A B两个数组大小 相同，逐像素相乘 即：D(I,J)=A(I,J)*B(I,J) * * 运算符 名称 Matlab函数 示例 + 数组和矩 阵加减 plus(A,B) - minus(A,B) .* 数组乘 times(A,B) * 矩阵乘 mtimes(A,B) ./ 数组右除 rdivide(A,B) C=A./B c(i,j)=A (i,j)/ B(i,j) / 矩阵右除 mrdivide(A,B) .\ 数组左除 ldiv