最近鄰插值與雙線性插值基本原理及Python實(shí)現(xiàn)

1 最近鄰插值

1.1 基本原理

假定

  • 源圖像A的尺寸為m \times n
  • 經(jīng)最近鄰插值處理后的圖像B的尺寸為a \times b

由此可以得到縮放倍數(shù)
K = \frac{a}{m} = \frac芽唇{n}

現(xiàn)取B圖像中的一點(diǎn)(x_0, y_0),可以知道對(duì)應(yīng)在A圖像中的理論計(jì)算位置為

(\frac{x_0 \times m}{a}, \frac{y_0 \times n}符匾)

而上式得出來的在A中的理論位置數(shù)值顯然可能為小數(shù)祈坠,這代表該點(diǎn)在A中無實(shí)際對(duì)應(yīng)點(diǎn)害碾,此時(shí)對(duì)其四舍五入,即是把A圖像中距離該理論點(diǎn)最近的一個(gè)的點(diǎn)當(dāng)作它赦拘,即為最近鄰插值慌随。

1.2 Code

import numpy as np
import math

def nearest_neighbor(input_signal, zoom_multiples):
    '''
    最近鄰插值(適用于灰度圖)
    :param input_signal: 輸入圖像
    :param zoom_multiples:  縮放倍數(shù)
    :return: 縮放后的圖像
    '''
    input_signal_cp = np.copy(input_signal)   # 輸入圖像的副本

    input_row, input_col = input_signal_cp.shape # 輸入圖像的尺寸(行、列)

    # 輸出圖像的尺寸
    output_row = int(input_row * zoom_multiples)
    output_col = int(input_col * zoom_multiples)

    output_signal = np.zeros((output_row, output_col)) # 輸出圖片

    for i in range(output_row):
        for j in range(output_col):
            # 輸出圖片中坐標(biāo) (i躺同,j)對(duì)應(yīng)至輸入圖片中的(m阁猜,n)
            m = round(i / output_row * input_row)
            n = round(j / output_col * input_col)
            # 防止四舍五入后越界
            if m >= input_row:
                m = input_row - 1
            if n >= input_col:
                n = input_col - 1
            # 插值
            output_signal[i, j] = input_signal_cp[m, n]

    return output_signal

1.3 對(duì)比

  • 放大倍數(shù) = 10
原圖
使用最近鄰插值放大十倍
  • 放大倍數(shù)=0.5


    原圖
使用最近鄰插值縮小一倍

2 雙線性插值

2.1 基本原理

實(shí)際上,雙線性插值的原理同最近鄰類似蹋艺,首先同樣的步驟根據(jù)放大倍數(shù)與縮放前后的圖像尺寸找到B圖像中(x_0, y_0)對(duì)應(yīng)A圖像中(\frac{x_0 \times m}{a}, \frac{y_0 \times n}剃袍),然后找到它最近的四個(gè)點(diǎn)捎谨,根據(jù)下列表達(dá)式計(jì)算該點(diǎn)的值民效,即完成雙線性插值。

y = (1-u)(1-v)A(i,j)+(1-u)vA(i,j+1)+u(1-v)A(i+1, j)+uvA(i+1,j+1)

其中

  • i=[\frac{x_0 \times m}{a}], j =[\frac{y_0 \times n}涛救], []為向下取整
  • u = \frac{x_0 \times m}{a} - i
  • v = \frac{y_0 \times n}畏邢- j

相對(duì)于最近鄰插值簡(jiǎn)單的四舍五入,雙線性插值的處理更為科學(xué)州叠,優(yōu)化了邊緣保護(hù)棵红。

2.2 Code

import numpy as np
import math

def double_linear(input_signal, zoom_multiples):
    '''
    雙線性插值
    :param input_signal: 輸入圖像
    :param zoom_multiples: 放大倍數(shù)
    :return: 雙線性插值后的圖像
    '''
    input_signal_cp = np.copy(input_signal)   # 輸入圖像的副本

    input_row, input_col = input_signal_cp.shape # 輸入圖像的尺寸(行凶赁、列)

    # 輸出圖像的尺寸
    output_row = int(input_row * zoom_multiples)
    output_col = int(input_col * zoom_multiples)

    output_signal = np.zeros((output_row, output_col)) # 輸出圖片

    for i in range(output_row):
        for j in range(output_col):
            # 輸出圖片中坐標(biāo) (i咧栗,j)對(duì)應(yīng)至輸入圖片中的最近的四個(gè)點(diǎn)點(diǎn)(x1逆甜,y1)(x2, y2),(x3致板, y3)交煞,(x4,y4)的均值
            temp_x = i / output_row * input_row
            temp_y = j / output_col * input_col

            x1 = int(temp_x)
            y1 = int(temp_y)

            x2 = x1
            y2 = y1 + 1

            x3 = x1 + 1
            y3 = y1

            x4 = x1 + 1
            y4 = y1 + 1

            u = temp_x - x1
            v = temp_y - y1

            # 防止越界
            if x4 >= input_row:
                x4 = input_row - 1
                x2 = x4
                x1 = x4 - 1
                x3 = x4 - 1
            if y4 >= input_col:
                y4 = input_col - 1
                y3 = y4
                y1 = y4 - 1
                y2 = y4 - 1

            # 插值
            output_signal[i, j] = (1-u)*(1-v)*int(input_signal_cp[x1, y1]) + (1-u)*v*int(input_signal_cp[x2, y2]) + u*(1-v)*int(input_signal_cp[x3, y3]) + u*v*int(input_signal_cp[x4, y4])
    return output_signal

