linearm3_temp.m

线性M3网络及其示意图

function Z=LinearM3(X,Y)
%%GZCM3 function
global M

NC1=size(M,1);
NC2=size(M{1},1);


%%% find min distance Point with border by means of Index
minD=pointDist(X,M{i}{1}(1:2));
for i=1:NC1
    pDist(:,:,i)=pointDist(X,Y,M{i}{1}(1:2));
    minD=min(minD,pDist(:,:,i));
end

minI=i*(minD==pDist(:,:,i));
for i=1:NC1
    mI(:,:,i)=i*(minD==pDist(:,:,i));
    minI=max(minI,mI(:,:,i)); 
end

%%%%%%%%%%





I=1;
for i=1:NC1
    
    AllZero=0;% 0 means all M{i}{j}(5) is zero for all j in i
    
    for j=1:NC2
       if M{i}{j}(5)==1 
           AllZero=1;
           break;
       end
    end
    
    if j==NC2 & AllZero==0
        continue;
    end
    
    TMIN=Linear(X,Y,i,j);  %TEMP MIN
      
    for k=j:NC2
        if M{i}{k}(5)==1 
            TT=Linear(X,Y,i,k);
            TMIN=min(TMIN,TT);
        end
    end
    
    
    MIN(:,:,I)=TMIN;
    I=I+1;
end


MAX=MIN(:,:,1);
if (I-1)>=2
  for i=2:(I-1)
    TMAX=MIN(:,:,i);    %TEMP MAX
    MAX=max(MAX,TMAX);
  end
end


Z=MAX;