Matlab讨论区 - 声振论坛 - 振动,动力学,声学,信号处理,故障诊断

eight 发表于 2007-4-25 20:17

原帖由程序255 于 2007-4-25 11:28 发表 http://forum.vibunion.com/forum/images/common/back.gif

利用形态学的知识做吧

eight 发表于 2007-4-25 20:18

原帖由程序255 于 2007-4-25 20:05 发表 http://forum.vibunion.com/forum/images/common/back.gif
Q(1:8)=(p2,p3,p4,p5,p6,p7,p8,p9);
Ｑ为一个一唯矩阵　后面是它的元素的值　　

Q(1:8)=;

程序255 发表于 2007-4-25 23:47

eight 看看我的这个表示对吗

比较R 和P 的元素是否完全一样
和把P用R代替

while R~=P
P=R;

eight 发表于 2007-4-26 10:02

原帖由程序255 于 2007-4-25 23:47 发表 http://forum.vibunion.com/forum/images/common/back.gif
比较R 和P 的元素是否完全一样
和把P用R代替

while R~=P
P=R;

help isequal

代替问题要两者维数相同才行

程序255 发表于 2007-4-26 10:42

eight 急死了

=size(P);
W=zeros(m,n);

for i=2:m-1
for j=2:n-1
count2=0;count3=0;count4=0;count5=0;count6=0;count7=0;count8=0;count9=0;
p1=P(i,j);
p2=P(i-1,j);
p3=P(i-1,j+1);
p4=P(i,j+1);
p5=P(i+1,j+1);
p6=P(i+1,j);
p7=P(i+1,j-1);
p8=P(i,j-1);
p9=P(i-1,j-1);
   if p1==255


   if p2==255
      count2=1;
   else count2=0;
   end;
   if p3==255
      count3=1;
   else count3=0;
   end;
   if p4==255
      count4=1;
   else count4=0;
   end;
   if p5==255
      count5=1;
   else count5=0;
   end;
   if p6==255
      count6=1;
   else count6=0;
   end;
   if p7==255
      count7=1;
   else count7=0;
   end;
   if p8==255
      count8=1;
   else count8=0;
   end;
   if p9==255
      count9=1;
   else count9=0;
   end;
   else W(i,j)=P(i,j);
end;
count=count2+count3+count4+count5+count6+count7+count8+count9;
if count<2|count>6
W(i,j)=P(i,j);
else
Q(1:8)=;
a=0;
for k=1:7
   if Q(k)<Q(k+1)
      a=a+1;
   end;
end;
if a~=1
   W(i,j)=P(i,j);
else
   if p2*p4*p6~=0|p4*p6*p8~=0

W(i,j)=P(i,j);
   else W(i,j)=0;
   end;
end;
end;
end;
end;
=size(W);
R=zeros(m,n);

for i=2:m-1
for j=2:n-1
p1=W(i,j);
p2=W(i-1,j);
p3=W(i-1,j+1);
p4=W(i,j+1);
p5=W(i+1,j+1);
p6=W(i+1,j);
p7=W(i+1,j-1);
p8=W(i,j-1);
p9=W(i-1,j-1);
if p1==255


   if p2==255
      count2=1;
   else count2=0;
   end;
   if p3==255
      count3=1;
   else count3=0;
   end;
   if p4==255
      count4=1;
   else count4=0;
   end;
   if p5==255
      count5=1;
   else count5=0;
   end;
   if p6==255
      count6=1;
   else count6=0;
   end;
   if p7==255
      count7=1;
   else count7=0;
   end;
   if p8==255
      count8=1;
   else count8=0;
   end;
   if p9==255
      count9=1;
   else count9=0;
   end;
else R(i,j)=W(i,j);
end;

count=count2+count3+count4+count5+count6+count7+count8+count9;
if count<2|count>6
R(i,j)=W(i,j);
else
Q(1:8)=;
a=0;
for k=1:7
   if Q(k)<Q(k+1)
      a=a+1;
   end;
end;
if a~=1
   R(i,j)=W(i,j);
else
   if p2*p4*p8~=0|p2*p6*p8~=0

R(i,j)=W(i,j);
   else R(i,j)=0;
   end;
