请教大神个问题。

waboer

clear all;
" Q1 C0 w" C7 W) D5 \I_Expan=imread('C:\Users\a\Desktop\112.png')
. x1 |! p! S+ W- Y+ P" g5 vimshow(I_Expan);
P=20;      %定义搜索范围
Q=75;
      %%  利用基元法检测卷缩轮轮廓
for i=1: P
) ~2 c) e, R1 }/ q5 ~8 J5 P    for j=3: Q;s=0;t=0;
q(1)=abs(I_Expan(i,j)-I_Expan(i,j-2));
q(2)=abs(I_Expan(i,j)-I_Expan(i,j-1));
4 `+ H3 k: n, H* }4 T4 Iq(3)=abs(I_Expan(i,j)-I_Expan (i,j+1));
q(4)=abs(I_Expan(i,j)-I_Expan(i,j+2));
for i=1:4
    if q(i)>9 p(i)=0;
  i. Y( P1 C4 }! r, `7 T  v    else if q(i)<4 p(i)=2;
  f3 Z0 k* X, [7 X' ~3 [; D        else  p(i)=1;         
        end
    end
end
O={p(1),p(2),p(3),p(4)};
if  O=={0,0,0,0}|{0,0,0,1}|{0,0,1,0}|{0,1,0,0}|{1,0,0,0}
        s=s+1;
    else if O=={2,2,2,2}|{2,2,2,1}|{2,2,1,2}|{2,1,2,2}|{1,2,2,2}
            t=t+1;        
, m  L6 \( M6 W3 d1 \        end
* T3 u4 {- i/ f- d( wend
" ]& v, i$ F; H% w4 Z" u% EV1=s/3^4;
V2=t/3^4;
V(j)=V1/V2;
    end
    for v=3: Q
    M(v)==max(V(v));         
# |" v. V# L  ^& p& E; L( Q    end
    for i=3: Q                                    
        [m n]=find(I_Expan(:,i)==M(v));            %最大值点所在坐标      
       [p,q]=size(m);
3 Y7 z0 F" D% @; p# s2 {# u& y- D    end
for j=1:p                                 
             if((m(j)>=P)&(m(j)<=Q))
  G: Q% h7 c6 q( `2 m3 H          a=m(j);
         end
9 G# |7 ^  c- X  }    end
     I_Expan(a,i)=255;                     %将极值点所在像素置为 255
/ q  E( d  t2 J9 \5 g7 {end                        
+ \; G: B& ]5 w+ g7 B5 E1 k1 O# |
  T; `. l% `) _8 U8 S$ }     I_grad=I_Expan;
1 w+ O3 B: V6 s1 Z9 I7 Y5 x! [' ufigure,imshow(I_grad);title('定位图像');


3 k6 c) Y: I5 s5 e要运行的图片

运行结果 出错的地方

麻烦大神帮忙看看，非常感谢！


) _7 z0 p3 d: x  w: U

Colbie

逻辑运算符没这样的比对功能，两种方法
; V/ b! ]" e- j% O- i4 L1 j3 g1. 自己写个子函数，实现你想要的胞组元素比对
2 {3 L6 w) e+ Y* h2. 简的变型一下，根据你提供的信息，可能考虑按以下这样来
, }6 y2 |2 i4 R$ XO=p(1)*1000+p(2)*100+p(3)*10+p(4);
# ^# ^) L" o. @* T. o+ fswicth O
case {0,1,10,100,1000}
3 V8 O0 r- y) A3 B: v  |( [  s=s+1;
" a/ Q# D. {4 z  ycase {2222,2221,2212,2122,1222}
   t=t+1;
- h, ~) P! e" qend