采用可编程逻辑器件实现Sobel边缘检测算法的研究和分析.docx

采用可编程逻辑器件实现Sobe1边缘检测算法的研究和分析引言随着设计复杂度的增加，使用IP核已经成为一种常用的设计方法。Quartus1I软件提供的Megaf1mt-巨OnS是基于A1tera底层硬件结构最合理的成熟应用模块，在代码中使用MegafImetiOnS这类IP资源，不但能将设计者从繁琐的代码编写中解脱出来，更重要的是，在大多数情况下MegafUnCtiOnS的综合和实现结果比用户编写的代码更优。而且只需要简单地设置选取宏功能模块的相关参数就可以在程序中调用它们，因此宏功能模块的使用也十分方便。QUartUSI1的Mega-functins中包含有算术运算(Arithmetic)、逻辑门(Gates).1/0.存贮器(St2Tage)等四个系列，可以根据系统设计需求灵活选用。在实时图像处理中，用Sobe1算子进行物体的边缘检遮是经常用到的篁法，由于对处理速度要求较高，因此用纯软件的方法很难达到要求。而FPGA对同时可完成的处理任务几乎没有限制，适合高速、并行信号处理,并且FPGA密度高、容量大，有内置存储器、容易实现，所以FPGA广泛用于实时图像处理系统中。在此通过调用基于RAM的移位寄存器a1tshifttaps、可编程乘加器a1tmu1tadd,可编程多路并行加法器ParaIIe1aCk1和参数化绝对值运算模块1pm_abs,实现了基于FPGA的SobeI边缘检测。最后给出设计系统的仿真结果，通过与MaHab仿真结果相比较,可以看出该设计获得了很好的边缘检测效果。1 Sobe1边缘检测算法图1给出了Sobe1边缘检测算法框图。从图中可以看出对一副图像进行Sobe1边缘检测时首先要利用Sobe1算子计算出水平梯度和垂直梯度，然后再把两个方向的梯度结合起来，最后应用门限处理模块判断图像边缘并输出边缘检测结果。国便输入用1Sobe1边缘检测算法框图图2(a)为一副图像的3X3区域，图2(b)和图2(C)分别为SObCI算子的X方向(垂直方向)梯度算子和y方向(水平方向)梯度算子。当采用Sobe1算子对图2(a)所示的3X3区域做梯度计算时，可得标记为z5的像素点X方向梯度和y方向的梯度分量分别为：然而实际执行时，为了运算方便可以采用公式Gx=+2=8+。)一(N1+2zt+zi)Gy=(zi+2z6+方)(ZI+2zi+z)f=血+G中梯度的计算嘉要Gx,Gy这两个分线公式联合使用.fGx1+1GyI对梯度进行近似。图2一副图像的3X3区域(Z是灰度值)和Sobe1算子该设计在门限处理时，采用基本全局门限：当某像素点(,y)的梯度值XXXXf(,y)大于或等于设定的门限T时，规定该点的灰度值为255,反之则为0。即：g(j,y)=2550V(xj)T/(W)VT2 Sobe1边缘检测的硬件实现根据图1所示的Sobe1边缘检测算法框图，可得FPGA硬件实现Sobe1边缘检测时应该包含梯度计算模块和门限处理模块。此外在硬件实现时还要采用图像数据缓冲模块对图像做缓冲处理，以便进一步对图像数据做模板处理。门限处理模块可以通过编写Veri1ogHD1代码实现。以下主要介绍图像数据缓冲模块和梯度计算模块。2. 1图像数据缓冲模块器在图像的空域滤波中，为了得到3X3的方形模板窗，常使用FIFO(FirstInFirstOut)模块作为图像数据的缓冲器。这里通过应用基于RAM的移位寄存器宏模块a1tshifttaps实现了同样的功能，而且还省去了一些控制信号,使用十分方便。a1tshift_taps宏功能模块是一*个可配置的、具有抽头(Taps)输出的移位寄存器，每个抽头在移位寄存器链的指定位置输出数据。图3(a)和图3(b)分别为定制的8位输入/8位输出、3抽头，且相邻两个抽头相距256个寄存器的a1tshiftjaps功能模块及其内部寄存器链结构图，图3(b)中的BufferO,Buffer1,BUffer2分别为由256个8位移位寄存器构成的寄存器链。当图像的第N行数据在像素四如同步下从Shiftin7：0端输入到a1tshift-taps的BUffero后,随着第N+1行图像数据输入到BUffCro中,第N行的图段数据依次存入BUffer1中，而当第N十2行图像数据存入BUfferO后，BUffer1和BUffer2中分别存放的是第N+1行和第N行的图像数据，从而实现缓冲图像数据的功能。这样在像素时钟的同步下，第N+2,N÷1,N行的同一列数据分别从tapx7：0,tapIx7：0,tap2x7：0端输出给梯度计算模块，进行梯度计算。图3ahshifttaps。功能模块及其内部存存器能结构图2.2梯度计算模块Sobe1边缘检测中，图像像素点梯度的计算可由So-be1算子与图像像素卷积运算的输出经梯度计算公式计算获得。图5是图4所示的3X3空间滤波模板与图2(a)所示的3X3图像区域卷积运算的原理图。从中可以看出，为了实现卷积运算需要做乘法和加法运算，之前的文献中大都采用分立的D触发器和加法器以及乘法器来完成卷积运算，它的结构复杂。在此采用可编程乘加器a1tmu1t_add模块和可编程多路并行加法器Pa1a1Ie1add模块实现卷积运算，大夫简化了设计。ffi43X3/波模板图5卷积运算原理府可编程乘加器a1tmu1t_and可以接收多组数据输入，各组数据相乘后相加或相减作为结果输出。而且a1tmu1tadd在使用时可以根据需要设置乘法器个数、输入/输出数据格式、流水线控制时钟等参数，同时它还支持输入数据内部移位功能。使用可编程多路并行加法器Para11e1add模块时，用户可以自由设计输入数据位宽，累加数据个数，定义累加输入薮据类型，模块最终自动生成适当位宽的数据输出。而且使用Para11e1add模块做加法运算时，可以通过指定时钟延时以实现流水线设计，从而改善电路的性能，提高整个系统的工作频率。图6为实现梯度计算而定制的a1tmu1t_addO模块，该乘加器包含3个乘法器，1个加法器而且为了改善电路的性能和提高系统的工作频率，在该模块中使用了寄存器。同时由于该模块使用了数据的内部移位功能，当某行的图像数据从dataaj)7：0依次输入时，经过3个时钟周期后，相邻的3个数据分别与固定的模板系数datab0,datab1,datab2相乘，并将结果送给加法器完成加法运算。当采用三个这样类似的a1tmu1t_add模块并联时，便可实现卷积运算，而当卷积模板采用Sobe1算子时，就可以获得像素点的水平梯度和垂直梯度。dtaa.0(70)d«tab_0(70)datab1(70)daiab_2(70AItnW1Iaddresu1t11170)图6a1tmu1tadd功能模块vg+gj,在获得GX和Gy后，还要通过公式，-'”:才可以计算出对应像素点的梯度值。设计中可以使用QUartUSI1提供的1pm_abs宏功能模块完成绝对值的计算。图8梯度计算模块的Quartus仿真结果图7为3X3图像区域，图8是基于a1tmu1t_add模块、ParaIIe1aCkI模块和1pm_abs模块构建的梯度计算模块对进行Sobe1算子处理的QUartUS11仿真结果。由于采用了流水线设计，在第7个时钟的上升沿从GX端，Gy端分别输出有效的X方向和y方向梯度值，并在第9个时钟的上升沿从Gf端输出图7的IG.I+IGvI中心像素点的梯度，这个仿真结果与利用公式j-'7计算出的梯度值完全相同。3仿真结果为了更加直观地验证该系统的边缘检测效果，在系统功能验证时采用Mat1ab和Mode1豆业进行混合仿真，仿真结果如图9所示。图9(a)为一幅256X256的原始图像，图9(b)为设定门限T为125的MatIab程序对原始图像的SobeI边缘检测结果。图9(c)为该设计的SobeI边缘检测结果。比较图9(b)和图9(c),可以看出该设计得到的结果几乎与Mat1ab的仿真结果完全一样，即该设计取得了很好的边缘检测效果。OiKKiff1ft(b)MatIabttA用9仿真结果(c)4工设计4结语该设计利用QuartusII软件提供的可编程乘加器a1tmu1t.add模块、可编程多路并行加法器para11e1_add模块和绝对值计算模块1pm_abs,及使用Veri1ogHD1设计的门限处理模块和其他相应的控制模块，完成基于FPGA的Sobe1边缘检测的硬件设计。该方法既避免了自己编写大量程序代码的繁琐,又获得很好的综合和实现结果。最后通过与Mat1ab的仿真结果相比较证明了本设计可以有效地实现Sobe1边缘检测。责任gt