系统实验(DSP)--图像的锐化处理、图像的边缘检测

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;y<lmageHeight;y+)for (x=O;x<lmageWidth;x+)(buffery*lmageWidth+x=2*

2、buffery*lmageWidth+x-(buffery*lmageWidth+x-1+b uffer(y-1 )*lmageWidth+x+buffery*lmageWidth+x+1 +buffer(y+1 )*lmageWidth+ x)/4;if(bufferlmageWidth*y+x>255)bufferlmageWidth*y+x=255;else if (bufferlmageWidth*y+x<0) bufferlmageWidth*y+x=O;运行结果：锐化前锐化后分析：从上面两幅图可以看出锐化后的图像轮廓变得明显，且噪声变得强 烈。3、 汇编程序及运行结果程

3、序：lmageAcuteSub(lmageWidth,lmageHeight,buffer_red);lmageAcuteSub(lmageWidthJmageHeight,buffer_green);lmageAcuteSub(lmageWidthJmageHeight,buffer_blue);.mmregs.def _lmageAcuteSub .textJmageAcuteSub:mov tO,brd ;IMAGE WIDTH movt1,brcO ;IMAGE HEIGHT rptb yjooprptb xjoop mov *ar0(#1 ),ac1 ;f(x+1 ,y)add *a

4、r0(#-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 ,ac0<0 sub #255,ac0,ac1 bcc branch2,ac1>0 mov acO/arO+ b x_loopbranchi: mov#O/arO+ b x_loopbranch2: mov #255,*ar0+ xjoop: nop yjoop: nopRET运行结果：锐化前锐化后

5、分析：可以看出汇编的结果和 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|WTg(x,y)- J 0 else实验中的阈值设为120o2、

6、 C程序及运行结果程序：Edge_Grey_lmage()(int x,y；int delta_x,delta_y;for (y=O;y<lmageHeight;y+)for (x=O;x<lmageWidth;x+) delta_x=abs(buffer_greylmageWidth*(y-1)+x+1+buffer_greylmageWidth*y+x +1 +buffer_greylmageWidth*(y+1 )+x+1 -buffer_greylmageWidth*(y-1 )+x-1 -b uffer_greylmageWidth*y+x-1-buffer_greylm

7、ageWidth*(y+1 )+x-1 )/3;delta_y=abs(buffer_greylmageWidth*(y-1)+x-1+buffer_greylmageWidth*(y-1 )+x+buffer_greylmageWidth*(y-1)+x+1-buffer_greylmageWidth*(y+1 )+x-1 -buffer_greylmageWidth*(y+1 )+x-buffer_greylmageWidth*(y+1 )+x+1 )/3;if(delta_x+delta_y)>=Threshhold) buffer_orglmageWidth*y+x=255; e

8、lse buffer_org(lmageWidth*y+x=0;运行结果：边缘处理前边缘处理后分析：从上面两幅图可以看出，原图像的边缘经过处理后可以很好的用白 线来勾勒出，没有边缘的地方都用黑色表示，可以很清楚的从处理后的图 片看出斑马的形状。3、 汇编程序及运行结果程序：lmageEdgeSub(lmageWidthJmageHeight5buffer_org);.mmregs.def JmageEdgeSub.bss TEMP,1.textJmageEdgeSub:mov tO,BRC1mov t1,BRCO mov xar0,ac3 mov #256«#11 ,ac2 mov

9、ac2,xar1RPTB Y_LOOPRPTB X_LOOPmov*arO+,*ar1+;将art)的内容复制到ar1里面X_LOOP: nopY_LOOP: nopmov tO,brc1 ;IMAGE WIDTH movt1,brcO ;IMAGE HEIGHT mov ac2,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

10、),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(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,ac0>0 mov #O,*arO+ b xjoopbranchi: mov#255,*ar0+ xjoop: nopyjoop:no PRETDIVION:amov #080000h,XAR7mov #3,*ar7mov #1