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

    function [im2] = window9
%%%%%%for indian data
%%%%%%for 3*3 window
load indian_pines
im1= indian_pines;



for i=1:147
    for j=1:147
        a=linspace(0,0,220);
        im3(i,j,:)=a;
        end
end


p=1;
for i=1:145
    p=p+1;s=2;
    for j=1:145
        
        im3(p,s,:)=im1(i,j,:);
        s=s+1;                                                                 
    end
end

p=1;
for i=1:145
    p=p+1;s=2;
    for j=1:145
        
        im2(i,j,:)=(im3(p,s,:)+im3(p-1,s-1,:)+im3(p-1,s-1,:)+im3(p-1,s,:)+im3(p-1,s+1,:)+im3(p,s-1,:)+im3(p,s+1,:)+im3(p+1,s-1,:)+im3(p+1,s,:)+im3(p+1,s+1,:))/9;
        s=s+1;                                                                 
    end
end