一种适于通用光学系统的亚像素级非线性标定方法与流程
本发明属于测量技术领域，尤其涉及一种适于通用光学系统的亚像素级非线性标定方法。
背景技术：
每当使用光学成像系统自动观测或控制现实世界中的对象时，需要通过它们在传感器上的投影推测观察对象的世界位置。确定投影系统的准确模型至关重要。标定将尺度从三维空间桥接到二维空间，从获得的二维图像估计相机的变换矩阵。因此，在视觉测量和光学检测中广泛研究并应用了标定算法，如对象重建、计算机断层摄影、姿态估计和机器臂定位。此外，标定对于生物学、材料、图像处理、光子成像、物理测量、物体检测的应用不可或缺。
当前已提出多种系统标定算法，其中自标定方法仅需建立图像间的对应，较为方便，但其准确性不高，不宜应用于精度要求高的场合。传统的相机标定算法基于特定的标定物，经过图像处理算法，并结合单应性矩阵，求取相机的内外参数。该类标定算法相比自标定方法精度较高，但没有考虑到光学系统的畸变等像差，还有提升的空间。
在传统的标定算法中，tsai提出基于径向排列约束的标定算法，但其仅考虑到径向畸变，对切向畸变较大的系统误差较大；张正友提出基于平面棋盘标定板的标定算法，方法简单，且精度较高，但对平面棋盘标定板的精度有较高的要求，当平面棋盘标定板的特征点坐标存在5％的误差，系统一阶径向畸变系数的误差会超过20％。
技术实现要素：
本发明的目的在于克服现有技术中存在的缺点，提供一种适于通用光学系统的亚像素级非线性标定方法，使测量整体仪器的标定误差达到亚像素量级，提高图像成像品质，可应用于视觉测量和光学检测领域。
为了实现上述目的，本发明提出了一种适于通用光学系统的亚像素级非线性标定方法，该方法包括以下步骤：
步骤(1)、调制lcd板晶格亮暗并拍摄图像；
将标定所需的图像显示在lcd板上，精确到lcd板上每个晶格，并拍摄图像；
步骤(2)、图像处理获得像面晶格形心坐标并计算像面晶格形心所对应的物面晶格形心坐标；
对步骤(1)拍摄的图像进行二值化，再筛选完整的连通区域求取像面晶格形心坐标，并根据已知lcd板每个晶格的间隔，由晶格的位置得到准确的物面晶格形心坐标；
步骤(3)、由计算得到的物面晶格形心坐标和其对应的像面晶格形心坐标进行非线性函数建模；
综合考虑畸变的三种模型——径向畸变、离心畸变、薄棱镜畸变，根据选取的物面晶格形心坐标和其对应的像面坐标求得结合三种畸变模型的非线性函数模型，根据非线性函数模型由拍摄到的图像反投影求得校正像差后的图像；
步骤(4)、迭代优化非线性函数模型中的系数；
将校正像差后的图像作为步骤(2)中拍摄的图像，再重复步骤(2)-步骤(3)，直到非线性函数模型的误差小于0.5像素，优化得到非线性函数模型中的系数。
进一步的，所述步骤(2)中二值化处理采用otsu最大类间方差算法。
进一步的，所述根据已知lcd板每个晶格的间隔，由晶格的位置得到准确的物面晶格形心坐标的步骤具体为：
centery＝codey*c(2)
其中，centerx为物面晶格形心x轴坐标，centery为物面晶格形心y轴坐标，codex为晶格在x轴上的位置，codey为晶格在y轴上的位置，a为一个像素(包含红、绿、蓝三个晶格)的横向间隔，b为单个晶格的横向间隔(a的1/3)，c为一个像素(包含红、绿、蓝三个晶格)的纵向间隔。
进一步的，所述步骤(3)具体为：
未考虑成像畸变的理想针孔模型用齐次坐标表示为：
其中，zc为一个比例系数，[xy]t为物面晶格形心坐标，[uv]t为像面晶格形心坐标，r1、r2是从物面坐标系到像面坐标系的旋转变换系数，t是从物面坐标系到像面坐标系的平移变换矩阵，k是相机内参矩阵，h3*3为从物面坐标系到像面坐标系的单应性矩阵,h11、h12、h13、h21、h22、h23、h31、h32、h33为单应性矩阵中的各系数，[u0v0]t为像面的主点坐标，α、β分别是图像在u轴和v轴的焦距参数，γ为图像扭曲参数，表示图像两坐标轴偏斜度；
以下是考虑到的三种畸变模型：
(1)径向畸变
δxr(x,y)为x轴的径向畸变，δyr(x,y)为y轴的径向畸变，k1为径向畸变系数；
(2)离心畸变
δxd(x,y)为x轴的离心畸变，δyd(x,y)为y轴的离心畸变，p1,p2为离心畸变系数；
(3)薄棱镜畸变
δxl(x,y)为x轴的薄棱镜畸变，δyl(x,y)为y轴的薄棱镜畸变，s1,p2薄棱镜畸变系数；
建立考虑畸变的非线性函数模型为：
给定n幅拍摄到的图像，每幅图像上有m个采样晶格，函数评价模型为：
其中ε为非线性函数模型的误差，mij是第i幅图像上第j个晶格的物面形心坐标，ri是第i幅图像对应的旋转矩阵，mj是第j个晶格的像面形心坐标，
根据非线性函数模型由拍摄到的图像反投影求得校正像差后的图像。
相对于现有技术，本发明的有益效果如下：通过本发明的方法处理后的图像质量高，解决了光学系统内在像差的影响；引入额外的非线性参数，矫正了整体光学仪器的畸变固有光学像差，精度更高；由于lcd板晶格的微小尺寸，在视场内可观察到较多数量的晶格，得到足够多的样本数量对非线性函数进行建模，且像面坐标和对应的物面坐标精确到微米数量级，大大减少标定的采样误差；在非线性函数建模的基础上进行迭代优化，进一步校正残余像差，实现对标定误差的有效抑制，获得整体系统的亚像素级标定误差。
附图说明
图1为本发明实施例中使用的lcd板放大图；
图2为本发明实施例中lcd板图像的晶格形心示意图；
图3为畸变的物理模型图。
具体实施形式
下面通过实施例并结合附图对本发明做进一步的说明。
实施例：
本实施例的具体过程和步骤包括：
步骤(1)、调制lcd板晶格亮暗并拍摄图像；
将标定所需的图像显示在lcd板上，精确到lcd板上每个晶格；根据标定所需的图像编程控制lcd板晶格亮暗，使用光学系统对不同角度和位置的lcd板进行成像，拍摄几幅图像，保持图像清晰。
步骤(2)、图像处理获得像面晶格形心坐标并计算像面晶格形心所对应的物面晶格形心坐标；
对步骤(1)拍摄的图像进行二值化，二值化采用otsu最大类间方差算法，该算法按照图像的灰度特性，将图像分为背景和目标两部分，背景和目标之间的类间方差越大，则图像分割的两部分的差距越大，当部分目标错分为背景或部分背景错分为目标将导致两部分的差距变小；再选取完整的连通区域求取像面晶格形心坐标，对二值化图像的每个晶格进行轮廓提取，晶格理论上为排列整齐的规律的小方块，并不由图像的灰度信息分布决定，采用形心法而非质心法筛选连通区域求取像面晶格形心坐标，如图2所示，每个晶格中间的黑点即为求得的图像lcd板晶格形心；在选取像面晶格形心坐标时，可根据实际要求灵活选择，如对于拼接边界部分的误差要求更高，则较多选取视场范围边界部分的像面晶格形心坐标，如对于视场中心的内容更为关注，则较多选取视场范围中间部分的像面晶格形心坐标；
根据已知lcd板每个晶格的间隔，由晶格的位置得到准确的物面晶格形心坐标：
centery＝codey*c(10)
其中，centerx为物面晶格形心x轴坐标，centery为物面晶格形心y轴坐标，codex为晶格在x轴上的位置，codey为晶格在y轴上的位置，a为一个像素(包含红、绿、蓝三个晶格)的横向间隔，b为单个晶格的横向间隔(a的1/3)，c为一个像素(包含红、绿、蓝三个晶格)的纵向间隔。
步骤(3)、由计算得到的物面晶格形心坐标和其对应的像面晶格形心坐标进行非线性函数建模；
未考虑成像畸变的理想针孔模型用齐次坐标表示为：
其中，zc为一个比例系数，[xy]t为物面晶格形心坐标，[uv]t为像面晶格形心坐标，r1、r2是从物面坐标系到像面坐标系的旋转变换系数，t是从物面坐标系到像面坐标系的平移变换矩阵，k是相机内参矩阵，h3*3为从物面坐标系到像面坐标系的单应性矩阵,h11、h12、h13、h21、h22、h23、h31、h32、h33为单应性矩阵中的各系数，[u0v0]t为像面的主点坐标，α、β分别是图像在u轴和v轴的焦距参数，γ为图像扭曲参数，表示图像两坐标轴偏斜度；
实际成像过程中光学系统存在一定的畸变，理想针孔模型不能准确描述成像几何关系，为使建立模型更精确描述系统的成像过程，引入反映畸变影响的修正系数，建立精确的非线性函数。以下是考虑到的三种畸变模型：
(1)径向畸变
δxr(x,y)为x轴的径向畸变，δyr(x,y)为y轴的径向畸变，k1为径向畸变系数；
(2)离心畸变
δxd(x,y)为x轴的离心畸变，δyd(x,y)为y轴的离心畸变，p1,p2为离心畸变系数；
(3)薄棱镜畸变
δxl(x,y)为x轴的薄棱镜畸变，δyl(x,y)为y轴的薄棱镜畸变，s1,p2薄棱镜畸变系数；
建立考虑畸变的非线性函数模型为：
给定n幅拍摄到的图像，每幅图像上有m个采样晶格，函数评价模型为：
其中ε为非线性函数模型的误差，mij是第i幅图像上第j个晶格的物面形心坐标，ri是第i幅图像对应的旋转矩阵，mj是第j个晶格的像面形心坐标，
根据非线性函数模型由拍摄到的图像反投影求得校正像差后的图像。
步骤(4)、迭代优化非线性函数模型中的系数；
将校正像差后的图像作为步骤(2)中拍摄的图像，再重复步骤(2)-步骤(3)，直到非线性函数模型的误差小于0.5像素，优化得到非线性函数模型中的系数。
表1线性函数模型的测试误差(单位：像素)
表2非线性函数模型的测试误差(单位：像素)
在保证表1和表2测试样本相同的情况下，表1为应用传统方法建立线性函数模型的误差，表2为应用本发明建立迭代优化的非线性函数模型的误差。由表1和表2的数据得：传统非线性函数模型的平均误差在2.08像素，迭代优化的非线性函数模型的平均误差在0.2像素，与传统线性函数模型相比，本发明误差提高了一个数量级，使测量整体仪器的标定误差达到亚像素量级，可应用于高精度视觉测量和光学检测领域。
- 还没有人留言评论。精彩留言会获得点赞！