% radiale Gauss-Filter unterschiedlicher Breite bauen

% gewünschte Bildschirmauflösung
X_Pixel = 1280;
Y_Pixel = 1024;
% X_Pixel = 1024;
% Y_Pixel = 768;
% X_Pixel = 800;
% Y_Pixel = 600;

for groesse=10:20:190
    maxx=X_Pixel;
    maxy=Y_Pixel;
    [xi,yi] = meshgrid(1:1:maxx,1:1:maxy);
    xc=xi-round(maxx/2);
    yc=yi-round(maxy/2);
    zi=exp(-(xc.*xc+yc.*yc)/(groesse*maxy));
%     [xn,yn] = meshgrid(1:20:X_Pixel,1:20:Y_Pixel);
%     zn=interp2(xi,yi,zi,xn,yn,'linear',0);
    fname=sprintf('GaussHole_%d_%dx%d.dat', (groesse/10+1)/2, X_Pixel, Y_Pixel);
    f=fopen(fname,'w','ieee-le');
    fwrite(f,zi','float32');
    fclose(f);
end
