XRay几何标定论文(5)阅读:Automated Flaw Detection in Aluminum Castings-编程实验室

是的，这篇论文详细介绍了X射线成像系统的几何校正（标定）方法，主要在第II章几何模型和附录A中详细说明。以下是校正方法的总结：

一、为什么要进行几何校正？

X射线成像系统存在多种畸变：

图像增强器输入屏的曲率（双曲面形状）
电磁场畸变
透视投影参数（焦距、位置等）

这些畸变会导致投影坐标与实际图像坐标不一致，影响多视图几何约束的准确性。

二、校正的核心模型

1.坐标系定义（图3）

论文定义了4个坐标系：

世界坐标系{X, Y, Z}：原点在X射线源
物体坐标系{X, Y, Z}：附着在被测物体上
X射线投影坐标系{x, y}：虚构平面（无畸变）
图像坐标系{u, v}：CCD相机实际图像

2.投影模型

理想投影（公式4）：
λpmp=PpM \lambda_p \mathbf{m}_p = \mathbf{P}_p \mathbf{M}λpmp=PpM
其中：
- M\mathbf{M}M：物体3D点
- mp\mathbf{m}_pmp：投影平面2D点
- Pp\mathbf{P}_pPp：透视投影矩阵（含7个参数）
- λp\lambda_pλp：比例因子
实际畸变模型（公式5, 14）：
wp=f(mp)=Ag(mp) \mathbf{w}_p = \mathbf{f}(\mathbf{m}_p) = \mathbf{A} \mathbf{g}(\mathbf{m}_p)wp=f(mp)=Ag(mp)
其中g\mathbf{g}g是非线性畸变函数，A\mathbf{A}A是仿射变换矩阵。

3.畸变具体形式

图像增强器畸变（双曲面模型）

输入屏是双曲面：
Z=F(X^,Y^)=f1+(X^a)2+(Y^b)2 \mathcal{Z} = F(\hat{X},\hat{Y}) = f\sqrt{1 + \left(\frac{\hat{X}}{a}\right)^2 + \left(\frac{\hat{Y}}{b}\right)^2}Z=F(X^,Y^)=f1+(aX^)2+(bY^)2

(f)：双曲面实半轴（等于焦距）
(a, b)：虚半轴

投影点(\mathbf{m}'_p)是射线与双曲面的交点：
mp′=g(mp)=[cxp,cyp,1]T \mathbf{m}'_p = \mathbf{g}(\mathbf{m}_p) = [cx_p, cy_p, 1]^Tmp′=g(mp)=[cxp,cyp,1]T
c=1/1−(xp/a)2−(yp/b)2 c = 1 / \sqrt{1 - (x_p/a)^2 - (y_p/b)^2}c=1/1−(xp/a)2−(yp/b)2

CCD成像仿射变换

wp=[+kxcos⁡θ+kysin⁡θu0−kxsin⁡θ+kycos⁡θv0001]mp′ \mathbf{w}_p = \begin{bmatrix} +k_x\cos\theta & +k_y\sin\theta & u_0 \\ -k_x\sin\theta & +k_y\cos\theta & v_0 \\ 0 & 0 & 1 \end{bmatrix} \mathbf{m}'_pwp=+kxcosθ−kxsinθ0+kysinθ+kycosθ0u0v01mp′

kx,kyk_x, k_ykx,ky：缩放因子
(u0,v0)(u_0, v_0)(u0,v0)：平移
θ\thetaθ：旋转角

三、标定过程

1.标定板

使用带规则网格孔的标定板（图4）
已知每个孔的三维坐标(Xi,Yi,Zi)(X_i, Y_i, Z_i)(Xi,Yi,Zi)

2.数据采集

在(N)个不同位置（绕Z轴旋转）拍摄标定板
测量每个孔在图像中的坐标(uip,vip)(u_{ip}, v_{ip})(uip,vip)

3.参数估计

待估参数（共13个）：

投影参数：tX,tY,tZ,ωX,ωY,ft_X, t_Y, t_Z, \omega_X, \omega_Y, ftX,tY,tZ,ωX,ωY,f
畸变参数：a,ba, ba,b
仿射参数：kx,ky,θ,u0,v0k_x, k_y, \theta, u_0, v_0kx,ky,θ,u0,v0

优化目标：最小化重投影误差
min⁡∑i,p∥(uip,vip)−(u^ip,v^ip)∥ \min \sum_{i,p} \| (u_{ip}, v_{ip}) - (\hat{u}_{ip}, \hat{v}_{ip}) \|mini,p∑∥(uip,vip)−(u^ip,v^ip)∥
其中(u^ip,v^ip)(\hat{u}_{ip}, \hat{v}_{ip})(u^ip,v^ip)是通过模型从(Xi,Yi,Zi)(X_i, Y_i, Z_i)(Xi,Yi,Zi)计算的理论投影。