%%
%	caller script for lrf extraction
%

% load data which loads tr_data , te_data, 
% te_labels, tr_labels variables on environment
load data.mat

% load 3d-positions of each voxel
pos = vXYZ';

% call for lpc with specified neighborhood set
[a_tr, a_te]=find_lpc(pos,2,tr_data,te_data);




