计算摄影学Lab2:图像滤波和傅里叶变换

任务清单

  • 实现盒装均值滤波
  • 实现高斯滤波
  • 实现中值滤波
  • 实现简单的双边滤波
  • 利用傅里叶变换完成图像的频域变换

空域滤波

空域滤波的基本思想是取像素的邻域,将这个邻域内的全部图像信息进行整合以得到对应像素的估计

均值滤波

均值滤波是一种最简单的滤波方式。它以像素邻域内的平均值代替原来的像素。
因为均值滤波是线性的,所以说我们可以引入一个卷积核,用空域卷积来表示它。卷积核和图像域坐标无关,所以均值滤波也是均匀的。
本次实验中,我们采用计算卷积核,利用卷积核进行卷积的方法实现均值滤波。
在OpenCV中,图像的卷积计算可以通过cv::filter2D来完成,而且边界已经处理所以我们只需要实现自己写出卷积核即可。
源代码放在了 BoxFilter.cpp 当中,可执行程序为 BoxFilter.exe。参数顺序同实验要求一致,实现的结果如下所示,其中,上面为原图像,下面为滤波过后的图像:

均值滤波前
均值滤波后

高斯滤波

通常我们认为图像像素之间的相关性随着距离增加应该不断减弱。然而均值滤波并没有体现这一性质。在对图像进行均值滤波时,如果图像中有一些很显著的亮点,滤波后它的周围会形成光斑。这正是因为均值滤波无视了距离,对很远处的像素依旧采用同样的权重导致的。一些场合,我们为了美感会需要这种效果。另一些场合,这种结果是不利的,特别是在做图像处理时。
高斯滤波是另外一种均匀的线性滤波器。在这部分中,我们依旧使用cv::filter2D函数,自己生成高斯卷积核,实现高斯滤波。
源代码放在了 Gaussian-Filter.cpp 当中, 可执行程序为 GaussianFilter.exe。参数顺序同实验要求一致,实现的结果如下所示,其中,上面为原图像,下面为滤波过后的图像:

高斯滤波前
高斯滤波后

中值滤波

高斯滤波适用于图像带有高斯噪声情况下的去噪。在图像被非高斯噪声污染的情况下,高斯滤波不一定能得到理想的去噪效果。
中值滤波是一种序统计滤波器(Order-Statistic Filter),序统计滤波器是依据邻域的值在统计上的次序关系来进行过滤的。这里,中值滤波器用邻域内像素亮度的中值来取代原本的像素值。所以中值滤波是非线性的。
冲击噪音(随机的01噪音)不是高斯的,利用前面的卷积滤波方法不能得到很好的结果。 中值滤波对于冲击噪音有比较好的抑制效果。直观上看,在邻域内的样本中,冲击噪音平均地分布在最大和最小端,通过取中值可以较好地回避冲击噪音。
在这部分中,我自己实现了中值滤波,由于不能用卷积核卷积的方式来实现,所以我们不能用之前的函数来处理。在这里我采用了一个点一个点的去计算附近的中值。对于边界的点,我选择不处理(用原来的像素点代替)。
源代码放在了 MedianFilter.cpp 当中,可执行程序为 MedianFilter.exe。参数顺序同实验要求一致,实现的结果如下所示,其中,上面为原图像,下面为滤波过后的图像:

中值滤波前
中值滤波后

双边滤波

前面的三种滤波器都会破坏图像的边界,在卷积核很大的时候,均值滤波和高斯滤波都会让边界变得模糊,在邻域很大时,中值滤波会减小边界的曲率。由于物体边界是物体的一个重要特征,很多任务里我们不希望图像边界被破坏。
双边滤波提供了一种降噪同时保持边界的方法。它的思路很简单:如果邻域内像素的亮度差异很大,它在加权平均时的贡献也应当小。我们可以在高斯滤波加权平均的基础上引入一个新的项,反应亮度差带来的加权。
在这部分内容中,我自己实现了双边滤波,实现的方法还是对每一个点进行计算,分别计算时间和空间上的系数,然后求得总和。由于需要大量的计算,这个程序的运行时间会较长。对于边界点,我同样选择不处理(用原来的像素点代替)。注意,在这个实验中,不应对高斯参数归一化,因为在计算结果最后,会对点进行归一化的处理。
源代码放在了 BilateralFilter.cpp 当中,可执行程序为 BilateralFilter.exe。参数顺序同实验要求一致,实现的结果如下所示,其中,上面为原图像,下面为滤波过后的图像:

双边滤波前
双边滤波后

傅里叶变换

本次实验要求对一个图像进行傅里叶变换,源代码放在了 dft.cpp 当中,可执行程序为 dft.exe。参数顺序同实验要求一致,实现的结果如下所示,其中,上面为原图像,下面为傅里叶变换过后的图像:

傅里叶变换前
傅里叶变换后

实验总结

本次实验中,我们采用不同的方法分别实现了各种滤波方式以及傅里叶变换,对于图像滤波的概念有了更加深刻的理解,同时对 OpenCV 处理图像的方式更加熟悉。

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 203,456评论 5 477
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 85,370评论 2 381
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 150,337评论 0 337
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 54,583评论 1 273
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 63,596评论 5 365
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 48,572评论 1 281
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 37,936评论 3 395
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 36,595评论 0 258
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 40,850评论 1 297
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 35,601评论 2 321
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 37,685评论 1 329
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 33,371评论 4 318
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 38,951评论 3 307
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 29,934评论 0 19
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 31,167评论 1 259
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 43,636评论 2 349
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 42,411评论 2 342

推荐阅读更多精彩内容

  • 不同图像灰度不同,边界处一般会有明显的边缘,利用此特征可以分割图像。需要说明的是:边缘和物体间的边界并不等同,边缘...
    大川无敌阅读 13,815评论 0 29
  • http://blog.csdn.net/x454045816/article/details/52153250 ...
    G风阅读 7,014评论 0 1
  • 参考资料: 图像卷积与滤波的一些知识点 图像处理基本概念——卷积,滤波,平滑 1.卷积的基本概念 首先,我们有一个...
    keloli阅读 9,965评论 0 26
  • 最新刚好遇到个需求是要求做高斯模糊的,虽然现有已经有一些框架可以提供调用,但关键还是要理解原理才行,思考的过程才是...
    Hohohong阅读 13,530评论 1 37
  • 经常玩唱吧都知道最近出来了一个直播app火星。主播都是唱吧里有点名气的人。 因为有个一直关注的女神也跑去直播了,本...
    粽子鸡阅读 2,621评论 0 2