function out=edf(in)

    [sizex,sizey]=size(in);
    in=double(in);
    
    blocksize=10;
    
    for (x=1:sizex-blocksize)
        for (y=1:sizey-blocksize)
            block=in(x:x+blocksize-1,y:y+blocksize);
            blockmean=mean2(block);
            blockstd=std2(block);
            out(x,y)=blockmean+blockstd*j;
        end
    end


end