2.3 對(duì)比

  • 放大倍數(shù) = 0.5


    原圖
使用最近鄰插值縮小一倍
  • 放大倍數(shù) = 10


    原圖
使用最近鄰插值放大十倍

3 GitHub

click me!

最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
  • 序言:七十年代末斟或,一起剝皮案震驚了整個(gè)濱河市素征,隨后出現(xiàn)的幾起案子,更是在濱河造成了極大的恐慌萝挤,老刑警劉巖御毅,帶你破解...
    沈念sama閱讀 217,734評(píng)論 6 505
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件,死亡現(xiàn)場(chǎng)離奇詭異怜珍,居然都是意外死亡端蛆,警方通過查閱死者的電腦和手機(jī),發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 92,931評(píng)論 3 394
  • 文/潘曉璐 我一進(jìn)店門酥泛,熙熙樓的掌柜王于貴愁眉苦臉地迎上來今豆,“玉大人,你說我怎么就攤上這事柔袁〈舳悖” “怎么了?”我有些...
    開封第一講書人閱讀 164,133評(píng)論 0 354
  • 文/不壞的土叔 我叫張陵捶索,是天一觀的道長(zhǎng)插掂。 經(jīng)常有香客問我,道長(zhǎng)情组,這世上最難降的妖魔是什么燥筷? 我笑而不...
    開封第一講書人閱讀 58,532評(píng)論 1 293
  • 正文 為了忘掉前任,我火速辦了婚禮院崇,結(jié)果婚禮上肆氓,老公的妹妹穿的比我還像新娘。我一直安慰自己底瓣,他們只是感情好谢揪,可當(dāng)我...
    茶點(diǎn)故事閱讀 67,585評(píng)論 6 392
  • 文/花漫 我一把揭開白布。 她就那樣靜靜地躺著捐凭,像睡著了一般拨扶。 火紅的嫁衣襯著肌膚如雪。 梳的紋絲不亂的頭發(fā)上茁肠,一...
    開封第一講書人閱讀 51,462評(píng)論 1 302
  • 那天患民,我揣著相機(jī)與錄音,去河邊找鬼垦梆。 笑死匹颤,一個(gè)胖子當(dāng)著我的面吹牛仅孩,可吹牛的內(nèi)容都是我干的。 我是一名探鬼主播印蓖,決...
    沈念sama閱讀 40,262評(píng)論 3 418
  • 文/蒼蘭香墨 我猛地睜開眼辽慕,長(zhǎng)吁一口氣:“原來是場(chǎng)噩夢(mèng)啊……” “哼!你這毒婦竟也來了赦肃?” 一聲冷哼從身側(cè)響起溅蛉,我...
    開封第一講書人閱讀 39,153評(píng)論 0 276
  • 序言:老撾萬榮一對(duì)情侶失蹤,失蹤者是張志新(化名)和其女友劉穎他宛,沒想到半個(gè)月后船侧,有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體,經(jīng)...
    沈念sama閱讀 45,587評(píng)論 1 314
  • 正文 獨(dú)居荒郊野嶺守林人離奇死亡厅各,尸身上長(zhǎng)有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點(diǎn)故事閱讀 37,792評(píng)論 3 336
  • 正文 我和宋清朗相戀三年勺爱,在試婚紗的時(shí)候發(fā)現(xiàn)自己被綠了。 大學(xué)時(shí)的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片讯检。...
    茶點(diǎn)故事閱讀 39,919評(píng)論 1 348
  • 序言:一個(gè)原本活蹦亂跳的男人離奇死亡琐鲁,死狀恐怖,靈堂內(nèi)的尸體忽然破棺而出人灼,到底是詐尸還是另有隱情围段,我是刑警寧澤,帶...
    沈念sama閱讀 35,635評(píng)論 5 345
  • 正文 年R本政府宣布投放,位于F島的核電站奈泪,受9級(jí)特大地震影響,放射性物質(zhì)發(fā)生泄漏灸芳。R本人自食惡果不足惜涝桅,卻給世界環(huán)境...
    茶點(diǎn)故事閱讀 41,237評(píng)論 3 329
  • 文/蒙蒙 一、第九天 我趴在偏房一處隱蔽的房頂上張望烙样。 院中可真熱鬧冯遂,春花似錦、人聲如沸谒获。這莊子的主人今日做“春日...
    開封第一講書人閱讀 31,855評(píng)論 0 22
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽批狱。三九已至裸准,卻和暖如春,著一層夾襖步出監(jiān)牢的瞬間赔硫,已是汗流浹背炒俱。 一陣腳步聲響...
    開封第一講書人閱讀 32,983評(píng)論 1 269
  • 我被黑心中介騙來泰國(guó)打工, 沒想到剛下飛機(jī)就差點(diǎn)兒被人妖公主榨干…… 1. 我叫王不留,地道東北人权悟。 一個(gè)月前我還...
    沈念sama閱讀 48,048評(píng)論 3 370
  • 正文 我出身青樓恼蓬,卻偏偏與公主長(zhǎng)得像,于是被迫代替她去往敵國(guó)和親僵芹。 傳聞我的和親對(duì)象是個(gè)殘疾皇子,可洞房花燭夜當(dāng)晚...
    茶點(diǎn)故事閱讀 44,864評(píng)論 2 354