OpenCV/Python：相机标定_opencv getperspectivetransform 角点-程序员宅基地

这也不是相机标定，只是了解相机标定之前的基础部分

import numpy as np
import cv2
import glob

'''
  在这里，我的棋盘格是8*8的，所以角点个数为7*7,当然棋盘格的行列个数可以不一样；
  如果想方便代码改变棋盘格数，是以定义两个变量w（列角点数）和h（行角点数），注意如果角点维数超出的话，标定的时候会报错。
'''
w = 7
h = 7
# glob是个文件名管理工具
images = glob.glob('test5/5test12.jpg')
print('...loading')
for fname in images:
    # 对每张图片，识别出角点，记录世界物体坐标和图像坐标
    print(f'processing img:{fname}')
    img = cv2.imread(fname)
    gray = cv2.cvtColor(img, cv2.COLOR_BGR2GRAY)  # 转灰度
    print('grayed')
    # 寻找角点，存入corners，ret是找到角点的flag
    ret, corners = cv2.findChessboardCornersSB(gray, (w, h),None)
    # 如果找到，添加对象点，图像点（精炼后）
    if ret == True:
        # 绘制并显示角点
        # img = cv2.drawChessboardCorners(img, (w,h), corners,ret)

        # cv2.getPerspectiveTransform()
        a = [10, 20, 30, 40, 50, 60, 70]
        dst_pts = np.zeros((49, 2))
        for i in range(7):
            for j in range(7):
                dst_pts[7 * i + j, :] = (a[j], a[i])
        dst_pts = dst_pts.reshape(-1, 1, 2)
        H, mask = cv2.findHomography(corners, dst_pts)
        # H, mask = cv2.findHomography(corners, dst_pts, cv2.RANSAC, 5.0)
        wrap = cv2.warpPerspective(img, H, (80, 80))

        img = cv2.drawChessboardCorners(img, (w,h), corners,ret)
        cv2.namedWindow('img', 0)
        cv2.resizeWindow('img', 500, 500)
        cv2.imshow('img', img)
        cv2.namedWindow('wrap', 0)
        cv2.resizeWindow('wrap', 500, 500)
        cv2.imshow('wrap', wrap)
        cv2.waitKey()
        cv2.destroyAllWindows()

先看一下结果

ret, corners = cv2.findChessboardCornersSB(gray, (w, h),None)
	input
			gray 		灰度图
			(w, h) 		内角点的行列数
	output
			ret 		True or False
			corners		角点坐标
	

img = cv2.drawChessboardCorners(img, (w,h), corners,ret)
	input
			img			待处理图像
			(w,h)		内角点的行列数
			corners		角点坐标
			ret 		True or False
	output
			img 		画了角点坐标后的图像

H, mask = cv2.findHomography(corners, dst_pts)
	input
			corners		角点坐标，shape=(49, 1, 2)
			dst_pts 	这些角点在原始图像中的坐标，shape=(49, 1, 2)
	output
			H 			变换矩阵

也可以不用findHomography，而使用getPerspectiveTransform测试一下

corners_2D=corners.reshape(-1,2)  # 将[49*1*2]降维到[49*2]
corners_2D_4=corners_2D[5:9,:]  # 取四个点
H = cv2.getPerspectiveTransform(np.array(corners_2D_4, np.float32), np.array(dst_pts_2D_4, np.float32))

要注意的是，getPerspectiveTransform的输入应该为二维的[4*2]的np.float32类型，且返回值只有一个

上面从49个角点中选择的连续四个点，很容易出现问题，因为连续的四个点很容易共线，构不成四边形，或者效果很差，因为是局部的四个点，可以像下面这样，取49个点的四个顶点，这四个点构成的四边形是面积最大的，误差就会小点

corners_2D_4=[corners_2D[0,:],corners_2D[6,:],corners_2D[42,:],corners_2D[48,:]]
dst_pts_2D_4 = [dst_pts_2D[0, :], dst_pts_2D[6, :], dst_pts_2D[42, :], dst_pts_2D[48, :]]

本文链接：https://blog.csdn.net/cd_yourheart/article/details/119423126

原作者删帖不实内容删帖广告或垃圾文章投诉

智能推荐

从零开始搭建Hadoop_创建一个hadoop项目-程序员宅基地

文章浏览阅读331次。第一部分：准备工作1 安装虚拟机2 安装centos73 安装JDK以上三步是准备工作，至此已经完成一台已安装JDK的主机第二部分：准备３台虚拟机以下所有工作最好都在root权限下操作1 克隆上面已经有一台虚拟机了,现在对master进行克隆,克隆出另外2台子机;1.1 进行克隆21.2 下一步1.3 下一步1.4 下一步1.5 根据子机需要,命名和安装路径1.6 ..._创建一个hadoop项目

