图像存储与显示•图像存储•图像显示1、有格式图像显示•对有格式图像，可用imread（）函数读入，用imshow（）函数显示。例：I＝imread(‘d:\image\x.bmp‘)；或I＝imread（’d:\image\x’,’bmp’）；路径格式figure(1)；imshow（I）；％显示图像I2、无格式图像显示（以二进制格式存储）•首先要用fopen（）函数打开某一路径下的文件，然后用fread（）函数读到一个变量里，该变量为矩阵变量。•例：fid＝fopen（‘d:\img\lena.img’,’r’）;%fid为文件句柄data＝（fread（fid，[256，256]，‘uint8’））’;figure(1);images(data,[255]);%图像为256级灰度colormap(gray)；％显示灰度图像axisimage％对图像加坐标。•若不能直接对data进行处理，可以进行转换后在处理；转换语句：data=uint8(data);3、同屏显示多个图像•可用subplot(m,n)将图形窗分为m*n个子窗口，然后取第一、第二…子窗口显示不同的图像，实现同屏显示多个图像。例如：•figure(1)；%取2×2个子屏中的第一个子屏subplot（2,2,1);……%显示第一个图像imshow(I1)；……%取2×2个子屏中的第四个子屏subplot(2,2,4)；……%显示第四个图像imshow(I4)；•Size（）函数•Zeros（）函数•Fft2(),ifft2()函数•Imhist()函数•Histeq()函数•Imrotate()函数•Imnoise()函数•Edge()函数•Title()函数•Xlable(),Ylable()函数•Size（）函数获取图像矩阵大小。一般是应用于有格式图像，因无格式（二进制）图像的大小在读入时已知。例：I＝imread(‘d:\img\radar’,’bmp’);%读入图像[x,y]＝size(I)；%获取图像大小得到x,y的值，该图像大小就是X×Y。•Zeros（）函数零矩阵函数。例：I=zeros(100,100)；%I为100×100的零矩阵，矩阵％中元素全为零。Imshow(I)；%显示一个100×100的黑方块。例如作业题1，生成一个外边黑中间一块是白的图像，可以先生成一个全黑的图像，然后在中间作一个双重循环，赋象素值为255或1。I1＝zeros(128,128)；%生成一个128×128的全黑图像forI=38:1:90forj=58:1:70I1（i,j）＝255；%或I1(i,j)=1；endendimshow（I1）；%I1即为所求图形。•fft2(),ifft2()函数fft2()函数为二维快速傅立叶变换函数；ifft2()函数为二维逆快速傅立叶变换函数。对一幅图像进行傅立叶变换后，得到它的频谱。例：I2=fft2(I1)；%对图像I1进行二维快速傅立叶变换Imshow(I2);%显示频谱图。I3=ifft2(I2)；%对图像I2进行二维逆快速傅立%叶变换,得到原图像•Imhist()函数图像直方图函数。例：imhist(I)；%显示图像I的直方图；•Histeq()函数直方图均衡化函数。例：I1=histeq(I)；%对图像I进行直方图均衡化。Imshow(I1);%显示均衡化后的图像I1•Imrotate()函数：旋转图像函数；格式：J=imrotate(I,angle,method)I：被旋转图像J：旋转后的图像angle：旋转角度method：可为'nearest‘、'bilinear‘、'bicubic'例：I=imread(‘ic.tif’);%读入图像ic.tifJ＝imrotate(I,45,’bilinear’);%对图像I旋转45度imshow(I)；%显示原图Ifigure,imshow(J)；%显示旋转后的图像J•Imnoise()函数：给图像增加噪声。格式：J=imnoise(I,’噪声类型’,参数);I:待加噪声图像；噪声类型：高斯噪声、盐噪声等；参数：噪声密度(0--1)；J：加入噪声的图像。例：J=imnoise(I,‘salt&pepper’，0.02)；%给图象I增加盐噪声•Edge()函数边缘检测函数。J＝edge(I,‘检测算子‘)；％J是对图像I用某算子%进行边缘检测后的边缘图；%J是二值图像，即黑、白二色。例：J＝edge(I,’Roberts’)；%用Roberts算子对图像I进行边缘检测imshow(J);%显示边缘图•Title()