数字图像处理——图像插值

网上有很多介绍插值算法的,但感觉收获都不大



简书的编译器可能有点问题tex公式在预览模式公式可以显示,从主页浏览就显示不出来,都换了图片

介绍三种图像插值算法:最近邻内插,双线性内插,双三次内插(双立方内插)


1.数学原理

1.1线性插值

linear

线性插值是指插值函数为一次多项式的插值方式,在已知原函数上两点的情况下,用线性函数拟合原函数(应该还有很多用途)。假设拟合的线性函数为

将已知两点坐标带入,满足
上式的解为
则线性拟合函数为

对两点的坐标做一个线性变化
得到的表达式为

1.2三次插值(立方插值)

三次插值即用三阶多项式拟合原函数(也应该有其他用途)。假设三次拟合函数为

若我们可以得到原函数的一阶导数信息,那么只需要两点函数以及导数信息就可求解拟合函数。当不能得到一阶导数信息时,需要至少用四点坐标值求解拟合函数。
为了方便计算,在采样时可以选择等间隔采样,进行线性变化后,将横坐标变为
将已知四点坐标带入,得到
不妨将上式表示为
其中
这里还有另一种三次内插的方式,就是利用两点之间的斜率来替代导数信息,根据两点的函数以及一阶导数值(指上一句话计算出的斜率)列一组方程解出拟合的方程。

2.程序设计思路

2.1坐标预处理

2.1.1坐标映射

在matlab中,图像被定义为一个三维向量,若不考虑图像的通道数,可以将图像看作一个二维矩阵处理。matlab图像矩阵中坐标值映射到二维坐标系中,每一个像素块对应的是一个点,但实际的像素块是有一定尺寸的。


mapping

为了保证图像缩放时的相对位置不变,选取图像左上角顶点为原点,每个像素块中心坐标为像素的坐标,那么matlab中像素点坐标映射到坐标系中时,需要减去一个偏移量
在将缩放后的照片矩阵坐标映射回原照片矩阵坐标时,需要满足如下关系
根据映射回原图片矩阵的坐标进行插值操作。

2.1.2局部坐标系

在进行双线性插值和双三次插值时,需要用坐标值拟合函数,为了简化计算,总是选取(\lfloor x \rfloor,\lfloor y \rfloor)作为局部坐标系原点,其中(x,y)为待插值坐标。

2.2边框处理

在进行双线性内插放大和双立方内插缩放是,需要用到的像素块超出了图像的边界,即

当出现这些情况时,补充这些像素的灰度值为图像内最相邻像素块的灰度值。

2.3最近邻内插

进行坐标变换后,选取与内插点(x,y)欧式距离最近的像素值进行插值。在程序中,使用将(x,y)按照四舍五入的舍入方式选取最近邻的像素块。

2.4双线性内插

双线性内插是线性内插的二维实现,在x维度先进行线性插值,再由得到的值对y维度进行插值。在局部坐标系中,选取(x,y)相邻的四个像素进行双线性内插。由在数学原理中的推导可知

将上式合并可得

bilinear

当然,也可以跳过分析的步骤,直接求解双线性拟合方程的系数。设方程为

将四点与相应的坐标对带入求解方程组,所得结果相同。

2.5双立方内插

双三次内插是三次插值的二维实现。选取与(x,y)相邻的16个像素进行双三次内插,局部坐标系中x与y坐标范围均为[-1,0,1,2]。由数学原理中的推到可知

其中
将上述几式合并,可得(有点复杂,不过比矩阵求逆快多了,这个结果的推导我是用matlab的)
与双线性插值相似,双三次插值也可以直接求解拟合方程的系数。设拟合的方程为
将16各坐标对与灰度值带入方程,求解方程组可以得到的系数解为
将系数带入,可得到与之前相同的结果
根据上述公式对图像进行双三次内插。

3.实战结果

lena

对比三种插值方式得到的图像,可以看出双立方内插效果最好,双线性内插次之,最近邻内插最差。lena图像中对比帽子中的纹理可以得出最近邻帽子纹理有许多马塞克和断层,双线性的纹理更平滑,双立方中帽子的纹理比起双线性更锐化。放大图像中比较瞳孔与眼白之间的区域可以看出,最近邻的马赛克现象很严重,双线性也有一些马赛克现象,双立方的边缘很平滑,在不放大的情况下几乎看不出马赛克现象。
cameraman

cameraman缩小图像中对比机架和人耳部分,可以看出双立方内插效果最好,双线性内插次之,最近邻内插最差。但是对于人脸部分和相机机身部分之类细节比较多的地方,由于图像缩小的原因,像素块减少,三种内插方法都比较模糊。放大图像中比较机架与人身体的边缘可以看出,最近邻的锯齿现象很严重,双线性也有一些锯齿现象,双立方的边缘很平滑,在不放大的情况下几乎看不出锯齿现象。
building

building缩小图像中对比屋顶纹理和牌楼图案部分,可以看出双立方内插效果最好并保留了更多的边缘和细节的部分,双线性内插次之,最近邻内插最差。但是整体上来看,双线性内插与双立方内插区别并不大,原图的细节过多,这几种方法都没能很好保留原图的细节。放大图像中可以很明显的看出双立方插值的优势,线条的轮廓边缘更细致。

4.结论

4.1最近邻插值法

最近邻插值法的优点是计算量很小,运算速度较快。但它仅使用离待测采样点最近的像素的灰度值作为该采样点的灰度值,而没考虑其他相邻像素点的影响,因而重新采样后灰度值有明显的不连续性,会产生明显的马赛克和锯齿现象。

4.2双线性插值法

双线性插值法效果要好于最近邻插值,计算量较大。缩放后图像质量高,基本克服了最近邻插值灰度值不连续的特点,因为它考虑了待测采样点周围四个直接邻点对该采样点的相关性影响。但是,此方法未考虑到各邻点间灰度值变化率的影响, 具有低通滤波器的性质, 从而导致缩放后图像的高频分量受到损失, 图像边缘在一定程度上变得较为模糊,丢失了一些细节信息。

4.3双立方插值法

双立方插值计算量最大,运算速度慢。双立方插值用三阶函数逼近,不仅考虑到周围四个直接相邻像素点灰度值的影响,还考虑到它们灰度值变化率的影响,能够产生比双线性插值更为平滑的边缘,计算精度很高,处理后的图像细节损失最少,效果最佳。

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

推荐阅读更多精彩内容

  • 在我看来笑容是我们每个人最美的一种表情,真诚的笑容不仅能够表现出自己愉快,更是一种分享快乐的方式。笑是是我们所有人...
    忆术品阅读 279评论 0 0
  • 2019.2.26星期②阴 读经人员:妈妈、芷晴、芷翔 读经第793天 读经内容:妈妈、芷晴 1、《伤寒论》 2、...
    欢欢成长足迹阅读 218评论 0 1
  • 3. 华灯初上,身影斑驳,他们两个人缓慢走在夜里显得尤其深黑的跑道上,操场上嬉闹声,奔跑声……不绝于耳,而他只想听...
    程景至阅读 251评论 0 4