主权项 |
一种切向和径向双近似圆综合的全景环带图像校正展开方法,其特征在于,具体步骤如下:步骤1、将全景环带图像视作切向近似圆,利用最小二乘法拟合圆心法计算切向近似圆圆心坐标参数:所述步骤1中全景环带图像切向近似圆圆心坐标通过最小二乘拟合圆心法计算得出,切向近似圆圆心(x<sub>tan0</sub>,y<sub>tan0</sub>)的计算方法具体为:<img file="FDA0000767099860000011.GIF" wi="1479" he="613" />其中,A<sub>i</sub>=2(x<sub>i+1</sub>‑x<sub>i</sub>),B<sub>i</sub>=2(y<sub>i+1</sub>‑y<sub>i</sub>),<img file="FDA0000767099860000012.GIF" wi="526" he="85" />i=1,2,…,m‑1;(x<sub>i</sub>,y<sub>i</sub>),i=1,2,…,m为切向近似圆边缘若干灰度值相近点的坐标测量数据;步骤2、将全景环带图像视作径向近似圆,确定近似圆上任意两条半径所在直线方程,继而得到径向近似圆圆心坐标参数:所述步骤2中全景环带图像径向近似圆圆心坐标通过计算全景图上两条半径所在直线的交点得到,径向近似圆心(x<sub>rad0</sub>,y<sub>rad0</sub>)的计算方法具体为:<img file="FDA0000767099860000013.GIF" wi="1358" he="173" />其中,K<sub>1</sub>、K<sub>2</sub>分别为径向近似圆上两条半径所在直线方程斜率,D<sub>1</sub>、D<sub>2</sub>分别为径向近似圆上两条半径所在直线方程截距;步骤3、利用步骤1和步骤2所获得的切向和径向近似圆圆心坐标参数,以及展开图几何特性,推导建立全景环带图像与柱面全景展开图像之间的仿射变换关系:所述步骤3中全景环带图像与柱面全景展开图像之间的仿射变换关系由切向和径向双近似圆参数以及全景环带图像展开前后坐标的几何特性推导建立,其仿射变换关系推导过程为:切向近似圆环上的点(x,y)满足二次方程:(x‑x<sub>tan0</sub>)<sup>2</sup>+(y‑y<sub>tan0</sub>)<sup>2</sup>=r<sup>2</sup> (3) 径向近似圆环上的点坐标(x,y)满足线性方程:y=kx+d (4) 在公式(4)中,径向近似圆一、三象限上点所在线性方程斜率k<sub>1,3</sub>=(y<sub>r1</sub>‑y<sub>rad0</sub>)/(x<sub>r1</sub>‑x<sub>rad0</sub>),截距d<sub>1,3</sub>=y<sub>r1</sub>‑k<sub>1,3</sub>x<sub>r1</sub>;二、四象限点所在线性方程斜率k<sub>2,4</sub>=(y<sub>r2</sub>‑y<sub>rad0</sub>)/(x<sub>r2</sub>‑x<sub>rad0</sub>),截距d<sub>2,4</sub>=y<sub>r2</sub>‑k<sub>2,4</sub>x<sub>r2</sub>;其中(x<sub>r1,</sub>y<sub>r1</sub>)为径向近似圆第一象限点坐标,(x<sub>r2</sub>,y<sub>r2</sub>)为近似圆第二象限点坐标,利用径向圆心坐标参数可求得x<sub>r1</sub>=j*cos(angle)+x<sub>rad0</sub>,y<sub>r1</sub>=y<sub>rad0</sub>‑j*sin(angle),angle=i/R<sub>c</sub>,x<sub>r2</sub>=2x<sub>rad0</sub>‑x<sub>r1</sub>,y<sub>r2</sub>=y<sub>r1</sub>,其中i、j分别为柱面全景展开图像中的行列变量,R<sub>c</sub>为设定的展开基准圆半径,联合公式(3)(4)可得以径向圆心为原点,向右为X轴正方向,向上为Y轴正方向的全景环带图像各象限像素坐标对应关系为:第一象限像素坐标对应关系为<img file="FDA0000767099860000021.GIF" wi="1246" he="269" />第二象限像素坐标对应关系为<img file="FDA0000767099860000022.GIF" wi="1230" he="269" />第三象限像素坐标对应关系为<img file="FDA0000767099860000023.GIF" wi="1229" he="275" />第四象限像素坐标对应关系为<img file="FDA0000767099860000024.GIF" wi="1237" he="270" />其中,(x<sub>1</sub>,y<sub>1</sub>)、(x<sub>2</sub>,y<sub>2</sub>)、(x<sub>3</sub>,y<sub>3</sub>)、(x<sub>4</sub>,y<sub>4</sub>)分别为全景环带图像径向近似圆上四个象限坐标,a<sub>1,3</sub>、b<sub>1,3</sub>、c<sub>1,3</sub>和a<sub>2,4</sub>、b<sub>2,4</sub>、c<sub>2,4</sub>为线性方程y=kx+d代入二次方程(x‑x<sub>tan0</sub>)<sup>2</sup>+(y‑y<sub>tan0</sub>)<sup>2</sup>=r<sup>2</sup>中化简所得二次方程(1+k<sup>2</sup>)x<sup>2</sup>‑2[x<sub>tan0</sub>+k(y<sub>tan0</sub>‑d)]x+(d‑y<sub>tan0</sub>)<sup>2</sup>‑r<sup>2</sup>=0的系数;步骤4、根据仿射变换关系,利用采样点周围四个邻点的灰度值作双线性插值,得到待采样点灰度,完成图像像素复制,展开全景图;步骤5、以方形网格图纸作为标定图像,包裹于检测装置上,获取其全景图像,并通过步骤1~4将其展开,建立图像修正前后行坐标反向映射关系式,修正步骤4中展开图像帧纵横比,实现全景环带图像的几何校正展开。 |