指针式仪表识别方法，确定指针式仪表图像分割阈值，为排除指针阴影的影响所设计的新方法

基于三次样条插值的指针式仪表图像分割方法，首先提取包含指针的圆形区域作为信息图，对信息图进行中值滤波、灰度拉伸；其次使用16个固定的阈值对信息图进行分割，并统计指针的长度和面积；最后使用三次样条插值法对面积序列插值并计算其一阶差分，以指针长度大于信息图半径的最小阈值作为起点，求得一阶差分的第一个波谷点，以该点对应的灰度值作为分割阈值对信息图进行图像分割。解决了指针阴影被分割为指针的问题。

function [erzhihuayuzhixinxitulashenminlashenmax]=newend（Lrr1）
tic %L为仪表灰度图  r为信息区域半径  r1为转到轴心半径
%%%%%%%%%%%%%%%%%%%灰度拉伸增大对比度%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
[counts0 ~] = imhist（L）;
for i=255:-1:1
    if counts0（i）~=0
        max0=i-1; % 直方图最后一个不为0的灰度值
        break;
    end
end
for i=2:255 
    if counts0（i）~=0
        min0=i-1; %直方图第一个不为0的灰度值
        break;
    end
end             %拉伸区域有原始图min0和max1决定
lashenL=imadjust（L[min0/255 max0/255][01]）;%对信息图拉伸
lashenL=uint8（lashenL）;%转换数据类型
lashenmin=min0;lashenmax=max0;
%%%%%%%%%%%%%%%%%%%%求带入的灰度值数组%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Tk=0:15:255;
Tk=Tk+1;%方便作图、拟合、给矩阵赋值 下标从1开始
k=18;
%%%%%%%计算模块，得到对应灰度值二值化后的指针所在连通域的长度，面积%%%%%%
changdu=zeros（1256）;
areanice=zeros（1256）;
for i=1:k
    if i~=k  %避免把圆形区域以外的填充色计算在内
        erzhihua=im2bw（lashenL（Tk（i）-1）/255）;
    else
        erzhihua=im2bw（lashenL（Tk（i）-2）/255）;%最后一个用254分割，排除圆形以外的区域
    end
    erzhihua=~erzhihua; %取反便于调用函数
    [zhizhentumaxarea]=maxandlabel（erzhihua）;%CV是只保留最大连通域的黑底图，maxarea为最大连通域的像素数
    chang=regionprops（zhizhentu‘BoundingBox‘）;
    changdu（Tk（i））=sqrt（chang.BoundingBox（13）^2+chang.BoundingBox（14）^2）;%取最大包络矩形的最大边长，非指针长度
    areanice（Tk（i））=maxarea;
    %imwrite（zhizhentu[‘F:\MATLABCODE\idtry\二值化系列图片\‘num2str（Tk（i）-1）‘.jpg‘]）;%保存到制定文件夹
end
 %%%%%%%%%%%%%%%%%%%%%%%%计算面积落差做大处作为拟合上限mianjiok%%%%%%%%%%%%
mianjidif=abs（areanice（Tk（1:k-1））-areanice（Tk（2:k）））;
[~ind]=sort（mianjidif‘descend‘）;%ind为面积降序排列的下角标k索引值
mianjiok=ind（11）;
%%%%%%%%%%%%%%%%%%计算指针所在连通域达到半径处作为zhizhenok下限%%%%%%%%%%%%
for i=1:k
    if changdu（Tk（i））>r+0.5*r1 %将转动轴心半径也考虑在内
        zhizhenok=i;
        break;
    end
end
%%%%%%%%%%%%%%%%%%%%%获取拟合横纵坐标，并作5阶曲线拟合%%%%%%%%%%%%%%%%%%%% 
dianshu=mianjiok-zhizhenok+1;%拟合数据点的个数
niheyuzhi0=zeros（1dianshu）;
nihearea0=zeros（1dianshu）;
for i=zhizhenok:mianjiok
    niheyuzhi0（i-zhizhenok+1）=Tk（i）-1;
    nihearea0（i-zhizhenok+1）=areanice（Tk（i））;
end
zeroshu=0;
for i=1:dianshu-1
    if nihearea0（i）==nihearea0（i+1）
        nihearea0（i）=0;
        zeroshu=zeroshu+1;
    end
end
ig=1;