function [e] = silopien(k,s,z)
%silopien(k,n,z)
%Silottaa ja pienentää matriisia z
%Pienentää osaan 1/k ja silottaa 2s+1-levyisen neliön keskiarvoksi

[n,m]=size(z);
N=floor(n/k);M=floor(m/k);
e=zeros(N,M);
for i=1:N
for j=1:M
e(i,j)=mean(mean(z(max(1,k*i-s):min(k*i+s,n),max(k*j-s,1):min(k*j+s,m))));
end
end


end

