Halcon相机标定

1，前言。

        相机的成像过程实质上是坐标系的转换。首先空间中的点由“世界坐标系”转换到“相机坐标系”，然后再将其投影到成像平面（图像物理坐标系），最后再将成像的平面上的数据转换为图像像素坐标系。但是由于透镜的制造精度以及组装工艺的偏差会出现不同的畸变，导致图像失真。

        镜头畸变分为径向畸变和切向畸变。径向畸变是由镜头形状缺陷造成，它又分为枕形畸变和桶形畸变。

        对于一般低精度的场合，只考虑径向畸变即可，因为切向畸变的影响远小于径向畸变。因此如果要进行图像畸变矫正就需要通过标定获取相机的参数。同时为了将图像的坐标系中的像素坐标系与世界坐标系中的坐标距离对应起来，就需要了解相机的外部参数信息，换算其在世界坐标系中的实际距离。相机标定就是获取相机内部参数和外部参数的过程。

2，相机标定的方法。

        传统标定法：标定精度较高，但在非线性相机模型下计算十分复杂，标定速度较低，通常用于精度要求较高的线性系统。

        自标定法：不需要借助参考物，操作简单且计算速度较快，但其标度精度及鲁棒性较低，只适用于速度要求较高精度较低的系统。

        张正友标定：介于前面两种标定法之间，具有操作简单，标定精度高，鲁棒性好的优点，较为广泛使用。

3，内部参数和外部参数。

        内部参数：相机的焦距、畸变系数、像素间距、中心点坐标、图像的宽和高等。这些参数虽然可以从制造商提供的说明书中查询，单其精度性不一定满足要求，因此只能作为参考，在实际中还需要对其进行标定，通过标定可以对这些参数进行修正。不同的相机，其内部参数也不相同，这取决于相机和镜头的类型，以及选择的畸变类型。

        外部参数：相机的外部参数描述相机坐标系与世界坐标系的关系，即相机在世界坐标系中的三位位置，如相机的X轴坐标，Y轴坐标，Z轴坐标，以及相机的朝向。

4，相机的畸变。

        径向畸变：矢量端点沿长度方向的变化，即矢量长度的变化。

        切向畸变：矢量端点沿切线方向发生的变化，即矢量角度的变化。

5，Halcon标定。

1，准备标定板。

用halcon代码生成用于打印标定板的描述文件及标定板图像的ps文件。

*生成标定板文件
 gen_caltab (7, 7, 0.0125, 0.5, 'caltab.descr', 'caltab.ps')

生成的.descr与caltab.ps文件。

备注：

caltab.descr：对标记点的半径，坐标，行，列进行描述的文件。

有.cpd与.descr两种扩展 名的标定板描述文件：

    .cpd描述文件：使用create_caltab函数生成，这种格式适用于六边形排列的标定板，通常在深色背景上有浅色的标记。

     .descr文件：使用gen_caltab函数生成，用于矩形排列的标定板，通常用于浅色背景下深色标记。

Halcon已提供的.cpd与.descr文件位于：D:\Program Files\MVTec\HALCON-12.0\calib目录下。

caltab.descr文件详情如下：

# Plate Description Version 2
# HALCON Version 12.0 --  Sun Dec 29 11:09:19 2024
# Description of the standard calibration plate
# used for the camera calibration in HALCON
# (generated by gen_caltab)
#
#

# 7 rows x 7 columns
# Width, height of calibration plate [meter]: 0.1, 0.1
# Distance between mark centers [meter]: 0.0125

# Number of marks in y-dimension (rows)
r 7

# Number of marks in x-dimension (columns)
c 7

#   offset of coordinate system in z-dimension [meter] (optional):
z 0

# Rectangular border (rim and black frame) of calibration plate
#   rim of the calibration plate (min x, max y, max x, min y) [meter]:
o -0.05125 0.05125 0.05125 -0.05125
#   outer border of the black frame (min x, max y, max x, min y) [meter]:
i -0.05 0.05 0.05 -0.05 
#   triangular corner mark given by two corner points (x,y, x,y) [meter]
#   (optional):
t -0.05 -0.0375 -0.0375 -0.05

