数字图像处理大作业

1、数字图像处理1 .图像工程的三个层次是指哪三个层次？各个层次对应的输入、输出对象分别是什么？图像处理特点：输入是图像，输出也是图像，即图像之间进行的变换。图像分割特点：输入是图像，输出是数据。图像识别特点：以客观世界为中心，借助知识、经验等来把握整个客观世界。“输入是数据，输出是理解。2 .常用的颜色模型有哪些（列举三种以上）？并分别说明颜色模型各分量代表的意义。RGB（红、绿、蓝）模型CMY（青、品红、黄）模型HSI（色调、饱和度、亮度）模型3 .什么是图像的采样？什么是图像的量化？1 .采样采样的实质就是要用多少点来描述一幅图像，采样结果质量的高低就是用前面所说的图像分辨率来衡量。简单来讲

2、，对二维空间上连续的图像在水平和垂直方向上等间距地分割成矩形网状结构，所形成的微小方格称为像素点。一副图像就被采样成有限个像素点构成的集合。例如：一副640*480分辨率的图像，表示这幅图像是由640*480=307200个像素点组成。2 .量化量化是指要使用多大范围的数值来表示图像采样之后的每一个点。量化的结果是图像能够容纳的颜色总数，它反映了采样的质量。针对数字图像而言：采样决定了图像的空间分辨率，换句话说，空间分辨率是图像中可分辨的最小细节。量化决定了图像的灰度级，即指在灰度级别中可分辨的最小变化。数字图像处理（第三次课）调用图像格式转换函数实现彩色图像、灰度图像、二值图像、索引图像之间

3、的转换。图像的类型转换：对于索引图像进行滤波时，必须把它转换为RGB图像，否则对图像的下标进行滤波,得到的结果是毫无意义的；灰度图像到索引图像X,map=gray2ind(I,n)将灰度图像I转化成索引图像真彩图像到二值图像BW=im2bw(RGB,level)基于level将彩色图像RGB转化为二值图像真彩色图像到灰度图像I=rgb2gray(RGB)将真彩图像RGB转换为灰度图像I真彩色图像到索引图像X,map=rgb2ind(RGB,n)X,map=rgb2ind(RGB,tol)直接将真彩图像RGB转化为索引图像X并返回颜色图图像抖动转换图像X=dither(RGB,map)BW=di

4、ther(I)通过颜色图将RGB图像转化为近似的索引图X;或把灰度图像I转化为一值图像2.用MATLAB完成灰度图像直方图统计代码设计。1functionladylenna()2a=imread('C:lenna.jpg');%W取图像至工作空间3imfinfo('C:lenna.jpg')%查看图像信息4subplot(2,2,1);5imshow(a);1庐改瘦如为身lad小伊士rtM6title('lady-lenna');7ifisrgb(a);8b=rgb2gray(a);%RGB转换为灰度图像9end10subplot(2,2,2);

