买男装最好的购物网站,龙岗网站建设技术,门户网站部署方案,深圳建筑工程公司排名import cv2
import numpy as np
import glob# 设置寻找亚像素角点的参数#xff0c;采用的停止准则是最大循环次数30和最大误差容限0.001
criteria (cv2.TERM_CRITERIA_MAX_ITER | cv2.TERM_CRITERIA_EPS, 30, 0.001)# 获取标定板角点的位置
objp np.zeros((6 * 7, 3), np.f…import cv2
import numpy as np
import glob# 设置寻找亚像素角点的参数采用的停止准则是最大循环次数30和最大误差容限0.001
criteria (cv2.TERM_CRITERIA_MAX_ITER | cv2.TERM_CRITERIA_EPS, 30, 0.001)# 获取标定板角点的位置
objp np.zeros((6 * 7, 3), np.float32)
objp[:, :2] np.mgrid[0:7, 0:6].T.reshape(-1, 2) # 将世界坐标系建在标定板上所有点的Z坐标全部为0所以只需要赋值x和yobj_points [] # 存储3D点
img_points [] # 存储2D点images glob.glob(C:\Users\Administrator\Desktop\cal_camera\cam\*.bmp)
for fname in images:img cv2.imread(fname)gray cv2.cvtColor(img, cv2.COLOR_BGR2GRAY)size gray.shape[::-1]ret, corners cv2.findChessboardCorners(gray, (7, 6), None)if ret:obj_points.append(objp)corners2 cv2.cornerSubPix(gray, corners, (5, 5), (-1, -1), criteria) # 在原角点的基础上寻找亚像素角点if True:img_points.append(corners2)else:img_points.append(corners)cv2.drawChessboardCorners(img, (7, 6), corners, ret) # 记住OpenCV的绘制函数一般无返回值cv2.imshow(img, img)cv2.waitKey(50)print len(img_points)
cv2.destroyAllWindows()# 标定
ret, mtx, dist, rvecs, tvecs cv2.calibrateCamera(obj_points, img_points, img.shape[:2], None, None)
print 标定
print ret:, ret
print mtx:\n, mtx # 内参数矩阵内参数的输出向量。输出顺序为:
# (fx,fy,cx,cy,
# k1,k2,p1,p2,k3,k4,k5,k6,s1,s2,s3,s4,τx,τy 如果不估计其中某一个参数值等于0
print dist:\n, dist # 畸变系数 distortion cofficients (k_1,k_2,p_1,p_2,k_3)
print rvecs:\n, rvecs # 旋转向量 # 外参数
print tvecs:\n, tvecs # 平移向量 # 外参数print(-----------------------------------------------------)
# 畸变校正
img cv2.imread(C:\Users\Administrator\Desktop\cal_camera\cam\cam2.bmp)
h, w img.shape[:2]#使用cv.getOptimalNewCameraMatrix()优化内参数和畸变系数
# 通过设定自由自由比例因子alpha。当alpha设为0的时候
# 将会返回一个剪裁过的将去畸变后不想要的像素去掉的内参数和畸变系数
# 当alpha设为1的时候将会返回一个包含额外黑色像素点的内参数和畸变系数并返回一个ROI用于将其剪裁掉。newcameramtx, roi cv2.getOptimalNewCameraMatrix(mtx, dist, (w, h), 1, (w, h))
print newcameramtx
print(------------------使用undistort函数-------------------)
#该功能可转换图像以补偿径向和切向镜头失真。
dst cv2.undistort(img, mtx, dist, None, newcameramtx)
x, y, w, h roi
dst1 dst[y:y h, x:x w]
cv2.imwrite(calibresult11.jpg, dst1)
print 方法一:dst的大小为:, dst1.shape# undistort方法二
print(-------------------使用重映射的方式-----------------------)该函数计算联合不变形和整流变换并代表。
以重映射的映射形式。未失真的图像看起来像原始图像就像它一样。
使用相机矩阵 newCameraMatrix和零失真用相机拍摄。如果是。
单目相机newCameraMatrix通常等于cameraMatrix或者它可以通过计算。cv :: getOptimalNewCameraMatrix用于更好地控制缩放。mapx, mapy cv2.initUndistortRectifyMap(mtx, dist, None, newcameramtx, (w, h), 5) # 获取映射方程
# dst cv2.remap(img,mapx,mapy,cv2.INTER_LINEAR) # 重映射
dst cv2.remap(img, mapx, mapy, cv2.INTER_CUBIC) # 重映射后图像变小了
x, y, w, h roi
dst2 dst[y:y h, x:x w]
cv2.imwrite(calibresult11_2.jpg, dst2)
print 方法二:dst的大小为:, dst2.shape # 图像比方法一的小print(-------------------计算反向投影误差-----------------------)
tot_error 0
for i in xrange(len(obj_points)):img_points2, _ cv2.projectPoints(obj_points[i], rvecs[i], tvecs[i], mtx, dist)error cv2.norm(img_points[i], img_points2, cv2.NORM_L2) / len(img_points2)tot_error errormean_error tot_error / len(obj_points)
print total error: , tot_error
print mean error: , mean_error
