预置位偏差计算方法
【专利摘要】本发明公开了一种预置位偏差计算方法，应用于监控系统中，用于测试电动云台、球机的预置位精度。本发明采用叠加标志点坐标及数学建模的方式，极大程度减小测试系统误差，从而大大提高预置位测试精度，以满足日益迫切的高精度预置位精度测试需求。
【专利说明】预置位偏差计算方法
【技术领域】
[0001]本发明涉及监控【技术领域】，尤其涉及一种预置位偏差计算方法。
【背景技术】
[0002]目前已有的预置位精度自动测试方法，主要描述了预置位精度测试原理及测试系统的组成，其测试原理是在被测设备上配置预置位映射点装置(如激光灯)，控制被测设备转动使标志点映射在数据采装置上，通过数据采集装置获取各个预置位映射点在画面中的像素坐标，将测试中所检测的预置位映射点像素坐标与初始映射点像素坐标对比，便能得到预置位偏差像素数，再以距整个画面每个像素所等效的实际长度，便能计算出预置位水平方向及垂直方向的偏差距离。
[0003]但是，上述测试方法中，无论数据采集装置中使用任何摄像装置，均会因光学成像而导致检测设备所获取的图像数据存在一定的畸变，而预置位偏差距离均是在畸变后的画面数据基础上计算的，故上述测试方法所计算的预置位精度存在较大的误差。为提高测试系统的测试精度，则要从提高预置位偏差距离的计算精度入手。
[0004]鉴于此，迫切需要提供一种预置位偏差计算方法及系统。
【发明内容】
[0005]本发明的目的是提供一种预置位偏差计算方法及系统，以解决上述现有技术中存在的预置位精度计算存在较大误差，从而影响预置位精度准确性的问题。
[0006]为实现本发明的目的，本发明提供了一种预置位偏差计算方法，应用于监控系统中，包括下述步骤:
[0007]a.利用标志点坐标纸，将标志点叠加到畸变画面上，并对所述畸变画面进行区域划分，所述标志点坐标纸与数据采集屏大小一致，其上相邻两标志点间的距离均相等；
[0008]b.获取所述畸变画面中标志点的坐标值(Xaj), Y(i，」))，其中，其中1、j分别表示所述标志点所在的行数和列数；
[0009]c.获取预置位映射点N的坐标，确定所述N所在的区域，获取所述N所在的区域的相邻四个标志点A、B、C、D，并以其建立单元四边形模型。
[0010]d.在所述单元四边形模型中，取线段AD、Be、AB、⑶的中点为L、K、P、Q，连接线段LK、PQ将四边形AB⑶分成左上、左下、右上、右下四个子区域，以所述预置位映射点N为基点作相应的四边形边的平行线与线段AD、BC、AB、⑶分别交于E、F、G、H，构成预置位映射点定位模型；
[0011]e.基于所述单元四边形模型和预置位映射点定位模型，计算线段EN及GN所代表的实际长度SeN和SEN，根据所述SeN和Sen计算所述N距离中心参考点的水平距离Sx和垂直距离Sy ；
[0012]f.基于步骤a-e中对所述N距离中心参考点的水平距离和垂直距离的计算方法，计算初始预置位映射点距中心参考点的实际水平距离Sxtl及垂直距离Sytl，从而计算出预置位偏差。
[0013]本发明的有益效果，与现有技术相比，以使用外界标志点坐标纸的方式，将预置位数据采集摄像机的图像畸变细化，并以微分的思想将整体图像畸变分割为多个小的单元四边形，以便于在尽可能减小畸变误差的前提下建立计算模型，反算出预置位标志点映射的实际距离。使用此方法能够极大减小测试误差，从而提高测试系统的测试精度，满足高精度预置位测试需求。
【专利附图】
【附图说明】
[0014]图1是本发明提供的标志点坐标纸示意图；
[0015]图2是本发明提供的畸变画面标志点示意图；
[0016]图3是本发明提供的单元四边形建模示意图；
[0017]图4是本发明提供的预置位映射点定位建模示意图；
[0018]图5是本发明的方法流程图。
【具体实施方式】
[0019]为了使本发明的目的、技术方案及有益效果更加清楚明白，以下结合实施例，对本发明进行进一步详细说明。应当理解为此处所描述的具体实施例仅仅用以解释本发明，并不用于限制本发明的保护范围。
[0020]监控系统中电动云台、球机的预置位精度一直是监控设备的重要性能指标。本发明提供了一种应用于上述电动云台、球机等监控系统中的高精度计算预置位偏差的方法。本发明采用叠加标志点坐标及数学建模的方式，极大程度减小测试系统误差，从而大大提高预置位测试精度，以满足日益迫切的高精度预置位精度测试需求。
[0021]如图5所示，本发明包括如下步骤:
[0022]步骤S501，通过叠加标志点的方式来对畸变图像分解为多个单元四边形，将如图1所示的标志点坐标纸与数据采集屏完全重合，通过屏幕摄像机获取此时标志点坐标画面，如图2所示，即为畸变后的标志点坐标，以此种方式便能将均匀等间距的标志点叠加到畸变画面上。
[0023]需要说明的是，该标志点坐标纸上的各标志点间的间距均一致，且从理论上看该间距越小计算精度越高，但此时计算工作量很大，而间距较大时计算误差比较大，达不到高精度计算的要求，故在制作标志点坐标纸时应选择适当的标志点间距。
[0024]步骤S502，再利用MATLAB读取畸变后的标志点坐标图中的各标志点的像素坐标(Xaj)jYajj)),其中X、Y的下标1、j分别表示该标志点所在的行数和列数。
[0025]步骤S503，以每四个相邻的标志点建立单元四边形模型，为下一步建立预置位映射点定位模型作准备。从第I行第I列的标志点开始，每相邻的四个标志点连线便构成一个单元四边形，设四边形位AB⑶,其各点坐标通式可表示为A(X(i，j), Y(i，j))、B(X(i+1，j), Y(i+1，j))、C (X(i+1, j+l)，Y(i+1, j+1))、D(X(i，j+1)，Y(i，j+1) )，并将此四边形命名为第i行第j列四边形；
[0026]步骤S504，线段 AD、BC、AB、CD 的中点分别为 L((X(i，j)+X(i，j+1))/2，(Ya,j)+Ya,J+i))/2)> K((Xa+1,J)+Xa+1,J+1))/2, (Υα+1, ^+Υα+1, J+1))/2), P ((Χα,^+Χα+1, /2,(Y(i, j) +Y(i+i, j)) /2) > Q ((X(i+1, J+1) +X(i, J+1)) /2, (Y(i+1，j+1)+Y(i，j+1))/2)。连接 LK、PQ 将单元四边形AB⑶划分成四个区域，分别称为左上、左下、右上、右下区域，如图3所示。
[0027]设检测到预置位映射点N的像素坐标为(X，Y)，线段AB、BC、⑶、AD所在的直线分别为a、b、C、d，通过点与直线的关系判断出包含点N的单元四边形及在该单元四边形中所在的区域。现假设N点包含于第i行第j列四边形，且位于该单元四边形的左上区域，则以此为基点作点N相邻的两四边形边的平行线，即过N点作直线a、d的平行形e、f,分别交直线a、b、c、d于点G、F、H、E，具体图像如图4所示，以此建立预置位映射点定位模型。同理，当点N位于该单元四边形的左下区域时，则是过点N作直线a、b的平行线；当点N位于该单元四边形的右上区域时，则是过点N作直线c、d的平行线；当点N位于该单元四边形的右下区域时，则是过点N作直线b、c的平行线。O
[0028]步骤S505，基于以上所建立的数学计算模式，便能计算出线段GN、EN所代表的实际距离SeN、Sen,具体计算步骤如下:
[0029]设直线a、b、C、d、e、f 的斜率分别为 Ka、Kb、Kc、Kd、Ke、Kf,则:
[0030]Ka= (Y(i; j) -Y(i+1, j)) / (X(i, j) (i+1, j))；
[0031]Kb — (Y(i+i，j+i)_Y(i+i，j)) / (X(i+i，j+i)_X(i+i，j));
[0032]Kc — (Y(i，j+i)_Y(i+i，j+i)) / (X(i，j+i)_X(i+i，j+i));
[0033]Kd= (Y(i; J+i) -Y(i; j)) / (X(i; J+i) (i, j));
[0034]由于直线a//e、d//f则有:
[0035]Ke = Ka ；
[0036]Kf = Kd ；
[0037]直线3、13、(3、(1、6、€的方程为:
[0038]Υ3 = Κ3(Χ3-Χ(^))+Υ(υ) ；......(I)
[0039]；......(2)
[0040]Yc = Kc(Xc-X(i+1,J+1))+Y(i+1,J+1) ；......(3)
[0041]Yd = Kd(Xd-X(i；J+1))+Y(i；J+1) ；......(4)
[0042]Ye = Ke (Xe-X) +Y ；......(5)
[0043]Yf = Kf (Xf-X) +Y ；......(6)
[0044]联立直线方程(4)与(5)、(2)与(5)、(I)与(6)、(3)与(6)，分别解出交点E、F的垂直像素坐标及交点G、H的水平像素坐标:
[0045]Ye = [KdYe (X(i,J+1)-X) -KeY(i, J+1)+KdY] / (Kd-Ke)；
[0046]Yf = [KbKe (X(i+1,j)-X) -KeY(i+1,^+KbY] / (Kb-Ke)；
[0047]Xg = [(YfX-Y)-(KaXaj)-Y(U)V(Kf-Ka)；
[0048]Xh = [ (YfX-Y) - (KcX(i+1, J+1)-Y(i+1, J+1)) ] / (Kf-Kc)；
[0049]设相邻标志点的实际间距为h，根据相似三角形原理，线段GN、EN的实际距离可表示为:
[0050]Sgn = I h (Xg-X)/(Xg-Xh) |
[0051]Sen = I h (Ye-Y)/(Ye-Yf) |
[0052]选定标志点坐标纸中心区域的标志点M(X(i(UCI),Y(i(UCI))为中心参考点，即以M为参考来计算预置位映射点距离点M的水平距离和垂直距离，则包含在第i行第j列单元四边形中的点N距离中心参考点的水平距离Sx和垂直距离Sy为:
[0053]Sx = I (j-j0)h+SGN I
[0054]Sy = I (1-10)h+SEN I
[0055]步骤S506，基于以上计算预置位映射点距离中心参考点的水平距离及垂直距离的方法，可计算出初始预置位映射点距中心参考点的实际水平距离Sxtl及垂直距离sy(l。则预置位水平方向偏差距离Λ Sx和垂直偏差距离Λ Sy分别为:
[0056]Λ Sx =丨 Sx-Sx0 I
[0057]Δ Sy = I Sy-Sy0 I
[0058]上述算法可计算畸变画面中各个区域的预置位实际偏差距离。
[0059]其中需要说明的是，此算法在进行单元四边形建模及预置位映射点定位建模中，分别将畸变后两标志点间的弧线近似为直线、将过N点所做的平行线EF和GH的实际长度均近似等于标志点间的间距h，则计算结果存在一定的误差。但经误差分析与计算，由于选取的标志点间距h远小于整个画面的实际大小，单元四边形建模及预置位映射点定位建模所引起的误差值远小与实际偏差距离值，故因建模所产生的误差值可忽略不计。另外，应用此方法计算预置位偏差实际距离，可消除N点所在的单元四边形与中心参考点间间隔的各个单元四边形的图像畸变误差。
[0060]综上所述，本发明通过叠加标志点的方式将畸变图像分解成多个单元四边形，并结合正确的数学建模方法减少计算误差，从而极大的减少了因图像畸变所引起的实际距离计算误差，很大程度提高了预置位实际偏差距离的计算准确度。该计算方法可应用在各类具备预置位功能设备的测试系统中，尤其是高精度预置位测试。
[0061]以上所述仅是本发明的优选实施方式，应当指出，对于本【技术领域】的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。
【权利要求】
1.一种预置位偏差计算方法，应用于监控系统中，其特征在于，包括下述步骤: a.利用标志点坐标纸，将标志点叠加到畸变画面上，并对所述畸变画面进行区域划分，所述标志点坐标纸与数据采集屏大小一致，其上相邻两标志点间的距离均相等； b.获取所述畸变画面中标志点的坐标值(X(ij，Y(ij)，其中，其中1、j分别表示所述标志点所在的行数和列数； c.获取预置位映射点N的坐标，确定所述N所在的区域，获取所述N所在的区域的相邻四个标志点A、B、C、D，并以其建立单元四边形模型； d.在所述单元四边形模型中，取线段AD、BC、AB、⑶的中点为1^、1(、？、0，连接线段1^、PQ将四边形AB⑶分成左上、左下、右上、右下四个子区域，以所述预置位映射点N为基点作相应的四边形边的平行线与线段AD、BC、AB、⑶分别交于E、F、G、H，构成预置位映射点定位模型； e.基于所述单元四边形模型和预置位映射点定位模型，计算线段EN及GN所代表的实际长度SeN和SEN，根据所述SeN和Sen计算所述N距离中心参考点的水平距离Sx和垂直距离Sy ； f.基于步骤a-e中对所述N距离中心参考点的水平距离和垂直距离的计算方法，计算初始预置位映射点距中心参考点的实际水平距离Sxtl及垂直距离Sytl，从而计算出预置位偏差。
2.根据权利要求1所述的预置位偏差计算方法，其特征在于，在步骤e中，计算所述SeN和Sen步骤包括: 设直线a、b、c、d、e、f分别对应线段AB、BC、CD、AD、EF、GH,其斜率分别为Ka、Kb、K。、Kd、Ke、Kf,则:
Ka= (Y(i, j) -Y(i+i, j)) / (X(i, j) -X(i+i, j))；
Kb 一 (Y(i+1, j+l)_Y(i+l, j)) / G(i+1, j+l)~^(i+l, j));
Kc (Y(i，j+1) Y(i+1，j+1)) / (X(i，j+1) X(i+l，j+l));
Kd = (y(i, j+D _y(i, j)) / (χ(?, j+D _x(i, j)); 由于直线a//e、d//f则有:
Ke = Ka ;
Kf = Kd ； 直线a、b、c、d、e、f的方程为: (1)Ya= Ka(H)HYaj)； (2)Yb= Κ“\-Χ_’Λ)+Υ_’」)；
(3)Yc = Kc(Xc-X(i+1；J+1))+Y
(i+l’j+1);
(4)Yd= Kd(Xd-Xa,J+1))+Y(i,J+1)；
(5)Ye = Ke (Xe-X)+Y ；
(6)Yf= Kf (Xf-X)+Y ； 联立直线方程⑷与(5)、(2)与(5)、⑴与(6)、(3)与(6)，分别解出交点E、F的垂直坐标及交点G、H的水平坐标:
Ye= [KdYe (Xa, J+1)-X) -KeYa,J+1)+KdY] / (Kd-Ke)；
Yf= [KbKe (Xa+1, j0-X) -KeY(i+1, ^+KJ] / (Kb-Ke)；
Xg= [(YfX-Y)-(KaX(U-Yaj))V(Kf-Ka)；
Xh= [ (YfX-Y) - (KcX(i+1, J+1)-Y(i+1, J+1)) ] / (Kf-Kc)； 设相邻标志点的实际间距为h，根据相似三角形原理，线段GN、EN的实际距离可表示为:
Sgn= I h (Xg-X)/(Xg-Xh) I
Sen= I h (Ye-Y)/(Ye-Yf) I。
3.根据权利要求2所述的预置位偏差计算方法，其特征在于，在步骤e中，根据所述SGN和SEN计算所述N距离中心参考点的水平距离Sx和垂直距离Sy，包括如下步骤: 选定标志点坐标纸中心区域的标志点M(X(iQ, j0), Y(i0, j。))为中心参考点，即以M为参考来计算预置位映射点距离点M的水平距离和垂直距离，则包含在第i行第j列单元四边形中的点N距离中心参考点的水平距离Sx和垂直距离Sy为:
Sx = I (j-j0) h+SGN I
Sy = I (1-10) h+SEN I。
4.根据权利要求3所述的预置位偏差计算方法，其特征在于，在步骤f中，计算预置位偏差包括如下步骤: 设预置位水平方向偏差距离Λ Sx和垂直偏差距离Λ Sy分别为:
ASx= I Sx-Sx0 I
Δ Sy = I Sy-Sy0 I。
【文档编号】G06F19/00GK104134005SQ201410378711
【公开日】2014年11月5日 申请日期:2014年7月31日 优先权日:2014年7月31日
【发明者】叶晨, 赵俊芳, 王杉, 刘鹏翔 申请人:天津市亚安科技股份有限公司