Fork me on GitHub

图像处理-深入探索插值操作

  • 最近邻:选择离它所映射到的位置最近的输入像素的灰度值为插值结果。

 

      图像的缩放很好理解,就是图像的放大和缩小。传统的绘画工具中,有一种叫做“放大尺”的绘画工具,画家常用它来放大图画。当然,在计算机上,我们不再需要用 放大尺去放大或缩小图像了,把这个工作交给程序来完成就可以了。下面就来讲讲计算机怎么来放大缩小图象;在本文中,我们所说的图像都是指点阵图,也就是用 一个像素矩阵来描述图像的方法,对于另一种图像:用函数来描述图像的矢量图,不在本文讨论之列。越是简单的模型越适合用来举例子,我们就举个简单 的图像:3X3 的256级灰度图,也就是高为3个象素,宽也是3个象素的图像,每个象素的取值可以是 0-255,代表该像素的亮度,255代表最亮,也就是白色,0代表最暗,即黑色 。假如图像的象素矩阵如下图所示(这个原始图把它叫做源图,Source):

234 38 22


67 44 12

 
89 65 63

      如果想把这副图放大为 4X4大小的图像,那么该怎么做呢?那么第一步肯定想到的是先把4X4的矩阵先画出来再说,好了矩阵画出来了,如下所示,当然,矩阵的每个像素都是未知数,等待着我们去填充(这个将要被填充的图的叫做目标图,Destination):

? ? ? ?


? ? ? ?
 

? ? ? ?


? ? ? ?

