source: trunk/src/get_file_series.m @ 1160

Last change on this file since 1160 was 1134, checked in by sommeria, 6 months ago

adapted to reading hdf5 data from fluidimage, extract_rdvision improved

File size: 7.8 KB
RevLine 
[757]1%'get_file_series': determine the list of input file names and file indices for functions called by 'series'.
[376]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
[809]11%
[376]12% INPUT:
13% Param: structure of input parameters as read from the GUI series (by the function read_GUI)
14
[809]15%=======================================================================
[1126]16% Copyright 2008-2024, LEGI UMR 5519 / CNRS UGA G-INP, Grenoble, France
[809]17%   http://www.legi.grenoble-inp.fr
[1127]18%   Joel.Sommeria - Joel.Sommeria (A) univ-grenoble-alpes.fr
[809]19%
20%     This file is part of the toolbox UVMAT.
21%
22%     UVMAT is free software; you can redistribute it and/or modify
23%     it under the terms of the GNU General Public License as published
24%     by the Free Software Foundation; either version 2 of the license,
25%     or (at your option) any later version.
26%
27%     UVMAT is distributed in the hope that it will be useful,
28%     but WITHOUT ANY WARRANTY; without even the implied warranty of
29%     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
30%     GNU General Public License (see LICENSE.txt) for more details.
31%=======================================================================
32
[373]33function [filecell,i1_series,i2_series,j1_series,j2_series]=get_file_series(Param)
[376]34
[757]35filecell={};
[373]36InputTable=Param.InputTable;
37first_i=Param.IndexRange.first_i;
[1097]38incr_i=Param.IndexRange.incr_i;
[373]39last_i=Param.IndexRange.last_i;
[635]40first_j=[];last_j=[];incr_j=1;%default
[754]41if isfield(Param.IndexRange,'first_j')&& isfield(Param.IndexRange,'last_j')
[373]42    first_j=Param.IndexRange.first_j;
43    last_j=Param.IndexRange.last_j;
44end
[754]45if isfield(Param.IndexRange,'incr_j')
46    incr_j=Param.IndexRange.incr_j;
47end
[391]48
[373]49%% determine the list of input file names
[635]50NbView=size(InputTable,1);
51i1_series=cell(NbView,1);% initiate index series with empty cells
52i2_series=cell(NbView,1);
53j1_series=cell(NbView,1);
54j2_series=cell(NbView,1);
[1134]55
56%LOOP ON INPUT FILE SERIES
[635]57for iview=1:NbView
[408]58    r.mode='';
59    if isfield (Param.IndexRange,'PairString')
[450]60        if ischar(Param.IndexRange.PairString)
61            Param.IndexRange.PairString={Param.IndexRange.PairString};
62        end
[871]63        if size(Param.IndexRange.PairString,1)>=iview && ~isempty(Param.IndexRange.PairString{iview,1})
64            r=regexp(Param.IndexRange.PairString{iview,1},'(?<mode>(Di=)|(Dj=)) -*(?<num1>\d+)\|(?<num2>\d+)','names');%look for mode=Dj or Di
65            if isempty(r)
66                r=regexp(Param.IndexRange.PairString{iview,1},'(?<num1>\d+)(?<mode>-)(?<num2>\d+)','names');%look for burst pairs
67            end
[537]68        end
[408]69        % TODO case of free pairs:
70        %r=regexp(pair_string,'.*\D(?<num1>[\d+|*])(?<delim>[-||])(?<num2>[\d+|*])','names');
71    end
[440]72    if isempty(r)||isempty(r.mode)
73        r(1).num1='';
74        r(1).num2='';
[871]75        if isfield (Param.IndexRange,'PairString') && size(Param.IndexRange.PairString,1)>=iview &&...
76                                                           strcmp(Param.IndexRange.PairString{iview,1},'j=*-*')
[537]77            r(1).mode='*-*';
78        else
[635]79            r(1).mode='';
[537]80        end
[408]81    end
[1134]82    % case of free pairs or increment
[639]83    if isempty(incr_i) || isempty(incr_j) || isequal(r(1).mode,'*-*')|| isequal(r(1).mode,'*|*')% free pairs or increment
[537]84        FilePath=fullfile(InputTable{iview,1},InputTable{iview,2});
85        fileinput=[InputTable{iview,3} InputTable{iview,4} InputTable{iview,5}];
[783]86        [tild,tild,tild,i1_series{iview},i2_series{iview},j1_series{iview},j2_series{iview},NomType,FileInfo,MovieObject,...
[635]87            i1_input,i2_input,j1_input,j2_input]=find_file_series(FilePath,fileinput);
88        i1_series{iview}=squeeze(i1_series{iview}(1,:,:)); %select first  pair index as ordered by find_file_series
[639]89        i2_series{iview}=squeeze(i2_series{iview}(1,:,:)); %select first  pair index as ordered by find_file_series
[537]90        j1_series{iview}=squeeze(j1_series{iview}(1,:,:)); %first  pair index
[635]91        j2_series{iview}=squeeze(j2_series{iview}(1,:,:)); %second  pair index
92        if isempty(incr_i)
[639]93            if isempty(first_j) || isempty(incr_j) % no j index or no defined increment for j
[643]94                ref_j=find(max(i1_series{iview},[],2));
95                ref_i=find(max(i1_series{iview},[],1));
[635]96                ref_i=ref_i-1;
97                ref_j=ref_j-1;
98                ref_i=ref_i(ref_i>=first_i & ref_i<=last_i);
99                ref_j=ref_j(ref_j>=first_j & ref_j<=last_j);
100            else
101                ref_j=first_j:incr_j:last_j;
[643]102                ref_i=find(max(i1_series{iview}(:,ref_j),[],1));
[635]103                ref_i=ref_i-1;
104                ref_i=ref_i(ref_i>=first_i & ref_i<=last_i);
105            end
106        else
107            ref_i=first_i:incr_i:last_i;%default
[639]108            if isempty(first_j) ||isempty(incr_j)% no j index or no defined increment for j
[643]109                ref_j=find(max(i1_series{iview},[],2));
[635]110                ref_j=ref_j-1;
111                ref_j=ref_j(ref_j>=first_j & ref_j<=last_j);
112            else
113                ref_j=first_j:incr_j:last_j;
114            end
115        end
[639]116        if isempty(ref_j)
117            i1_series{iview}=i1_series{iview}(2,ref_i+1);
118            if ~isempty(i2_series{iview})
119                i2_series{iview}=i2_series{iview}(2,ref_i+1);
120            end
121        else
122            i1_series{iview}=i1_series{iview}(ref_j+1,ref_i+1);
123            if ~isempty(i2_series{iview})
124                i2_series{iview}=i2_series{iview}(ref_j+1,ref_i+1);
125            end
[635]126        end
[639]127        if ~isempty(j1_series{iview})
128            j1_series{iview}=j1_series{iview}(ref_j+1,ref_i+1);
129            if ~isempty(j2_series{iview})
130                j2_series{iview}=j2_series{iview}(ref_j+1,ref_i+1);
131            end
132        end
[1134]133    % case of imposed file index increment
[537]134    else
[635]135        ref_i=first_i:incr_i:last_i;%default
136        ref_j=first_j:incr_j:last_j;%default
[537]137        [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);
138    end
[639]139    %     if ~isequal(r(1).mode,'*-*')% imposed pairs or single i and/or j index
140    %         [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);
141    %     end
[635]142   
[639]143    %list of files
[373]144    i2=[];j1=[];j2=[];
[453]145    for ifile=1:numel(i1_series{iview})
146        i1=i1_series{iview}(ifile);
147        if ~isempty(i2_series{iview})
148            i2=i2_series{iview}(ifile);
[373]149        end
[453]150        if ~isempty(j1_series{iview})
151            j1=j1_series{iview}(ifile);
[373]152        end
[453]153        if ~isempty(j2_series{iview})
154            j2=j2_series{iview}(ifile);
[373]155        end
[635]156        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]157    end
158end
159
160
161function [i1_series,i2_series,j1_series,j2_series]=find_file_indices(ref_i,ref_j,num1,num2,mode)
162i1_series=ref_i;%default
163j1_series=[];
164if ~isempty(ref_j)
[400]165    [i1_series,j1_series]=meshgrid(ref_i,ref_j);
[373]166end
[393]167i2_series=i1_series;
168j2_series=j1_series;
[391]169
[373]170switch mode
[393]171    case 'Di='  %  case 'series(Di)')
[373]172        i1_series=i1_series-num1;
173        i2_series=i2_series+num2;
[393]174    case 'Dj='  %  case 'series(Dj)'
175        j1_series=j1_series-num1;
176        j2_series=j2_series+num2;
177    case '-'  % case 'bursts'
178        j1_series=num1*ones(size(i1_series));
179        j2_series=num2*ones(size(i1_series));
[809]180end
Note: See TracBrowser for help on using the repository browser.