www.gusucode.com > 关于空间信息的处理，能较好的对空间信息进行窗口处理 > 关于空间信息的处理，能较好的对空间信息进行窗口处理/window9/hsi_read.m

    function [C,T,B,B17] = hsi_read
load Indian_pines_gt
im3=indian_pines_gt;
[im4]=window9;

%%%%%训练样本 1(46) 2(1428) 3(830) 4(237) 5(483) 6(730) 7(28) 8(478) 9(20) 10(972) 11(2455) 12(593) 13(205) 14(1265) 15(386) 16(93) 0(10776) 
%%%%%%%%%%%%%%%%%
a=1;b=1;c=1;d=1;e=1;f=1;g=1;h=1;r=1;s=1;k=1;l=1;m=1;n=1;o=1;p=1;q=1;
for i=1:145
    for j=1:145
        if im3(i,j)<1;
         %%%%%%%%%%%%%已知每一类别的总个数这里假设为k 第16类（k=93）
       A17(:,q)=im4(i,j,:);
       B17(q,:)=[i j];
        q=q+1;
        else if im3(i,j)==1;
        %%%%%%%%%%%%%已知每一类别的总个数这里假设为k 第一类（k=46）
       A1(:,a)=im4(i,j,:);
       B1(a,:)=[i j];
       a=a+1; 
%figure(2);imshow(imb,[])
   else if im3(i,j)==2;
        %%%%%%%%%%%%%已知每一类别的总个数这里假设为k 第2类（k=1428）
       A2(:,b)=im4(i,j,:);
       B2(b,:)=[i j];
       b=b+1;
%figure(2);imshow(imb,[])
   else if im3(i,j)==3;
          %%%%%%%%%%%%%已知每一类别的总个数这里假设为k 第3类（k=830）
       A3(:,c)=im4(i,j,:);
       B3(c,:)=[i j];
       c=c+1;
%figure(2);imshow(imb,[])
   else if im3(i,j)==4;
       %%%%%%%%%%%%%已知每一类别的总个数这里假设为k 第4类（k=237）
       A4(:,d)=im4(i,j,:);
       B4(d,:)=[i j];
       d=d+1;     
%figure(2);imshow(imb,[])
   else if im3(i,j)==5;
        %%%%%%%%%%%%%已知每一类别的总个数这里假设为k 第5类（k=483）
       A5(:,e)=im4(i,j,:);
       B5(e,:)=[i j];
       e=e+1;       
%figure(2);imshow(imb,[])
       else if im3(i,j)==6;
        %%%%%%%%%%%%%已知每一类别的总个数这里假设为k 第6类（k=730）
       A6(:,f)=im4(i,j,:);
       B6(f,:)=[i j];
      f=f+1;
%figure(2);imshow(imb,[])
     else if im3(i,j)==7;
         %%%%%%%%%%%%%已知每一类别的总个数这里假设为k 第7类（k=28）
       A7(:,g)=im4(i,j,:);
       B7(g,:)=[i j];
       g=g+1; 
%figure(2);imshow(imb,[])

         else if im3(i,j)==8;
       %%%%%%%%%%%%%已知每一类别的总个数这里假设为k 第8类（k=478）
       A8(:,h)=im4(i,j,:);
       B8(h,:)=[i j];
       h=h+1; 
%figure(2);imshow(imb,[])
 
             else if im3(i,j)==9;
        %%%%%%%%%%%%%已知每一类别的总个数这里假设为k 第9类（k=20）
       A9(:,r)=im4(i,j,:);
       B9(r,:)=[i j];
       r=r+1;
%figure(2);imshow(imb,[])
  
                 else if im3(i,j)==10;
          %%%%%%%%%%%%%已知每一类别的总个数这里假设为k 第10类（k=972）
       A10(:,s)=im4(i,j,:);
       B10(s,:)=[i j];
       s=s+1;  
%figure(2);imshow(imb,[])
  
                     else if im3(i,j)==11;
         %%%%%%%%%%%%%已知每一类别的总个数这里假设为k 第11类（k=2455）
       A11(:,k)=im4(i,j,:);
       B11(k,:)=[i j];
       k=k+1;
%figure(2);imshow(imb,[])
 
                         else if im3(i,j)==12;
         %%%%%%%%%%%%%已知每一类别的总个数这里假设为k 第12类（k=205）
       A12(:,l)=im4(i,j,:);
       B12(l,:)=[i j];
       l=l+1; 
%figure(2);imshow(imb,[])
  
                             else if im3(i,j)==13;
         %%%%%%%%%%%%%已知每一类别的总个数这里假设为k 第13类（k=205）
       A13(:,m)=im4(i,j,:);
       B13(m,:)=[i j];
       m=m+1; 
%figure(2);imshow(imb,[])
 
                                 else if im3(i,j)==14;
          %%%%%%%%%%%%%已知每一类别的总个数这里假设为k 第14类（k=1265）
       A14(:,n)=im4(i,j,:);
       B14(n,:)=[i j];
      n=n+1;  
%figure(2);imshow(imb,[])
 
                                     else if im3(i,j)==15;
          %%%%%%%%%%%%%已知每一类别的总个数这里假设为k 第15类（k=386）
       A15(:,o)=im4(i,j,:);
       B15(o,:)=[i j];
       o=o+1;
%figure(2);imshow(imb,[])
  
                                         else if im3(i,j)==16;
         %%%%%%%%%%%%%已知每一类别的总个数这里假设为k 第16类（k=93）
       A16(:,p)=im4(i,j,:);
       B16(p,:)=[i j];
        p=p+1; 
%figure(2);imshow(imb,[])
          
                                         end
                                             end
                                         end
                                     end
                                 end
                             end
                         end
                     end
                 end
             end
         end
           end
       end
       end
       end
       end
   end
    end
end
T=[a-1 b-1 c-1 d-1 e-1 f-1 g-1 h-1 r-1 s-1 k-1 l-1 m-1 n-1 o-1 p-1];
C=[A1';A2';A3';A4';A5';A6';A7';A8';A9';A10';A11';A12';A13';A14';A15';A16'];
B=[B1;B2;B3;B4;B5;B6;B7;B8;B9;B10;B11;B12;B13;B14;B15;B16];