-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathcamera_calibration.py
65 lines (54 loc) · 1.97 KB
/
camera_calibration.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
import numpy as np
import cv2 as cv
# sub pixel find termination criteria
criteria = (cv.TERM_CRITERIA_EPS + cv.TERM_CRITERIA_MAX_ITER, 30, 0.001)
numPointsY = 9
numPointsX = 6
chessboardSize = (numPointsY,numPointsX)
# prepare object points, like (0,0,0), (1,0,0), (2,0,0) ....,(6,5,0)
objp = np.zeros((numPointsY * numPointsX, 3), np.float32)
objp[:,:2] = (np.mgrid[0:numPointsY, 0:numPointsX].T.reshape(-1, 2)) * 26
# Arrays to store object points and image points from all the images.
objpoints = [] # 3d point in real world space
imgpoints = [] # 2d points in image plane.
captured = 0
imageShape = None
cap = cv.VideoCapture(0)
if not cap.isOpened():
print("Cannot open camera")
exit()
else:
print("capturing images")
while captured < 50:
# Capture frame-by-frame
ret, img = cap.read()
# if frame is read correctly ret is True
if not ret:
print("Can't receive frame (stream end?). Exiting ...")
break
gray = cv.cvtColor(img, cv.COLOR_BGR2GRAY)
imageShape = gray.shape
# Find the chess board corners
ret, corners = cv.findChessboardCorners(gray, chessboardSize, None, cv.CALIB_CB_ADAPTIVE_THRESH)
# If found, draw image points (after refining them)
if ret == True:
objpoints.append(objp)
corners2 = cv.cornerSubPix(gray, corners, (11,11), (-1,-1), criteria)
imgpoints.append(corners)
captured = captured + 1
# Draw and display the corners
cv.drawChessboardCorners(img, chessboardSize, corners2, ret)
print(f"frame {captured}")
#cv.imwrite('./image.png', img)
print("Frame")
cv.imshow('img', img)
cv.waitKey(500)
#cv.destroyAllWindows()
print("images captured")
print("calibrating")
ret, mtx, dist, rvecs, tvecs = cv.calibrateCamera(objpoints, imgpoints, imageShape[::-1], None, None)
if ret:
print("calibration successful, saving values")
np.savez('laptop_calibration.npz', mtx=mtx, dist=dist, rvecs=rvecs, tvecs=tvecs)
else:
print("calibration failed")