#   width of the black frame [meter]:
w 0.003125

# calibration marks:  x y radius [meter]

# calibration marks at y = -0.0375 m
-0.0375 -0.0375 0.003125
-0.025 -0.0375 0.003125
-0.0125 -0.0375 0.003125
0 -0.0375 0.003125
0.0125 -0.0375 0.003125
0.025 -0.0375 0.003125
0.0375 -0.0375 0.003125

# calibration marks at y = -0.025 m
-0.0375 -0.025 0.003125
-0.025 -0.025 0.003125
-0.0125 -0.025 0.003125
0 -0.025 0.003125
0.0125 -0.025 0.003125
0.025 -0.025 0.003125
0.0375 -0.025 0.003125

# calibration marks at y = -0.0125 m
-0.0375 -0.0125 0.003125
-0.025 -0.0125 0.003125
-0.0125 -0.0125 0.003125
0 -0.0125 0.003125
0.0125 -0.0125 0.003125
0.025 -0.0125 0.003125
0.0375 -0.0125 0.003125

# calibration marks at y = 0 m
-0.0375 0 0.003125
-0.025 0 0.003125
-0.0125 0 0.003125
0 0 0.003125
0.0125 0 0.003125
0.025 0 0.003125
0.0375 0 0.003125

# calibration marks at y = 0.0125 m
-0.0375 0.0125 0.003125
-0.025 0.0125 0.003125
-0.0125 0.0125 0.003125
0 0.0125 0.003125
0.0125 0.0125 0.003125
0.025 0.0125 0.003125
0.0375 0.0125 0.003125

# calibration marks at y = 0.025 m
-0.0375 0.025 0.003125
-0.025 0.025 0.003125
-0.0125 0.025 0.003125
0 0.025 0.003125
0.0125 0.025 0.003125
0.025 0.025 0.003125
0.0375 0.025 0.003125

# calibration marks at y = 0.0375 m
-0.0375 0.0375 0.003125
-0.025 0.0375 0.003125
-0.0125 0.0375 0.003125
0 0.0375 0.003125
0.0125 0.0375 0.003125
0.025 0.0375 0.003125
0.0375 0.0375 0.003125

caltab.ps：可直接用于打印的标定板图像文件。

打印出的标定板如下：

2，内部参数与外部参数获取
1，使用Halcon算子获取

相机内部参数标定算子所在位置：

算子如下：

*生成标定板文件
 *gen_caltab (7, 7, 0.0125, 0.5, 'caltab.descr', 'caltab.ps')
ImgPath := '3d_machine_vision/calib/'
dev_close_window ()
dev_open_window (0, 0, 652, 494, 'black', WindowHandle)
dev_update_off ()
dev_set_draw ('margin')
dev_set_line_width (3)
OpSystem := environment('OS')
set_display_font (WindowHandle, 14, 'mono', 'true', 'false')
* 
* Calibrate the camera.
* StartCamPar := [Focus,Kappa,Sx,Sy,Cx,Cy,ImageWidth,ImageHeight]
*Focus:相机的焦距，如果是远心相机，则焦距为0
*Kappa:畸变系数，初始值可以设置为0
*Sx、Sy:单个像元的宽，高（可从相机说明书中获取）
*Cx、Cy:图像的原点坐标，初始值可认为是图像的中心点，即坐标分为为图像的宽度和高度的一半
*ImageWidth,ImageHeight:采集图像的宽和高
StartCamPar := [0.016,0,0.0000074,0.0000074,326,247,652,494]
create_calib_data ('calibration_object', 1, 1, CalibDataID)
*设置为面阵模式
set_calib_data_cam_param (CalibDataID, 0, 'area_scan_division', StartCamPar)
*设置标定描述文件
set_calib_data_calib_object (CalibDataID, 0, 'caltab.descr')
NumImages := 10
* Note, we do not use the image from which the pose of the measurement plane can be derived
for I := 1 to NumImages by 1
    read_image (Image, ImgPath + 'calib_' + I$'02d')
    dev_display (Image)
    find_calib_object (Image, CalibDataID, 0, 0, I, [], [])
    get_calib_data_observ_contours (Caltab, CalibDataID, 'caltab', 0, 0, I)
    dev_set_color ('green')
    dev_display (Caltab)
