本发明涉及一种适用于电力作业机器人的颜色型状态指示器识别方法，所属技术领域为数字图像处理、模式识别和机器学习领域。
背景技术：
近年来中国经济飞速发展，电力是各行各业发展的基石。在电力调度过程中主要依据断路器辅助结点判断开关的状态和位置。由于腐蚀、磨损、老化等原因，致使辅助开关切换不到位，有时不能正确判断开关的真实位置，为调度指挥提供了错误信息。然而主开关设备均联动一个合分指示器，通过观测指示器的合分状态能识别出电力高压断路器开关的合分状态和故障告警。
常见指示器主要分为两种：颜色型指示器和指向型指示器。
当前颜色型状态指示器主要依靠变电站值班人员观测。依靠人力观测有若干问题：
第一、对观测人员专业技能要求高，需要专门培训；
第二、工作环境不安全，变电站存在高压和放电问题，可能对观测检人员造成伤害；
第三、除去观测任务，还有日常巡检任务。巡检人员的巡检周期较长，不能及时发现问题。
随着技术发展，图像处理技术在各领域得到了广泛应用。在电力巡检方向图像处理技术也取得了喜人的成绩。机器人技术今年发展迅猛，逐渐代替了值班人员观测和巡检。通过图像处理技术和机器人巡检相结合不但节约了人力成本，避免人在危险的环境下工作，而且可以及时的发现异常，做出处理。
现有的技术依旧存在一些问题：
第一、机器人每次停靠的位置有差异，不能保证每次拍摄的视角完全一致；
第二、外界环境干扰问题，如光照过强，光照过暗，大雾等；
技术实现要素：
本发明目的就是为了解决上述问题，提供一种适用于电力作业机器人的颜色型状态指示器识别方法。它能快速准确的识别颜色型状态指示器及其状态，可使调度员能真实地了解开关的位置，调度指挥更加准确及时。
本发明采用以下技术方案：
一种适用于电力作业机器人的颜色型状态指示器识别方法，包括如下步骤：
(1-1)采集并标定指示器并保存相关信息，保存的相关信息包括指示器位置、指示器图像、指示器图像特征、指示器检测方法和指示器当前状态；
(1-2)采集指示器当前时刻的图像为待检测图像；
(1-3)根据检测方法，对步骤(1-2)获取的图像做不同预处理；
(1-4)根据指示器检测方法，检测步骤(1-3)预处理后图像中的指示器；
(1-5)根据颜色型状态指示器特征，判定指示器状态。
所述步骤(1-1)中指示器检测方法有4种：矩形指示器检测方法、椭圆指示器检测方法、颜色地图指示器检测方法和特征匹配指示器检测方法，
矩形指示器检测方法适用于指示器为矩形且无边框；
椭圆指示器检测方法适用于指示器为圆形或者椭圆且无边框；
颜色地图指示器检测方法适用于指示器形状无要求，指示器形状为矩阵、圆形或异形，指示器颜色与背景色不接近；
特征匹配指示器检测方法适用于有边框的指示器。
所述步骤(1-3)预处理是为了消除采集待检测图像噪声和后续计算加速；指示器检测方法为矩形指示器检测方法或椭圆指示器检测方法时，预处理的步骤为：
步骤(1-3-1)将采集的待检测图像灰度化，得到灰度图像；
步骤(1-3-2)统计灰度图像距离参考亮度的灰度平均偏移值：
其中：E代表平均偏移值；mean代表参考偏移值，一般取128；G(i,j)代表图像在(i,j)处灰度值；w代表图像宽度；h代表图像高度；
统计图像加权偏移：
其中D代表加权偏移；k代表灰度值，取值范围0-255；E代表灰度平均偏移值；mean代表参考偏移值，取128；Hist(k)代表图像中灰度值为k的点个数；w代表图像宽度；h代表图像高度；如果|E|＞D；图像亮度存在异常，E>0表示过量，E<0表示过暗，根据E的值调整gamma校正的变换参数，得到gamma校正图像；
步骤(1-3-3)计算gamma校正图像的双边滤波去噪图像；
步骤(1-3-4)对双边滤波图像做自适应尺度变换，得到待检测预处理后图像；
自适应尺度变换具体为：确定一个合适的变换因子scale，变换因子scale计算公式如下：
scale＝max(min(1，scale_X)，min(1，scale_Y))
scale_X＝sw/w，scale_Y＝sh/h
其中scale代表变换因子；scale_X代表X方向变换因子；scale_Y代表Y方向变换因子；w代表待检图像宽度，h代表待检图像高度；sw代表参考图像宽度，取1920，sh代表参考图像高度，取1080。
步骤(1-3-5)对自适应尺度变换图像做自适应canny变换。自适应canny变换会自动计算canny变换的高低阈值，比普通canny变换具有更好的鲁棒性。具体步骤为：
首先按照以下公式计算x方向和y方向的一阶导数。
其中：G表示原始图像，Gx表示G在x方向求一阶导，Gy表示G在y方向求一阶导。
然后计算融合梯度图像；
magGrad(i，j)＝Gx(i，j)+Gy(i，j)
mag＝max(magGrad(i，j)，0)
其中：magGrad为融合梯度图；mag为图像梯度峰值。
再后压缩融合梯度图像灰度级到0-mag，统计灰度直方图；
pmagGrad(i，j)＝magGrad(i，j)/bin_size
bin_size＝mag/NUM_BINS
其中：pmagGrad为压缩后的梯度图，bin_size位压缩比例因子，NUM_BINS为压缩的灰度级，一般取64。
接着累加直方图，当能量大于阈值时映射回原始灰度级作为高阈值；
high_thresh＝high_thresh*bin_size
其中：hist代表灰度直方图；high_thresh代表高阈值；w代表图像宽度，h代表图像高度not_edge代表非边界占比，一般取0.95；
通过高阈值计算出低阈值：
low_thresh＝high_thresh*ratio
其中low_thresh代表低阈值；ratio代表高低阈值比例，一般取0.3；
最后计算出高低阈值后再做canny变换，得到边界图像。
所述步骤(1-3)检测方法为颜色地图时，预处理涉及对待检测图像做自适应尺度变换、图像由RGB空间转换到HSI空间，针对I分量做双边滤波，针对S分量做直方图均衡化和图像由HSI空间转换到RGB空间，具体步骤为：
(1-3-6)自适应尺度变换得到尺度变换图像：确定一个合适的变换因子scale，变换因子scale计算公式如下：
scale＝max(min(1，scale_X)，mim(1，scale_Y))
scale_X＝sw/w，scale_Y＝sh/h
其中scale代表变换因子；scale_X代表X方向变换因子；scale_Y代表Y方向变换因子；w代表待检图像宽度，h代表待检图像高度；sw代表参考图像宽度，取1920，sh代表参考图像高度，取1080；
(1-3-7)把尺度变换图像转换到HSI色彩空间：
转换公式如下：
其中为HSI色彩空间色调分量的角度值；R、G、B分别为红、绿、蓝三种颜色分量；H为色调分量，S为饱和度分量，I为亮度分量；
(1-3-8)针对步骤(1-3-7)中的I分量进行双边滤波，针对步骤(1-3-7)中的S分量进行直方图均衡化增强；
(1-3-9)针对步骤(1-3-8)的图像做转换回RGB色彩空间处理；
时，B＝I*(1-S)，G＝3*I-B-R
时，R＝I*(1-S)，B＝3*I-G-R
时，G＝I*(1-S)，R＝3*I-G-B
其中：R、G、B分别为红、绿、蓝三种颜色分量；H为色调分量，S为饱和度分量，I为亮度分量。
所述步骤(1-4)检测方法为矩形时，指示器检测涉及到形态学膨胀处理和矩形检测。具体为：
(1-4-1)膨胀步骤(1-3)中预处理后图像的边界图像；
(1-4-2)统计膨胀后图像轮廓，去除不符合条件的部分包括：去除尺寸过小或者过大的轮廓；去除长宽比不符合条件的轮廓；去除面积和周长不符合条件的轮廓和去除直线分布不符合要求的轮廓；
(1-4-3)计算符合条件的轮廓包围矩形；
(1-4-4)非极大值抑制合并矩形，具体步骤为：计算矩形重叠面积比例，重叠面积大于阈值时，标记当前面积和面积重叠比例；当矩形重叠面积为局部极值时合并矩形，成为新的矩形加入计算。
所述步骤(1-4)检测方法为椭圆时，指示器检测涉及到椭圆检测。根据椭圆的几何特性检测，速度和精度比传统的hough圆检测高。特别是针对大离心率椭圆，hough检测几乎无法胜任，本发明的方法可以有很高的检测精度。
(1-4-5)检测步骤(1-3)预处理后图像的边界图像弧段的凹凸性：
其中：图像x方向一阶导数Gx；图像y方向一阶导数Gy，其中(i,j)表示图像点的坐标；
(1-4-6)抛弃过短和包含矩形过小的弧；
(1-4-7)将凸弧分为第一和第三象限，凹弧分为第二和第四象限；具体步骤为：
其中：area_b代表弧线下方面积，area_p代表弧线上方面积。
(1-4-8)三个象限上分别取弧a，b，c；
(1-4-9)连接a的起点和b的中点作为起始参考弦，利用夹逼方法求出一组平行弦，两两连接平行弦中点得到直线。分别针对a，b，c执行上述操作得到一系列直线l1，l2，l3...ln；
(1-4-10)计算出不同组直线的交点c1，c2，c3...cm，当点落在一个较小邻域内，可认为这组弧能构成椭圆，否则转到步骤(1-4-8)；
(1-4-11)如果构成椭圆，按照以下公式估算椭圆方程参数。
γ＝q1q2-q3q4
β＝(q3q4+1)(q1+q2)-(q1q2+1)(q3+q4)
θ＝cos-1Kp
B＝A*Np
其一般方程为：
x′＝x-x0
y′＝y-y0
其中：q1为a,b弧平行弦斜率，q2为a,b弧平行弦终点连线斜率，q3为b,c弧平行弦斜率，q4为b,c弧平行弦终点连线斜率，γ，β，Kp，Np，Ax，x′，y′为计算过程中间变量，θ为椭圆倾斜角度，(x0,y0)为椭圆中心点，A为椭圆长轴，B为椭圆短轴；
(1-4-12)a,b,c弧段上的点(xi,yi)满足下式则认为落在椭圆上，统计占比为socre1，
x′＝xi-x0
y′＝yi-y0
其中：d为点和椭圆的距离，一般取1；
(1-4-13)计算弧段与1/4椭圆相似度score2；
(1-4-14)按照下式计算椭圆得分：
score＝w1*score1+w2*score2
w1+w2＝1
其中：score为椭圆得分，w1，w2为分别权重，一般都取0.5；
(1-4-15)按照得分由高到低给椭圆排序；
(1-4-16)非极大值抑制合并椭圆，具体步骤为：计算所有椭圆包含矩形；将得分最高椭圆的包含矩形作为初始值，做矩形非极大值抑制合并；在剩余集合中按照做矩形非极大值抑制合并，直到集合为空。
所述步骤(1-4)检测方法为数字地图时，指示器检测涉及到颜色地图搜索。
(1-4-17)利用标定时建立的目标区域数字地图，先按照x方向搜索整个数字地图，得到n个条状图像bar1，bar2，……barn；
其中：min_num为最小满足条件点数，一般取45；G为彩色图像；(i,j)为图像坐标；R为图像行数；F为关系函数；
(1-4-18)针对每个条状图像按照y方向搜索，得到一系列小的子区域图像；
其中：min_num为最小满足条件点数，一般取45；bar为切割后的条状图像；(i,j)为图像坐标；C为图像列数；F为关系函数；
(1-4-19)搜索整个数字地图具体步骤：搜索落在指定RGB区间的连续像素；搜索R，G，B分量关系与模板图像一致的连续像素；统计这些像素的起点和终点，波峰和波谷补齐区域；
(1-4-20)然后按照先y后x方向搜索整个待检测图像，得到一系列子区域图像；
(1-4-21)利用子区域坐标关系确认指示器区域。
所述步骤(1-4)检测方法为特征匹配时，指示器检测涉及到尺度不变(SIFT)特征匹配。
SIFT特征是图像的局部特征，其对旋转、尺度缩放、亮度变化保持不变性，对视角变化、仿射变换、噪声也保持一定程度的稳定性；独特性好，信息量丰富，适用于在海量特征数据库中进行快速、准确的匹配；多量性，即使少数的几个物体也可以产生大量的SIFT特征向量；高速性，经优化的SIFT匹配算法甚至可以达到实时的要求。
检测步骤为：
(1-4-22)提取图像特征为尺度不变特征变换(SIFT)特征；
(1-4-23)利用模板图像的SIFT特征建立kd树；
(1-4-24)待匹配图像的SIFT特征在kd树上做二叉搜索，在搜索过程中利用最小优先队列建立回溯搜索的索引，最小优先队列键值为对应维度特征值差值的绝对值；按照最小优先队列顺序回溯搜索，当达到最小优先队列为空或者达到搜索次数上限时停止搜索；当模板图像特征点对应多个待匹配特征点时，仅保留搜索过程中的最优值和次优值；搜索完毕后，根据匹配特征点与最优值和次优值的欧氏距离关系筛选出最终匹配结果；参考方法是min_Dis＜max_Dis*0.6；其中min_Dis是特征点与最优值的欧氏距离，max_Dis是特征点与次优值的欧氏距离；
(1-4-25)根据待匹配图像和匹配图像匹配点的坐标关系计算对应变换矩阵H；
(1-4-26)利用变换矩阵H的逆变换矩阵H_inv将待匹配图像映射到模板图像同一视角。
所述步骤(1-5)判定指示器状态，涉及到根据颜色型状态指示器特征判定。具体步骤为：检测落在指定RGB区间的像素占比；检测R，G，B分量关系与模板图像一致的像素占比；利用两个占比关系推导当前指示器状态。
判定规则：与标定时一致，否则与标定相反；
其中：G代表彩色图像，H代表图像高度，W代表图像宽度，F为关系函数，s代表阈值，一般取0.6。
附图说明
图1为本发明的算法流程图；
具体实施方式
下面结合附图与实施例对比发明做进一步说明。
如图1所示。本发明的方法步骤如下：
第一步：采集并标定指示器并保存相关信息，保存的相关信息包括指示器位置、指示器图像、指示器图像特征、指示器检测方法和指示器当前状态；
第二步：采集指示器当前时刻的图像为待检测图像，对图像做预处理。根据不同的检测方法，自适应不同的预处理方式。
矩形指示器检测方法适用于指示器为矩形且无边框；
椭圆指示器检测方法适用于指示器为圆形或者椭圆且无边框；
颜色地图指示器检测方法适用于指示器形状无要求，指示器形状为矩阵、圆形或异形，指示器颜色与背景色不接近；
特征匹配指示器检测方法适用于有边框的指示器。
当指示器检测方法为指示器检测方法为矩形指示器检测方法或椭圆指示器检测方法时，预处理的步骤为：(1)图像灰度化；(2)自适应亮度矫正；(3)双边滤波去噪；(4)自适应尺度变换，(5)自适应canny滤波检测边界。
步骤(1-3-1)将采集的待检测图像灰度化，得到灰度图像；
灰度化公式为：Gray＝0.299*R+0.587*G+0.114*B+0.5，其中Gray为灰度值，R、G、B分别为红、绿、蓝三种颜色分量。
步骤(1-3-2)自适应亮度矫正：灰度图像距离参考亮度的灰度平均偏移值，统计图像加权偏移，其中：E代表平均偏移值；mean代表参考偏移值，一般取128；G(i,j)代表图像在(i,j)处灰度值；w代表图像宽度；h代表图像高度；其中D代表加权偏移；k代表灰度值，取值范围0-255；E代表灰度平均偏移值；mean代表参考偏移值，取128；Hist(k)代表图像中灰度值为k的点个数；w代表图像宽度；h代表图像高度；如果|E|＞D，图像亮度存在异常，E>0表示过量，E<0表示过暗。根据E的值调整gamma校正的变换参数。
步骤(1-3-3)计算gamma校正图像的双边滤波去噪图像；同时考虑了灰度值和空间位置关系，去噪同时不改变边界位置。
步骤(1-3-4)对双边滤波图像做自适应尺度变换，得到待检测预处理后图像：这里的核心是确定一个合适的变换因子，保证后续计算加速且不会影响精度。计算公式如下：
scale＝max(min(1，scale_X)，min(1，scale_Y))
scale_X＝sw/w，scale_Y＝sh/h
其中scale代表变换因子，scale_X代表X方向变换因子，scale_Y代表Y方向变换因子；w代表图像宽度，h代表图像高度；sw代表参考图像宽度，一般取1920，sh代表参考图像高度，一般取1080；
步骤(1-3-5)对自适应尺度变换图像做自适应canny变换。自适应canny变换会自动计算canny变换的高低阈值，比普通canny变换具有更好的鲁棒性。具体步骤为：
首先按照以下公式计算x方向和y方向的一阶导数。
其中：G表示原始图像，Gx表示G在x方向求一阶导，Gy表示G在y方向求一阶导。
然后计算融合梯度图像；
magGrad(i，j)＝Gx(i，j)+Gy(i，j)
mag＝max(magGrad(i，j)，0)
其中：magGrad为融合梯度图；mag为图像梯度峰值。
再后压缩融合梯度图像灰度级到0-mag，统计灰度直方图；
pmagGrad(i，j)＝magGrad(i，j)/bin_size
bin_size＝mag/NUM_BINS
其中：pmagGrad为压缩后的梯度图，bin_size位压缩比例因子，NUM_BINS为压缩的灰度级，一般取64。
接着累加直方图，当能量大于阈值时映射回原始灰度级作为高阈值；
high_thresh＝high_thresh*bin_size
其中：hist代表灰度直方图；high_thresh代表高阈值；w代表图像宽度，h代表图像高度not_edge代表非边界占比，一般取0.95；
通过高阈值计算出低阈值：
low_thresh＝high_thresh*ratio
其中low_thresh代表低阈值；ratio代表高低阈值比例，一般取0.3；
最后计算出高低阈值后再做canny变换，得到边界图像。
当指示器检测方法为颜色地图指示器检测方法时，预处理的步骤为：(1)自适应尺度变换，(2)图像由RGB空间转换到HSI空间，(3)针对I分量做双边滤波，(4)针对S分量做直方图均衡化，(5)图像由HSI空间转换到RGB空间。
步骤(1-3-6)自适应尺度变换得到尺度变换图像：确定一个合适的变换因子scale，变换因子scale计算公式如下：
scale＝max(min(1，scale_X)，min(1，scale_Y))
scale_X＝sw/w，scale_Y＝sh/h
其中scale代表变换因子；scale_X代表X方向变换因子；scale_Y代表Y方向变换因子；w代表待检图像宽度，h代表待检图像高度；sw代表参考图像宽度，取1920，sh代表参考图像高度，取1080；
步骤(1-3-7)把尺度变换图像转换到HSI色彩空间：
转换公式如下：
其中为HSI色彩空间色调分量的角度值；R、G、B分别为红、绿、蓝三种颜色分量；H为色调分量，S为饱和度分量，I为亮度分量；
步骤(1-3-8)针对步骤(1-3-7)中的I分量进行双边滤波，针对步骤(1-3-7)中的S分量进行直方图均衡化增强；
步骤(1-3-9)针对步骤(1-3-8)的图像做转换回RGB色彩空间处理；
时，B＝I*(1-S)，G＝3*I-B-R
时，R＝I*(1-S)，B＝3*I-G-R
时，G＝I*(1-S)，R＝3*I-G-B
其中：R、G、B分别为红、绿、蓝三种颜色分量；H为色调分量，S为饱和度分量，I为亮度分量。
当指示器检测方法为特征匹配检测方法时，不做预处理。
第三步：检测指示器所在位置。根据不同的检测方法自适应不同的检测方式。
指示器检测方法为矩形检测方法时，指示器检测步骤为：(1)针对步骤(1-3-5)的边界图像做膨胀处理，(2)计算矩形。
步骤(1-4-1)膨胀步骤(1-3)中预处理后图像的边界图像；
步骤(1-4-2)统计膨胀后图像轮廓，去除不符合条件的部分包括：去除尺寸过小或者过大的轮廓；去除长宽比不符合条件的轮廓；去除面积和周长不符合条件的轮廓和去除直线分布不符合要求的轮廓；
步骤(1-4-3)计算符合条件的轮廓包围矩形；
步骤(1-4-4)非极大值抑制合并矩形，具体步骤为：计算矩形重叠面积比例，重叠面积大于阈值时，标记当前面积和面积重叠比例；当矩形重叠面积为局部极值时合并矩形，成为新的矩形加入计算。
指示器检测方法为椭圆检测方法时，指示器检测步骤为：(1)针对步骤(1-3-5)的边界图像做椭圆检测。根据椭圆的几何特性检测，速度和精度比传统的hough圆检测高。特别是针对大离心率椭圆，hough检测几乎无法胜任，本发明的方法可以有很高的检测精度。
步骤(1-4-5)检测步骤(1-3-5)预处理后图像的边界图像弧段的凹凸性：
其中：图像x方向一阶导数Gx；图像y方向一阶导数Gy，其中(i,j)表示图像点的坐标；
步骤(1-4-6)抛弃过短和包含矩形过小的弧；
步骤(1-4-7)将凸弧分为第一和第三象限，凹弧分为第二和第四象限；具体步骤为：
其中：area_b代表弧线下方面积，area_p代表弧线上方面积。
步骤(1-4-8)三个象限上取弧a，b，c；
步骤(1-4-9)连接a的起点和b的中点作为起始参考弦，利用夹逼方法求出一组平行弦，两两连接平行弦中点得到直线。分别针对a，b，c执行上述操作得到一系列直线l1，l2，l3...ln；
步骤(1-4-10)计算出不同组直线的交点c1，c2，c3...cm，当点落在一个较小邻域内，可认为这组弧能构成椭圆，否则转到步骤(1-4-8)；
步骤(1-4-11)如果构成椭圆，按照以下公式估算椭圆方程参数。
γ＝q1q2-q3q4
β＝(q3q4+1)(q1+q2)-(q1q2+1)(q3+q4)
θ＝cos-1Kp
B＝A*Np
其一般方程为：
x′＝x-x0
y′＝y-y0
其中：q1为a,b弧平行弦斜率，q2为a,b弧平行弦终点连线斜率，q3为b,c弧平行弦斜率，q4为b,c弧平行弦终点连线斜率，γ，β，Kp，Np，Ax，x′，y′为计算过程中间变量，θ为椭圆倾斜角度，(x0,y0)为椭圆中心点，A为椭圆长轴，B为椭圆短轴；
步骤(1-4-12)a,b,c弧段上的点(xi,yi)满足下式则认为落在椭圆上，统计占比为socre1，
x′＝xi-x0
y′＝yi-y0
其中：d为点和椭圆的距离，一般取1；
步骤(1-4-13)计算弧段与1/4椭圆相似度score2；
步骤(1-4-14)按照下式计算椭圆得分：
score＝w1*score1+w2*score2
w1+w2＝1
其中：score为椭圆得分，w1，w2为权重，一般都取0.5；
步骤(1-4-15)按照得分由高到低给椭圆排序；
步骤(1-4-16)非极大值抑制合并椭圆，具体步骤为：计算所有椭圆包含矩形；将得分最高椭圆的包含矩形作为初始值，做矩形非极大值抑制合并；在剩余集合中按照做矩形非极大值抑制合并，直到集合为空。
指示器检测方法为数字地图检测方法时，指示器检测涉及到颜色地图搜索。
步骤(1-4-17)利用标定时建立的目标区域数字地图，先按照x方向搜索整个数字地图，得到n个条状图像bar1，bar2，……barn；
其中：min_num为最小满足条件点数，一般取45；G为彩色图像；(i,j)为图像坐标；R为图像行数；F为关系函数；
步骤(1-4-18)针对每个条状图像按照y方向搜索，得到一系列小的子区域图像；
其中：min_num为最小满足条件点数，一般取45；bar为切割后的条状图像；(i,j)为图像坐标；C为图像列数；F为关系函数；
步骤(1-4-19)搜索整个数字地图具体步骤：搜索落在指定RGB区间的连续像素；搜索R，G，B分量关系与模板图像一致的连续像素；统计这些像素的起点和终点，波峰和波谷补齐区域；
步骤(1-4-20)然后按照先y后x方向搜索整个待检测图像，得到一系列子区域图像；
步骤(1-4-21)利用子区域坐标关系确认指示器区域。
指示器检测方法为特征匹配检测方法时，指示器检测涉及到尺度不变(SIFT)特征匹配。
SIFT特征是图像的局部特征，其对旋转、尺度缩放、亮度变化保持不变性，对视角变化、仿射变换、噪声也保持一定程度的稳定性；独特性好，信息量丰富，适用于在海量特征数据库中进行快速、准确的匹配；多量性，即使少数的几个物体也可以产生大量的SIFT特征向量；高速性，经优化的SIFT匹配算法甚至可以达到实时的要求。
检测步骤为：
步骤(1-4-22)提取图像特征为尺度不变特征变换(SIFT)特征；
步骤(1-4-23)利用模板图像的SIFT特征建立kd树；
步骤(1-4-24)待匹配图像的SIFT特征在kd树上做二叉搜索，在搜索过程中利用最小优先队列建立回溯搜索的索引，最小优先队列键值为对应维度特征值差值的绝对值；按照最小优先队列顺序回溯搜索，当达到最小优先队列为空或者达到搜索次数上限时停止搜索；当模板图像特征点对应多个待匹配特征点时，仅保留搜索过程中的最优值和次优值；搜索完毕后，根据匹配特征点与最优值和次优值的欧氏距离关系筛选出最终匹配结果；参考方法是min_Dis＜max_Dis*0.6；其中min_Dis是特征点与最优值的欧氏距离，max_Dis是特征点与次优值的欧氏距离；
步骤(1-4-25)根据待匹配图像和匹配图像匹配点的坐标关系计算对应变换矩阵H；
步骤(1-4-26)利用变换矩阵H的逆变换矩阵H_inv将待匹配图像映射到模板图像同一视角。
第四步：检测落在指定RGB区间的像素占比；检测R，G，B分量关系与模板图像一致的像素占比；利用两个占比关系推导当前指示器状态。
判定规则：与标定时一致，否则与标定相反；
其中：G代表彩色图像，H代表图像高度，W代表图像宽度，F为关系函数，s代表阈值，一般取0.6。统计目标区域的R分量，G分量和B分量落在指定区间的百分比r1；统计R分量，G分量和B分量间满足关系函数的百分比r2。根据统计值判定当前指示器状态。
根据下式判定指示器当前状态：
r1≥s或r2≥s，当前指示器状态可认定为与标定时一致；
r1＜s且r2＜s，当前指示器状态可认定为与标定时相反；
上述虽然结合附图对本发明的具体实施方式进行了描述，但非对本发明保护范围的先知，所属领域计时人员应该明白，在本发明的技术方案基础上，本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围内。