使用python opencv對畸變圖像進(jìn)行矯正的實(shí)現(xiàn)
代碼:
__Author__ = "Shliang"
__Email__ = "shliang0603@gmail.com"
import os
import cv2
import numpy as np
from tqdm import tqdm
def undistort(frame):
fx = 685.646752
cx = 649.107905
fy = 676.658033
cy = 338.054431
k1, k2, p1, p2, k3 = -0.363219, 0.093818, 0.006178, -0.003714, 0.0
# 相機(jī)坐標(biāo)系到像素坐標(biāo)系的轉(zhuǎn)換矩陣
k = np.array([
[fx, 0, cx],
[0, fy, cy],
[0, 0, 1]
])
# 畸變系數(shù)
d = np.array([
k1, k2, p1, p2, k3
])
h, w = frame.shape[:2]
mapx, mapy = cv2.initUndistortRectifyMap(k, d, None, k, (w, h), 5)
return cv2.remap(frame, mapx, mapy, cv2.INTER_LINEAR)
# 對攝像頭實(shí)時(shí)視頻流做畸變矯正
def distortion_correction_cam():
cap = cv2.VideoCapture(0)
while (cap.isOpened()):
ret, frame = cap.read()
undistort_frame = undistort(frame)
compare = np.hstack((frame, undistort_frame))
cv2.imshow('frame', compare)
if cv2.waitKey(1) & 0xFF == ord('q'):
break
cap.release()
cv2.destroyAllWindows()
# 對目錄下的所有圖片做畸變矯正,并把畸變矯正后的圖片保存下來
def distortion_correction_imgs(input_dir, output_dir):
in_imgs = os.listdir(input_dir)
for img_name in tqdm(in_imgs):
image = cv2.imread(os.path.join(input_dir, img_name))
distroted_img = undistort(image)
cv2.imwrite(os.path.join(output_dir, img_name), distroted_img)
if __name__ == '__main__':
input_dir = "/home/shl/extract_rosbag_data/0324_bags/plycal_calib/root/images"
output_dir = "/home/shl/extract_rosbag_data/0324_bags/plycal_calib/root/distro_imgs"
# distortion_correction_imgs(input_dir, output_dir)
distortion_correction_cam()對圖片進(jìn)行矯正效果:
原圖:

矯正后的圖片:

采集的攝像頭畫面矯正效果:

從上面的換面可以看到,左邊是未矯正的畫面,右邊是矯正后的畫面:
矯正后的畫面會被裁切,明顯可以看到畫面中的信息是有裁切的,例如左邊的椅子已經(jīng)被裁切掉了- 矯正后的畫面和原畫面是保持相同的分辨率:
640x480,但是,為什么會看到畫面會出現(xiàn)橫向的拉伸,這是因?yàn)?code>標(biāo)定相機(jī)內(nèi)參的時(shí)候畫面的分辨率設(shè)置為1280x720=16:9,但是opencv讀取攝像頭默認(rèn)的分辨率卻是640x480=4:3,兩者的比例都不一樣,所以肯定會出現(xiàn)拉伸
解決拉伸的方式,就是把讀取攝像頭的時(shí)候,把攝像頭的分辨率設(shè)置成和標(biāo)定的時(shí)候一樣的分辨率,設(shè)置為1280x720,下面是如何在opencv讀取攝像頭的時(shí)候設(shè)置攝像頭分辨率:
# 對攝像頭實(shí)時(shí)視頻流做畸變矯正
def distortion_correction_cam():
cap = cv2.VideoCapture(0)
# 獲取攝像頭讀取畫面的寬和高
width = cap.get(3)
height = cap.get(4)
fps = cap.get(5)
print(width, height, fps) # 640.0 480.0 30.0
# 在這里把攝像頭的分辨率修改為和我們標(biāo)定時(shí)使用的一樣的分辨率 1280x720
cap.set(3, 1280)
cap.set(4, 720)
width = cap.get(3)
height = cap.get(4)
print(width, height, fps) # 1280.0 720.0 30.0
while (cap.isOpened()):
ret, frame = cap.read()
print(frame.shape)
undistort_frame = undistort(frame)
compare = np.hstack((frame, undistort_frame))
cv2.imshow('frame', compare)
if cv2.waitKey(1) & 0xFF == ord('q'):
break
cap.release()
cv2.destroyAllWindows()重新設(shè)置分辨率后,矯正前后畫面對比,可以看到幾乎是沒有橫向或縱向拉伸的!

參考:
https://blog.csdn.net/weixin_40516558/article/details/103494029
https://blog.csdn.net/guaiderzhu1314/article/details/96306509
https://www.codenong.com/cs110623399/
到此這篇關(guān)于使用python opencv對畸變圖像進(jìn)行矯正的實(shí)現(xiàn)的文章就介紹到這了,更多相關(guān)python opencv畸變圖像矯正內(nèi)容請搜索腳本之家以前的文章或繼續(xù)瀏覽下面的相關(guān)文章希望大家以后多多支持腳本之家!
相關(guān)文章
python神經(jīng)網(wǎng)絡(luò)MobileNetV2模型的復(fù)現(xiàn)詳解
這篇文章主要為大家介紹了python神經(jīng)網(wǎng)絡(luò)MobileNetV2模型的復(fù)現(xiàn)詳解,有需要的朋友可以借鑒參考下,希望能夠有所幫助,祝大家多多進(jìn)步,早日升職加薪2022-05-05
Python基于字典實(shí)現(xiàn)switch case函數(shù)調(diào)用
這篇文章主要介紹了Python基于字典實(shí)現(xiàn)switch case函數(shù)調(diào)用,文中通過示例代碼介紹的非常詳細(xì),對大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價(jià)值,需要的朋友可以參考下2020-07-07
教你如何把Python CSV 合并到多個(gè)sheet工作表
今天我們要實(shí)現(xiàn)將多個(gè)CSV文件,合并到一個(gè)Excel文件中的,多個(gè)sheet工作表,網(wǎng)上大多方法都是將csv直接合并在一起,也不分別創(chuàng)建sheet表。接下來通過本文給大家詳細(xì)介紹下,需要的朋友參考下吧2021-11-11
pytorch 如何實(shí)現(xiàn)HWC轉(zhuǎn)CHW
這篇文章主要介紹了pytorch HWC轉(zhuǎn)CHW的實(shí)現(xiàn)方式,具有很好的參考價(jià)值,希望對大家有所幫助。如有錯(cuò)誤或未考慮完全的地方,望不吝賜教2021-05-05
Matplotlib實(shí)戰(zhàn)之折線圖繪制詳解
折線圖是一種用于可視化數(shù)據(jù)變化趨勢的圖表,它可以用于表示任何數(shù)值隨著時(shí)間或類別的變化,本文主要介紹了如何利用Matplotlib實(shí)現(xiàn)折線圖的繪制,感興趣的可以了解下2023-08-08
Django在pycharm下修改默認(rèn)啟動端口的方法
今天小編就為大家分享一篇Django在pycharm下修改默認(rèn)啟動端口的方法,具有很好的參考價(jià)值,希望對大家有所幫助。一起跟隨小編過來看看吧2019-07-07
Pycharm 如何連接遠(yuǎn)程服務(wù)器并debug調(diào)試
本文主要介紹了Pycharm 如何連接遠(yuǎn)程服務(wù)器并debug調(diào)試,文中通過示例代碼介紹的非常詳細(xì),對大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價(jià)值,需要的朋友們下面隨著小編來一起學(xué)習(xí)學(xué)習(xí)吧2022-06-06