endfor
*标定相机内参
calibrate_cameras (CalibDataID, Error)
*获取标定后的相机内部参数
get_calib_data (CalibDataID, 'camera', 0, 'params', CamParam)
* 将相机内部参数写入到文件
write_cam_par (CamParam, 'camera_parameters.dat')
Message := 'Interior camera parameters have'
Message[1] := 'been written to file'
disp_message (WindowHandle, Message, 'window', 12, 12, 'red', 'false')
clear_calib_data (CalibDataID)

相机外部参数标定算子所在位置：

算子如下：

* Attention:
* This program reads the internal camera parameters from the file
* 'camera_parameters.dat', which, e.g., could be generated by the program
* 'camera_calibration_internal.hdev'
* 
ImgPath := '3d_machine_vision/calib/'
dev_close_window ()
dev_open_window (0, 0, 652, 494, 'black', WindowHandle)
dev_update_off ()
dev_set_draw ('margin')
dev_set_line_width (1)
set_display_font (WindowHandle, 14, 'mono', 'true', 'false')
* 从文件中读取内部参数
try
    read_cam_par ('camera_parameters.dat', CamParam)
catch (Exception)
    * run 'camera_calibration_internal.hdev' first to generate camera
    * parameter file 'camera_parameters.dat'
    stop ()
endtry
* 
* Determine the external camera parameters and world coodinates from image points
* 
* The external camera parameters can be determined from an image, where the
* calibration plate is positioned directly on the measurement plane
read_image (Image, ImgPath + 'calib_11')
dev_display (Image)
*标定描述文件
CaltabName := 'caltab.descr'
create_calib_data ('calibration_object', 1, 1, CalibDataID)
* Here, the final camera parameters are already known and can be used instead
* of the starting values used in the program 'camera_calibration_internal.hdev'
set_calib_data_cam_param (CalibDataID, 0, 'area_scan_division', CamParam)
set_calib_data_calib_object (CalibDataID, 0, CaltabName)
find_calib_object (Image, CalibDataID, 0, 0, 1, [], [])
get_calib_data_observ_contours (Caltab, CalibDataID, 'caltab', 0, 0, 1)
*获取粗略估计的姿态（外部参数）
get_calib_data_observ_points (CalibDataID, 0, 0, 1, RCoord, CCoord, Index, PoseForCalibrationPlate)
dev_set_color ('green')
dev_display (Caltab)
dev_set_color ('red')
disp_caltab (WindowHandle, CaltabName, CamParam, PoseForCalibrationPlate, 1)
dev_set_line_width (3)
disp_circle (WindowHandle, RCoord, CCoord, gen_tuple_const(|RCoord|,1.5))
* caltab_points (CaltabName, X, Y, Z)
* calibrate_cameras (CalibDataID, Error)
* To take the thickness of the calibration plate into account, the z-value
* of the origin given by the camera pose has to be translated by the
* thickness of the calibration plate.
* Deactivate the following line if you do not want to add the correction.
*设置标定板的厚度，如果没有则省略
set_origin_pose (PoseForCalibrationPlate, 0, 0, 0.00075, PoseForCalibrationPlate)
disp_continue_message (WindowHandle, 'black', 'true')
stop ()
* Alternatively, the external camera parameters can  be determined from
* at least three point correspondances between the WCS and the pixel coordinate system
read_image (Image, ImgPath + 'caliper_01')
dev_display (Image)
* Set the world coordinates of three points on the rule
*在尺上选择点设置世界坐标
X := [0,50,100,80]
Y := [5,0,5,0]
Z := [0,0,0,0]
* Set the respective image plane coordinates of the three points
*设置选择点在图像上对应的像素坐标
RCoord := [414,227,85,128]
CCoord := [119,318,550,448]
* 
disp_cross (WindowHandle, RCoord, CCoord, 6, 0)
* create_pose (-50, 25, 400, 0, 0, -30, 'Rp+T', 'gba', 'point', InitialPose)
*生成外部参数（姿态参数）
vector_to_pose (X, Y, Z, RCoord, CCoord, CamParam, 'iterative', 'error', FinalPose, Errors)
*将外部参数写入到文件
write_pose (FinalPose, 'pose_from_three_points.dat')

