多媒体应用基础_第五节_MATLAB.pptVIP

* * * * 信息管理与信息系统 延边大学 MATLAB图像处理 延边大学 信息管理与信息系统 多媒体信息处理系统设计 1 数字图像的基本概念 图像分辨率：图像的采样点数N×M 256*256 128*128 …… 量化级Q：Q＝28＝256 不同分辨率下的图像 1024 ×1024 → 512 × 512 → 256 × 256 → 128 × 128 → 64 × 64 → 32 × 32 不同灰度级的图像 二值图像与灰度图像 2 图像文件的读写/显示 将图像文件读入内存 imread() 如果图像是彩色图，内存数据为3维矩阵，第3维分别表示红、绿、蓝空间 如果图像是灰度图，内存数据为2维矩阵 数据类型为uint8 将内存中的数据已图片形式保存 imwrite() 内存数据的图像显示 imshow() 可根据图像数据显示灰度图或彩色图 数据类型必须为uint8，如果数据是double，可用uint8()函数转换 另：image()函数可以将矩阵的数据作为图像显示，可接收double、uint16、uint8类型，同时需通过colormap函数进行颜色配置 3 图像的直方图 定义：灰度级的函数，描述的是图像中每种灰度级像素的个数，反映图像中每种灰度出现的频率。横坐标是灰度级，纵坐标是灰度级出现的频率。 直方图计算 依据定义，若图像具有L（通常L=256,即8位灰度级）级灰度，则大小为MxN的灰度图像f(x,y)的灰度直方图hist[0…L-1]可用如下计算获得： 初始化 hist[k]=0; k=0,…,L-1 统计 hist[f(x,y)]++; x=0,…,M-1, y =0,…,N-1 归一化 hist[f(x,y)]/=M*N Matlab函数imhist() 直方图的用处 1）数字化参数 一般一幅数字图像应该利用全部或几乎全部可能的灰度级； 对直方图做快速检查。 2 ）边界阈值选择 使用轮廓线确定简单物体的边界的方法，称为阈值化； 对物体与背景有较强对比的景物的分割特别有用； 基于直方图的图像矫正、分割…… 直方图均衡化 使一输入图像转换为在每一灰度级上都有相同的像素点数（即输出的直方图是平的） 进一步的作用在于图像比较和分割 实际是完成直方图变换：DB＝f(DA) 累积分布密度为 变换 例 %初始化及得到lenna的灰度 lenna=zeros(128,128); lenna_equ=zeros(128,128); histgram=zeros(256); cdf=zeros(256); [lenna,map]=imread(c:\temp\lenna.bmp,bmp); %get histogram for i=1:128 for j=1:128 k=lenna(i,j); histgram(k)=histgram(k)+1; end end %get cdf cdf(1)=histgram(1); for i=2:256 cdf(i)=cdf(i-1)+histgram(i); end %run point operation for i=1:128 for j=1:128 k=lenna(i,j); lenna_equ(i,j)=cdf(k)*256/(128*128); end end %生成直方图均衡化后的lenna图 imwrite(lenna_equ,map,c:\temp\lenna_equ.bmp); 4 图像操作的基本函数 imread() / imwrite() /imshow imresize(A,[mrows ncols],method) 图像缩放，A为缩放因子，method为nearest(默认)(最近邻插值)bilinear(双线性插值)bicubic(三线性插值) imrotate(A,angle,method) 图像逆时针旋转，angle为角度 imcrop(A,rect) 图像剪切，其中rect为[x y width height] imhist(): 图像直方图计算和显示 histeq(): 直方图均衡化 imnoise(): 图像中添加噪声 …… 5 图像处理 图像变换： fft2(傅里叶变换)、dct2(离散余弦变换)、radon(Radon变换)、ifft2/idct2/iradon…. 图像类型转换 rgb2gray(彩色转灰度)、im2bw(转为2值图像) 空域滤波 filter2(线性平滑滤波)、wiener2(维纳滤波)、 Medfilt2(中值滤波)…… 边缘检测edge()：提供6种检测方法 Sobel