final.m

这几项M文件提供了扫描仪如何识别电子标签的方法

clear;
P1=load('200540-1.txt')'
PP=P1(351:450);
PS=interp(PP,10);%10内插
P=medfilt1(PS,9);%中值滤波
np=min(P); %信号最小值
mp=max(P); %信号最大值
dp=np+(mp-np)./4; %信号下四分之一
up=np+2.*(mp-np)./3;%信号上三分之二
%P=medfilt1(Pp,3);
N=[];
for n=2:999;
    pc1=P(n)-P(n-1);
    pc2=P(n)-P(n+1);
    pcc=pc1.*pc2;
    if pcc>=0;
        N=[N,n];%把符合条件的峰谷像素值放进数组N中
    end  
end
s=length(N);%计算数组N的长度
%N1=N(2:s);% 寻访从第2个元素到最后一个元素，并将其放入数组N1中
px=0; %存放不独立的峰谷位置标志
for i=2:s
    xp=N(i)-N(i-1);
    pp=P(N(i))-P(N(i-1));
    if (xp==1)&(abs(pp)==0) %该峰谷位置比前一个大一个像素，且灰度值相等，则该峰（谷）位置不独立
        px(i)=1; %不独立的峰谷位置标志为1
    elseif (abs(pp)<=2)&((P(N(i))<dp)|(P(N(i))>up))%灰度值大于up或小于dp的位置，应为峰谷区，此时灰度值有波动或峰谷位置不相邻也可能是不独立的峰谷;(灰度值处于up和dp之间的波动应视为峰谷波动）
         px(i)=1;
    else
        px(i)=0;
    end
end
for j=1:s-1
    dpx(j)=px(j)-px(j+1);
end
apx=find(dpx==-1);
bpx=find(dpx==1);
lapx=length(apx);
lbpx=length(bpx);
lpx=min(lapx,lbpx);
N2=[];
tp=0;
for ii=1:s-1
    if (px(ii)==0)&(dpx(ii)==0)
       N2=[N2,N(ii)];
    elseif (px(ii)==0)&(dpx(ii)==-1)&(tp<=lpx)
       uua=N(ii);%不独立峰谷位置的起始
    elseif (px(ii)==1)&(dpx(ii)==1)&(tp<=lpx)
       uub=N(ii);
       tp=tp+1; %不独立峰谷位置的结束
       N0=fix((uua+uub)./2);
       N2=[N2,N0];
    end
end

lg=length(N2);%N2存储的峰谷位置像素所组成的数组的长度

%下面找各峰谷中心位置坐标
w=P(N2);
M=[];
MM=[];
for i=2:lg-1
    T1=(w(i)+w(i-1))./2;%计算左阈值
    T2=(w(i)+w(i+1))./2;%计算右阈值
    t1=N2(i-1)+1;
    tm1=abs(P(t1)-T1);
    a=t1; %左阈值坐标
    while t1<N2(i)
          if abs(P(t1)-T1)<tm1;
              tm1=abs(P(t1)-T1);
              a=t1;
          else 
              t1=t1+1;
          end
    end 
    t2= N2(i)+1;
    tm2=abs(P(t2)-T2);
    b=t2; %右阈值坐标
    while t2<N2(i+1)
          if abs(P(t2)-T2)<tm2;
              tm2=abs(P(t2)-T2);
              b=t2;
          else 
              t2=t2+1;
          end
    end 
     S=0;
     T=0;
    for k=a:b
        S=S+k.*P(k);
        T=T+P(k);
    end
    m=S./T;%做加权平均的计算，得出中心位置
    MM=[MM,m];
    M=[M,fix(m)];%把中心位置存放入数组M中
end
%M=N2;
%plot(P),hold on
%plot(M,P(M),'*'),hold on 
%plot(N2,P(N2),'ro'),hold on
l=length(M);
ds=0;
% 下面计算条码宽度
for i=2:1000
    ds(i)=P(i)-P(i-1); %信号求导
    %if ds(i)==0
        %ds(i)=1;
    %end
end
cds=find(ds==max(ds)); %导数最大值坐标
for ia=1:l-1
    if MM(ia)<cds(1)&cds(1)<MM(ia+1)
       as=ia;   %边缘左边峰（谷）位置
       bs=ia+1; %边缘右边谷（峰）位置
    end  
end  
S=0;
T=0;
for ks=M(as):M(bs)
        S=S+ks.*ds(ks);
        T=T+ds(ks);
end
e0=S./T; %找条码的独立边缘e0
%计算各边缘位置
E=[];%存储各边缘位置
E(as)=e0;
iq=as-1;
while iq>=1
     ae(iq)=2.*MM(iq+1);
     be(iq)=E(iq+1);
     E(iq)=ae(iq)-be(iq);
     iq=iq-1;
end
for ip=(as+1):l-1
    E(ip)=2.*MM(ip)-E(ip-1);
end 
FE=fix(E);
%计算各条码宽度值
ll=l-1;%（l-1）个边缘
wid=[];
if ds(FE(1))<0
     lc=0
     if rem(ll,2)==0
       la=ll./2;
     else 
       la=(ll-1)./2;
     end
else
     lc=1;
     if rem((ll-1),2)==0
       la=(ll-1)./2;
     else 
       la=(ll-2)./2;
     end
end
   
for lb=1:la
    lbb=2.*lb;
    wid(lb)=E(lbb+lc)-E(lbb-1+lc); %存储各条码的宽度值
end
MM=fix(M);
%subplot(2,1,1)
subplot(2,1,1),plot(P),hold on
plot(FE,P(FE),'ro')
subplot(2,1,2), plot(P),hold on
plot(M,P(M),'black*'),hold on
%plot(N2,P(N2),'r*'),hold on