* Now, transform a point measured interactively into the WCS
dev_update_window ('on')
dev_display (Image)
while (1)
    disp_message (WindowHandle, 'Measure one point: left mouse button', 'window', 12, 12, 'red', 'false')
    disp_message (WindowHandle, 'Exit measure mode: right mouse button', 'window', 36, 12, 'red', 'false')
    get_mbutton (WindowHandle, Row, Column, Button)
    if (Button == 4)
        break
    endif
    dev_display (Image)
    dev_set_color ('green')
    disp_cross (WindowHandle, Row, Column, 6, 0)
    image_points_to_world_plane (CamParam, FinalPose, Row, Column, 1, X1, Y1)
   
    disp_message (WindowHandle, 'X = ' + X1, 'window', 320, 400, 'red', 'false')
    disp_message (WindowHandle, 'Y = ' + Y1, 'window', 340, 400, 'red', 'false')
endwhile
2，使用助手获取相机内部参数与外部参数。

选择标定助手

设置参数

选择拍摄的各个角度标定板图像

不同角度的标定板图像

查看标定结果

6，应用。

* gen_caltab (7, 7, 0.0125, 0.5, 'caltab.descr', 'caltab.ps')
dev_close_window ()
read_image (Image, '3d_machine_vision/calib/calib_10')
get_image_size(Image, Width, Height)
dev_open_window (0, 0, Width, Height, 'black', WindowHandle)
dev_display (Image)
*测量两个点之间的世界距离
*因为标定板相邻的两个圆远心之间的距离是12.5mm，所以用这个进行测试
*读取相机内参
read_cam_par ('camera_parameters.dat', CameraParam)
*读取相机的外参
read_pose ('pose_from_three_points.dat', Pose)
*获取两个相邻圆的圆形图像距离
threshold (Image, Region, 0, 100)
connection (Region, ConnectedRegions)
select_shape (ConnectedRegions, SelectedRegions, ['area','row','column'], 'and', [0,305.96,360.09], [131651,345.41,500])
*获取该两点的相机坐标
gen_contour_region_xld (SelectedRegions, Contours, 'border')
dev_set_color ('green')
dev_set_line_width (2)

dev_display (Image)
dev_display (Contours)
*获取两个选定区域的坐标
area_center (SelectedRegions, Area, Row, Column)

*映射成世界坐标
image_points_to_world_plane (CameraParam, Pose, Row, Column, 'mm', X, Y)
*计算两点之间的距离
distance_pp (X[1], Y[1], X[0], Y[0], Distance)
dev_set_color ('red')
disp_line (WindowHandle, Row[0], Column[0], Row[1], Column[1])
disp_message (WindowHandle,'Distance:'+ Distance$'.2f'+'mm', 'window', 20, 20, 'green', 'false')

 

7，代码文件资源链接。

https://download.csdn.net/download/lingxiao16888/90384839

本文来自互联网用户投稿，该文观点仅代表作者本人，不代表本站立场。本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如若转载，请注明出处：/a/969698.html

如若内容造成侵权/违法违规/事实不符，请联系我们进行投诉反馈qq邮箱809451989@qq.com，一经查实，立即删除！

相关文章

学习星开源在线考试教育系统

学习星开源在线考试教育系统

学习星开源考试系统 项目介绍 项目概述： 学习星在线考试系统是一款基于Java和Vue.js构建的前后端分离的在线考试解决方案。它旨在为教育机构、企业和个人提供一个高效、便捷的在线测试平台，支持多种题型，包括但不限于单选题、多选题、判断…
阅读更多...
趣味魔法项目 LinuxPDF —— 在 PDF 中启动一个 Linux 操作系统

