-
Notifications
You must be signed in to change notification settings - Fork 0
/
load_centData_oneBased.m
89 lines (74 loc) · 3.07 KB
/
load_centData_oneBased.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
%Function to load tracking structure data
%*it was moved here to make sure it works properly
function varargout = load_centData(varargin)
% charge les donn�es associ�e au jeu dataset
dataset = varargin{1};
if iscell(dataset)
prefix = dataset{1};
dataset = dataset{2};
else
prefix = '';
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% d�finition de datapath reprise dans les autres programmes %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
fsep = filesep; %
datapath = ['../../MechanicsData/']; %
%datapath='/home/public/MechanicsData';
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if strcmp(dataset,'datapath'); varargout{1}=datapath; return; end
path2 = [datapath prefix dataset '_t/'];
infofile = [path2 dataset '_info.mat'];
if exist(infofile,'file')
a = load(infofile);
a= a.dataset_info;
t0 = a.t0; % en heure
dt = a.dt; % en heure
extent = a.extent; % en pixel le long de chaque dimension
spac = a.spac; % dimension d'un voxel en micron dans chaque dimension
type = a.type; % nomenclature de Nadine pour l'orientation des donn�es
if isfield(a,'manualSph'), manualSph = a.manualSph; end % 4 valeurs, centre de la sph�re fittant les donn�es et sont rayon en deriere position (en microns)
if isfield(a,'manualRef'), manualRef = a.manualRef; end % matrice de rotation
if isfield(a,'axisRange'), axisRange = a.axisRange; end % 6 valeur xmin,xmax,... en microns
if isfield(a,'viewAngle'), viewAngle = a.viewAngle; end % azimuth elevation (voir 'view' function)
if isfield(a,'nbstepValid'), nbstepValid = a.nbstepValid; end % nombre de pas de temps, � partir de 1, estimer ok pour les mesures
else
error(['pas de donn�e correspondant � ' dataset]);
end
timeline = t0:dt:t0+(nbstepValid-1)*dt ; % temps en heure de chaque timestep
%% output preparation
nbstepValid
pause
for i = 2:length(varargin)
if strcmp(varargin{i},'cent')
varargin{i} = 'XYZ';
end
end
filelist = ls(path2);
filelist = [ filelist char(32*ones(size(filelist,1),1)) ]; % on rajoute un espace � la fin de chaque ligne (necessaire � la compatibilit� windows)
filelist = reshape(filelist',[1,numel(filelist)]);
filelist(strfind(filelist,char(10))) = ' '; % on remplace les retour chariot par des espace (compatibilit� mac)
filelist = strsplit(' ',filelist);
for j = 2:length(varargin)
varargout{j-1} = -1;
i=1;
[dataset '_' varargin{j} '.mat']
while i <= length(filelist)
name = filelist{i};
%strcmp([dataset '_' varargin{j} '.mat'],name)
if strcmp([dataset '_' varargin{j} '.mat'],name)
[path2 name]
a = load( [path2 name] );
a = struct2cell(a);
varargout{j-1} = a{1};
size(a{1})
%else
% 'no match'
end
i=i+1;
end
if exist(varargin{j},'var')==1
eval(['varargout{j-1}=' varargin{j} ';']);
end
end
varargout