インデックスが範囲外のmatlab

Question

私はそれを実行する方法を得ることができません、誰でもこの助けてください。

私はこのエラーが発生します：ライン212と213で

K>> [minutiae_valid]=postprocess_TICO(B,im,~K) 
??? Attempted to access wind(273,386); index out of bounds because size(wind)=[368,385]. 

Error in ==> postprocess_TICO at 214 
if (thinned(p,q)==1)&&(wind(p,q)==0) 

postprocess_Tico

function [minutiae_valid]=postprocess_TICO(minutiae_valid,im1,thinned) 
% clc; 
% close all; 
% clear all; 
cell=16;%total no. of cells needed in quantization; 
% im1= imread('E:\FINGERPRINT DATABASE\FVC2002\Dbs\Db1_a\16_6.tif'); 
wl=(33-1)/2; 
% no=20; 
% no2=20; 
% nox = 0; 
% noy = 0; %to choose border area 
% 
m=size(im1,1); 
n=size(im1,2); 
%m=(round(m/cell))*cell; 
%n=(round(n/cell))*cell; 
% im1=imresize(im1,[m n]); 
% im=double(im1); 
% figure; 
% subplot(1,2,1); 
% imshow(uint8(im)); 
% title('original input image'); 
% 
% %Identify ridge-like regions and normalise image-------------------------- 
% blksze = 10; 
% thresh =30; 
% M0=100; 
% STD0=150; 
% [thinim, mask] = imgenhance(im, blksze, M0, STD0); 
% subplot(1,2,2); subimage((thinim)); 
% title('Improved thinned input image', 'FontSize', 12); 
% INPUT_IMG=thinim; 
% 
% %minutiae extraction 
% thinned= INPUT_IMG; 
% img = INPUT_IMG; 
% [minutiae, minutiae_img, combined] = findminutiae(thinned, img); 
% totminu_I = size(minutiae,1) 
% figure, 
% subplot(1,2,1), subimage(combined), title('minutiae before postprocess.') 
% %______________________________________________________________ 
% %TO FIND BOUNDARY OF INPUT IMAGE AS WELL AS REGION OF FALSE MINU 
% mask_input=mask; 
% border=uint8(zeros(m,n)); 
% 
% [m,n]=size(mask); 
%  validr=uint8(zeros(m,n)); 
%  thin_temp=thinned; 
%  for i=1:m 
%   for j=2:n-1 
%    if((mask_input(i,j-1)==0)&&(mask_input(i,j)==1)&&(j<=n-no)) 
%      border(i,j)=1; 
%     validr(i,j-nox:j+no)=1; 
%      thin_temp(i,j-nox:j+no)=1; 
%    end 
%     
%    if ((mask_input(i,j+1)==0)&&(mask_input(i,j)==1)&&(j>no)) 
%      border(i,j)=1; 
%     validr(i,j+nox:-1:j-no)=1; 
%     thin_temp(i,j+nox:-1:j-no)=1; 
%    end 
%     
%   end 
%  end 
% 
%   for i=1:m 
%   for j=1:n-1:n 
%    if (mask_input(i,j)==1) 
%     if (j==1) 
%      border(i,j)=1; 
%     validr(i,j:j+no)=1; 
%     thin_temp(i,j:j+no)=1; 
%     end 
%     if (j==n) 
%      border(i,j)=1; 
%     validr(i,j:-1:j-no)=1; 
%     thin_temp(i,j:-1:j-no)=1; 
%     end 
%    end 
%   end 
%   end 
%  
%   
%   
% for j=1:n 
%   for i=2:m-1 
%    if((mask_input(i-1,j)==0)&&(mask_input(i,j)==1)&&(i<=m-no)) 
%      border(i,j)=1; 
%     validr(i-noy:i+no,j)=1; 
%      thin_temp(i-noy:i+no,j)=1; 
%    end 
%    if((mask_input(i+1,j)==0)&&(mask_input(i,j)==1)&&(i>no)) 
%      border(i,j)=1; 
%      validr(i+noy:-1:i-no,j)=1; 
%      thin_temp(i+noy:-1:i-no,j)=1; 
%    end 
%     
%   end 
% end 
% 
% for j=1:n 
%   for i=1:m-1:m 
%    if (mask_input(i,j)==1) 
%    if (i==1) 
%      border(i,j)=1; 
%     validr(i:i+no,j)=1; 
%     thin_temp(i:i+no,j)=1; 
%    else validr(i:-1:i-no,j)=1; 
%     thin_temp(i:-1:i-no,j)=1; 
%    end 
%    end 
%   end 
% end 
% 
% %___________________________________________________________________ 
% % removing minu at the border of roi 
% for x=1:size(thinned,1) 
% for y=1:size(thinned,2) 
%   if ((validr(x,y)==1)&&(thinned(x,y)==1)) 
%    combined(x,y,:)=[255,255,255]; 
%    minutiae_img(x,y,:)=[0,0,0]; 
%    
%   end 
% end 
% end 
%  
% for i=1:totminu_I 
%  x=minutiae(i,1); 
%  y=minutiae(i,2); 
%  if ((validr(x,y)==1)&&(thinned(x,y)==1)) 
%   minutiae(i,:)=0; 
%  end 
% end 
% %__________________________ 
% 
% subplot(1,2,2), subimage(combined), title('after boundary effect ') 
% minu_count=1; 
% minutiae_valid(minu_count, :) = [0,0,0,0]; %to make x, y, CN, theta value of minutiae not in border region. 
% for i=1:totminu_I 
%  CN=minutiae(i,3); 
%  if (CN~=0) 
%   minutiae_valid(minu_count, :) = minutiae(i,:); 
%   minu_count = minu_count + 1; 
%  end 
% end 

