[f1, f2] = freqspace(512, 'meshgrid');
D1=VALKOINEN PISTE #1
D2=VALKOINEN PISTE #2

D = sqrt(f1.^2 + f2.^2);
Hlp((D>D2 & D<D1)) = 0;