EDA365电子论坛网

标题: 请教大神个问题。 [打印本页]

---

作者: waboer    时间: 2020-3-18 15:34
标题: 请教大神个问题。
clear all;
I_Expan=imread('C:\Users\a\Desktop\112.png')
) t" Q- k+ T3 {9 E; \imshow(I_Expan);
4 G. b+ ?, C, P7 F& O$ DP=20;      %定义搜索范围
Q=75;
      %%  利用基元法检测卷缩轮轮廓
for i=1: 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));
* ?" i( U( A: o5 Rq(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;
    else if q(i)<4 p(i)=2;
        else  p(i)=1;         
        end
8 p3 N/ s3 c' N& t7 N& h    end
end
: K" g" o2 X1 ~- @8 c+ eO={p(1),p(2),p(3),p(4)};
) k( X$ O1 g4 t/ n! q( b( {3 T) a6 o7 `+ Cif  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;        
        end
end
. I7 v4 c8 X# C  K% m$ m; SV1=s/3^4;
V2=t/3^4;
V(j)=V1/V2;
    end
8 K, m, _# G* `; p    for v=3: Q
) G, K* W' W! `, B5 x% i0 v    M(v)==max(V(v));         
9 u* w5 `4 Q& E/ z9 R# K    end
    for i=3: Q                                    
$ x# {6 O' k" x- [  O        [m n]=find(I_Expan(:,i)==M(v));            %最大值点所在坐标      
7 ?( H! v4 B1 [+ U       [p,q]=size(m);
    end
for j=1:p                                 
, p: ?' L6 j+ `             if((m(j)>=P)&(m(j)<=Q))
) F2 ?! T' `6 H4 N          a=m(j);
         end
& M& Q; B( ~  z" }    end
+ z5 }6 I9 X& T( V) N/ `     I_Expan(a,i)=255;                     %将极值点所在像素置为 255
end                        
% H. Q+ A9 ]6 K9 y6 n/ ^
- @; Z5 i8 a7 E+ W) F( I     I_grad=I_Expan;
* H8 c! Y" ?: K. \1 d7 `+ B( h% cfigure,imshow(I_grad);title('定位图像');
9 w, E+ I& E9 D
; m1 q* i& [+ G- y( X  g- y3 E6 ]8 V
要运行的图片
- V; R. p3 V$ Z4 B
运行结果 出错的地方

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


# k' o4 e- r6 g7 `

---

作者: Colbie    时间: 2020-3-18 18:23

/ u, K1 ~0 z4 Z9 ?( a" }7 ^/ L逻辑运算符没这样的比对功能，两种方法
1. 自己写个子函数，实现你想要的胞组元素比对
1 A2 L% N2 A8 d2. 简的变型一下，根据你提供的信息，可能考虑按以下这样来
O=p(1)*1000+p(2)*100+p(3)*10+p(4);
swicth O
- b% d0 ^2 U- O1 Scase {0,1,10,100,1000}
  s=s+1;
case {2222,2221,2212,2122,1222}
   t=t+1;
end

---

欢迎光临 EDA365电子论坛网 (https://bbs.eda365.com/)

Powered by Discuz! X3.2