【万用表识别】数字仪表识别【含Matlab源码 693期】

✅博主简介：热爱科研的Matlab仿真开发者，修心和技术同步精进，Matlab项目合作可私信。
🍎个人主页：海神之光
🏆代码获取方式：
海神之光Matlab王者学习之路—代码获取方式
⛳️座右铭：行百里者，半于九十。

更多Matlab仿真内容点击👇
Matlab图像处理（进阶版）
路径规划（Matlab）
神经网络预测与分类（Matlab）
优化求解（Matlab）
语音处理（Matlab）
信号处理（Matlab）
车间调度（Matlab）

⛄一、简介

1 灰度化 （grayscale）

将彩色图像转化为灰度图像的过程称为图像灰度化。彩色图像中的像素值由RGB三个分量决定，每个分量都有0-255（256种）选择，这样一个像素点的像素值可以有1600万种可能（256256256），而灰度图的像素点的像素值是RGB三个分量值相同的一种特殊的彩色图像， 只有256种可能。所以在图像处理中，往往将各种图像首先灰度化成灰度图像以便后续处理，降低计算量。灰度是指只含亮度信息，不含色彩信息的图像。黑白照片就是灰度图，特点是亮度由暗到明，变化是连续的。灰度图像的描述与彩色图像一样仍然反映了整幅图像的整体和局部的色度和亮度等级的分布和特征，

使用灰度图的好处：
① RGB的值都一样。
② 图像数据即调色板索引值，就是实际的RGB值，也就是亮度值。
③ 因为是256色调色板，所以图像数据中一个字节代表一个像素，很整齐。
所以，做图像处理时一般都采用灰度图。
要表示灰度图，就需要把亮度值进行量化，有四种方法：

（1）分量法
将彩色图像中的三分量的亮度作为三个灰度图像的灰度值，可根据应用需要选取一种灰度图像。
（2）最大值法
将彩色图像中的三分量亮度的最大值作为灰度图的灰度值。
（3）均值法
将彩色图像中的三分量亮度求平均得到灰度图的灰度值。
（4）加权平均法
根据重要性及其它指标，将三个分量以不同的权值进行加权平均。由于人眼对绿色的敏感最高，对蓝色敏感最低，因此，按下式对RGB三分量进行加权平均能得到较合理的灰度图像，f(i,j)=0.30R(i,j)+0.59G(i,j)+0.11B(i,j))。

2 二值化（binaryzation）

图像的二值化是将图像上的像素点的灰度值设置为0或255，也就是将整个图像呈现出明显的黑白效果。将256个亮度等级的灰度图像通过适当的阀值选取而获得仍然可以反映图像整体和局部特征的二值化图像。所有灰度大于或等于阀值的像素被判定为属于特定物体，其灰度值为255，否则这些像素点被排除在物体区域以外，灰度值为0，表示背景或者例外的物体区域。在数字图像处理中，二值图像占有非常重要的地位，首先，图像的二值化有利于图像的进一步处理，使图像变得简单，而且数据量减小，能凸显出感兴趣的目标的轮廓。其次，要进行二值图像的处理与分析，首先要把灰度图像二值化，得到二值化图像。
二值化的常用算法有：

全局二值化： 一幅图像包括目标物体、背景还有噪声，要想从多值的数字图像中直接提取出目标物体，最常用的方法就是设定一个全局的阈值T，用T将图像的数据分成两部分：大于T的像素群和小于T的像素群。将大于T的像素群的像素值设定为白色（或者黑色），小于T的像素群的像素值设定为黑色（或者白色）。全局二值化，在表现图像细节方面存在很大缺陷。为了弥补这个缺陷，出现了局部二值化方法。

局部二值化：按照一定的规则将整幅图像划分为N个窗口，对这N个窗口中的每一个窗口再按照一个统一的阈值T将该窗口内的像素划分为两部分，进行二值化处理。局部二值化也有一个缺陷。这个缺陷存在于那个统一阈值的选定。这个阈值是没有经过合理的运算得来，一般是取该窗口的平局值。这就导致在每一个窗口内仍然出现的是全局二值化的缺陷。为了解决这个问题，就出现了局部自适应二值化方法。

局部自适应二值化：在局部二值化的基础之上，将阈值的设定更加合理化。该方法的阈值是通过对该窗口像素的平均值E，像素之间的差平方P，像素之间的均方根值Q等各种局部特征，设定一个参数方程进行阈值的计算，例如：T=aE+bP+c*Q，其中a,b,c是自由参数。这样得出来的二值化图像就更能表现出二值化图像中的细节。