% totminu_I= minu_count-1; 


k=size(minutiae_valid); 
totminu_I=k(1,1); 
%FALSE RIDGE BIFURCATION 
for i=1:totminu_I 
if minutiae_valid(i,3)==3 
      wind=ones(m,n); 
     xi= minutiae_valid(i,1); 
     yi=minutiae_valid(i,2); 
     x=xi; 
     y=yi; 
     wind(xi-wl:xi+wl,yi-wl:yi+wl)=0; %window of size 25 cross 25 
     wind(xi,yi)=-1; 

     %mark '1', '2', '3'. 
     value=0; 

     r=x-1; 
     for c=y-1:y+1 
      if (thinned(r,c)==1)&&(wind(r,c)==0) 
       value=value+1; 
       wind(r,c)=value; 
      end 
     end 

     c=y+1; 
     for r=x-1:x+1 
      if (thinned(r,c)==1)&&(wind(r,c)==0) 
       value=value+1; 
       wind(r,c)=value; 
      end 
     end 

     r=x+1; 
     for c=y+1:-1:y-1 
      if (thinned(r,c)==1)&&(wind(r,c)==0) 
       value=value+1; 
       wind(r,c)=value; 
      end 
     end 

     c=y-1; 
     for r=x+1:-1:x-1 
      if (thinned(r,c)==1)&&(wind(r,c)==0) 
       value=value+1; 
       wind(r,c)=value; 
      end 
     end 


for value=1:3 %search '1'/'2/'3' value pixel arond the minutia point in the window. 
for p=x-1:x+1 
    for q=y-1:y+1 
     if wind(p,q)==value 
      xi=p; 
      yi=q; 
     end 
    end 
end 

cnt=0; 
     while(xi>x-wl)&&(xi<x+wl)&&(yi>y-wl)&&(yi<y+wl) %mark '1'/'2'/'3' to the adjacent connected pixels 
      k=0; 
      for p=xi-1:xi+1 
       for q=yi-1:yi+1 
        if (thinned(p,q)==1)&&(wind(p,q)==0) 
         wind(p,q)=value; 
          k=k+1; 
         if k==1 
          x1=p; 
          y1=q; 
         end 
         if k==2 %more than one connected pixel 
          cnt=cnt+1; 
          x2=p; 
          y2=q; 

          t1=x1; 
          x1=x2; 
          x2=t1; 

          t2=y1; 
          y1=y2; 
          y2=t2; 
          %___________ 
          xi=x2; 
          yi=y2; 

         while(xi>x-wl)&&(xi<x+wl)&&(yi>y-wl)&&(yi<y+wl) 
            k=0; 
            flag_see=0; 
          for p=xi-1:xi+1 

           for q=yi-1:yi+1 

            if (thinned(p,q)==1)&&(wind(p,q)==0) 
              wind(p,q)=value; 
              flag_see=1; 
              x3=p; 
              y3=q; 
              k=k+1;    
            end 
           end 
          end 
          if(flag_see==1) 
          xi=x3; 
          yi=y3; 
          else 
           xi=p; 
           yi=q; 

          end 

           if k==0 
           break 
           end 
          end 
