[求助] 一段图像处理的代码和应用

yg=gradient(yrpe);
ygg=ones([1 length(yrpe)]); ygg(abs(yg)>20)=0; ygl=bwlabel(ygg);
figure; imshow(mat2gray(Lbinrpe*0.5+Lmed));hold on;
palett=jet(max(ygl));
for iiih=1:max(ygl( :))
plot(x(ygl==iiih), yrpe(ygl==iiih),'Color',palett(iiih,:),'LineWidth',4);
end
pam_dl=[];
figure; imshow(mat2gray(Lbinrpe*0.5+Lmed)); hold on
for iiik=1:max(ygl(:))
for iiikk=iiik:max(ygl(:))
if iiik<=iiikk
ygk=[yrpe(ygl==iiik),yrpe(ygl==iiikk)];
xgk=[x(ygl==iiik);x(ygl==iiikk)];
else
ygk=[yrpe(ygl==iiikk),yrpe(ygl==iiik)];
xgk=[x(ygl==iiikk);x(ygl==iiik)];
end
if length(ygk)>10
P = POLYFIT(xgk',ygk,2); yrpes =
round(POLYVAL(P,x));
plot(yrpes,'g*-')
pam_dl=[pam_dl;[iiik iiikk sum(abs(yrpeyrpes')<
20)]];
end
end
end