end;
end;
end;
end;
while R~=P
P=R;
=size(P);

for i=2:m-1
for j=2:n-1
count2=0;count3=0;count4=0;count5=0;count6=0;count7=0;count8=0;count9=0;
p1=P(i,j);
p2=P(i-1,j);
p3=P(i-1,j+1);
p4=P(i,j+1);
p5=P(i+1,j+1);
p6=P(i+1,j);
p7=P(i+1,j-1);
p8=P(i,j-1);
p9=P(i-1,j-1);
   if p1==255


   if p2==255
      count2=1;
   else count2=0;
   end;
   if p3==255
      count3=1;
   else count3=0;
   end;
   if p4==255
      count4=1;
   else count4=0;
   end;
   if p5==255
      count5=1;
   else count5=0;
   end;
   if p6==255
      count6=1;
   else count6=0;
   end;
   if p7==255
      count7=1;
   else count7=0;
   end;
   if p8==255
      count8=1;
   else count8=0;
   end;
   if p9==255
      count9=1;
   else count9=0;
   end;
   else W(i,j)=P(i,j);
end;
count=count2+count3+count4+count5+count6+count7+count8+count9;
if count<2|count>6
W(i,j)=P(i,j);
else
Q(1:8)=;
a=0;
for k=1:7
   if Q(k)<Q(k+1)
      a=a+1;
   end;
end;
if a~=1
   W(i,j)=P(i,j);
else
   if p2*p4*p6~=0|p4*p6*p8~=0

W(i,j)=P(i,j);
   else W(i,j)=0;
   end;
end;
end;
end;
end;
=size(W);

for i=2:m-1
for j=2:n-1
p1=W(i,j);
p2=W(i-1,j);
p3=W(i-1,j+1);
p4=W(i,j+1);
p5=W(i+1,j+1);
p6=W(i+1,j);
p7=W(i+1,j-1);
p8=W(i,j-1);
p9=W(i-1,j-1);
if p1==255


   if p2==255
      count2=1;
   else count2=0;
   end;
   if p3==255
      count3=1;
   else count3=0;
   end;
   if p4==255
      count4=1;
   else count4=0;
   end;
   if p5==255
      count5=1;
   else count5=0;
   end;
   if p6==255
      count6=1;
   else count6=0;
   end;
   if p7==255
      count7=1;
   else count7=0;
   end;
   if p8==255
      count8=1;
   else count8=0;
   end;
   if p9==255
      count9=1;
   else count9=0;
   end;
else R(i,j)=W(i,j);
end;

count=count2+count3+count4+count5+count6+count7+count8+count9;
if count<2|count>6
R(i,j)=W(i,j);
else
Q(1:8)=;
a=0;
for k=1:7
   if Q(k)<Q(k+1)
      a=a+1;
   end;
end;
if a~=1
   R(i,j)=W(i,j);
else
   if p2*p4*p8~=0|p2*p6*p8~=0

R(i,j)=W(i,j);
   else R(i,j)=0;
   end;
end;
end;
end;
end;
end;
figure,imshow(R);

我的那个循环为什么不执行啊我估计问题出在条件的表示上可是弄不好上面程序是对P进行边缘细化　，需要重复那个循环　直到不再有变化为止　　关键给我看看怎么让那个循环执行　　:'( 谢了

程序255 发表于 2007-4-26 10:50

我手工又细化了一次是这个效果　可那循环的就没再处理

eight 发表于 2007-4-26 10:53

原帖由程序255 于 2007-4-26 10:42 发表 http://forum.vibunion.com/forum/images/common/back.gif
=size(P);
W=zeros(m,n);

for i=2:m-1
for j=2:n-1
count2=0;count3=0;count4=0;count5=0;count6=0;count7=0;count8=0;count9=0;
p1=P(i,j);
p2=P(i-1,j);
p3=P(i-1,j+1);
p4=P(i,j+1);
p5 ...

你至少得说明是哪个循环吧，我可不是天天有空坐在这里等你们来提问的

程序255 发表于 2007-4-26 11:09

这个循环

while R~=P
P=R;

eight 发表于 2007-4-26 11:17

