Python實現RGB等圖片的圖像插值算法
前言
問題:我們在放大圖片的過程中,放大的目標圖像和原圖圖像之間會存在不同。
放大的基本思想:
第一步:
將目標圖像進行縮小到原圖像的尺寸,虛擬的將二者進行對重疊在一起,這樣目標圖像上面的像素點就和原圖像上面的像素點並不是一一對應的。
第二步:
將目標圖像與原圖像的像素點進行建立一個映射關系,這個像素點的映射關系不是單一的像素映射,因為我們把圖片存放在三維數組裡面,所以我們在得到圖像的像素點時,我們是通過索引數組的高和寬,在坐標系中對應的就是x坐標、y坐標.
第三步:
怎樣建立映射關系呢?
- 最臨近插值法
- 線性插值法
- 雙線性插值法
- 三線性插值法
本文中用到瞭前面三種插值算法,並且進行瞭代碼實現,代碼僅供參考,大佬隨便看一下小菜雞的就行瞭,還希望大佬能指出其中的錯誤或不足之處,跪謝。
RGB彩色圖像和數組理解
對於這個圖片讀取到數組中形成的三維數組,我剛理解瞭很久,在網上找大佬的資料自己看等等,然後慢慢的才形成瞭自己的理解,大佬們都是縱說紛紜,自己走瞭很多彎路,這裡吧我的理解寫下來,第一,方便自以後學習查看;第二,方便和我一樣的初學者理解。
先話不多說,直接上一個圖:
這裡我直接把彩色圖片的三個通道數畫成瞭三位圖片顯示出來方面自己理解。彩色圖片是三通道的,分別為R、G、B,這三個通道的重疊,通過調節每個通道數灰度值的亮度,從而構成瞭五顏六色的彩色世界!!
紅色區域為第0通道(R),以此類推,第1通道(G)、第2通道(B)。讀取回來的數組是一個三維數組,這個三維數組又分為瞭很多二維矩陣,每個二維矩陣有三個列(三個通道)。上圖中有多少個i(高h),就表示三維數組中有多少個二維數組,有多少個j(寬h),就表示二維數組有多少個行,通道數k表示通道數,彩色圖像有三個通道,所以 k 是個固定值 3。
數組中的值怎麼理解欸?
下面插入我這個靈魂畫師的一張圖示來說明吧哈哈哈哈哈哈(不禁默默的流下瞭眼淚)
上圖中左邊給出瞭一個圖片讀入到數組(假設)中的樣子,不同顏色的線條表示組成的不同通道圖片,然後再把這三個通道裡面的圖片按照上上圖片一樣堆疊起來,就構成瞭這個三維空間圖瞭。
RGB圖像的理解就說到這裡吧,希望能幫助到和我一樣的初學者少走彎路。
圖片坐標對其
要問這公式怎麼來的,第一個可以根據放大後像素點的位置成比例算出來,第二個公式暫時還沒想出來咋個算出來的,要是有大佬給我指出來就好瞭。
在代碼中有傳入參數指定使用哪一種對齊方式align = left,該參數默認是居中對其center。
左對齊
src_X = dst_X*(src_Width/dst_Width) src_Y = dst_Y*(src_Height/dst_Height)
這裡的src_X 就是目標圖像上的點映射到原圖像上的x坐標點,同理src_Y 就是映射到原圖像上的y坐標點。
dst_X表示目標圖像的x坐標,dst_Y表示目標圖像的y坐標。
中心對齊
src_X = (dst_X+0.5)*(src_Width/dst_Width) - 0.5 src_Y = (dst_Y+0.5)*(src_Height/dst_Height) - 0.5
這裡的src_X 就是目標圖像上的點映射到原圖像上的x坐標點,同理src_Y 就是映射到原圖像上的y坐 標點。
dst_X表示目標圖像的x坐標,dst_Y表示目標圖像的y坐標。
臨近插值算法
最鄰近插值算法是最簡單的,我們算出來的坐標點:i ,j,使用round函數對其進行四舍五入,就可以得到img[i,j]這個像素點的臨近值瞭(是一個像素值)。
這個圖片中就能看出,A點是目標圖像映射到原圖像上面的位置,整個背景是原圖像,圖中的Q1、Q2、Q3、Q4四個點就是A點的最鄰近的四個像素點。
最鄰近插值算法的代碼是最簡單的,但是放大後的圖片效果是最差的,下面看代碼實現:
def nearest_inter(self): """最鄰近插值法""" new_img = np.zeros((self.goal_h, self.goal_w, self.goal_channel), dtype=np.uint8) for i in range(0, new_img.shape[0]): for j in range(0, new_img.shape[1]): src_i, src_j = self.convert2src_axes(i, j) new_img[i, j] = self.img[round(src_i), round(src_j)] return new_img
線性插值法
線性插值公式說白瞭就是,咱們數學中的直線的參數方程形式。我們知道兩點可以確定一條直線,但是在確定好的這條直線中我們怎麼確定直線中的某個點呢,這樣我們就可以根據兩條平行線之間的關系來建立該兩點確定的直線方程。看下面的圖(圖是我在網上搬的,勿噴,侵刪)
所以這樣就能得該直線方程瞭。
這裡給出線性插值算法的代碼:
def linear_inter(self): """線性插值算法""" new_img = np.zeros((self.goal_h, self.goal_w, self.goal_channel), dtype=np.uint8) for i in range(0, new_img.shape[0]): for j in range(0, new_img.shape[1]): src_i, src_j = self.convert2src_axes(i, j) if ((src_j - int(src_j)) == 0 and (src_i - int(src_i)) == 0) \ or ((src_i - int(src_i)) == 0) \ or ( (src_j - int(src_j)) == 0): # 表明border_src_j是一個整數 如果是整數,直接將原圖的灰度值付給目標圖像即可 new_img[i, j] = self.img[round(src_i), round(src_j)] # 直接將原圖的灰度值賦值給目標圖像即可 else: """否則進行向上和向下取整,然後進行(一維)線性插值算法""" src_i, src_j = self.convert2src_axes(i, j) j1 = int(src_j) # 向下取整 j2 = math.ceil(src_j) # 向上取整 new_img[i, j] = (j2 - src_j)*self.img[round(src_i), j1] + (src_j-j1)*self.img[round(src_i), j2] return new_img
雙線性插值
雙線性插值算法實現起來就比線性插值算法要難一些,本質上還是和線性插值算法一樣,都是為瞭去找更目標圖像映射到原圖像上的點,把這個點周圍盡像素盡可能多的信息傳遞到這個點中。線性插值使用到瞭周圍的兩個點,這裡的雙線性就是使用到瞭周圍的四個點的像素值信息,所以理論上來說,采用雙線線性插值算法的效果會明顯優於線性插值算法(單線性插值算法)。更具體的圖示我這裡就不畫出來瞭,就使用我自己手繪的簡化版(因為懶),網上其他帖子講得更加明白,這裡隻貼出我自己的理解,這樣以後來看的時候也能一目瞭然。
圖中的x就是圖片的高h,y就是圖片的寬度w,因為我們讀出來的圖片是彩色圖片,所以這樣操作的同時,圖片的三個通道也會同步更新這樣的操作。
下面給出雙線性插值算法的python代碼:
def double_linear_inter(self): new_img = np.zeros((self.goal_h-2, self.goal_w-2, self.goal_channel), dtype=np.uint8) # 這裡減的目的就是為瞭可以進行向上向下取整 for i in range(0, new_img.shape[0]): # 減1的目的就是為瞭可以進行向上向下取整數 for j in range(0, new_img.shape[1]): inner_src_i, inner_src_j = self.convert2src_axes(i, j) # 將取得到的變量參數坐標映射到原圖中,並且返回映射到原圖中的坐標 inner_i1 = int(inner_src_i) # 對行進行向上向下取整數 inner_i2 = math.ceil(inner_src_i) inner_j1 = int(inner_src_j) # 對列進行向上向下取整數 inner_j2 = math.ceil(inner_src_j) Q1 = (inner_j2 - inner_src_j) * self.img[inner_i1, inner_j1] + \ (inner_src_j - inner_j1) * self.img[inner_i1, inner_j2] Q2 = (inner_j2 - inner_src_j) * self.img[inner_i2, inner_j1] + \ (inner_src_j - inner_j1) * self.img[inner_i2, inner_j2] new_img[i, j] = (inner_i2 - inner_src_i) * Q1 + (inner_src_i - inner_i1) * Q2 return new_img
三種插值算法的綜合使用
在進行圖片處理的時候,對於圖片的四個角,沒有四個或者兩個鄰域,所以四個角就采用最鄰近插值算法實現,二上下,左右這幾行,沒有四個鄰域,所以采用線性插值算法,其餘中心部分每個點都有四個鄰域,所以采用雙線性插值算法來實現,這樣的圖片效果會更好,以下是這三種算法的具體實現:
def all_transform(self): # source_h, source_w, source_channel = self.img.shape # 獲得原圖像的高度,寬度和通道量 # goal_h, goal_w = round(source_h*self.h_rate), round(source_w*self.w_rate) # 將原圖像的size進行按照傳入進來的rate參數等比的放大 # goal_channel = source_channel """進行圖像轉換瞭""" new_img = np.zeros((self.goal_h-1, self.goal_w-1, self.goal_channel), dtype=np.uint8) # 得到一個空的數組用來存放轉換後的值,即為新的圖片 """邊界使用線性插值算法""" temp_row = [0, new_img.shape[0]-1] # 上下兩行進行線性插值 for i in temp_row: # i -> h -> x # j -> w -> y for j in range(0, new_img.shape[1]): """邊界線(除瞭四個角落)采用線性插值法""" t_border_src_i, t_border_src_j = self.convert2src_axes(i, j) if ((t_border_src_j - int(t_border_src_j)) == 0 and (t_border_src_i - int(t_border_src_i)) == 0) \ or (t_border_src_i - int(t_border_src_i)) == 0 \ or (t_border_src_j - int(t_border_src_j)) == 0: # 表明t_border_src_j是一個整數 如果是整數,直接將原圖的灰度值付給目標圖像即可 new_img[i, j] = self.img[round(t_border_src_i), round(t_border_src_j)] # 直接將原圖的灰度值賦值給目標圖像即可 else: """否則進行向上和向下取整,然後進行(一維)線性插值算法""" t_border_src_i, t_border_src_j = self.convert2src_axes(i, j) j1 = int(t_border_src_j) # 向下取整 j2 = math.ceil(t_border_src_j) # 向上取整 new_img[i, j] = self.img[round(t_border_src_i), j1] + (t_border_src_j - j1) * \ (self.img[round(t_border_src_i), j2] - self.img[round(t_border_src_i), j1]) # 左右兩列進行線性插值 temp_col = [0, new_img.shape[1]-1] for i in temp_col: # i -> w -> y # j -> h -> x for j in range(0, new_img.shape[0]): """邊界線(除瞭四個角落)采用線性插值法""" t_border_src_i, t_border_src_j = self.convert2src_axes(i, j) if ((t_border_src_j - int(t_border_src_j)) == 0 and (t_border_src_i - int(t_border_src_i)) == 0) \ or (t_border_src_i - int(t_border_src_i)) == 0 \ or (t_border_src_j - int(t_border_src_j)) == 0: # 表明border_src_j是一個整數 如果是整數,直接將原圖的灰度值付給目標圖像即可 new_img[j, i] = self.img[round(t_border_src_i), round(t_border_src_j)] # 直接將原圖的灰度值賦值給目標圖像即可 else: """否則進行向上和向下取整,然後進行(一維)線性插值算法""" t_border_src_i, t_border_src_j = self.convert2src_axes(j, i) j1 = int(t_border_src_i) # 向下取整 j2 = math.ceil(t_border_src_i) # 向上取整 new_img[j, i] = self.img[j1, round(t_border_src_j)] + (t_border_src_i - j1) * \ (self.img[j2, round(t_border_src_j)] - self.img[j1, round(t_border_src_j)]) """四個角落(頂點)使用最臨近插值算法""" corner_low = [0, new_img.shape[0]-1] corner_height = [0, new_img.shape[1]-1] for i in corner_low: for j in corner_height: src_i, src_j = self.convert2src_axes(i, j) new_img[i, j] = self.img[round(src_i), round(src_j)] """中間的使用雙線性插值法""" for i in range(1, new_img.shape[0] - 1): # 減1的目的就是為瞭可以進行向上向下取整數 for j in range(1, new_img.shape[1] - 1): inner_src_i, inner_src_j = self.convert2src_axes(i, j) # 將取得到的變量參數坐標映射到原圖中,並且返回映射到原圖中的坐標 inner_i1 = int(inner_src_i) # 對行進行向上向下取整數 inner_i2 = math.ceil(inner_src_i) inner_j1 = int(inner_src_j) # 對列進行向上向下取整數 inner_j2 = math.ceil(inner_src_j) Q1 = (inner_j2 - inner_src_j) * self.img[inner_i1, inner_j1] + \ (inner_src_j - inner_j1) * self.img[inner_i1, inner_j2] Q2 = (inner_j2 - inner_src_j) * self.img[inner_i2, inner_j1] + \ (inner_src_j - inner_j1) * self.img[inner_i2, inner_j2] new_img[i, j] = (inner_i2 - inner_src_i) * Q1 + (inner_src_i - inner_i1) * Q2 return new_img
附件
下面附上我這個插值算法所有的代碼
import numpy as np
import cv2
import math
import logging
class Image_inter_lines(object):
"""設置傳入參數"""
def __init__(self, img, modify_size=(3, 3), *, align='center'):
self.img = img
self.h_rate = modify_size[0] # 高度的縮放率
self.w_rate = modify_size[1] # 寬度的縮放率
self.align = align # 設置居中模式,進而進行判斷其對齊方式
self.source_h = img.shape[0] # 對應 i 列 -> x
self.source_w = img.shape[1] # 對飲 j 列 -> y
self.goal_channel = img.shape[2] # 通道數
self.goal_h = round(self.source_h * self.h_rate) # 將原圖像的size進行按照傳入進來的rate參數等比的放大
self.goal_w = round(self.source_w * self.w_rate)
if self.align not in ['center', 'left']:
logging.exception(f'{self.align} is not a valid align parameter')
self.align = 'center' # 如果傳入的參數不是居中或者居左,則強制將其置為居中
pass
def set_rate(self, new_modify_size=(None, None)):
self.h_rate = new_modify_size[0]
self.w_rate = new_modify_size[1]
def convert2src_axes(self, des_x, des_y):
if self.align == 'left': # 左對齊
src_x = float(des_x * (self.source_w / self.goal_w))
src_y = float(des_y * (self.source_h / self.goal_h))
src_x = min((self.source_h - 1), src_x)
src_y = min((self.source_w - 1), src_y)
else: # 幾何中心對齊
src_x = float(des_x * (self.source_w / self.goal_w) + 0.5 * (self.source_w / self.goal_w))
src_y = float(des_y * (self.source_h / self.goal_h) + 0.5 * (self.source_h / self.goal_h))
src_x = min((self.source_h - 1), src_x)
src_y = min((self.source_w - 1), src_y)
return src_x, src_y # 這裡返回的數值可以是小數,也可能是整數例如:23.00,但是這個數仍然是小數
def nearest_inter(self):
"""最鄰近插值法"""
new_img = np.zeros((self.goal_h, self.goal_w, self.goal_channel), dtype=np.uint8)
for i in range(0, new_img.shape[0]):
for j in range(0, new_img.shape[1]):
src_i, src_j = self.convert2src_axes(i, j)
new_img[i, j] = self.img[round(src_i), round(src_j)]
return new_img
def linear_inter(self):
"""線性插值算法"""
new_img = np.zeros((self.goal_h, self.goal_w, self.goal_channel), dtype=np.uint8)
for i in range(0, new_img.shape[0]):
for j in range(0, new_img.shape[1]):
src_i, src_j = self.convert2src_axes(i, j)
if ((src_j - int(src_j)) == 0 and (src_i - int(src_i)) == 0) \
or ((src_i - int(src_i)) == 0) \
or (
(src_j - int(src_j)) == 0): # 表明t_border_src_j是一個整數 如果是整數,直接將原圖的灰度值付給目標圖像即可
new_img[i, j] = self.img[round(src_i), round(src_j)] # 直接將原圖的灰度值賦值給目標圖像即可
else:
"""否則進行向上和向下取整,然後進行(一維)線性插值算法"""
src_i, src_j = self.convert2src_axes(i, j)
j1 = int(src_j) # 向下取整
j2 = math.ceil(src_j) # 向上取整
new_img[i, j] = (j2 - src_j)*self.img[round(src_i), j1] + (src_j-j1)*self.img[round(src_i), j2]
return new_img
def double_linear_inter(self):
new_img = np.zeros((self.goal_h-2, self.goal_w-2, self.goal_channel), dtype=np.uint8) # 這裡減的目的就是為瞭可以進行向上向下取整
for i in range(0, new_img.shape[0]): # 減1的目的就是為瞭可以進行向上向下取整數
for j in range(0, new_img.shape[1]):
inner_src_i, inner_src_j = self.convert2src_axes(i, j) # 將取得到的變量參數坐標映射到原圖中,並且返回映射到原圖中的坐標
inner_i1 = int(inner_src_i) # 對行進行向上向下取整數
inner_i2 = math.ceil(inner_src_i)
inner_j1 = int(inner_src_j) # 對列進行向上向下取整數
inner_j2 = math.ceil(inner_src_j)
Q1 = (inner_j2 - inner_src_j) * self.img[inner_i1, inner_j1] + \
(inner_src_j - inner_j1) * self.img[inner_i1, inner_j2]
Q2 = (inner_j2 - inner_src_j) * self.img[inner_i2, inner_j1] + \
(inner_src_j - inner_j1) * self.img[inner_i2, inner_j2]
new_img[i, j] = (inner_i2 - inner_src_i) * Q1 + (inner_src_i - inner_i1) * Q2
return new_img
def all_transform(self):
# source_h, source_w, source_channel = self.img.shape # 獲得原圖像的高度,寬度和通道量
# goal_h, goal_w = round(source_h*self.h_rate), round(source_w*self.w_rate) # 將原圖像的size進行按照傳入進來的rate參數等比的放大
# goal_channel = source_channel
"""進行圖像轉換瞭"""
new_img = np.zeros((self.goal_h-1, self.goal_w-1, self.goal_channel), dtype=np.uint8) # 得到一個空的數組用來存放轉換後的值,即為新的圖片
"""邊界使用線性插值算法"""
temp_row = [0, new_img.shape[0]-1]
# 上下兩行進行線性插值
for i in temp_row:
# i -> h -> x
# j -> w -> y
for j in range(0, new_img.shape[1]):
"""邊界線(除瞭四個角落)采用線性插值法"""
t_border_src_i, t_border_src_j = self.convert2src_axes(i, j)
if ((t_border_src_j - int(t_border_src_j)) == 0 and (t_border_src_i - int(t_border_src_i)) == 0) \
or (t_border_src_i - int(t_border_src_i)) == 0 \
or (t_border_src_j - int(t_border_src_j)) == 0: # 表明t_border_src_j是一個整數 如果是整數,直接將原圖的灰度值付給目標圖像即可
new_img[i, j] = self.img[round(t_border_src_i), round(t_border_src_j)] # 直接將原圖的灰度值賦值給目標圖像即可
else:
"""否則進行向上和向下取整,然後進行(一維)線性插值算法"""
t_border_src_i, t_border_src_j = self.convert2src_axes(i, j)
j1 = int(t_border_src_j) # 向下取整
j2 = math.ceil(t_border_src_j) # 向上取整
new_img[i, j] = self.img[round(t_border_src_i), j1] + (t_border_src_j - j1) * \
(self.img[round(t_border_src_i), j2] - self.img[round(t_border_src_i), j1])
# 左右兩列進行線性插值
temp_col = [0, new_img.shape[1]-1]
for i in temp_col:
# i -> w -> y
# j -> h -> x
for j in range(0, new_img.shape[0]):
"""邊界線(除瞭四個角落)采用線性插值法"""
t_border_src_i, t_border_src_j = self.convert2src_axes(i, j)
if ((t_border_src_j - int(t_border_src_j)) == 0 and (t_border_src_i - int(t_border_src_i)) == 0) \
or (t_border_src_i - int(t_border_src_i)) == 0 \
or (t_border_src_j - int(t_border_src_j)) == 0: # 表明border_src_j是一個整數 如果是整數,直接將原圖的灰度值付給目標圖像即可
new_img[j, i] = self.img[round(t_border_src_i), round(t_border_src_j)] # 直接將原圖的灰度值賦值給目標圖像即可
else:
"""否則進行向上和向下取整,然後進行(一維)線性插值算法"""
t_border_src_i, t_border_src_j = self.convert2src_axes(j, i)
j1 = int(t_border_src_i) # 向下取整
j2 = math.ceil(t_border_src_i) # 向上取整
new_img[j, i] = self.img[j1, round(t_border_src_j)] + (t_border_src_i - j1) * \
(self.img[j2, round(t_border_src_j)] - self.img[j1, round(t_border_src_j)])
"""四個角落(頂點)使用最臨近插值算法"""
corner_low = [0, new_img.shape[0]-1]
corner_height = [0, new_img.shape[1]-1]
for i in corner_low:
for j in corner_height:
src_i, src_j = self.convert2src_axes(i, j)
new_img[i, j] = self.img[round(src_i), round(src_j)]
"""中間的使用雙線性插值法"""
for i in range(1, new_img.shape[0] - 1): # 減1的目的就是為瞭可以進行向上向下取整數
for j in range(1, new_img.shape[1] - 1):
inner_src_i, inner_src_j = self.convert2src_axes(i, j) # 將取得到的變量參數坐標映射到原圖中,並且返回映射到原圖中的坐標
inner_i1 = int(inner_src_i) # 對行進行向上向下取整數
inner_i2 = math.ceil(inner_src_i)
inner_j1 = int(inner_src_j) # 對列進行向上向下取整數
inner_j2 = math.ceil(inner_src_j)
Q1 = (inner_j2 - inner_src_j) * self.img[inner_i1, inner_j1] + \
(inner_src_j - inner_j1) * self.img[inner_i1, inner_j2]
Q2 = (inner_j2 - inner_src_j) * self.img[inner_i2, inner_j1] + \
(inner_src_j - inner_j1) * self.img[inner_i2, inner_j2]
new_img[i, j] = (inner_i2 - inner_src_i) * Q1 + (inner_src_i - inner_i1) * Q2
return new_img
if __name__ == '__main__':
pic1 = cv2.imread(r'C:\Users\heshijie\Desktop\Py_study\image processing\Carmen.jpg')
pic2 = cv2.imread(r'C:\Users\heshijie\Desktop\Py_study\image processing\girl.jpg')
pic3 = cv2.imread(r'C:\Users\heshijie\Desktop\Py_study\image processing\architecture.jpg')
Obj_pic1 = Image_inter_lines(pic1, modify_size=(2, 2), align='center')
new_pic1 = Obj_pic1.nearest_inter()
cv2.imshow('origin', pic1)
cv2.imshow('nearest_inter', new_pic1)
new_pic2 = Obj_pic1.linear_inter()
cv2.imshow('liner_inter', new_pic2)
new_pic3 = Obj_pic1.all_transform()
cv2.imshow('double_liner_inter', new_pic3)
cv2.waitKey()
到此這篇關於Python實現RGB等圖片的圖像插值算法的文章就介紹到這瞭,更多相關Python 圖像插值算法內容請搜索WalkonNet以前的文章或繼續瀏覽下面的相關文章希望大傢以後多多支持WalkonNet!
推薦閱讀:
- Python中八大圖像特效算法的示例詳解
- 基於Opencv制作的美顏相機帶你領略美顏特效的效果
- 詳解Python實現圖像分割增強的兩種方法
- 深入瞭解Python Opencv數據增強
- Python實現邊緣提取的示例代碼