趣味魔法项目 LinuxPDF —— 在 PDF 中启动一个 Linux 操作系统

最近，一位开源爱好者开发了一个LinuxPDF 项目（ading2210/linuxpdf: Linux running inside a PDF file via a RISC-V emulator），它的核心功能是在一个 PDF 文件中启动并运行 Linux 操作系统。它通过巧妙地使用 PDF 文件格式中的 Ja…
阅读更多...
k8s的安装

k8s的安装

1. k8s的安装 192.168.48.6 master01 192.168.481.6 node01 192.168.48.26 node02 三台机器一起操作 1.swapoff -a ：关闭交换分区 2. iptables -F && iptables -t nat -F && iptables -t mangle -F && iptables -X 3. cat > /etc/sy…
阅读更多...
Sentinel 持久化配置

Sentinel 持久化配置

前言 在微服务架构中，Sentinel 是一个非常流行的流量控制和熔断组件，它可以帮助我们保护系统免受高流量的冲击。然而，Sentinel 的配置在默认情况下是不持久化的，这意味着一旦服务重启，所有配置就会丢失。为了解决这个…
阅读更多...
不到1M的工具，使用起来非常丝滑！

不到1M的工具，使用起来非常丝滑！

今天给大家推荐两款电脑上超实用的小软件，分别是倒计时工具和关机助手，用起来特别顺手，希望能帮到大家。 关机助手 帮你自动关机 先说说关机助手。这款软件简直太方便了！它是一款免安装的小工具，体积超小，…
阅读更多...
day9手机创意软件

day9手机创意软件

趣味类 in:记录趣味生活（通用） 魔漫相机：真人变漫画（通用） 活照片：让照片活过来（通用） 画中画相机：与众不同的艺术 年龄检测仪：比一比谁更年轻&#xf…
阅读更多...
【前端 DevOps】GitHub Actions 与 GitLab CI 实战：实现前端项目的自动化测试与部署

【前端 DevOps】GitHub Actions 与 GitLab CI 实战：实现前端项目的自动化测试与部署

网罗开发 （小红书、快手、视频号同名） 大家好，我是 展菲，目前在上市企业从事人工智能项目研发管理工作，平时热衷于分享各种编程领域的软硬技能知识以及前沿技术，包括iOS、前端、Harmony OS、Java、Python等…
阅读更多...
无耳科技 Solon v3.0.8 发布，Java 企业级应用开发框架

无耳科技 Solon v3.0.8 发布，Java 企业级应用开发框架

Solon 框架！ Solon 是新一代，Java 企业级应用开发框架。是杭州无耳科技有限公司的“根级”开源项目（最近“杭州六小龙”很火啊，我们也是杭州的哦）。从零开始构建（No Spring、No Java-EE、No Servlet&#…
阅读更多...
Flutter 异步编程利器：Future 与 Stream 深度解析

Flutter 异步编程利器：Future 与 Stream 深度解析

目录 一、Future：处理单次异步操作 1. 概念解读 2. 使用场景 3. 基本用法 3.1 创建 Future 3.2 使用 then 消费 Future 3.3 特性 二、Stream：处理连续异步事件流 1. 概念解读 2. 使用场景 3. 基本用法 3.1 创建 Stream 3.2 监听 Stream 3.…
阅读更多...
Agents Go Deep 智能体深入探索

Agents Go Deep 智能体深入探索

Agents Go Deep 智能体深入探索 核心事件 OpenAI发布了一款先进的智能体“深度研究”，它能借助网络搜索和推理生成研究报告。 最新进展 功能特性：该智能体依据数百个在线资源生成详细报告，目前仅支持文本输出，不过很快会增加对图…
阅读更多...
STM32单片机芯片与内部85 RS232 RS485 UART ISP下载硬件选择 电路设计 IO分配

STM32单片机芯片与内部85 RS232 RS485 UART ISP下载硬件选择 电路设计 IO分配

