数字图像处理学2

• 当用imshow显示另一副图像g时，为保持第 一副图像并同时显示第二副图像，可以使用 figure函数 >>imshow(f),figure,imshow(g)
2.4 图像保存
• imwrite函数
语法：imwrite(f,‘filename’)
>>imwrite(f,‘patient10_run1.tif’)
parameter: none,packbits,ccitt [colres rowres]:行列分辨率 >>res=round(200*2.25/1.5); >>imwrite(f,‘sf.tif’,‘compression’,… ‘none’,‘resolution’,res)
• print函数
imwrite只适用于JPEG图像保存的形式：
imwrite(f,‘filename.jpg’,‘quality’,q)
质量参数q是一个在0到100之间的整数
• imfinfo函数 给出压缩图像的详细信息
语法：imfinfo filename
>>imfinfo bubbles25.jpg
压缩比的计算：
2.6 数据类与图像类型间的转换
• 数据类之间的转换 B=data_class_name(A) >>B=double(A); >>D=uint8(C); • 图像类与类型间的转换 转换函数包括： im2uint8,im2uint16,mat2gray,im2double, im2bw
f= -0.5 0.5 0.75 1.5 >>g=im2uint8(f) g= 0 128 191 255 >>h=uint8([25,50;128,200]); >>g=im2double(h) g= 0.0980 0.1961 0.4706 0.7843
>>C3=A(:,3) --C3=A(1:3,3) C3=3 6 9 >>R2=A(2,:) R2=4 5 6 >>T2=A(1:2,1:3) T2=1 2 3 4 5 6
>>B=A; >>B(:,3)=0 B=1 2 0 4 5 0 7 8 0
• 使用end操作符选取一个二维数据块:
>>A(end,end)
• size函数
给出一副图像的行数和列数。
>>size(f)
ans=1024
1024
>>[M,N]= size(f);
%返回图像的行数和列数
• whos函数 显示数组的信息 >>whos(f) Name Size f 1024*1024
Byte 1048576
Class uint8 array
2.3 图像显示 • imshow函数 语法：imshow(f,G) G为显示该图像的灰度级数。 imshow(f) imshow(f,[low high]) imshow(f,[ ])
• 线性分隔函数linspace
语法: x=linspace(a,b,n)
产生一个含有n个元素的行向量x，这n个
元素之间线性地间隔开并且包含a和b。 >>x=linspace(1,10,4) x=
1
4
ห้องสมุดไป่ตู้
7
10
• 矩阵索引 矩阵的表示： A=[1 2 3;4 5 6;7 8 9]; 从矩阵中选取元素: >>A(2,3) ans=6 使用冒号操作符选取一个二维数据块:
• 把矩阵索引转化为列向量: >>T2=[1 2 3; 4 5 6]; >>V=T2(:) V =1 4 2 5 3 6
• 使用SUM函数对矩阵求和: >>A=[1 2 3;4 5 6;7 8 9]; >>s=sum(A(:)) s=45 >>s=sum(sum(A)) s=45
• 例2.5 使用数组索引进行简单的图像操作 >>f=imread('rose.tif');
>>comressed_ratio=image_bytes/compress ed_bytes
imwrite只适用于tif图像保存的形式：
imwrite(g,‘filename.tif’,‘compression’,…
‘parameter’,‘resolution’,[colres rowres])
>>f=[1,2;3,4]; >>g=mat2gray(f) g= 0 0.3333 0.6667 1.0000 >>gb=im2bw(g,0.6) gb= 0 0 1 1 >>gbd=im2double(gb) gbd= 0 0 1 1
MATLAB支持嵌套语句，上述操作可用语 句表示为： >>gbd=im2double(im2bw(mat2gray(f),0.6));
char：字符，每像素2B logical：逻辑数据，每像素1B
2.6 图像类型
• • • • • MATLAB和IPT支持的图像类型： 亮度图像 二值图像 索引图像 彩色图像 二值图像：一个取值只有0和1的逻辑数组 >>B=logical(A) %创建一个逻辑数组B >>islogical(B) %测试数组B是否为逻辑数组
保存图像窗口的内容
>>print –f1 –dtiff –r300 hi_res_rose
2.5 数据类
• MATLAB和IPT支持的数据类： double：双精度浮点数，每像素8B unit8 :无符号8Bits整数，每像素1B uint16：无符号16Bits整数，每像素2B uint32：无符号32Bits整数，每像素4B int8： 有符号8Bits整数，每像素1B int16：有符号16Bits整数，每像素2B int32：有符号32Bits整数，每像素4B single：单精度浮点数，每像素4B
• 选择数组的维数 operation(A,dim) >>k=size(A,1); >>plot(f(size(f,1)/2,:)) d=ndims(A)
2.9 一些重要的标准数组 • • • • • • • zeros(M,N) ones(M,N) ture(M,N) false(M,N) magic(M,N) rand(M,N) randn(M,N)
ans=9
>>A(2:end,end:-2:1)
ans=6 4
9 7
• 使用向量作为矩阵的索引:
>>A=[1 2 3;4 5 6;7 8 9]
>>E=A([1 3],[2 3])
E =2
8
3
9
• 使用逻辑数组寻址: >>D=logical([1 0 0;0 0 1;0 0 0] ); >>A(D) ans =1 6
2.8 数组索引 • 向量索引 向量即一维数组 >>v=[1,3,5,7,9]; >>v(2) ans=3 >>w=v.'; >>v(1:3) ans=1 3 5
>>v(3:end) ans=5 7 9 >>v(:); %产生一个列向量 >>v(1:end); %产生一个行向量 >>v(1:2:end) ans=1 5 9 >>v(end:-2:1) ans=9 5 1
（Width*Height*BitDepth/8）/压缩图像字节数
• 压缩比另一种计算方法： >>k=imfinfo(‘bubbles25.jpg’); %将函数imfinfo显示的信息域捕获到结构变 量k中 >>image_bytes=k.Width*k.Height*k.BitDept h/8; >>compressed_bytes=k.FileSize;
第2章 MATLAB图像及操作 2.1 MATLAB图像 • MATLAB图像的表示 MATLAB图像使用(r,c)来表示行与列； MATLAB图像坐标系统的原点位于(1,1)处。 MATLAB图像的矩阵表示：
f (1,1) f (1, N ) f (2,1) f (2, N ) f f ( M ,1) f ( M , N )
• MATLAB的图像文件格式 JPEG格式（.jpg)：联合图像专家组 TIFF格式(.tif)：加标识的图像文件格式 BMP格式(.bmp)：Windows 位图 GIF格式(.gif)：图像交换格式 PNG格式(.png)：可移植的网络图像 XWD格式(.xwd)：X Window转储
2.2 图像的读取 • imread函数 语法：imread(‘filename’) >>f=imread(‘chestxray.jpg’); %将当前目录下的chestxray.jpg图像读入 到数组f中。 >>f=imread(‘D:\myimages\chestxray.jpg’); %将D盘上的chestxray.jpg图像读入到数 组f中。
• >>fa=f(end:-1:1,:); • >>imshow(fa)
• >>fb=f(257:758,257:758); >>figure,imshow(fb)
• >> fc=f(1:2:end,1:2:end); >> figure,imshow(fc)
• >> plot(f(512,:))