5、11imshow(b);%显示图像12title('ladygaga-lenna');13m,n=size(a);%返回图像大小14e=zeros(1,256);15fork=0:25516fori=1:m17forj=1:n18ifa(i,j)=k19e(k+1)=e(k+1)+1;%灰度值相同的进行累加20end21end22end23end24subplot(2,2,4);25bar(e);%画图像的灰度直方图26title('灰度直方图)；27c=imrotate(a,20);28%图像的旋转29subplot(2,2,3);30imshow(c);时也巴，鞘S

6、LMJJLlbTOlCii旧理Rrm数字图像处理(第四次课)编写matlab函数，实现在医学图像中数字减影血管造影。即实现两幅图像的差值图像。1实验代码2I=imread(图像处理1.jpg');%原始图像3imshow(I);4BK=imread(图像处理2.jpg');%背景5figure;6imshow(BK);7Ip=imsubtract(I,BK);%图像减背景8figure;9imshow(Ip,)实验结果：*的LFE3文件CD“ES嚏八(D工MIL算面©虫匚阳菱先M|口自心电堂31图像剪切:12345678910任意选择一幅图片，实现图像的镜像。这题其实

7、就是一个函数imcrop()的使用而已closeall;clearall;clc;img=imread(图像处理3.jpg');imshow(img);imgn=imcrop();%这里把要裁剪的图像框出来figure;imshow(imgn)%显示裁剪的图像实验结果:充I岫Al实现图像的镜像:12实验代码clearall;3closeall;4clc;5A=imread(图像处理4.jpg');6height,width,dim=size(A);7tform=maketform('affine',-18B=imtransform(A,tform,'ne

8、arest'9像10tform2=maketform('affine',111C=imtransform(A,tform2,'nearest12像13subplot(1,3,1),imshow(A);14title('原图像');15subplot(1,3,2),imshow(B);16title('水平镜像');17subplot(1,3,3),imshow(C);18title('坚直镜像);19subplot(2,2,3);20imshow(c);000;);0;010;width01);%B中存储的是经过水平镜像变

9、换后的图0-10;0height1);)；%C中存储的是经过竖直镜像变换后的图实验结果:/Figure1X文件（f）编号CD®#（V）SAW工具莫国窗口的智助（Hi、f：口亡|»|04仃!8夏/，|国|国|口原图像水平镜像竖直憧像数字图像处理（第五次课）数字图像增强有哪些方法？忖像增强可分成两大类：频率域法和空间域法。前者把图像看成一种二维信号，对其进行基于二维傅里叶变换的信号增强。采用低通滤波（即只让低频信号通过）法，可去掉图中的噪声；采用高通滤波法，则可增强边缘等高频信号，使模糊的图片变得清晰。后者空间域法中具有代表性的算法有局部求平均值法和中值滤波（取局部邻域中的中间

10、像素值）法等，它们可用于去除或减弱噪声。数字设一幅图像有如下图所示直方图，若对其进行直方图均衡处理，求出增强后的灰度分布,并给出处理后的直方图。0.0.0.0.0.1740.0880.0860.080.0680.0580.0620.3840.01234567实验代码:1clc2clearall3prob0=.174.088.086.08.068.058.062.3844zero0=zeros(1,8);5fori=1:86ifi=17zero0(1)=prob0(1);8else9zero0(1,i)=zero0(1,i-1)+prob0(1,i);10end11end12fori=1:813

11、ifi=114zero0(1)=prob0(1);15else16zero0(1,i)=zero0(1,i-1)+prob0(1,i);17end18end1920zero1=zeros(1,8);21fori=1:822forj=0:723a=abs(zero0(i)-j/7);24ifa<1/1425zero1(j+1)=zero1(j+1)+prob0(i);26end27end28end29b=find(zero1=0);30zero1(b)=;31prob1=zero1;32subplot(121)33bar(prob0)34title（'均衡化前）35subplot(

12、122)36bar(prob1)37title（'均衡化后）38title（'均衡化前）39subplot(122)40bar(prob1)41title（'均衡化后）实验结果:匚X.Figure1数字图像处理(第六次课)利用help,学习imfilter和fspecial函数的用法，编写程序创建3*3的均值滤波器和半径为5的圆形滤波器对图像进行平滑处理；MATLAB提供了基于卷积的图像滤波函数imfilter进行平滑滤波，该函数用指定的滤波器模板对图像进行运算。其调用格式为：B=imfilter(h,A)B=imfilter(h,A,shape)其中B=imfilte

13、r(h,A)返回图像A经算子h滤波后的结果，参数shape指定滤波的计算范围，即：shape='full'时，作边界补零shape='same'时，返回图像B与输入图像A大小相同shape='valid'时，不考虑边界补零，只计算有效输出部分。此外，也可以使用imfilter函数进行卷积或相关性滤波来实现图像平滑。该函数的调用格式为：B=imfilter(A,H)B=imfilter(A,H,option1,option2,)其中B=imfilter(A,H)返回图像A经算子h滤波后的结果；B=imfilter(A,H,option1,optio

14、n2,)是根扼琳危的参数实现图像滤波。option参数可以有下列取值：边界填充选项对一幅图像进行滤波时，零填充可能会导致图像被一个黑框围绕，为了消除零填充的人工痕迹，imfilter函数支持三种可选的边界填充方法：'symmetric边界对称'replicate边界复制，缺省值'circular边界循环输出尺寸选项其选项与imfilter函数的shape参数相同。滤波选项corr使用相关性来进行滤波，缺省使用此方法conv使用卷积方法实现滤波学习ordfilt2函数的用法，编程调用该函数利用不同的统计排序滤波器实现图像滤波。中文理解函数名就是顺序滤波函数，后面的2表示的

15、是2维滤波，常用的调用形式如下：1.B=ordfilt2(A,order,domain)这种形式中，A为被滤波的矩阵，order表示选取由domain确定的邻域中第order个元素替换A中的元素，存入B中。示例如下：Y=ordfilt2(X,5,ones(3);相当于3*3窗口的中值滤波Y=ordfilt2(X,1,ones(3);相当于3*3窗口的最小值滤波需要注意的是，domain矩阵中可能会有某些位置为零(不为零的位置为1),这种为零的位置不算在顺序统计的领域之中。也就是说：Y=ordfilt2(X,1,010;101;010)表示用窗口表示的邻域中不为零的位置的像素值中选择最小的替代原

16、来的像素值。2.B=ordfilt2(A,order,domain,S)暂时没搞明白S是做啥用的，好像一时半会儿也用不到，先记录着，以后想明白了再更新利用help,学习edge函数的用法，并思考如何利用edge函数使用不同锐化算子对图像进行锐化。1)利用matlab提供的edge函数，选择三种边缘检测算子，分别对图像Lena.tif和baboon.tif进行边缘检测显示检测结果。2)对Lena.tif和baboon.tif分别添加高斯和椒盐噪声，然后对有噪声的图像进行边缘检测，显示检测出的边缘图像1clearall%青空工作空间2I=imread('lena.jpg');%W入

17、图像3BW=edge(I,'sobel');%sobel算#析4BW2=edge(I,'log');%&普拉斯算子分析5BW3=edge(I,'prewitt');%W度算子分析6J1=imnoise(I,'gaussian',0,0.02);播加高斯噪声7J2=imnoise(I,'salt&pepper',0.02);9加椒盐噪声8BW4=edge(J1,'sobel');%高斯噪声后sobel算9子分析10BW5=edge(J2,'log');%M盐噪声后lap

18、lace11算子分析12BW6=edge(J1,'log');%高斯噪声后laplace13算子分析14BW7=edge(J2,'sobel');%W盐噪声后sobel算15转析16171819202122232425262728293031323334BW8=edge(J1,'prewitt');%高斯噪声后梯度算子分析BW9=edge(J2,'prewitt');%椒盐噪声后梯度算子分析figure;%显示所需图像subplot(2,2,1);imshow(I);xlabel('原图')subplot(2,2,2);imshow(BW);x