📄 test1.m
字号:
clear;
[P1]=bargene();%调用原始条码
P=P1(401:500);
N=-1;
for n=2:99
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:end);% 寻访从第2个元素到最后一个元素,并将其放入数组N1中
N2=N1(1);%先把N1的第1个元素放入N2,再对后面的进行筛选
for i=2:s-1
if (N1(i)-N1(i-1))~=1 %排除同一峰谷中多的相临像素点,仅保留一个
N2=[N2,N1(i)]; %从N1的第2个元素开始筛选,存储符合条件的峰谷坐标位置(象素)
end
end
%plot(P)
%hold on;
%plot(N2,P(N2),'black*')
l=length(N2);%N2存储的峰谷位置像素所组成的数组的长度
E=[];%存储各边缘位置
for ii=2:l
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
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -