edge.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;