一种定焦、变焦镜头通光量强度不均校正方法及装置的制造方法
【技术领域】
[0001] 本发明涉及原始图像信号处理领域，尤其是涉及一种定焦、变焦镜头通光量强度 不均校正方法及装置。
【背景技术】
[0002] 镜头，尤其是长焦镜头，由于其光学物理特性和镜片组的设计局限可能会出现成 像中间区域亮边缘区域暗的亮度不均匀问题。特别是在户外环境下使用的长焦镜头，无法 保证在理想环境下工作。受温度影响镜头热胀冷缩，通光量强度不均匀的程度也会随之变 化。
[0003] 现有一些技术在相机出厂前对镜头进行校正，并将校正数据存储在相机中。但是 这要求相机与镜头固定配对。一旦配对关系改变，相机记录的校正数据将不再适用新镜头， 需要专业的设备和方法对镜头进行重新校正。除此以外，即便镜头在出厂前进行过亮度不 均匀校正，也难以保证镜头在所有工作环境下（尤其是温差大的环境下）不会出现亮度不 均匀的问题。
【发明内容】
[0004] 本发明所要解决的技术问题是：针对定焦镜头自身缺陷或者由于镜头形变引起的 成像亮度不均匀问题，本发明提出一种自动对通光量强度区域分布不均匀问题进行校正的 方法及装置。其能适应定焦镜头，进一步的在两次自动校正方法运行完成后，可以手动设置 镜头形变系数G，快速校正定焦镜头通光量强度区域分布不均匀问题，进一步的，通过选取 焦距样本点计算变焦镜头的平滑曲面参数从而到达变焦镜头通光量强度分布不均匀问题 进行校正。
[0005] 本发明采用的技术方案如下：
[0006] -种定焦镜头通光量强度不均校正的方法，其特征在于包括：
[0007] 步骤1 :采集大量图像生成平均图像IA;
[0008] 步骤2 :对平均图像进行平滑曲面最小误差拟合，计算平滑曲面参数C;
[0009] 步骤3 :通过平滑曲面参数c生成校正模板Ιτ;
[0010] 步骤4 :在计算出校正模板、后，每次拍照生成图像I，对I执行
[0011] Γ = I. /Ιτ
[0012] 其中./为点除运算，即让对应行数和列数的像素做除法；然后将Γ像素数值的分 布归一化处理至〇到255区间，对图像亮度不均的像素进行校正。
[0013] IR=sJ' +s2
[0014] 其中
[0015] 进一步的，所述步骤1具体包括：
[0016] 步骤11 :以视频或者连续快速拍照的方式采集图像，图像数据应尽量保持原始 性，不做数据处理；同时保持相机转动使拍摄的照片不重复；在采集100到1000张图像后， 将所有图像转化为灰度图像，然后计算平均图像IA;
[0017] 步骤12 :将^用Sobel边缘检测算子计算边缘强度，即分别使月
和
:寸IA做卷积运算，如果生成的结果中所有像素的值均小于Ti，则完成平均图 像生成，否则执行步骤13 ;1\为20到50之间；
[0018] 步骤13 :再采集100到500张图像后，重新计算包含原有图像和新采集图像的所 有图像的灰度转化平均图像IA，然后重新执行步骤12。
[0019] 进一步的，所述步骤2具体包括：
[0020] 步骤21 :对平均图像^的所有行列按照5到100等像素间隔进行亚采样，样本记 为匕^^⑴……^…^^…^丄其中"^为像素坐标^为像素对应的像素灰度 值，角标为采样点序号，&为亚采样样本点个数；然后选择步骤22或者步骤23对亚采样的 平均图像进行最小误差拟合，计算计算平滑曲面参数C;
[0021] 步骤22 :根据二阶多项式，构造矩阵X和V;
[0022]
[0023] 计算C= (XTX)￥V，其中矩阵右上角的T表示转置，右上角的-1表示求逆矩阵；
[0024] 步骤23 :根据加权二阶多项式，构造矩阵X和V
[0025]
[0026] 其中w是对每一个采样点设定的权重值，角标为采样点序号，一般图片中心样本 的的权重值较大，图像边缘样本点的权重值偏小；计算C= (XTX)YV，其中矩阵右上角的T 表示转置，右上角的-1表示求逆矩阵；
[0027] 进一步的，所述步骤2具体包括：
[0028] 步骤21 :对平均图像^的所有行列按照5到100等像素间隔进行亚采样，样本记 为{xdyi,v1;x2,y2,v2;. . .;xN1,yN1,vN1}，其中（X，y)为像素坐标，v为像素灰度值，角标为采 样点序号，&为亚采样样本点个数；然后选择步骤22或者步骤23对亚采样的平均图像进 行最小误差拟合，计算计算平滑曲面参数C;
[0029] 步骤22 :根据二维高斯，让〇7的取值从0到IA中的最小值min(IA)，对每一个取值 构造X和V
[0030]
[0031] 计算C' = (XTX)YV，其中矩阵右上角的T表示转置，右上角的-1表示求逆矩阵， 然后计算err= (XC' -V)T(XC' -V);在〇7的所有取值中选取err最小的一组C'和c7共同构 造C，
[0032]
[0033] 步骤23 :根据加权二维高斯，让c7的取值从0到IA中的最小值min(IA)，对每一 个取值构造X和V
[0034]
[0035]
[0036] 其中w是对每一个采样点设定的权重值，角标为采样点序号，一般图片中心样本 的权重值较大，图像边缘样本点的权重值偏小；计算C' = (XTX)YV，其中矩阵右上角的T 表示转置，右上角的-1表示求逆矩阵，然后计算err= (XC' -V)T (XC' -V);在c7的所有取值 中选取err最小的一组C'和c7共同构造C;
[0037]
[0038] 进一步的，所述步骤3具体包括
[0039] 步骤31 :在二阶多项式曲面及加权二阶多项式曲面时，对、中的每一个像素，令 其行数为X，列数为y，则该像素的值z为：z=cy+CjjXy+cY+qx+c^y+Cj；
[0040] 步骤32:在二维高斯曲面及加权二维高斯曲面时，对、中的每一个像素，令其行 数为X，列数为y，则该像素的值Z为：z=exp(cy+CjjXy+c^^+qx+Csy+Ce)+c7。
[0041] 进一步的，在执行N2次所述步骤1和步骤2得到多组曲面参数之后，不再执行步 骤1及步骤2,在此基础上使用步骤24能够快速得到近似的曲面参数适应镜头形变，具体包 括：
[0042] 步骤241:当得到曲面参数为Cj={cu，c2」，c3」，c4」，c5」，c6」}的二阶多项式曲面 参数或者加权二阶多项式曲面参数，则执行步骤242 ;当得到参数为Cj={ci」，c2」，c3」，c4_ 彡c5」，c6」，c7」}的二维高斯曲面参数或者加权二维高斯曲面参数，则执行步骤245 ;其中 表示第j次自动校正操作结果，1 <j<N2;
[0043] 步骤242 :对所有的Cl」，1彡j彡N2，提取最近两次自动校正算法结果c1N2i和c N2，令ki=cu-Cu；；i，执行步骤 243 ;
[0044] 步骤243:对c2」，c3」，c4」，c5」，c6」，重复步骤242,计算k2,k3,k4,k5,k6，执行步骤 244；
[0045] 步骤244:设置镜头形变参数G，则新的二阶多项式曲面参数Cnew= {clnew，c2new，c3_ new?。4-new，。5-new，。6-new^
[0046] clnew= c^2+kfi
[0047] c2-new= c2-N2+k2G
[0048] c3 new= c3N2+k3G
[0049] c4 new= c4N2+k4G
[0050] c5-new= c5-N2+k5G
[0051] c6-new= c6-N2+k6G;
[0052] 其中G的范围是任意正数；
[0053] 步骤245 :对所有的Cl」，1彡j彡N2，提取最近两次自动校正算法结果c1N2i和c N2，令ki=cu-Cu；；i，执行步骤 246 ;
[0054] 步骤246:对c2」，c3」，c4」，c5」，c6」，c7」，重复步骤245,计算k2,k3,k4