KernelVisual.html

KernelVisual.m

% KernelVisual :: Script to visualize each of the various spatial kernels
% that are employed in the MODE model
%
%% Reference
% Grossberg, S. and Pilly, P. K. (2008). Temporal dyanamics of decision-making during motion perception in the visual cortex. Vision Research, 48(12), 1345-1373.
%
%% Author
% Praveen K. Pilly ([e-mail])
%
%% License policy
% Written by Praveen K. Pilly, Department of Cognitive and Neural Systems, Boston University
% Copyright 2009, Trustees of Boston University
%
% Permission to use, copy, modify, distribute, and sell this software and its documentation for any purpose is hereby granted
% without fee, provided that the above copyright notice and this permission notice appear in all copies, derivative works and
% associated documentation, and that neither the name of Boston University nor that of the author(s) be used in advertising or
% publicity pertaining to the distribution or sale of the software without specific, prior written permission. Neither Boston
% University nor its agents make any representations about the suitability of this software for any purpose. It is provided "as
% is" without warranty of any kind, either express or implied. Neither Boston University nor the author indemnify any
% infringement of copyright, patent, trademark, or trade secret resulting from the use, modification, distribution or sale of
% this software.
%
%% Last modified
% June 25, 2009

%%
clear all

numdir=8;

%%% Various kernel generations
% Kernel G in Stage 5
sigG1=2;
sigG2=0.5;
Ga=100;

% Kernel J in Stage 6
sigJ1=3;
sigJ2=1;
Ja=75;

% Kernel K in Stage 6
sigK=4;
Ka=75;

% Kernel L in Stage 7
sigL1=10;
sigL2=3;
La=10;

% Kernel P in Stages 7 and 8
sigP=8;
Pa=1;

% Anisotropic kernels [G, J, and L] and offset isotropic kernel [K]
for d=1:numdir
    theta=(d-1)*2*pi/numdir;

    G(:,:,d)=anisot(Ga,theta,sigG1,sigG2);
    J(:,:,d)=anisot(Ja,theta,sigJ1,sigJ2);
    L(:,:,d)=anisot(La,theta,sigL1,sigL2);

    K(:,:,d)=isot(Ka,sigK,2,theta);
end

% Isotropic kernel
P=isot(Pa,sigP,1,0);


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% To visualize, choose any kernel and the direction number (if applicable)
%
% e.g., J and d=2 (pi/4)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

figure
imagesc(J(:,:,2)), colormap(jet), colorbar