1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
|
function [image] = MakeSincImage(freq,i0,j0,nRowPixels,nColPixels)
% image = MakeSincImage(freq,i0,j0,nRowPixels,[nColPixels])
%
% Computes a two-dimensional sinc function image.
% The sinc function has cut-off frequency freq
% (in cycles/image) and is positioned at (i0,j0).
%
% The image has dimensions nRowPixels by nColPixels.
% If nColPixels is omitted, a square image is returned.
%
% 8/15/94 dhb Added optional column pixels argument.
% 10/28/98 dhb Change to call Matlab sinc rather than my Sinc.
usefreq = 2*freq;
x = ((1:nColPixels) - j0)/nColPixels;
y = ((1:nRowPixels) - i0)/nRowPixels;
sincx = sinc(usefreq*x);
sincy = sinc(usefreq*y);
image = sincy'*sincx;
|
