在圖像通信,、遙感圖像分析、醫(yī)學(xué)成像診斷等應(yīng)用領(lǐng)域,，為了便于顯示,、觀察或進(jìn)行進(jìn)一步的處理，常常需要對(duì)原始的數(shù)字圖像進(jìn)行特征提取(如邊緣檢測(cè),、邊緣銳化),、噪聲平滑濾波、幾何校正等處理,，這類(lèi)圖像處理技術(shù)稱(chēng)為圖像的預(yù)處理,。在實(shí)際應(yīng)用中，空域?yàn)V波算法被廣泛地應(yīng)用于圖像的預(yù)處理技術(shù)中,。
空域?yàn)V波算法是圖像增強(qiáng)技術(shù)的一種,，直接對(duì)圖像的象素進(jìn)行處理，不需要進(jìn)行變換,。常見(jiàn)的濾波算子如銳化算子,、高通算子、平滑算子等，可以完成圖像的邊緣提取,、噪聲去除等處理,。這些濾波算子盡管功能不同，實(shí)現(xiàn)方法卻都是類(lèi)似的,，都是通過(guò)模板卷積的方法來(lái)實(shí)現(xiàn)的,。
VLSI技術(shù)的迅猛發(fā)展為數(shù)字圖像實(shí)時(shí)處理技術(shù)提供了硬件基礎(chǔ)，其中FPGA(現(xiàn)場(chǎng)可編程門(mén)陣列)的特點(diǎn)使其非常適用于數(shù)字圖像處理,。本文研究的就是在FPGA設(shè)計(jì)平臺(tái)上設(shè)計(jì)硬件電路,，實(shí)現(xiàn)數(shù)字圖像的空域?yàn)V波算法。
1 數(shù)字圖像空域?yàn)V波算法
數(shù)字圖像空域?yàn)V波算法的實(shí)現(xiàn)步驟如圖1所示,，左邊的部分是要處理的圖像的某一部分,，中間是對(duì)圖像進(jìn)行處理的3×3模板。
具體的處理步驟是：
將模板在圖像上漫游,，并將模板中心與圖中某個(gè)象素位置重合,；
將模板上的系數(shù)與模板下對(duì)應(yīng)的象素相乘；
將所有的乘積相加,。
把和(模板的輸出響應(yīng))賦給圖像中對(duì)應(yīng)模板中心位置的象素,。圖1中所示是圖像的一部分，S0～S8是象素點(diǎn)的灰度值,，K0～K8是3×3的模板系數(shù),。用這個(gè)3×3模板來(lái)進(jìn)行空域?yàn)V波的過(guò)程是：將模板中心點(diǎn)，即模板系數(shù)為K0的點(diǎn)所在位置,，與圖像中灰度值為S0的點(diǎn)重合,，模板的輸出響應(yīng)R為：
R=K0*S0+K1*S1+…+K8*S8 (1)
這樣增強(qiáng)后的圖像在原來(lái)位置為(x，y)處的象素點(diǎn)的灰度值就由S0變?yōu)镽,。如果對(duì)圖像中的每個(gè)象素點(diǎn)都這樣進(jìn)行模板操作,，就可以得到增強(qiáng)后的圖像在所有位置的新灰度值。如果在設(shè)計(jì)濾波器時(shí)給模板系數(shù)賦予不同的值,，就可以得到不同的高通,、低通效果。
本文采用的圖像是256×256大小的灰度圖像,，濾波模板3×3大小,。如何設(shè)計(jì)硬件電路來(lái)完成上述空域?yàn)V波算法,，分析上述算法實(shí)現(xiàn)過(guò)程,，可以得出結(jié)論，實(shí)現(xiàn)空域?yàn)V波算法可采用3個(gè)三階的FIR濾波器+延時(shí)單元來(lái)描述,。
2 FIR數(shù)字濾波器的FPGA設(shè)計(jì)
在設(shè)計(jì)實(shí)現(xiàn)空域?yàn)V波算法的3個(gè)三階的FIR濾波器+延時(shí)單元的電路時(shí),，要考慮的主要問(wèn)題是：如何縮短硬件電路設(shè)計(jì)時(shí)的關(guān)鍵路徑以及提高系統(tǒng)的數(shù)據(jù)吞吐率。為解決這些實(shí)際的FPGA設(shè)計(jì)中的關(guān)鍵問(wèn)題，在具體設(shè)計(jì)電路時(shí),，主要從以下幾個(gè)方面進(jìn)行考慮：
2.1 FIR數(shù)字濾波器與流水線結(jié)構(gòu)
現(xiàn)代微處理器,、數(shù)字信號(hào)處理器、高速數(shù)字系統(tǒng)設(shè)計(jì)中都廣泛應(yīng)用了流水線(Pipelining)技術(shù),，其核心設(shè)計(jì)思想是把一個(gè)周期內(nèi)執(zhí)行的邏輯操作分成幾步較小的操作,，在多個(gè)高速的時(shí)鐘周期內(nèi)完成。每一次邏輯小操作的結(jié)果都存儲(chǔ)在寄存器中,，被高速時(shí)鐘同步,，在下一流水線單元使用，因而是速度優(yōu)化中最常用的技術(shù)之一,，可以大大地提高數(shù)字系統(tǒng)的總體運(yùn)行速度,。
下面分析一下三階的FIR濾波器的基本結(jié)構(gòu)和采用了流水線技術(shù)后的FIR結(jié)構(gòu)，以及FIR濾波器的數(shù)據(jù)廣播結(jié)構(gòu),。
三階有限沖擊響應(yīng)(FIR)數(shù)字濾波器可以表示如下：
y(n)=ax(n)+bx(n-1)+cx(n-2) (2)
這個(gè)三階的FIR濾波器實(shí)現(xiàn)的結(jié)構(gòu),，如圖2所示。
圖2中,，這種結(jié)構(gòu)的FIR濾波器的關(guān)鍵路徑(處理一個(gè)新樣點(diǎn)的最小時(shí)間)由1個(gè)乘法器和2個(gè)加法器的時(shí)間來(lái)限定,。如果采樣周期小于這個(gè)最小時(shí)間的話(huà)，那么這種結(jié)構(gòu)的FIR濾波器就不能滿(mǎn)足要求,。這時(shí)就要考慮流水技術(shù),。采用流水線技術(shù)可以縮短關(guān)鍵路徑，如圖3所示,。
在流水線結(jié)構(gòu)的FIR濾波器中,，啟動(dòng)當(dāng)前的迭代計(jì)算時(shí)，節(jié)點(diǎn)2的加法器正在完成前次迭代結(jié)果的計(jì)算,。因此,，這時(shí)的關(guān)鍵路徑由1個(gè)乘法器和2個(gè)加法器的時(shí)間縮短為1個(gè)乘法器和1個(gè)加法器的時(shí)間。
采用在結(jié)構(gòu)中適當(dāng)插入流水線鎖存器來(lái)減小關(guān)鍵路徑長(zhǎng)度的流水技術(shù)時(shí),，鎖存器的插入并不是隨意的,。當(dāng)把數(shù)據(jù)流圖切開(kāi)時(shí)，數(shù)據(jù)方向要一致向前,。這樣加入的流水線,，才不會(huì)影響功能。圖3中,，插入流水線鎖存器時(shí),，就是沿著結(jié)構(gòu)中數(shù)據(jù)流的正向，在上下兩個(gè)路徑均加入了鎖存器,，這樣FIR濾波器的邏輯才不會(huì)混亂,。一種結(jié)構(gòu)的速度(時(shí)鐘周期),，通常由任意兩個(gè)鎖存器之間、一個(gè)輸入與一個(gè)鎖存器間,、一個(gè)鎖存器與一個(gè)輸出間或輸入與輸出間路徑中最長(zhǎng)的路徑限定,。流水線鎖存器可以有效地縮短最長(zhǎng)路徑。
除了以上兩種FIR濾波結(jié)構(gòu)外,，還有一種數(shù)據(jù)廣播結(jié)構(gòu)的FIR數(shù)字濾波器,，這種結(jié)構(gòu)通過(guò)轉(zhuǎn)置結(jié)構(gòu)來(lái)縮短關(guān)鍵路徑，同時(shí)不需要引入任何的流水線鎖存器,。具體的改造方法是：改變輸入和輸出,；顛倒信號(hào)流程的方向；由一個(gè)分支來(lái)取代加法器,，反之亦然,。數(shù)據(jù)廣播結(jié)構(gòu)的FIR數(shù)字濾波器如圖4所示。
在這種結(jié)構(gòu)中,，數(shù)據(jù)不是存儲(chǔ)下來(lái),，而是同時(shí)廣播到所有的乘法器。這種結(jié)構(gòu)的關(guān)鍵路徑和圖3中插人流水線鎖存器的FIR濾波器結(jié)構(gòu)的關(guān)鍵路徑一樣,。但是對(duì)輸入不需要額外的移位寄存器,，對(duì)部分積的和也不需要額外的流水線就能達(dá)到很高的通過(guò)率。這就是FIR濾波器數(shù)據(jù)廣播結(jié)構(gòu)的優(yōu)點(diǎn),。
本文在設(shè)計(jì)空域?yàn)V波算法電路時(shí),，就是分別采取了以上介紹的3種不同的FIR數(shù)字濾波器結(jié)構(gòu)形式。
2.2 乘法器模塊的硬件設(shè)計(jì)
從式(1)的模板運(yùn)算表達(dá)式和式(2)的FIR濾波器表達(dá)式可以看出,，完成模板運(yùn)算和實(shí)現(xiàn)FIR數(shù)字濾波還有一個(gè)重要的環(huán)節(jié),，就是乘法運(yùn)算。乘法器模塊是影響空域?yàn)V波算法的運(yùn)算速度的關(guān)鍵模塊之一,。
乘法運(yùn)算基本上可以分為兩步：一是求出所有的基本乘積項(xiàng),，二是將所有的基本乘積項(xiàng)相加。因此,，要設(shè)計(jì)快速的乘法器電路模塊,，就要針對(duì)這兩步進(jìn)行改進(jìn)，一方面要減少部分積的數(shù)目,，同時(shí)另一方面要提高部分積求和陣列的累加速度,。因此，為了加快乘法器模塊的運(yùn)算速度,，在設(shè)計(jì)乘法器電路時(shí),，特別考慮到了采用基4-BOOTH算法來(lái)減少部分和的數(shù)目，同時(shí)采用Wallace Tree減少陣列乘法器中部分積加法陣列的進(jìn)位傳輸延遲,，加快整個(gè)加法陣列的運(yùn)算速度,。
基4-BOOTH算法的基本原理是對(duì)乘數(shù)進(jìn)行編碼,，根據(jù)編碼表來(lái)產(chǎn)生部分積,，一次只考慮3位：本位,、相鄰高位、相鄰低位,。Wallace Tree比較規(guī)則,，易于布局布線，這種方法不是直接將所有的部分積完全的一對(duì)一地相加,，而是采用將各個(gè)部分積中具有相同權(quán)重的數(shù)據(jù)位相加合并,。通常采用全加器(Full Adder)來(lái)完成相同權(quán)重的位相加。采用一位全加器,，那么Wallace樹(shù)的每一層,，就可以將部分積的向量數(shù)目按照3：2的比例縮減。也可以采用2個(gè)全加器,，來(lái)獲得4：2的縮減比例,。本文中采用3：2計(jì)數(shù)器(全加器)來(lái)進(jìn)行部分和縮減，這樣當(dāng)部分積的數(shù)目很多時(shí),，采用Wallace Tree乘法器來(lái)縮減部分和速度很快,。
3 仿真綜合結(jié)果
第2節(jié)介紹了本文中在FPGA平臺(tái)上設(shè)計(jì)實(shí)現(xiàn)數(shù)字圖像空域?yàn)V波算法的高速數(shù)字濾波器FIR的過(guò)程，其中主要考慮的是縮短關(guān)鍵路徑,、提高數(shù)據(jù)吞吐率,。本節(jié)給出仿真和綜合后的結(jié)果比較表。本文中,，測(cè)試圖像選用的是256×256大小,、8 b的灰度圖像，設(shè)計(jì)軟件是XILINX公司的ISE集成綜合開(kāi)發(fā)環(huán)境,，仿真工具是Modelsim SE 5.8b,，綜合工具是ISE自帶的綜合軟件XST，實(shí)現(xiàn)芯片是XILINX公司的XC2V1000,。仿真,、綜合結(jié)果表明，設(shè)計(jì)電路完全符合要求,。
下面給出根據(jù)上述不同的FIR濾波結(jié)構(gòu)設(shè)計(jì)的,，用于實(shí)現(xiàn)空域?yàn)V波算法的3種電路結(jié)構(gòu)的綜合結(jié)果比較表。
從表1可以看出：從資源占用角度看,，結(jié)構(gòu)三的等效門(mén)數(shù)最多,，結(jié)構(gòu)二的最少。從延時(shí)／最大頻率可以看出,，結(jié)構(gòu)一最好,。
3種結(jié)構(gòu)的結(jié)構(gòu)時(shí)延比較數(shù)據(jù),，如表2所示。
4 結(jié) 語(yǔ)
本文論述了數(shù)字圖像空域?yàn)V波算法以及FIR濾波器的基本設(shè)計(jì)方法,，在對(duì)關(guān)鍵路徑分析的基礎(chǔ)上,，引入流水線設(shè)計(jì)提高運(yùn)算速度，提出了濾波器的3種設(shè)計(jì)結(jié)構(gòu),，給出了濾波器的設(shè)計(jì)過(guò)程,，通過(guò)仿真和綜合結(jié)果可以看出，有效地節(jié)省了硬件資源,，大大減小了硬件體積,，增加了系統(tǒng)的可靠性。
