最近邻:选择离它所映射到的位置最近的输入像素的灰度值为插值结果。
最临近插值
3X3 的256级灰度图,也就是高为3个象素,宽也是3个象素的图像,每个象素的取值可以是 0-255,代表该像素的亮度,255代表最亮,也就是白色,0代表最暗,即黑色 。假如图像的象素矩阵如下图所示(这个原始图把它叫做源图,Source):
234 38 22
67 44 12
89 65 63
如果想把这副图放大为 4X4大小的图像,那么该怎么做呢?那么第一步肯定想到的是先把4X4的矩阵先画出来,然后要往这个空的矩阵里面填值了,要填的值从哪里来来呢?是从源图中来,好,先填写目标图最左上角的象素,坐标为(0,0),那么该坐标对应源图中的坐标可以由如下公式得出:
srcX=dstX* (srcWidth/dstWidth) , srcY = dstY * (srcHeight/dstHeight)
套用公式,就可以找到对应的原图的坐标了(0*(3/4),0*(3/4))=>(0*0.75,0*0.75)=>(0,0),找到了源图的对应坐标,就可以把源图中坐标为(0,0)处的234象素值填进去目标图的(0,0)这个位置了。
寻找目标图中坐标为(1,0)的象素对应源图中的坐标,套用公式:(1*0.75,0*0.75)=>(0.75,0)
那么按照四舍五入的方法就得到坐标(1,0),完整的 运算过程就是这样的:(1*0.75,0*0.75)=>(0.75,0)=>(1,0)
那么就可以再填一个象素到目标矩阵中了,同样是把源图中坐标为(1,0)处的像素值38填入目标图中的坐标。
234 38 22 22
67 44 12 12
89 65 63 63
89 65 63 63
这种放大图像的方法叫做最临近插值算法,这是一种最基本、最简单的图像缩放算法,效果也是最不好的,放大后的图像有很严重的马赛克,缩小后的图像有很严重的 失真;效果不好的根源就是其简单的最临近插值方法引入了严重的图像失真,比如,当由目标图的坐标反推得到的源图的的坐标是一个浮点数的时候,采用了四舍五 入的方法,直接采用了和这个浮点数最接近的象素的值,这种方法是很不科学的,当推得坐标值为 0.75的时候,不应该就简单的取为1,既然是0.75,比1要小0.25 ,比0要大0.75 ,那么目标象素值其实应该根据这个源图中虚拟的点四周的四个真实的点来按照一定的规律计算出来的,这样才能达到更好的缩放效果。双线型内插值算法就是一种 比较好的图像缩放算法,它充分的利用了源图中虚拟点四周的四个真实存在的像素值来共同决定目标图中的一个像素值,因此缩放效果比简单的最邻近插值要好很多。
scipy.ndimage.zoom上采样与下采样
双线性插值将是order = 1,
最临近插值的是order = 0,
立方体是默认值(顺序= 3)。
Examples:
import numpy as np
import scipy.ndimage
x = np.arange(64).reshape(8,8)
print 'Original array:'
print x
print 'Resampled by a factor of 2 with nearest interpolation:'
print scipy.ndimage.zoom(x, 2, order=0)
print 'Resampled by a factor of 2 with bilinear interpolation:'
print scipy.ndimage.zoom(x, 2, order=1)
print 'Resampled by a factor of 2 with cubic interpolation:'
print scipy.ndimage.zoom(x, 2, order=3)
print 'Downsampled by a factor of 0.5 with default interpolation:'
print(scipy.ndimage.zoom(x, 0.5))
Results:
Original array:
array([[ 0, 1, 2, 3, 4, 5, 6, 7],
[ 8, 9, 10, 11, 12, 13, 14, 15],
[16, 17, 18, 19, 20, 21, 22, 23],
[24, 25, 26, 27, 28, 29, 30, 31],
[32, 33, 34, 35, 36, 37, 38, 39],
[40, 41, 42, 43, 44, 45, 46, 47],
[48, 49, 50, 51, 52, 53, 54, 55],
[56, 57, 58, 59, 60, 61, 62, 63]])
Resampled by a factor of 2 with nearest interpolation:
[[ 0 0 1 1 2 2 3 3 4 4 5 5 6 6 7 7]
[ 0 0 1 1 2 2 3 3 4 4 5 5 6 6 7 7]
[ 8 8 9 9 10 10 11 11 12 12 13 13 14 14 15 15]
[ 8 8 9 9 10 10 11 11 12 12 13 13 14 14 15 15]
[16 16 17 17 18 18 19 19 20 20 21 21 22 22 23 23]
[16 16 17 17 18 18 19 19 20 20 21 21 22 22 23 23]
[24 24 25 25 26 26 27 27 28 28 29 29 30 30 31 31]
[24 24 25 25 26 26 27 27 28 28 29 29 30 30 31 31]
[32 32 33 33 34 34 35 35 36 36 37 37 38 38 39 39]
[32 32 33 33 34 34 35 35 36 36 37 37 38 38 39 39]
[40 40 41 41 42 42 43 43 44 44 45 45 46 46 47 47]
[40 40 41 41 42 42 43 43 44 44 45 45 46 46 47 47]
[48 48 49 49 50 50 51 51 52 52 53 53 54 54 55 55]
[48 48 49 49 50 50 51 51 52 52 53 53 54 54 55 55]
[56 56 57 57 58 58 59 59 60 60 61 61 62 62 63 63]
[56 56 57 57 58 58 59 59 60 60 61 61 62 62 63 63]]
Resampled by a factor of 2 with bilinear interpolation:
[[ 0 0 1 1 2 2 3 3 4 4 5 5 6 6 7 7]
[ 4 4 5 5 6 6 7 7 7 8 8 9 9 10 10 11]
[ 7 8 8 9 9 10 10 11 11 12 12 13 13 14 14 14]
[11 12 12 13 13 14 14 14 15 15 16 16 17 17 18 18]
[15 15 16 16 17 17 18 18 19 19 20 20 21 21 21 22]
[19 19 20 20 21 21 21 22 22 23 23 24 24 25 25 26]
[22 23 23 24 24 25 25 26 26 27 27 28 28 28 29 29]
[26 27 27 28 28 28 29 29 30 30 31 31 32 32 33 33]
[30 30 31 31 32 32 33 33 34 34 35 35 35 36 36 37]
[34 34 35 35 35 36 36 37 37 38 38 39 39 40 40 41]
[37 38 38 39 39 40 40 41 41 42 42 42 43 43 44 44]
[41 42 42 42 43 43 44 44 45 45 46 46 47 47 48 48]
[45 45 46 46 47 47 48 48 49 49 49 50 50 51 51 52]
[49 49 49 50 50 51 51 52 52 53 53 54 54 55 55 56]
[52 53 53 54 54 55 55 56 56 56 57 57 58 58 59 59]
[56 56 57 57 58 58 59 59 60 60 61 61 62 62 63 63]]
Resampled by a factor of 2 with cubic interpolation:
[[ 0 0 1 1 2 2 3 3 4 4 5 5 6 6 7 7]
[ 2 3 3 4 4 5 5 6 6 7 7 8 8 9 9 9]
[ 7 8 8 9 9 10 10 11 11 12 12 12 13 13 14 14]
[12 12 12 13 13 14 14 15 15 16 16 17 17 18 18 19]
[15 15 16 16 17 17 18 18 19 19 20 20 21 21 22 22]
[19 19 19 20 20 21 21 22 22 23 23 24 24 25 25 26]
[22 23 23 24 24 25 25 26 26 27 27 27 28 28 29 29]
[26 26 27 28 28 28 29 29 30 30 31 31 32 32 33 33]
[30 30 31 31 32 32 33 33 34 34 35 35 35 36 37 37]
[34 34 35 35 36 36 36 37 37 38 38 39 39 40 40 41]
[37 38 38 39 39 40 40 41 41 42 42 43 43 44 44 44]
[41 41 42 42 43 43 44 44 45 45 46 46 47 47 48 48]
[44 45 45 46 46 47 47 48 48 49 49 50 50 51 51 51]
[49 49 50 50 51 51 51 52 52 53 53 54 54 55 55 56]
[54 54 54 55 55 56 56 57 57 58 58 59 59 60 60 61]
[56 56 57 57 58 58 59 59 60 60 61 61 62 62 63 63]]
Downsampled by a factor of 0.5 with default interpolation:
[[ 0 2 5 7]
[19 21 23 26]
[37 40 42 44]
[56 58 61 63]]
本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)