1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
|
function lsf = GeislerLSFMinutes(distance)
%GeislerLSFMinutes Compute Geisler estimate of 2mm pupil human LSF
% lsf = GeislerLSFMinutes(distance)
%
% Compute the 2mm pupil LSF from Geisler (1984, JOSA A, 1, pp. 775-782),
% Figure 1. This is also used in Banks, Geisler & Bennett (1987,Vision
% Research, 27, 1915-1924).
%
% They give the parameters of a sum of Gaussians LSF as a function of
% passed distance. Distance passed in minutes of arc.
%
% Return is normalized to a maximum of 1.
%
% See also WestPSFMinutes, WestLSFMinutes, LsfToPsf, PsfToLsf
weight1 = 0.684;
weight2 = 0.587;
sigma1 = 0.443;
sigma2 = 2.035;
lsf = weight1*normpdf(distance,0,sigma1) + ...
weight2*normpdf(distance,0,sigma2);
lsf = lsf/max(lsf(:));
|
