发明名称 切向和径向双近似圆综合的全景环带图像校正展开方法
摘要 本发明公开了一种切向和径向双近似圆综合的全景环带图像校正展开方法,本方法首先将全景环带图像分别视作切向与径向近似圆,利用最小二乘拟合圆心法计算切向近似圆圆心坐标参数,根据图像特征确定径向近似圆上两条半径所在的直线方程计算径向近似圆心坐标参数;再根据所得到的双近似圆圆心坐标参数及展开图几何特性,推导建立全景环带图像与柱面全景展开图像之间的仿射变换关系,并利用双线性插值法完成图像像素复制,展开全景图;最后将展开后图像的帧纵横比通过标定方法建立的图像映射关系式进行修正,从而获得高精度的全景环带几何校正展开图像。本发明使图像展开结果更真实,同时也提高了后期图像处理检测精度。<sup />
申请公布号 CN103208120B 申请公布日期 2016.01.20
申请号 CN201310112319.9 申请日期 2013.04.01
申请人 南京理工大学 发明人 卜雄洙;王新征;杨波;李桂娟;于靖
分类号 G06T7/00(2006.01)I 主分类号 G06T7/00(2006.01)I
代理机构 南京理工大学专利中心 32203 代理人 朱显国
主权项 一种切向和径向双近似圆综合的全景环带图像校正展开方法,其特征在于,具体步骤如下:步骤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中展开图像帧纵横比,实现全景环带图像的几何校正展开。
地址 210094 江苏省南京市玄武区孝陵卫200号