圖像的幾何變換從原理上看主要包括兩種:基于2×3矩陣的仿射變換(平移、縮放、旋轉和翻轉等)、基于3×3矩陣的透視變換。
- ?仿射變換
基本的圖像變換就是二維坐標的變換:從一種二維坐標(x,y)到另一種二維坐標(u,v)的線性變換:
如果寫成矩陣的形式,那就是:
作如下定義:
矩陣T(2×3)就稱為仿射變換的變換矩陣,R為線性變換矩陣,t為平移矩陣,簡單來說,仿射變換就是線性變換+平移。變換后直線依然是直線,平行線依然是平行線,直線間的相對位置關系不變,因此非共線的三個對應點便可確定唯一的一個仿射變換,線性變換4個自由度+平移2個自由度→仿射變換自由度為6。
仿射變換在OpenCV中的實現如下:
import cv2
import numpy as np
import matplotlib.pyplot as pltimg = cv2.imread('drawing.jpg')
rows, cols = img.shape[:2]# 變換前的三個點
pts1 = np.float32([[50, 65], [150, 65], [210, 210]])
# 變換后的三個點
pts2 = np.float32([[50, 100], [150, 65], [100, 250]])# 生成變換矩陣,維數:2*3
M = cv2.getAffineTransform(pts1, pts2)
dst = cv2.warpAffine(img, M, (cols, rows))plt.subplot(121), plt.imshow(img), plt.title('input')
plt.subplot(122), plt.imshow(dst), plt.title('output')
plt.show()
?運行結果:
其實平移、旋轉、縮放和翻轉等變換就是對應了不同的仿射變換矩陣,下面分別來看:
(1)平移
平移就是x和y方向上的直接移動,可以上下/左右移動,自由度為2,變換矩陣可以表示為:
# 平移圖片
import numpy as nprows, cols = img.shape[:2]# 定義平移矩陣,需要是numpy的float32類型
# x軸平移100,y軸平移50
M = np.float32([[1, 0, 100], [0, 1, 50]])
dst = cv2.warpAffine(img, M, (cols, rows))cv2.imshow('shift', dst)
cv2.waitKey(0)
?
?
?
(1)通過比例進行縮放
import cv2 as cv
import numpy as np# 圖片縮放
img = cv.imread('images/animal.jpg', flags=1) # flags=1讀取為彩色,flags=0讀取為灰度
cv.imshow('i', img)
h, w, channel = img.shape # 以行列形式存儲, 第幾行到第幾行為圖像高度
dst_h = int(h*0.5)
dst_w = int(w*0.5)
# 最近鄰域差值 雙線性插值 像素關系重采樣 立方差值
dst = cv.resize(img, (dst_w, dst_h)) # 默認雙線性差值
cv.imshow('img', dst)
cv.waitKey(0)
OpenCV提供了resize函數來改變圖像的大小,C++中的函數原型如下:
void resize(InputArray src, OutputArray dst, Size dsize, double fx=0, double fy=0, int interpolation=INTER_LINEAR );
函數參數說明:
src:輸入,原圖像,即待改變大小的圖像;
dst:輸出,改變大小之后的圖像,這個圖像和原圖像具有相同的內容,只是大小和原圖像不一樣而已;
dsize:輸出圖像的大小。如果這個參數不為0,那么就代表將原圖像縮放到這個Size(width,height)指定的大小;如果這個參數為0,那么原圖像縮放之后的大小就要通過下面的公式來計算:
dsize = Size(round(fx*src.cols), round(fy*src.rows))
其中,fx和fy就是下面要說的兩個參數,是圖像width方向和height方向的縮放比例。
fx:width方向的縮放比例,如果它是0,那么它就會按照(double)dsize.width/src.cols來計算;
fy:height方向的縮放比例,如果它是0,那么它就會按照(double)dsize.height/src.rows來計算;
interpolation:這個是指定插值的方式,圖像縮放之后,肯定像素要進行重新計算的,就靠這個參數來指定重新計算像素的方式,有以下幾種:
- INTER_NEAREST?- 最鄰近插值
- INTER_LINEAR?- 雙線性插值,如果最后一個參數你不指定,默認使用這種方法
- INTER_AREA?- resampling using pixel area relation. It may be a preferred method for image decimation, as it gives moire’-free results. But when the image is zoomed, it is similar to the?INTER_NEAREST?method.
- INTER_CUBIC?-?4x4像素鄰域內的雙立方插值
- INTER_LANCZOS4?-?8x8像素鄰域內的Lanczos插值
函數使用說明:
- dsize和fx/fy不能同時為0,要么你就指定好dsize的值,讓fx和fy空置直接使用默認值,就像resize(img, imgDst, Size(30,30));?要么你就讓dsize為0,指定好fx和fy的值,比如fx=fy=0.5,那么就相當于把原圖兩個方向縮小一倍!
- 至于最后的插值方法,正常情況下使用默認的雙線性插值就夠用了。幾種常用方法的效率是:最鄰近插值>雙線性插值>雙立方插值>Lanczos插值;但是效率和效果成反比,所以根據自己的情況酌情使用。
- 正常情況下,在使用之前dst圖像的大小和類型都是不知道的,類型從src圖像繼承而來,大小也是從原圖像根據參數計算出來。但是如果你事先已經指定好dst圖像的大小,那么你可以通過下面這種方式來調用函數:
resize(src, dst, dst.size(), 0, 0, interpolation);
?
(2)通過矩陣變換進行縮放
import cv2 as cv
import numpy as np# 圖片縮放
img = cv.imread('../images/moon.jpg', flags=1) # flags=1讀取為彩色,flags=0讀取為灰度
h, w = img.shape[:2]
mat_shift = np.float32([[0.5, 0, 0], [0, 0.5, 0]]) # 縮放矩陣
dst = cv.warpAffine(img, mat_shift, (int(w/2), int(h/2)))
cv.imshow('img1', img)
cv.imshow('img2', dst)
cv.waitKey(0)
?