然后要往这个空的矩阵里面填值了,要填的值从哪里来来呢?是从源图中来,好,先填写目标图最左上角的象素,坐标为(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 ,那么目标象素值其实应该根据这个源图中虚拟的点四周的四个真实的点来按照一定的规律计算出来的,这样才能达到更好的缩放效果。双线型内插值算法就是一种 比较好的图像缩放算法,它充分的利用了源图中虚拟点四周的四个真实存在的像素值来共同决定目标图中的一个像素值,因此缩放效果比简单的最邻近插值要好很多。

import cv2
import numpy as np
def function(img):
    height,width,channels =img.shape
    emptyImage=np.zeros((512,512,channels),np.uint8)
    sh=512/height
    sw=512/width
    for i in range(512):
        for j in range(512):
            x=int(i/sh)
            y=int(j/sw)
            emptyImage[i,j]=img[x,y]
    return emptyImage

img=cv2.imread("./a.jpg")
zoom=function(img)
cv2.imshow("nearest neighbor",zoom)
cv2.imshow("image",img)
cv2.waitKey(0)
  • 双线性插值

双线性内插值算法描述如下: 对于一个目的像素,设置坐标通过反向变换得到的浮点坐标为(i+u,j+v) (其中i、j均为浮点坐标的整数部分,u、v为浮点坐标的小数部分,是取值[0,1)区间的浮点数),则这个像素得值 f(i+u,j+v) 可由原图像中坐标为 (i,j)、(i+1,j)、(i,j+1)、(i+1,j+1)所对应的周围四个像素的值决定,即:

  • f(i+u,j+v) = (1-u)(1-v)f(i,j) + (1-u)vf(i,j+1) + u(1-v)f(i+1,j) + uvf(i+1,j+1)    公式1    其中f(i,j)表示源图像(i,j)处的的像素值,以此类推。

      比 如象刚才的例子,现在假如目标图的象素坐标为(1,1),那么反推得到的对应于源图的坐标是(0.75 , 0.75), 这其实只是一个概念上的虚拟象素,实际在源图中并不存在这样一个象素,那么目标图的象素(1,1)的取值不能够由这个虚拟象素来决定,而只能由源图的这四 个象素共同决定:(0,0)(0,1)(1,0)(1,1),而由于(0.75,0.75)离(1,1)要更近一些,那么(1,1)所起的决定作用更大一 些,这从公式1中的系数uv=0.75×0.75就可以体现出来,而(0.75,0.75)离(0,0)最远,所以(0,0)所起的决定作用就要小一些, 公式中系数为(1-u)(1-v)=0.25×0.25也体现出了这一特点;

import cv2
import numpy as np
import math
def function(img,m,n):
    height,width,channels =img.shape
    emptyImage=np.zeros((m,n,channels),np.uint8)
    value=[0,0,0]
    sh=m/height
    sw=n/width
    for i in range(m):
        for j in range(n):
            x = i/sh
            y = j/sw
            p=(i+0.0)/sh-x
            q=(j+0.0)/sw-y
            x=int(x)-1
            y=int(y)-1
            for k in range(3):
                if x+1<m and y+1<n:
                    value[k]=int(img[x,y][k]*(1-p)*(1-q)+img[x,y+1][k]*q*(1-p)+img[x+1,y][k]*(1-q)*p+img[x+1,y+1][k]*p*q)
            emptyImage[i, j] = (value[0], value[1], value[2])
    return emptyImage

img=cv2.imread("./a.jpg")
zoom=function(img,512,512)
cv2.imshow("Bilinear Interpolation",zoom)
cv2.imshow("image",img)
cv2.waitKey(0)

C实现线性和双线性插值

int WarpAffineMoneyInt(IplImage_ *_src_ir_img, u8 *p_src_grey, IplImage_ *_new_grey_img, ZoneInfo_ *zone_info, int *Mat)
{
    int i, j;
    int i1, j1;
    int i0, j0;
    int decimal_i, decimal_j;
    int f12, f34;
    int coordinate1;

    int width = _src_ir_img->width;   // 源图像的宽度和高度
    int height = _src_ir_img->height;
    int widthstep = _src_ir_img->width_step;
    int new_width = _new_grey_img->width;       // 旋转后图像的宽度和高度
    int new_height = _new_grey_img->height;

    int y_start = zone_info->ystart;
    int x_start = zone_info->xstart;

    unsigned char *p_new_grey = (unsigned char*)_new_grey_img->image_data;

    zone_info->type = TYPE_ROTATE_LINE;

    int sjtmp = 0;
    int sitmp = 0;
    int srcpos = 0;
    if (zone_info->type == TYPE_ROTATE_LINE)    //线性插值
    {
        int flag = 1;

        i = y_start;
        j = x_start;
        j1 = (j * Mat[0] + i * Mat[1] + Mat[2] + 32768) >> 16;
        i1 = (j * Mat[3] + i * Mat[4] + Mat[5] + 32768) >> 16;

        if (!((j1 >= 0) && (j1 < width - 1) && (i1 >= 0) && (i1 < height - 1)))   // 判断是否在源图范围内
        {
            flag = 0;
        }

        i = y_start + new_height - 1;
        j = x_start;
        j1 = (j * Mat[0] + i * Mat[1] + Mat[2] + 32768) >> 16;
        i1 = (j * Mat[3] + i * Mat[4] + Mat[5] + 32768) >> 16;

        if (!((j1 >= 0) && (j1 < width - 1) && (i1 >= 0) && (i1 < height - 1)))   // 判断是否在源图范围内
        {
            flag = 0;
        }

        i = y_start;
        j = x_start + new_width - 1;
        j1 = (j * Mat[0] + i * Mat[1] + Mat[2] + 32768) >> 16;
        i1 = (j * Mat[3] + i * Mat[4] + Mat[5] + 32768) >> 16;

        if (!((j1 >= 0) && (j1 < width - 1) && (i1 >= 0) && (i1 < height - 1)))   // 判断是否在源图范围内
        {
            flag = 0;
        }

        i = y_start + new_height - 1;
        j = x_start + new_width - 1;
        j1 = (j * Mat[0] + i * Mat[1] + Mat[2] + 32768) >> 16;
        i1 = (j * Mat[3] + i * Mat[4] + Mat[5] + 32768) >> 16;

        if (!((j1 >= 0) && (j1 < width - 1) && (i1 >= 0) && (i1 < height - 1)))   // 判断是否在源图范围内
        {
            flag = 0;
        }

        if (flag == 1)
        {

            for (i = y_start; i < y_start + new_height; i++)
            {
                sjtmp = i * Mat[1] + Mat[2] + 32768 + (x_start - 1)*Mat[0];
                sitmp = i * Mat[4] + Mat[5] + 32768 + (x_start - 1)*Mat[3];

                for (j = x_start; j < x_start + new_width; j++)     // 针对图像每列进行操作
                {
                    sjtmp += Mat[0];
                    sitmp += Mat[3];
                    j1 = sjtmp >> 16;
                    i1 = sitmp >> 16;

                    srcpos = i1*widthstep;
                    *p_new_grey = p_src_grey[srcpos + j1];

                    p_new_grey++;
                }
                p_new_grey += _new_grey_img->width_step - new_width;
            }
        }
        else
        {

            for (i = y_start; i < y_start + new_height; i++)
            {
                sjtmp = i * Mat[1] + Mat[2] + 32768 + (x_start - 1)*Mat[0];
                sitmp = i * Mat[4] + Mat[5] + 32768 + (x_start - 1)*Mat[3];

                for (j = x_start; j < x_start + new_width; j++)     // 针对图像每列进行操作
                {
                    sjtmp += Mat[0];
                    sitmp += Mat[3];
                    j1 = sjtmp >> 16;
                    i1 = sitmp >> 16;

                    if ((j1 >= 0) && (j1 < width - 1) && (i1 >= 0) && (i1 < height - 1))   // 判断是否在源图范围内
                    {
                        srcpos = i1*widthstep;
                        *p_new_grey = p_src_grey[srcpos + j1];
                    }
                    else
                    {
                        *p_new_grey = 0;
                    }
                    p_new_grey++;
                }
                p_new_grey += _new_grey_img->width_step - new_width;
            }
        }

    }
    else     //双线性插值
    {
        for (i = y_start; i < y_start + new_height; i++)
        {
            for (j = x_start; j < x_start + new_width; j++)    // 针对图像每列进行操作
            {
                j0 = j * Mat[0] + i * Mat[1] + Mat[2];
                i0 = j * Mat[3] + i * Mat[4] + Mat[5];

                decimal_i = i0 & 0xffff;
                decimal_j = j0 & 0xffff;
                i1 = i0 >> 16;
                j1 = j0 >> 16;
                coordinate1 = i1*widthstep + j1;

                if ((j1 >= 0) && (j1 < width - 1) && (i1 >= 0) && (i1 < height - 1))    // 判断是否在源图范围内
                {
                    f12 = (p_src_grey[coordinate1] << 16) + decimal_j*(p_src_grey[coordinate1 + 1] - p_src_grey[coordinate1]);
                    f34 = (p_src_grey[coordinate1 + widthstep] << 16) + decimal_j*(p_src_grey[coordinate1 + widthstep + 1] - p_src_grey[coordinate1 + widthstep]);
                    *p_new_grey = abs(f12 + decimal_i* ((f34 - f12) >> 16)) >> 16;
                }
                else
                {
                    *p_new_grey = 0;
                }
                p_new_grey++;
            }
            p_new_grey += _new_grey_img->width_step - new_width;
        }
    }
    return 0;
}
  •  双三次插值

import cv2
import numpy as np
import math


def S(x):
    x = np.abs(x)
    if 0 <= x < 1:
        return 1 - 2 * x * x + x * x * x
    if 1 <= x < 2:
        return 4 - 8 * x + 5 * x * x - x * x * x
    else:
        return 0
def function(img,m,n):
    height,width,channels =img.shape
    emptyImage=np.zeros((m,n,channels),np.uint8)
    sh=m/height
    sw=n/width
    for i in range(m):
        for j in range(n):
            x = i/sh
            y = j/sw
            p=(i+0.0)/sh-x
            q=(j+0.0)/sw-y
            x=int(x)-2
            y=int(y)-2
            A = np.array([
                [S(1 + p), S(p), S(1 - p), S(2 - p)]
            ])
            if x>=m-3:
                m-1
            if y>=n-3:
                n-1
            if x>=1 and x<=(m-3) and y>=1 and y<=(n-3):
                B = np.array([
                    [img[x-1, y-1], img[x-1, y],
                     img[x-1, y+1],
                     img[x-1, y+1]],
                    [img[x, y-1], img[x, y],
                     img[x, y+1], img[x, y+2]],
                    [img[x+1, y-1], img[x+1, y],
                     img[x+1, y+1], img[x+1, y+2]],
                    [img[x+2, y-1], img[x+2, y],
                     img[x+2, y+1], img[x+2, y+1]],

                    ])
                C = np.array([
                    [S(1 + q)],
                    [S(q)],
                    [S(1 - q)],
                    [S(2 - q)]
                ])
                blue = np.dot(np.dot(A, B[:, :, 0]), C)[0, 0]
                green = np.dot(np.dot(A, B[:, :, 1]), C)[0, 0]
                red = np.dot(np.dot(A, B[:, :, 2]), C)[0, 0]

                # ajust the value to be in [0,255]
                def adjust(value):
                    if value > 255:
                        value = 255
                    elif value < 0:
                        value = 0
                    return value

                blue = adjust(blue)
                green = adjust(green)
                red = adjust(red)
                emptyImage[i, j] = np.array([blue, green, red], dtype=np.uint8)


    return emptyImage

img=cv2.imread("./a.jpg")
zoom=function(img,512,512)
cv2.imshow("cubic",zoom)
cv2.imshow("image",img)
cv2.waitKey(0)

 

posted @ 2018-05-16 19:32  ranjiewen  阅读(576)  评论(0编辑  收藏  举报