3 反色（inverse）
反色的实际含义是将R、G、B值反转。若颜色的量化级别是256，则新图的R、G、B值为255减去原图的R、G、B值。这里针对的是所有图，包括真彩图、带调色板的彩色图(又称为伪彩色图)、和灰度图。真彩图不带调色板，每个象素用3个字节，表示R、G、B三个分量。所以处理很简单，把反转后的R、G、B值写入新图即可，比如一个点的颜色为(0,0,0)，反色后为(255,255,255)。带调色板的彩色图，其位图中的数据只是对应调色板中的一个索引值，我们只需要将调色板中的颜色反转，形成新调色板，而位图数据不用动，就能够实现反转。

⛄二、部分源代码

%读数%20140102
close all
clear
clc

plantype = 2;% 1是读指针，2是读数字
file = ‘C:\Users\lenovo\Desktop\111186778Plan_Recognise\数字指针表盘识别总结\4.tif’;%指针
% file = ‘H:\picture\digital\28-20.tif’;%数字
I0 = imread(file);
figure;imshow(I0);title(‘待读数图片’);
scale = 1000/length(I0);
Irsz = imresize(I0,scale);
clear I0;
%%%%%%2.1：读指针%%%%%%
if plantype == 1
%%%%%%2.1.1:剪切图片，与模板图片相减%%%%%%
load plan1;
Ic = imcrop(Irsz,Rect);
Isub = Istd - Ic ;
level = graythresh(Isub);
Ib = im2bw(Isub,level);
% figure;imshow(Ib);title(‘bwimage’);%----------------------------------

