《系统实验(DSP)--图像的锐化处理、图像的边缘检测.docx》由会员分享，可在线阅读，更多相关《系统实验(DSP)--图像的锐化处理、图像的边缘检测.docx（6页珍藏版）》请在三一文库上搜索。
1、DSP实验报告一、图像的锐化处理(高通滤波处理)1、 实验原理处理模板如下：f y 1M u -a t -a |a =0.250.a0 J对应数学表达式：g(x, y) =(! *4a)f ( x, y)-a f ( x-1, y)* f (x, y -1) + f ( x*1, y)* f ( x, y*i)2、 C程序及运行结果程序：Acute_RGB_lmage(int *buffer)(int x,y;for (y=O;ylmageHeight;y+)for (x=O;x255)bufferlmageWidth*y+x=255;else if (bufferlmageWidth*y+x
2、0) bufferlmageWidth*y+x=O;运行结果：锐化前锐化后分析：从上面两幅图可以看出锐化后的图像轮廓变得明显，且噪声变得强 烈。3、 汇编程序及运行结果程序：lmageAcuteSub(lmageWidth,lmageHeight,buffer_red);lmageAcuteSub(lmageWidthJmageHeight,buffer_green);lmageAcuteSub(lmageWidthJmageHeight,buffer_blue);.mmregs.def _lmageAcuteSub .textJmageAcuteSub:mov tO,brd ;IMAGE W
3、IDTH movt1,brcO ;IMAGE HEIGHT rptb yjooprptb xjoop mov *ar0(#1 ),ac1 ;f(x+1 ,y)add *ar0(#-1),ac1 ;f(x-1,y)add *ar0(#-250),ad ;f(x,y-1) add *ar0(#250),ac1 ;f(x,y+1) sfts ac1 ,#-2mov *arO#1 ,acO;2f(x,y)sub ac1,acObcc branchi ,ac00 mov acO/arO+ b x_loopbranchi: mov#O/arO+ b x_loopbranch2: mov #255,*ar0
4、+ xjoop: nop yjoop: nopRET运行结果：锐化前锐化后分析：可以看出汇编的结果和 C程序的结果是一致的。二、图像的边缘检测1、 实验原理X方向和y方向的梯度分别为：1 0 1F 11111| Ji1I- 1 0 1y - I 000hod卜 1 - 1 r总梯度幅度|冬|葬 f (x*-1, y-1)+ f ( x, y)+ f (x -*1, y 1) -f (xT,y 3) - f (x T, y) - f (x T, y *1)l%l=Uf(Z,yrf(x,yr)+f(x1y 1) -f(xlyl) - f(x,y*1)-f(x*1,y*1)2f 1 if |Af|W
5、Tg(x,y)- J 0 else实验中的阈值设为120o2、 C程序及运行结果程序：Edge_Grey_lmage()(int x,y；int delta_x,delta_y;for (y=O;ylmageHeight;y+)for (x=O;x=Threshhold) buffer_orglmageWidth*y+x=255; else buffer_org(lmageWidth*y+x=0;运行结果：边缘处理前边缘处理后分析：从上面两幅图可以看出，原图像的边缘经过处理后可以很好的用白 线来勾勒出，没有边缘的地方都用黑色表示，可以很清楚的从处理后的图 片看出斑马的形状。3、 汇编程序及运行
6、结果程序：lmageEdgeSub(lmageWidthJmageHeight5buffer_org);.mmregs.def JmageEdgeSub.bss TEMP,1.textJmageEdgeSub:mov tO,BRC1mov t1,BRCO mov xar0,ac3 mov #256#11 ,ac2 mov ac2,xar1RPTB Y_LOOPRPTB X_LOOPmov*arO+,*ar1+;将art)的内容复制到ar1里面X_LOOP: nopY_LOOP: nopmov tO,brc1 ;IMAGE WIDTH movt1,brcO ;IMAGE HEIGHT mov a
7、c2,xar1 mov ac3,xar0 rptb yjoop rptb x_loopmov *ar1+,ac0mov *ar1(#-149),ac0;f(x+1 ,y-1) add *ar1(#1 ),acO ;f(x+1 ,y) add *ar1(#151),ac0 ;f(x+1,y+1) sub *ar1(#-151),ac0;f(x-1,y-1) sub *ar1(#-1),ac0 ;f(x-1,y) sub *ar1(#149),ac0;f(x-1 ,y+1) abs acOmov *ar1(#-151),ac1 ;f(x-1,y-1) add *ar1(#-150),ac1 ;f(
8、x,y-1) add *ar1(#-149),ac1 ;f(x+1,y-1) sub *ar1(#149),ac1 ;f(x-1 ,y+1) sub *ar1(#150),ac1 ;f(x,y+1) sub *ar1(#151 ),ac1 ;f(x+1,y+1) abs adadd ad,acO call DIVIONsub#120,ac0,ac0 bcc branchi,ac00 mov #O,*arO+ b xjoopbranchi: mov#255,*ar0+ xjoop: nopyjoop:no PRETDIVION:amov #080000h,XAR7mov #3,*ar7mov #16,tOmov ac0,ac1 branch2: subc *ar7,ad,ac2mov ac2,ac1sub #1,t0bcc branch2,t00mov ac1,t1mov t1 ,acORET运行结果:边缘处理前E Graphical D.o sfl H边缘处理后分析：可以看出汇编的结果和 C程序的结果是一致的。