source: trunk/src/get_file_series.m @ 756

Last change on this file since 756 was 754, checked in by sommeria, 11 years ago
File size: 6.7 KB
RevLine 
[376]1%'get_file_series': determine the list of file names and file indices for functions called by 'series'.
2%------------------------------------------------------------------------
3% [filecell,i1_series,i2_series,j1_series,j2_series]=get_file_series(Param)
4%
5% OUTPUT:
[446]6% filecell{iview,fileindex}: cell array representing the list of file names
7%        iview: line in the table corresponding to a given file series
8%        fileindex: file index within  the file series,
[453]9% i1_series{iview}(ref_j,ref_i)... are the corresponding arrays of indices i1,i2,j1,j2, depending on the input line iview and the two reference indices ref_i,ref_j
10% i1_series{iview}(fileindex) expresses the same indices as a 1D array in file indices
[446]11
[376]12% INPUT:
13% Param: structure of input parameters as read from the GUI series (by the function read_GUI)
14
[373]15function [filecell,i1_series,i2_series,j1_series,j2_series]=get_file_series(Param)
[376]16
[591]17filecell={};i1_series={};i2_series={};j1_series={};j2_series={};
[373]18InputTable=Param.InputTable;
19first_i=Param.IndexRange.first_i;
20incr_i=Param.IndexRange.incr_i;
21last_i=Param.IndexRange.last_i;
[635]22% ref_i=first_i:incr_i:last_i;
23%ref_j=[];
24first_j=[];last_j=[];incr_j=1;%default
[754]25if isfield(Param.IndexRange,'first_j')&& isfield(Param.IndexRange,'last_j')
[373]26    first_j=Param.IndexRange.first_j;
27    last_j=Param.IndexRange.last_j;
28end
[754]29if isfield(Param.IndexRange,'incr_j')
30    incr_j=Param.IndexRange.incr_j;
31end
[391]32
[373]33%% determine the list of input file names
34nbmissing=0;
[635]35NbView=size(InputTable,1);
36i1_series=cell(NbView,1);% initiate index series with empty cells
37i2_series=cell(NbView,1);
38j1_series=cell(NbView,1);
39j2_series=cell(NbView,1);
40for iview=1:NbView
[408]41    r.mode='';
42    if isfield (Param.IndexRange,'PairString')
[450]43        if ischar(Param.IndexRange.PairString)
44            Param.IndexRange.PairString={Param.IndexRange.PairString};
45        end
[635]46        r=regexp(Param.IndexRange.PairString{iview,1},'(?<mode>(Di=)|(Dj=)) -*(?<num1>\d+)\|(?<num2>\d+)','names');%look for mode=Dj or Di
[408]47        if isempty(r)
[635]48            r=regexp(Param.IndexRange.PairString{iview,1},'(?<num1>\d+)(?<mode>-)(?<num2>\d+)','names');%look for burst pairs
[537]49        end
[408]50        % TODO case of free pairs:
51        %r=regexp(pair_string,'.*\D(?<num1>[\d+|*])(?<delim>[-||])(?<num2>[\d+|*])','names');
52    end
[440]53    if isempty(r)||isempty(r.mode)
54        r(1).num1='';
55        r(1).num2='';
[540]56        if isfield (Param.IndexRange,'PairString') && strcmp(Param.IndexRange.PairString{iview,1},'j=*-*')
[537]57            r(1).mode='*-*';
58        else
[635]59            r(1).mode='';
[537]60        end
[408]61    end
[639]62   
63    if isempty(incr_i) || isempty(incr_j) || isequal(r(1).mode,'*-*')|| isequal(r(1).mode,'*|*')% free pairs or increment
[537]64        FilePath=fullfile(InputTable{iview,1},InputTable{iview,2});
65        fileinput=[InputTable{iview,3} InputTable{iview,4} InputTable{iview,5}];
[635]66        [tild,tild,tild,i1_series{iview},i2_series{iview},j1_series{iview},j2_series{iview},NomType,FileType,FileInfo,MovieObject,...
67            i1_input,i2_input,j1_input,j2_input]=find_file_series(FilePath,fileinput);
68        i1_series{iview}=squeeze(i1_series{iview}(1,:,:)); %select first  pair index as ordered by find_file_series
[639]69        i2_series{iview}=squeeze(i2_series{iview}(1,:,:)); %select first  pair index as ordered by find_file_series
[537]70        j1_series{iview}=squeeze(j1_series{iview}(1,:,:)); %first  pair index
[635]71        j2_series{iview}=squeeze(j2_series{iview}(1,:,:)); %second  pair index
72        if isempty(incr_i)
[639]73            if isempty(first_j) || isempty(incr_j) % no j index or no defined increment for j
[643]74                ref_j=find(max(i1_series{iview},[],2));
75                ref_i=find(max(i1_series{iview},[],1));
[635]76                ref_i=ref_i-1;
77                ref_j=ref_j-1;
78                ref_i=ref_i(ref_i>=first_i & ref_i<=last_i);
79                ref_j=ref_j(ref_j>=first_j & ref_j<=last_j);
80            else
81                ref_j=first_j:incr_j:last_j;
[643]82                ref_i=find(max(i1_series{iview}(:,ref_j),[],1));
[635]83                ref_i=ref_i-1;
84                ref_i=ref_i(ref_i>=first_i & ref_i<=last_i);
85            end
86        else
87            ref_i=first_i:incr_i:last_i;%default
[639]88            if isempty(first_j) ||isempty(incr_j)% no j index or no defined increment for j
[643]89                ref_j=find(max(i1_series{iview},[],2));
[635]90                ref_j=ref_j-1;
91                ref_j=ref_j(ref_j>=first_j & ref_j<=last_j);
92            else
93                ref_j=first_j:incr_j:last_j;
94            end
95        end
[639]96        if isempty(ref_j)
97            i1_series{iview}=i1_series{iview}(2,ref_i+1);
98            if ~isempty(i2_series{iview})
99                i2_series{iview}=i2_series{iview}(2,ref_i+1);
100            end
101        else
102            i1_series{iview}=i1_series{iview}(ref_j+1,ref_i+1);
103            if ~isempty(i2_series{iview})
104                i2_series{iview}=i2_series{iview}(ref_j+1,ref_i+1);
105            end
[635]106        end
[639]107        if ~isempty(j1_series{iview})
108            j1_series{iview}=j1_series{iview}(ref_j+1,ref_i+1);
109            if ~isempty(j2_series{iview})
110                j2_series{iview}=j2_series{iview}(ref_j+1,ref_i+1);
111            end
112        end
[537]113    else
[635]114        ref_i=first_i:incr_i:last_i;%default
115        ref_j=first_j:incr_j:last_j;%default
[537]116        [i1_series{iview},i2_series{iview},j1_series{iview},j2_series{iview}]=find_file_indices(ref_i,ref_j,str2num(r.num1),str2num(r.num2),r.mode);
117    end
[639]118    %     if ~isequal(r(1).mode,'*-*')% imposed pairs or single i and/or j index
119    %         [i1_series{iview},i2_series{iview},j1_series{iview},j2_series{iview}]=find_file_indices(ref_i,ref_j,str2num(r.num1),str2num(r.num2),r.mode);
120    %     end
[635]121   
[639]122    %list of files
[373]123    i2=[];j1=[];j2=[];
[453]124    for ifile=1:numel(i1_series{iview})
125        i1=i1_series{iview}(ifile);
126        if ~isempty(i2_series{iview})
127            i2=i2_series{iview}(ifile);
[373]128        end
[453]129        if ~isempty(j1_series{iview})
130            j1=j1_series{iview}(ifile);
[373]131        end
[453]132        if ~isempty(j2_series{iview})
133            j2=j2_series{iview}(ifile);
[373]134        end
[635]135        filecell{iview,ifile}=fullfile_uvmat(InputTable{iview,1},InputTable{iview,2},InputTable{iview,3},InputTable{iview,5},InputTable{iview,4},i1,i2,j1,j2);
[373]136    end
137end
138
139
140function [i1_series,i2_series,j1_series,j2_series]=find_file_indices(ref_i,ref_j,num1,num2,mode)
141i1_series=ref_i;%default
142j1_series=[];
143if ~isempty(ref_j)
[400]144    [i1_series,j1_series]=meshgrid(ref_i,ref_j);
[373]145end
[393]146i2_series=i1_series;
147j2_series=j1_series;
[391]148
[373]149switch mode
[393]150    case 'Di='  %  case 'series(Di)')
[373]151        i1_series=i1_series-num1;
152        i2_series=i2_series+num2;
[393]153    case 'Dj='  %  case 'series(Dj)'
154        j1_series=j1_series-num1;
155        j2_series=j2_series+num2;
156    case '-'  % case 'bursts'
157        j1_series=num1*ones(size(i1_series));
158        j2_series=num2*ones(size(i1_series));
[373]159end
Note: See TracBrowser for help on using the repository browser.