[Matlab] Define 2D sinc() function

%Define Cartesian Coordinates

radiusSZ = floor(imSize/2);

[X, Y] = meshgrid(-radiusSZ:radiusSZ, -radiusSZ:radiusSZ);

[theta, rho] = cart2pol(X, Y);

rho = round(rho);

%2d sinc function

sincFrq = 1;

sincFilter= sin(rho*sincFrq+0.00001)./(rho*sincFrq+0.00001);

sincFilter_clamped = sincFilter; %init

%If you want to clamp the kernel,

clampDist = 4; %variable

cen = ceil(imSize/2); %variable

sincFilter_clamped = sincFilter_clamped(cen-clampDist:cen+clampDist,cen-clampDist:cen+clampDist);