%%%%%%2.1.2：二值化，找指针%%%%%%
topy = find(sum(Ib')>0,1, 'first');
topx = find(Ib(topy,:)>0,1,'first');
% figure(1);hold on;
% plot([O(1),topx],[O(2),topy]);
%%%%%%2.1.3：读数%%%%%%
zz = [topx,topy];
thetazz = -atan( (zz(2)-O(2))/(zz(1)-O(1)) );
if thetazz<0
    thetazz = thetazz + pi;
end
thetakd = -atan( (Kd(2,:)-O(2)) ./ (Kd(1,:)-O(1)) );
nidx = find(thetakd<0);
thetakd(nidx) = thetakd(nidx) + pi;

stepidx = find( thetakd<=thetazz,1,'first');
para1 = Fsv / (length(Kd)-1);
para2 = para1*(stepidx-2);
dushu = para2 + ( thetazz - thetakd(stepidx-1) ) /( thetakd(stepidx)-thetakd(stepidx-1) )*para1;
msgbox(sprintf('%0.3f',dushu),'Result');

%%%%%%2.2：读数字%%%%%%
else if plantype ==2
if length(size(Irsz))==3
Ig=rgb2gray(Irsz);
else
Ig=Irsz;
end
clear Irsz;
If = imfill(Ig);
Ib = im2bw(If,graythresh(If));
clear If;
% figure;imshow(Ib);%-----------------------------------

    SE = strel('square',25);
    Ied = imerode(Ib,SE);
    clear Ib;
    L = bwlabel(Ied);
    stats = regionprops(L,'Area');
    idx = find([stats.Area] == max([stats.Area]));
    [r,c] = find( L==idx );
    clear L;
    r1=min(r);r2=max(r);c1=min(c);c2=max(c);
    rect = [c1,r1,c2-c1,r2-r1];
    Ic = imcrop(Ig,rect);clear Ig;
    % figure;imshow(Ic);%-------------------------------------

    level = graythresh(Ic)/2;   %剪切之后的图片二值化
    Ib=im2bw(Ic,level);clear Ic;
    mp = mean(mean(Ib(1,:))+mean(Ib(end,:))...
         +mean(Ib(:,1))+mean(Ib(:,end)));
    if mp>0.5
        Ib = ~Ib;
    end
    % figure;imshow(Ib);%-------------------

    [m,n] = size(Ib);
    SE = strel('square',round(m/50)); %初步去噪，图片开运算（先腐蚀后膨胀），去除噪声
    Iop = imopen(Ib,SE);
     %  figure;imshow(Iop);%------------------   

    SE = strel('square',round(m/25));    
    Icl = imclose(Iop,SE);clear Iop; %闭运算，希望将组成数字的各LED段连接在一起
     %   figure;imshow(Icl);title('Icl')%-------------------
    Ith = bwmorph(Icl,'thin',inf);%图片细化，将数字变为由细线组成
     %  figure; imshow(Ith);%----------------- 

    sRI = sum(Ith');    %进一步去噪
    upidx = find(sRI>0, 1, 'first' );
    dnidx = find(sRI>0, 1, 'last');
    pse = round((dnidx-upidx)/40);
    SE = strel('square',pse);
    Ioc = imopen(Icl,SE);
    Ith = bwmorph(Ioc,'thin',inf);clear Iop;
    %  figure; imshow(Ith);title('Ith')%-----------------

    %%%%%%%%%%%%%%%%%%%%%%%%%对字符的处理识别%%%%%%%%%%%%%%%%%%%%%%%%
    %%%cut(认为读数是单排的)
    %%横切
    sRI = sum(Ith');   
    upidx = find(sRI>0, 1, 'first' );
    dnidx = find(sRI>0, 1, 'last');
    udidx = upidx:dnidx;
    ud = length(udidx);
    thrlen = ud/2; %判定LED段是否点亮的初步阈值
    Ith = Ith(udidx,:);
    % figure; imshow(Ith);%-----------------
    %%竖切
    sCI=sum(Ith);
    BsCI = (sCI>0); %将x轴投影矩阵scI变为只有0和1的矩阵
    dBsCI = diff(BsCI); %找BscI中为1的元素的范围的第一步，dBscI中的元素只有-1，0，1三种
    idxp = find(dBsCI==1)+1; %记录dBscI中为1的元素的位置，加1是对应BscI中1元素群开始位置
    idxn = find(dBsCI==-1);%记录dBscI中为-1的元素位置，对应BscI中1元素群结束位置，与idxp长度一样
    if length(idxp)~=length(idxn)
        idxn(find(idxn<idxp(1)))=[];
    end
    Cidxchar = [idxp;idxn]; %一个2行多列矩阵，每列对应一个字符，列的第一行为起始位置，第二行为结束位置
    dCidxchar = diff(Cidxchar);
    tdCidx = find(dCidxchar < thrlen/2);
    %%处理切割得到的符号，将“'”型字符归到上一个字符去
    if ~isempty(tdCidx)
        ltdCidx = length(tdCidx);
        delemt0 = [];
        delemt1 = [];
        delemt2 = [];
        for i=1:ltdCidx
            lridx = Cidxchar(1,tdCidx(i)):Cidxchar(2,tdCidx(i));   
            tcharI = Ith(:,lridx)';
            [m,n]=size(tcharI);
            if m>1
                sRtc = sum(tcharI);
            else
                sRtc = tcharI;
            end
            BsRtc = (sRtc>0);
            dBsRtc = diff(BsRtc);
            cupidx = find(dBsRtc==1,1,'first')+1;
            cdnidx = find(dBsRtc==-1, 1, 'last' );
            hc = cdnidx-cupidx;
            if ~isempty(hc)&&hc<thrlen/2&&cdnidx<ud*0.3
                delemt0 = [tdCidx(i),delemt0];
            end
            if ~isempty(hc)&&hc>thrlen/2 && hc<thrlen && cdnidx<ud*2/3
                delemt1 = [tdCidx(i),delemt1];
            end
            if ~isempty(hc)&&hc>thrlen/2 && hc<thrlen && cupidx<ud*2/3
                delemt2 = [tdCidx(i),delemt2];
            end
        end

        if ~isempty(delemt1)
            Cidxchar(2,delemt1-1) = Cidxchar(2,delemt1);
        end
        if ~isempty(delemt2)
            Cidxchar(1,delemt2+1) = Cidxchar(1,delemt2);
        end
        delemt = sort([delemt0,delemt1,delemt2]);
        if ~isempty(delemt) 
            Cidxchar(:,delemt) = [];
        end
   end

    nchar = length(Cidxchar);%读数包含的字符个数（包括小数点）
    %     %%%%显示切割效果%---------------------------------
    %     for i=1:nchar
    %         figure;
    %         imshow(Ith(:,Cidxchar(1,i):Cidxchar(2,i)));
    %     end

    %%%确定各LED段点亮与否    
    ntab = [ 
             1 1 1 1 1 1 0  %0
             0 1 1 0 0 0 0  %1
             1 1 0 1 1 0 1  %2
             1 1 1 1 0 0 1  %3
             0 1 1 0 0 1 1  %4
             1 0 1 1 0 1 1  %5
             1 0 1 1 1 1 1  %6
             1 1 1 0 0 0 0  %7
             1 1 1 1 1 1 1  %8
             1 1 1 1 0 1 1  %9
             0 0 0 0 0 0 0  %.
             0 0 0 0 0 0 1  %-
            ];
    result = [];  
    for i=1:nchar
        lridx = Cidxchar(1,i):Cidxchar(2,i); 
         if length(lridx) < thrlen/2
            charI = [zeros(ud,int16(thrlen)-length(lridx)),Ith(:,lridx)];
         else
            charI = Ith(:,lridx);
         end
         [m,n] = size(charI);
    %               figure;imshow(charI)%-------------------
        %%%%%agd段LED，由于字符可能有倾斜，取四分之一为判断标准
        xthrlen = thrlen/4;
        ayidx = 1:int16(m/8);      %a段y轴方向长度取为总长的八分之一
        asCchar = sum( charI( ayidx , :) ); %a段LED的x轴投影
        BasCchar = asCchar(find(asCchar<2));
        sBa = sum(BasCchar);
        fa = (sBa > xthrlen); %a段LED是否点亮 
        gyidx = int16(7*m/16):int16(9*m/16);
        gsCchar = sum(charI(gyidx,:));%g段LED的x轴投影
        BgsCchar = gsCchar(find(gsCchar<2));
        sBg = sum(BgsCchar);
        fg = (sBg > xthrlen); %g段LED是否点亮
        dyidx = int16(7*m/8):m;
        dsCchar = sum(charI(dyidx,:));    %d段LED的x轴投影
        BdsCchar = dsCchar(find(dsCchar<2));
        sBd = sum(BdsCchar);
        fd = (sBd > xthrlen); %d段LED是否点亮
        %%%%%fedc段LED，投影到y轴一般不倾斜，取二分之一为判断标准
        ythrlen = thrlen/2;
        fyidx = 1:int16(m/2); fxidx = 1:int16(n/2);
        fsRchar = sum(charI(fyidx,fxidx)');%f段LED的y轴投影
        BfsRchar = (fsRchar>0);
        sBf = sum(BfsRchar);
        ff = (sBf>ythrlen);               %f段LED是否点亮
        eyidx = int16(m/2):m; exidx = 1:int16(n/2);
        esRchar = sum(charI(eyidx,exidx)');%e段LED的y轴投影
        BesRchar = (esRchar>0);
        sBe = sum(BesRchar);
        fe = (sBe>ythrlen);               %e段LED是否点亮
        cyidx = int16(m/2):m; cxidx = int16(n/2):n;
        csRchar = sum(charI(cyidx,cxidx)');%c段LED的y轴投影
        BcsRchar = (csRchar>0);
        sBc = sum(BcsRchar);
        fc = (sBc>ythrlen);               %c段LED是否点亮
        byidx = 1:int16(m/2); bxidx = int16(n/2):n;
        bsRchar = sum(charI(byidx,bxidx)');%b段LED的y轴投影
        BbsRchar = (bsRchar>0);
        sBb = sum(BbsRchar);
        fb = (sBb>ythrlen);               %b段LED是否点亮

⛄三、运行结果

⛄四、matlab版本及参考文献

1 matlab版本
2014a

2 参考文献
[1]曲超然,陈立伟,王建生,王水根.基于深度学习的工业数字仪表识别算法研究[J].应用科技. 2022,49(02)

3 备注
简介此部分摘自互联网，仅供参考，若侵权，联系删除

🍅 仿真咨询
1 各类智能优化算法改进及应用
生产调度、经济调度、装配线调度、充电优化、车间调度、发车优化、水库调度、三维装箱、物流选址、货位优化、公交排班优化、充电桩布局优化、车间布局优化、集装箱船配载优化、水泵组合优化、解医疗资源分配优化、设施布局优化、可视域基站和无人机选址优化

2 机器学习和深度学习方面
卷积神经网络（CNN）、LSTM、支持向量机（SVM）、最小二乘支持向量机（LSSVM）、极限学习机（ELM）、核极限学习机（KELM）、BP、RBF、宽度学习、DBN、RF、RBF、DELM、XGBOOST、TCN实现风电预测、光伏预测、电池寿命预测、辐射源识别、交通流预测、负荷预测、股价预测、PM2.5浓度预测、电池健康状态预测、水体光学参数反演、NLOS信号识别、地铁停车精准预测、变压器故障诊断

3 图像处理方面
图像识别、图像分割、图像检测、图像隐藏、图像配准、图像拼接、图像融合、图像增强、图像压缩感知

4 路径规划方面
旅行商问题（TSP）、车辆路径问题（VRP、MVRP、CVRP、VRPTW等）、无人机三维路径规划、无人机协同、无人机编队、机器人路径规划、栅格地图路径规划、多式联运运输问题、车辆协同无人机路径规划、天线线性阵列分布优化、车间布局优化

5 无人机应用方面
无人机路径规划、无人机控制、无人机编队、无人机协同、无人机任务分配

6 无线传感器定位及布局方面
传感器部署优化、通信协议优化、路由优化、目标定位优化、Dv-Hop定位优化、Leach协议优化、WSN覆盖优化、组播优化、RSSI定位优化

7 信号处理方面
信号识别、信号加密、信号去噪、信号增强、雷达信号处理、信号水印嵌入提取、肌电信号、脑电信号、信号配时优化

8 电力系统方面
微电网优化、无功优化、配电网重构、储能配置

9 元胞自动机方面
交通流 人群疏散 病毒扩散 晶体生长

10 雷达方面
卡尔曼滤波跟踪、航迹关联、航迹融合