%   xi=x1; 
%   yi=y1; 

         end 
        end 
       end 
      end 
      xi=x1; 
      yi=y1; 
      if k==0 
       break 
      end 
     end 

end 

% figure,subimage(wind); 
%to count 0-1, 0-2, 0-3 transition around the boundary clkwise 
T=1; 
for v=1:3 
     T0v=0; 
     r=x-wl; 
     for c=y-wl:y+(wl-1) 
      if (wind(r,c)==0)&&(wind(r,c+1)==v) 
       T0v=T0v+1; 
      end 
     end 

     c=y+wl; 
     for r=x-wl:x+(wl-1) 
      if (wind(r,c)==0)&&(wind(r+1,c)==v) 
       T0v=T0v+1; 
      end 
     end 

     r=x+wl; 
     for c=y+wl:-1:y-(wl-1) 
      if (wind(r,c)==0)&&(wind(r,c-1)==v) 
       T0v=T0v+1; 
      end 
     end 

     c=y-wl; 
     for r=x+wl:-1:x-(wl-1) 
      if (wind(r,c)==0)&&(wind(r-1,c)==v) 
       T0v=T0v+1; 
      end 
     end 

T1=T0v==1; 
    T=T & T1; 
end 
     if T~=1 
      minutiae_valid(i,:)=[0 0 0 0]; 
      combined(x,y,:)=[255,255,255]; 
     end 

end 
end 


minu_count=1; 
minutiae_valid_final(minu_count, :) = [0,0,0,0]; %to make x, y, CN, theta value of minutiae not in border region. 
for i=1:totminu_I 
    CN = minutiae_valid(i,3); 
    if (CN~=0) 
     minutiae_valid_final(minu_count, :) = minutiae_valid(i,:); 
      minu_count = minu_count + 1; 
    end 
end 

% totminu_I= minu_count-1; 
minutiae_valid = minutiae_valid_final; 

% figure, subimage(combined); 

Answer 1

あなたは風の9つの要素をループは、条件のいくつかの並べ替えを一致させる要素を探していることをネストされたループを作成します。見つけたら、9つの要素ブロックの中心を変更します。

ブロックの中心が風の端にあるときに問題が発生します。エッジの中央にある9つの要素を見ると、同じ方法でインデックスにアクセスしてはいけません。さもなければ、上のようにエラーが発生します。

a = rand(4,4); 
for p=2:4 
    for q = 1:3 
    tmp += a(p,q) 
    end 
end 

は、（3,2）を中心とする（正方形の配列の）9要素を合計します。配列のこの要素が存在しないため、

a = rand(4,4); 
for p=3:5 
    for q = 1:3 
    tmp += a(p,q) 
    end 
end 

は、それが（5,1）にアクセスしようとしたとき、あなたが持っている同じエラーが発生します。これを修正する正しい方法は、あなたがしようとしていることに依存しますが、それはあなたの質問の範囲を超えていることがわかります。

こちらがお役に立てば幸いです。

いくつかのコメントを耳にして、私はあなたにアドバイスをしたいと思います。将来的には、可能な限り最短のテストケースまで質問を沸かせることができれば、質の高い回答を得ることができます。これは、独自のコードをデバッグするための貴重な一歩です。テストケースを構築する方法を理解することは、しばしばあなたの問題に対する答えを明らかにします。明らかでないときは、Stackoverflowは優れたリソースです。