【博文精选】基于FPGA的图像FFT滤波处理

1关于傅里叶变换

关于傅里叶变换,这么一个神奇的变换,其基本原理和应用在教科书、网络上漫天飞舞,这里就不赘述了,以免有凑字数的嫌疑。前面的例子我们已经使用Matlab和Vivado的FFT IP核进行了初步的验证,掌握的FFT/IFFT IP核的脾气,那么接下来我们要玩点真的了,基于我们STAR/SF-AT7板采集到的MT9V034图像,我们要进行每个行的FFT和IFFT变换,当然,生成的FFT结果我们可以进行必要的滤波,然后再进行IFFT查看滤波效果。

2基于Matlab的FFT滤波

使用at7_img_ex06\matlab文件夹下的Matlab源码image_1D_fft_ifft.m或L1024_of_image_1D_fft_ifft.m(将640个点扩展为1024个点进行FFT变换,扩展的点以0填充,模拟FPGA的FFT IP核实际工作状况),对测试图像test进行FFT变换,进行必要的滤波,然后IFFT逆变换。
测试图像为彩色图像,原始图像如下。
首先进行彩色转灰度的变换,灰度图像如下。
提取出其中1行进行FFT变换后的图像频谱如下。很明显,大部分高频分量集中在前面几个点,而后面的点几乎频率都很小。
放大频谱图,看到细节如下。这里绘制了一条取值为300的直线,有将近50%的频谱集中在这条线以下。若是做图像压缩,其实我们可以把这些低频分量忽略了,那么数据量可能会大大降低,当然了,副作用是图像可能会有一定程度的失真,有失必有得嘛。滤除这些低频分量,也会使图像更锐一些。话说做FFT变换的目的可远不止这些,在一些特殊的应用场景中,我们总是希望从原始图像中提取出一些和应用直接相关的特征信息,那么做了FFT后的图像常常非常有益于这些操作。为了演示,这里我们的代码里面就将这些低于300的点都滤除,即取0。
从频谱图上看,如图所示,右侧的滤波后明显图像偏黑(很多值取0了)了。

我们重新把原图放到这里,和FFT滤波并IFFT以后的图像做比对,图像整体仍然保持不变,但是查看细节,可以发现处理后的图像明显锐了一些。

Matlab源码如下:

