基于FPGA的快速DCT医学图像的压缩算法_临床医学论文.docVIP

基于FPGA的快速DCT医学图像的压缩算法_临床医学论文 基于FPGA的快速DCT医学图像的压缩算法_临床医学论文 【摘要】 介绍了一种适用于医学图像压缩的二维DCT快速算法的FPGA实现结构，采用行列分解法来实现该算法，首先把8×8的二维DCT变换分解为两个一维DCT变换，通过对变换的系数矩阵进行化简，使加法器和乘法器数量减少到最少，并采用了FPGA特有的并行流水线技术，明显节省了计算时间，提高了图像处理速度。 【关键词】 FPGA;二维DCT;医学图像压缩;行列分解;并行流水线 　　Abstract：A FPGA design for fast Discrete Cosine transform(DCT) implementation architecture used in medical image compression is presented. The architecture is achieved using row-column decomposition. First, the 8×8 two-dimensional DCT is decomposed into two one-dimensional DCT, simplifying the coefficient matrix, we can reduce the number of adder and multiplier to minimum, and using the parallel pipeline technology of FPGA, the computation time is greatly saved and image processing speed is improved. 　　Key words：FPGA;2D DCT;Medical image compression;Row-column decomposition;Parallel pipeline 　　1 引 言 　　DCT变换是视频压缩编解码器中很重要的一部分，被广泛应用于各种视频格式的编码算法中，例如：JPEG，MPEG1，MPEG2，H.264等[1]。DCT变换虽然不能直接降低数据量,但是它可以利用图像的统计特性使得各种降低数据量的方法更为有效地工作。它能把图像的能量集中到少数的几个数据上，在很大程度上消除数据间的冗余性和相关性。 　　一般来讲，医学图像都比较大，例如腹部横断面的CT图像(512×512的bmp图像)，对于此类图像，如果要在整体上进行DCT变换，将耗费大量的时间，因此我们需要将整幅图像切割成若干的子块，对子块进行DCT变换，在本研究中采用的是8×8pixels的字块。 　　二维DCT变换需要进行大量的运算，大量的乘法及加法运算严重影响了变换速度，为减少运算次数，缩短运算时间，人们作了不懈的努力,并提出了多种快速算法，但处理过程中数学运算仍然相当复杂。我们在前人工作的基础上，对用行列分解法实现DCT变换的方法进行了优化，并充分利用了FPGA器件的嵌入式乘法器及存储器资源，明显缩短了变换时间。 　　2 DCT算法原理 　　2.1 一维DCT算法原理 　　设{x(n)}表示N个有限的一维实数信号序列集合，n=0,1,2,…,N-1，则一维DCT定义为： 　　y(k)=2Nc(k)∑N-1n=0x(n)cos((2n+1)kπ2N), 　　其中c(k)=12 k=0 　　1 k=1,2,…,N-1 　　2.2 二维DCT算法 　　对于一个N×N的图像，x(i,j)表示图像样值，其二维DCT变换公式如下： 　　y(u,v)=2Ne(u)e(v)∑N-1i=0∑N-1j=0x(i,j) 　　cos(2i+1)uπ2Ncos(2j+1)vπ2N， 　　其中 　　e(u),e(v)=12 u,v=0 　　1 u,v≠0(u,v=0,1,2,…,N 　　生 物 医 学 工 程 研 究 第28卷第3期 王 宁，等：基于FPGA的快速DCT医学图像的压缩算法 二维DCT变换具有很高的复杂度，从上式我们可以看出，二维DCT可以分解为两个一维的DCT变换： 　　y(u,v)=2Ne(v)∑N-1j=0 　　2Ne(u)∑N-1i=0x(i,j)cos(2i+1)uπ2Ncos(2j+1)vπ2N， 　　首先对二维DCT进行一维列变换： 　　x(u,j)=2Ne(u)∑N-1i=0x(i,j)cos(2i+1)uπ2N 　　转置之后再次进行一维列变换： 　　y(u,v)=2Ne(v)∑N-1j=0x(u,j)cos(2i+1)vπ2N 　　至此，我们就完成了一次完整的二维DCT变换。 3 DCT变换在FPGA上的实现 　　3.1 一维DCT的实现 　　在图像压缩