心脏滴血漏洞HeartBleed CVE-2014-0160深入代码层面的分析_heartbleed代码分析-程序员宅基地

文章浏览阅读1.7k次。心脏滴血漏洞HeartBleed CVE-2014-0160 是由heartbeat功能引入的，本文从深入码层面的分析该漏洞产生的原因_heartbleed代码分析

java读取ofd文档内容_ofd电子文档内容分析工具（分析文档、签章和证书）-程序员宅基地

文章浏览阅读1.4k次。前言ofd是国家文档标准，其对标的文档格式是pdf。ofd文档是容器格式文件，ofd其实就是压缩包。将ofd文件后缀改为.zip，解压后可看到文件包含的内容。ofd文件分析工具下载：点我下载。ofd文件解压后，可以看到如下内容：对于xml文件，可以用文本工具查看。但是对于印章文件(Seal.esl)、签名文件(SignedValue.dat)就无法查看其内容了。本人开发一款ofd内容查看器，..._signedvalue.dat

基于FPGA的数据采集系统（一）_基于fpga的信息采集-程序员宅基地

文章浏览阅读1.8w次，点赞29次，收藏313次。整体系统设计本设计主要是对ADC和DAC的使用，主要实现功能流程为：首先通过串口向FPGA发送控制信号，控制DAC芯片tlv5618进行DA装换，转换的数据存在ROM中，转换开始时读取ROM中数据进行读取转换。其次用按键控制adc128s052进行模数转换100次，模数转换数据存储到FIFO中，再从FIFO中读取数据通过串口输出显示在pc上。其整体系统框图如下：图1：FPGA数据采集系统框图从图中可以看出，该系统主要包括9个模块：串口接收模块、按键消抖模块、按键控制模块、ROM模块、D.._基于fpga的信息采集

微服务 spring cloud zuul com.netflix.zuul.exception.ZuulException GENERAL-程序员宅基地

文章浏览阅读2.5w次。1.背景错误信息：-- [http-nio-9904-exec-5] o.s.c.n.z.filters.post.SendErrorFilter : Error during filteringcom.netflix.zuul.exception.ZuulException: Forwarding error at org.springframework.cloud..._com.netflix.zuul.exception.zuulexception

邻接矩阵-建立图-程序员宅基地

文章浏览阅读358次。1.介绍图的相关概念　　图是由顶点的有穷非空集和一个描述顶点之间关系-边（或者弧）的集合组成。通常，图中的数据元素被称为顶点，顶点间的关系用边表示，图通常用字母G表示，图的顶点通常用字母V表示，所以图可以定义为:　　G=(V,E)其中，V(G)是图中顶点的有穷非空集合，E(G)是V(G)中顶点的边的有穷集合1.1 无向图：图中任意两个顶点构成的边是没有方向的1.2 有向图：图中..._给定一个邻接矩阵未必能够造出一个图

随便推点

MDT2012部署系列之11 WDS安装与配置-程序员宅基地

文章浏览阅读321次。（十二）、WDS服务器安装通过前面的测试我们会发现，每次安装的时候需要加域光盘映像，这是一个比较麻烦的事情，试想一个上万个的公司，你天天带着一个光盘与光驱去给别人装系统，这将是一个多么痛苦的事情啊，有什么方法可以解决这个问题了？答案是肯定的，下面我们就来简单说一下。WDS服务器，它是Windows自带的一个免费的基于系统本身角色的一个功能，它主要提供一种简单、安全的通过网络快速、远程将Window..._doc server2012上通过wds+mdt无人值守部署win11系统.doc

python--xlrd/xlwt/xlutils_xlutils模块可以读xlsx吗-程序员宅基地

文章浏览阅读219次。python–xlrd/xlwt/xlutilsxlrd只能读取，不能改,支持 xlsx和xls 格式xlwt只能改，不能读xlwt只能保存为.xls格式xlutils能将xlrd.Book转为xlwt.Workbook，从而得以在现有xls的基础上修改数据，并创建一个新的xls，实现修改xlrd打开文件import xlrdexcel=xlrd.open_workbook('E:/test.xlsx') 返回值为xlrd.book.Book对象,不能修改获取sheett_xlutils模块可以读xlsx吗

关于新版本selenium定位元素报错：‘WebDriver‘ object has no attribute ‘find_element_by_id‘等问题_unresolved attribute reference 'find_element_by_id-程序员宅基地

文章浏览阅读8.2w次，点赞267次，收藏656次。运行Selenium出现'WebDriver' object has no attribute 'find_element_by_id'或AttributeError: 'WebDriver' object has no attribute 'find_element_by_xpath'等定位元素代码错误，是因为selenium更新到了新的版本，以前的一些语法经过改动。..............._unresolved attribute reference 'find_element_by_id' for class 'webdriver