基于K-means算法的遥感图像分类的matlab实现

基于K-means算法的遥感图像分类的matlab实现，基本没什么用

clear all; 
  
s=imread（‘L1.bmp‘）;
 
s=double（s）;
 
 
X1=20;
X2=45;
X3=90; 
 
num1=0;
num2=0;
num3=0;
[MN]=size（s）;
 
R=zeros（MN）;    
 
a=10;
b=20;
counter=0;      
format long;
 
 
while（abs（b-a）>0.001）   
    counter=counter+1 
    a
    b=a; 
    a=0; 
    for i=1:M 
    for j=1:N 
            
    %-------------------------------------------------------
            p1=（s（ij）-X1）.^2;     
            p2=（s（ij）-X2）.^2; 
            p3=（s（ij）-X3）.^2;       
            min_p=min（[p1p2p3]）;        
            if （p1==min_p）  
                R（ij）=1;        
                num1=num1+1;   
            elseif （p2==min_p） 
                R（ij）=2;       
                num2=num2+1; 
            elseif （p3==min_p） 
                R（ij）=3;        
                num3=num3+1; 
            end        
        a=a+sqrt（min_p）;                     
        end 
    end 
    
    
   X1=X1*0;
   X2=X2*0;
   X3=X3*0;                             
    a=a/M/N; 
    
    %---------------------------------------------------
    for i=1:M 
    for j=1:N 
            if （R（ij）==1） X1=X1+s（ij）; 
            elseif （R（ij）==2） X2=X2+s（ij）; 
            elseif （R（ij）==3） X3=X3+s（ij）;
            end 
        end 
    end 
    
    X1=X1/num1; num1=0;     
    X2=X2/num2; num2=0; 
    X3=X3/num3; num3=0;     
end 
 
 
r=zeros（MN）; 
for  i=1:M 
for j=1:N                
        if （R（ij）==1）     r（ij）=30;       
        elseif（R（ij）==2）  r（ij）=100;      
        elseif（R（ij）==3）  r（ij）=200; 
        end 
    end 
end  
r=uint8（r）; 
figure; 
imshow（r）;                           
imwrite（r‘r1.bmp‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件    1542558  2014-01-16 21:36  图像分析期末设计\2012_图像分析作业_刘蕾_SC12023170 .pdf

     文件       1718  2014-01-17 21:11  图像分析期末设计\chengxu1.m

     文件      66616  2007-03-05 08:43  图像分析期末设计\L1.bmp

     文件      66616  2007-03-05 08:43  图像分析期末设计\L2.bmp

     文件      66616  2007-03-05 08:43  图像分析期末设计\L3.bmp

     文件      66616  2007-03-05 08:43  图像分析期末设计\L4.bmp

     文件      66616  2007-03-05 08:43  图像分析期末设计\L5.bmp

     文件      66616  2007-03-05 08:43  图像分析期末设计\L6.bmp

     文件      66616  2007-03-05 08:43  图像分析期末设计\L7.bmp

     文件     196662  2014-01-17 21:07  图像分析期末设计\r.bmp

     文件      66614  2014-01-17 20:57  图像分析期末设计\r1.bmp

     文件      66614  2014-01-15 20:39  图像分析期末设计\r2.bmp

     文件      66614  2014-01-15 20:40  图像分析期末设计\r3.bmp

     文件      66614  2014-01-15 20:41  图像分析期末设计\r4.bmp

     文件      66614  2014-01-15 20:44  图像分析期末设计\r5.bmp

     文件      66614  2014-01-15 20:44  图像分析期末设计\r6.bmp

     文件      66614  2014-01-15 20:45  图像分析期末设计\r7.bmp

     文件       2689  2014-01-17 21:08  图像分析期末设计\zongchengxu.m

    ..A.SH.       162  2014-01-16 09:43  图像分析期末设计\~$鑫SC13023015期末实验报告.doc

     文件    2582528  2013-11-19 10:00  图像分析期末设计\图像作业\DIA课程实验2013-2014.ppt

     文件      24576  2014-01-15 15:02  图像分析期末设计\图像作业\K-means图像分割.doc

     文件       2042  2014-01-15 15:52  图像分析期末设计\图像作业\k_means.m

     文件      19065  2014-01-15 16:03  图像分析期末设计\图像作业\K_means_Image.cpp

     文件      65651  2014-01-15 14:45  图像分析期末设计\图像作业\K_means算法在遥感分类中的应用.caj

     文件      66616  2007-03-05 08:43  图像分析期末设计\图像作业\L1.bmp

     文件      66616  2007-03-05 08:43  图像分析期末设计\图像作业\L2.bmp

     文件      66616  2007-03-05 08:43  图像分析期末设计\图像作业\L3.bmp

     文件      66616  2007-03-05 08:43  图像分析期末设计\图像作业\L4.bmp

     文件      66616  2007-03-05 08:43  图像分析期末设计\图像作业\L5.bmp

     文件      66616  2007-03-05 08:43  图像分析期末设计\图像作业\L6.bmp

............此处省略46个文件信息