Gateway routine to various human space-time sensitivity curves
[sens,spaceSamples,timeSamples] = humanSpaceTime(spaceVariable,timeVariable)
Current calculations includes
* Kelly's 1979 space time surface
* Watson's temporal impulse response
* Poirson/Wandell spatial color models
These functions are implemented in other functions and this routine
simply formulates the call.
The space and time variables are either in the time/space domain or
frequency domain, depending on the call. See below for examples.
Example:
[sens,fs,ft] = humanSpaceTime('kelly79',logspace(-0.2,log10(30),20),logspace(-0.2,log10(60),20));
surf(ft,fs,sens); set(gca,'xscale','log','yscale','log');
set(gca,'xlim',[0 60],'ylim',[0 30]);
xlabel('Temporal freq (Hz)');ylabel('Spatial freq (cpd)')
[tMTF,junk,ft] = humanSpaceTime('watsonTMTF',[],[0:60]);
plot(ft,tMTF,'r-o');
t = [0.001:0.002:0.200];
[impResp,junk,t] = humanSpaceTime('watsonImpulseResponse',[],t);
plot(t,impResp)
[spatialFilters,positions] = humanSpaceTime('poirsoncolor');
mesh(positions,positions,spatialFilters.by); xlabel('Position (deg)')
Copyright ImagEval Consultants, LLC, 2005.