edge_detection.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存储的峰谷位置像素所组成的数组的长度
E=[];%存储各边缘位置
for ii=2:lg
    clear mds;
    clear ds;
    mm=(N2(ii-1)+1):(N2(ii)-1);
    for j=(N2(ii-1)+1):(N2(ii)-1)
        ds(j)=P(j)-P(j-1);
    end
    ads=abs(ds);
    mds=find(ads==max(ads));  
    f=fix(mds);
    E=[E,f];
end
subplot(2,1,1),plot(P),hold on
plot(E,P(E),'ro')
subplot(2,1,2), plot(P),hold on
plot(N2,P(N2),'black*'),hold on