请教大神个问题。

waboer

clear all;
2 D7 X  s% O, r5 A" G& h) |5 K2 II_Expan=imread('C:\Users\a\Desktop\112.png')
imshow(I_Expan);
P=20;      %定义搜索范围
Q=75;
) U; B8 V( E. u7 v/ s; q      %%  利用基元法检测卷缩轮轮廓
for i=1: P
    for j=3: Q;s=0;t=0;
q(1)=abs(I_Expan(i,j)-I_Expan(i,j-2));
1 e3 i4 s" O+ x2 T5 q+ p- M1 ~q(2)=abs(I_Expan(i,j)-I_Expan(i,j-1));
q(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
3 y" A! O# M; y( u' \# W    if q(i)>9 p(i)=0;
" k8 {* R! u: X# w) i% a3 H    else if q(i)<4 p(i)=2;
        else  p(i)=1;         
        end
. X' o3 }6 `7 J    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;
# E# [7 h  Z* k: z) ?: W; e    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;        
        end
end
V1=s/3^4;
V2=t/3^4;
3 ~1 Y: I1 Y# v: E8 yV(j)=V1/V2;
" `! b$ r1 C- o    end
9 P! W' V7 F6 |0 {( R. f8 y    for v=3: Q
; A" u  F! L" V7 D' r; U    M(v)==max(V(v));         
    end
/ v. T) R. W% J6 v9 ?/ C: f4 }5 \0 Q: e    for i=3: Q                                    
        [m n]=find(I_Expan(:,i)==M(v));            %最大值点所在坐标      
. U: u( V9 [/ h       [p,q]=size(m);
    end
for j=1:p                                 
             if((m(j)>=P)&(m(j)<=Q))
, q* |: z9 b) i; N6 e, O! R5 D          a=m(j);
         end
* H5 A1 m" I2 P+ B; }2 N    end
7 p4 m. O5 J, Y2 N) [) G. @     I_Expan(a,i)=255;                     %将极值点所在像素置为 255
* T9 o! j3 s, G6 l, e  `3 m. fend                        
, h* q, X  r  Z) n6 A
     I_grad=I_Expan;
figure,imshow(I_grad);title('定位图像');
  o! ~7 \2 F4 ^( \8 ^% p

; E4 v/ N3 N0 Z要运行的图片
. m5 E. c; v- c1 ]$ n3 W
0 V, B2 z3 z, V运行结果 出错的地方
8 A, J$ Y7 f0 r8 a, a4 P: C9 s9 A
$ V& u: G8 P6 [2 z, l' F麻烦大神帮忙看看，非常感谢！
) v. q( ?1 ^5 ~) y
9 P/ c% V) a9 |- V
& H1 e: x6 w: P& d% x6 ]! W/ r

Colbie

" u( X) N" S: i9 {' h! p逻辑运算符没这样的比对功能，两种方法
0 X! d3 i, j9 `9 Q1. 自己写个子函数，实现你想要的胞组元素比对
2. 简的变型一下，根据你提供的信息，可能考虑按以下这样来
O=p(1)*1000+p(2)*100+p(3)*10+p(4);
7 A, o& A" Q: n$ u7 w# aswicth O
2 y( h- _6 w1 K1 G5 acase {0,1,10,100,1000}
  s=s+1;
% H0 o! e4 d3 \! C, @" I5 r- \! S- f3 ^case {2222,2221,2212,2122,1222}
3 S, i, F2 c5 Y9 g3 T8 T   t=t+1;
end