clear;
close all;

graphics_toolkit fltk;
pkg load signal;
pkg load image;

inImg = imread("pattern1.bmp");

tRadius = 20;
pRadius = 2;

curPatch = zeros((2*pRadius+1),(2*pRadius+1));
curSurrPatch = zeros((2*pRadius+1),(2*pRadius+1));
curDiff = zeros((2*pRadius+1),(2*pRadius+1));

distSurface = zeros((2*tRadius+1),(2*tRadius+1));

outImg = zeros((rows(inImg)-40)/(2*pRadius+1),(columns(inImg)-40)/(2*pRadius+1));

size(outImg)
pause

oI = 0;
oJ = 0;

for i=(1+pRadius+tRadius):(2*pRadius+1):(rows(inImg)-pRadius-tRadius)

    oI = oI+1;
    oJ = 0;

    for j=(1+pRadius+tRadius):(2*pRadius+1):(columns(inImg)-pRadius-tRadius)
         oJ = oJ + 1;

        curPatch(:,:) = inImg(i-pRadius:i+pRadius,j-pRadius:j+pRadius);

        for u=(i-tRadius):1:(i+tRadius)
            for v=(j-tRadius):1:(j+tRadius)
                curSurrPatch(:,:) = inImg(u-pRadius:u+pRadius,v-pRadius:v+pRadius);
                curDiff(:,:) = curPatch .- curSurrPatch;
                distSurface(u+1-(i-tRadius),v+1-(j-tRadius)) = sqrt(sum(sum(curDiff.^2)));
            endfor
        endfor

        outImg(oI,oJ) = max(max(distSurface));
    endfor

    oI
    oJ
    fflush(stdout);
endfor