clc;clear `all;close all;

IMAGE_WIDTH = 640;

IMAGE_HIGHT = 480;

%load origin image

%I = imread('Lena_gray_niose.bmp');

I = imread('test.bmp');

I = rgb2gray(I);

%fclose(fid1);

%% output image data in hex file

raw_image = reshape(I, IMAGE_HIGHT, IMAGE_WIDTH);

raw_image = raw_image';

fid2 = fopen('image_in_hex.txt', 'wt');

fprintf(fid2, '%04x\n', raw_image);

fid2 = fclose(fid2);

%show origin image

figure,imshow(I);

title('Original image');

%1D fft base on every image line

II = zeros(IMAGE_HIGHT,1024);

J = zeros(IMAGE_HIGHT,1024);

for i = 1:IMAGE_HIGHT

for j = 1:IMAGE_WIDTH

II(i,j) = I(i,j);

end

J(i,:) = fft(II(i,:));%fft(I(i,:));

end

%show 1 linefft result

t1 = (0:IMAGE_WIDTH);                % Time vector

line = ones(IMAGE_WIDTH) * 200;

figure;

plot(t1(1:IMAGE_WIDTH),abs(J(50,1:IMAGE_WIDTH)),t1(1:IMAGE_WIDTH),line(1:IMAGE_WIDTH))

title(['1 line image in the Frequency Domain'])

%show fft of origin image

figure,imshow(log(abs(J)),[]);

title('1D fft image base on every image line');

%colormap(jet(64)),colorbar;

%fftfiter

J(abs(J) < 300) = 0;

%J(abs(J) > 1000) = 1000;

%show fft of fft filter image

figure,imshow(log(abs(J)),[]);

title('1D fft image after filter');

%1D ifft base on every image line

K = zeros(IMAGE_HIGHT,1024);

for i = 1:IMAGE_HIGHT

K(i,:) = real(ifft(J(i,:)));

end

KK = zeros(IMAGE_HIGHT,IMAGE_WIDTH);

for i = 1:IMAGE_HIGHT

for j = 1:IMAGE_WIDTH

KK(i,j) = K(i,j);

end

end

%show ifft image

figure,imshow(KK,[])

title('1D ifft image');

3FPGA仿真

在Sources面板中,展开Simulation Sources à sim_1,将sim_fft.v文件设置为top module。同样是对前面的测试图像,经过FFT和IFFT变换后存储在image_view0.txt文本中(仿真测试结果位于at7_img_ex06\at7.sim\sim_1\behav文件夹下)。为了确认FFT和IFFT IP核运算的精度和效果,这里没有做任何的滤波处理。
使用draw_image_from_FPGA_result.m脚本(at7_img_ex06\matlab文件夹下)导入image_view0.txt文本的图像,和原始图像比对如下所示。看到图像几乎没有任何失真。

4 基于FPGA的图像平滑处理

工程文件夹at7_img_ex06\zstar.srcs\sources_1\new下的image_fft_filter.v模块以及3个子模块image_fft_controller.v、image_filter.v和image_ifft_controller.v实现了图像的FFT变换、滤波和IFFT变换处理。FPGA设计的功能框图如下。
image_fft_controller.v模块例化FFT IP核,将采集的图像留以行为单位输入到FFT IP核,输出FFT频域数据。
image_filter.v模块对FFT频域数据计算绝对值并进行必要的滤波处理,假设FFT结果的实部值为a,虚部值为b,那么其绝对值abs =sqrt(a^2+b^2)。如下代码,注释部分可以滤除低频分量,当前例程中为了验证FFT和IFFT变换后精度没有损失,未作滤波。

always @(posedgeclk or negedgerst_n)

if(!rst_n) begin

o_image_filter_data_image<= 20'd0;

o_image_filter_data_real<= 20'd0;

end

/*else if(sqrt_fft[19:0] < 20'd300) begin       //此处可以做必要的高频或低频滤波处理

o_image_filter_data_image<= 20'd0;

o_image_filter_data_real<= 20'd0;

end*/

else begin

o_image_filter_data_image<= r_image_fft_data_image[TOTAL_LATENCY-1];

o_image_filter_data_real<= r_image_fft_data_real[TOTAL_LATENCY-1];

end

image_ifft_controller.v模块将滤波处理后的FFT结果进行IFFT变换,图像转回时域值,供后续模块缓存DDR3并显示。

AT7_Xilinx开发板(USB3.0+LVDS)资料共享

腾讯链接:https://share.weiyun.com/5GQyKKc

百度网盘链接:https://pan.baidu.com/s/1M7PLzRs-yMJv7NFJE8GVAw   提取码:qr0t

 

(0)

相关推荐

  • 使用 Verilog HDL 在 FPGA 上进行图像处理

    来源:Hack电子该FPGA项目旨在详细展示如何使用Verilog处理图像,从Verilog中读取输入位图图像(.bmp),处理并将处理结果写入Verilog中的输出位图图像.提供了用于读取图像.图像 ...

  • 傅里叶级数

    一.理论基础   任何周期函数都可以看成是不同振幅,不同相位的正弦波的叠加. 1. 三角函数系的正交性:    三角函数系:          三角函数系中任意两个不同的函数之积在[-π,π]上的积分 ...

  • [Halcon] 图像增强相关算子及其分类

    图像增强就是指通过某种图像处理方法对退化的某些图像特征,如边缘.轮廓.对比度等进行处理,以改善图像的视觉效果,提高图像的清晰度,或是突出图像中的某些"有用",压缩其他"无 ...

  • 【学术论文】基于FPGA的图像自适应加权均值滤波设计

    摘要 针对CMOS图像传感器采集图像过程中的噪声预处理问题,提出一种在FPGA中实现的可配置的自适应加权均值滤波模块设计方案.该模块通过检测滤波窗口内不同方向的方差来确定纹理方向,从而自动生成相应的加 ...

  • 系统设计精选 | 基于FPGA的CAN总线控制器的设计(附代码)

    导读 CAN 总线(Controller Area Network)是控制器局域网的简称,是 20 世纪 80 年代初德国 BOSCH 公司为解决现代汽车中众多的控制与测试仪器之间的数据交换而开发的一 ...

  • 【原创博文】基于FPGA的异步FIFO验证

    现在开始对上一篇博文介绍的异步FIFO进行功能验证,上一篇博文地址:http://blog.chinaaet.com/crazybird/p/5100000872 .对异步FIFO验证的平台如图1所示 ...

  • 【原创博文】基于FPGA的异步FIFO设计

    今天要介绍的异步FIFO,可以有不同的读写时钟,即不同的时钟域.由于异步FIFO没有外部地址端口,因此内部采用读写指针并顺序读写,即先写进FIFO的数据先读取(简称先进先出).这里的读写指针是异步的, ...

  • 【精品博文】基于FPGA的VGA测试程序

    赢一个双肩背包 有多难? 戳一下试试看! →_→ 长摁识别 由于自己先买的黑金的开发板,后买的CB的书,最近在看基于SDRAM的VGA显示,一直有点问题,今天上午本来想写个验证程序看一下问题在哪,但是 ...

  • 【精品博文】基于FPGA的原理图PCB设计 | 电源供电方案选择

    由于FPGA的可编程性,同一芯片在不同应用需求下,使用的逻辑资源不一样时,功率可能会有很大差别,因此在不同应用下可根据功率等来选择不同电源方案. 之前看到看过大牛crazybingo写的一篇关于电源方 ...

  • 【精品博文】基于FPGA的原理图PCB设计

    我的开发工具:Cadence 16.6,所用 FPGA为Altera家的. 画FPGA原理图需准备以下资料: (1) 原理图库文件 (2) 芯片手册 (3) Pin out file 原理图库文件可以 ...

  • 【精品博文】基于FPGA的串口通信时序设计

    这篇博文主要总结一下串口通信的收发时序.上周基于IIC,SPI,RS232串口写了几个简单的小实验,感觉对于这种常用协议的写法大体上都差不多,主要是读懂时序图,然后根据时序图,用HDL语言描述出来就可 ...

  • 【精品博文】基于FPGA的计算器设计——(数码管消“0”)

    我在上一周做了一个数字钟的实验,当时也用到了数码管,我的开发板上有6个数码管,正好2个给小时,2个给分钟,2个给秒.每一个数码管都有自己的任务,所以就不用考虑消"0"的问题了.但这 ...