load SampleMatrix32;

I=imread('goldhill.bmp');

imshow(I)
size=256;

% Sampling
f=zeros(N*N,1);
g=zeros(size/2,size/2);

num1=size/N; 
num2=size/N;

for t1=1:1:num1
    for t2=1:1:num2
        for i=1:1:N
            for j=1:1:N
                f((i-1)*N+j)=I((t1-1)*N+i,(t2-1)*N+j);
            end
        end
        temp=SM*f;        
        for i=1:1:M
            for j=1:1:M
                g((t1-1)*M+i,(t2-1)*M+j)=temp((i-1)*M+j);
            end
        end
    end
end

I=g;
save goldhill-sample I;
figure;imshow(g,[])