原帖由程序255 于 2007-4-26 11:09 发表 http://forum.vibunion.com/forum/images/common/back.gif
while R~=P
P=R;

自己设置断点调试一下吧，你的程序太复杂，我没有那么多时间细看

你试试把 R~=P 改为 any(R~=P)

[ 本帖最后由 eight 于 2007-4-26 11:20 编辑 ]

程序255 发表于 2007-4-26 11:38

好的　

程序255 发表于 2007-4-26 13:11

大虾看看我的循环为什么只执行一次

K=1;
while K>0

=size(R);
for i=1:m
for j=1:n
   P(i,j)=R(i,j);
end;
end;

for i=2:m-1
for j=2:n-1
count2=0;count3=0;count4=0;count5=0;count6=0;count7=0;count8=0;count9=0;
p1=P(i,j);
p2=P(i-1,j);
p3=P(i-1,j+1);
p4=P(i,j+1);
p5=P(i+1,j+1);
p6=P(i+1,j);
p7=P(i+1,j-1);
p8=P(i,j-1);
p9=P(i-1,j-1);
   if p1==255


   if p2==255
      count2=1;
   else count2=0;
   end;
   if p3==255
      count3=1;
   else count3=0;
   end;
   if p4==255
      count4=1;
   else count4=0;
   end;
   if p5==255
      count5=1;
   else count5=0;
   end;
   if p6==255
      count6=1;
   else count6=0;
   end;
   if p7==255
      count7=1;
   else count7=0;
   end;
   if p8==255
      count8=1;
   else count8=0;
   end;
   if p9==255
      count9=1;
   else count9=0;
   end;
   else W(i,j)=P(i,j);
end;
count=count2+count3+count4+count5+count6+count7+count8+count9;
if count<2|count>6
W(i,j)=P(i,j);
else
Q(1:8)=;
a=0;
for k=1:7
   if Q(k)<Q(k+1)
      a=a+1;
   end;
end;
if a~=1
   W(i,j)=P(i,j);
else
   if p2*p4*p6~=0|p4*p6*p8~=0

W(i,j)=P(i,j);
   else W(i,j)=0;
   end;
end;
end;
end;
end;
=size(W);

for i=2:m-1
for j=2:n-1
p1=W(i,j);
p2=W(i-1,j);
p3=W(i-1,j+1);
p4=W(i,j+1);
p5=W(i+1,j+1);
p6=W(i+1,j);
p7=W(i+1,j-1);
p8=W(i,j-1);
p9=W(i-1,j-1);
if p1==255


   if p2==255
      count2=1;
   else count2=0;
   end;
   if p3==255
      count3=1;
   else count3=0;
   end;
   if p4==255
      count4=1;
   else count4=0;
   end;
   if p5==255
      count5=1;
   else count5=0;
   end;
   if p6==255
      count6=1;
   else count6=0;
   end;
   if p7==255
      count7=1;
   else count7=0;
   end;
   if p8==255
      count8=1;
   else count8=0;
   end;
   if p9==255
      count9=1;
   else count9=0;
   end;
else R(i,j)=W(i,j);
end;

count=count2+count3+count4+count5+count6+count7+count8+count9;
if count<2|count>6
R(i,j)=W(i,j);
else
Q(1:8)=;
a=0;
for k=1:7
   if Q(k)<Q(k+1)
      a=a+1;
   end;
end;
if a~=1
   R(i,j)=W(i,j);
else
   if p2*p4*p8~=0|p2*p6*p8~=0

R(i,j)=W(i,j);
   else R(i,j)=0;
   end;
end;
end;
end;
end;
for i=1:m
for j=1:n
   K=0;
   if R(i,j)~=P(i,j)
         K=K+1;
   else K=0;
   end;
end;
end;

end;
figure,imshow(R);

页: 1 2 [3]

声振论坛's Archiver

eight 看看我的这个表示对吗

eight 急死了

我手工又细化了一次是这个效果 可那循环的就没再处理

这个循环

好的

大虾 看看我的循环为什么只执行一次

我手工又细化了一次是这个效果　可那循环的就没再处理

好的　

大虾看看我的循环为什么只执行一次