目录 一、UART 1、硬件选择 2、电路设计 3、IO分配 4、其他设计 二、RS232 1、硬件选择 2、电路设计 3、IO分配 4、其他设计 三、RS485 1、硬件选择 2、电路设计 3、IO分配 4、其他设计 四、ISP下载 一、UART 1、硬件选择 一般选择CH340完成STM32的IO电平与US…
阅读更多...
期权帮 | 场外个股期权可以做吗，风险高吗?

期权帮 | 场外个股期权可以做吗，风险高吗?

锦鲤三三每日分享期权知识，帮助期权新手及时有效地掌握即市趋势与新资讯！ 场外个股期权可以做吗，风险高吗? 场外个股期权，就是在正式的交易所之外进行交易的个股期权。 注：这里的“场外”指的是这类交易不在像沪深…
阅读更多...
【DeepSeek】deepseek可视化部署

【DeepSeek】deepseek可视化部署

目录 1 -> 前文 2 -> 部署可视化界面 1 -> 前文 【DeepSeek】DeepSeek概述 | 本地部署deepseek 通过前文可以将deepseek部署到本地使用，可是每次都需要winR输入cmd调出命令行进入到命令模式，输入命令ollama run deepseek-r1:latest。体验很…
阅读更多...
USART串口协议

USART串口协议

USART串口协议 文章目录 USART串口协议1. 通信接口2.串口通信2.1硬件电路2.2电平标准2.3串口参数及时序（软件部分） 3.USART串口外设3.1串口外设3.2USART框图3.3USART基本结构3.4数据帧 4.输入电路4.1起始位侦测4.2数据采样 5.波特率发生器6.相关函数介绍…
阅读更多...
2025 西湖论剑wp

2025 西湖论剑wp

web Rank-l 打开题目环境： 发现一个输入框，看一下他是用上面语言写的 发现是python，很容易想到ssti 密码随便输，发现没有回显 但是输入其他字符会报错 确定为ssti注入 开始构造payload， {{(lipsum|attr(‘global…
阅读更多...
twisted实现MMORPG 游戏数据库操作封装设计与实现

twisted实现MMORPG 游戏数据库操作封装设计与实现

在设计 MMORPG（大规模多人在线角色扮演游戏）时，数据库系统是游戏架构中至关重要的一部分。数据库不仅承担了游戏中各种数据（如玩家数据、物品数据、游戏世界状态等）的存储和管理任务，还必须高效地支持并发访…
阅读更多...
PyCharm 批量替换

PyCharm 批量替换

选择替换的内容 1. 打开全局替换窗口 有两种方式可以打开全局替换窗口： 快捷键方式： 在 Windows 或 Linux 系统下，按下 Ctrl Shift R。在 Mac 系统下，按下 Command Shift R。菜单操作方式：点击菜单栏中的 Edit&…
阅读更多...
LabVIEW用户界面设计原则

LabVIEW用户界面设计原则

在LabVIEW开发中，用户界面（UI）设计不仅仅是为了美观，它直接关系到用户的操作效率和体验。一个直观、简洁、易于使用的界面能够大大提升软件的可用性，尤其是在复杂的实验或工业应用中。设计良好的UI能够减少操作错误&am…
阅读更多...
网络防御高级-第8章及之前综合作业

网络防御高级-第8章及之前综合作业

标准版 接口ip配置 r2 [r2]interface GigabitEthernet 0/0/0 [r2-GigabitEthernet0/0/0]ip address 13.0.0.3 24 [r2-GigabitEthernet0/0/0]interface GigabitEthernet 0/0/1 [r2-GigabitEthernet0/0/1]ip address 100.1.1.254 24 [r2-GigabitEthernet0/0/1]interface Gigab…
阅读更多...
若依系统环境搭建记录

若依系统环境搭建记录

开源若依系统网上资料也很全的，本篇博文记录下自己搭建环境过程中遇到的一些问题。 配置Maven和编辑器选择 我懒得配置Eclipse了，直接用vscode作为编辑器，后面构建运行都用命令行。 配置数据库连接 按照mysql5.7按网上教程即可&#xff1…
阅读更多...
最新文章

Copyright @ 2022~2023