Ignore:
Timestamp:
Mar 3, 2010, 1:47:56 PM (11 years ago)
Author:
sommeria
Message:

-subfunction find_file_indices put back in series (it had been extracted by mistake to check_files
-aver_stat: name of output file rationalized
-read_get_field : pb with vector fields, yet to solve

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/src/series/check_files.m

    r27 r32  
    150150    hh=uicontrol('Style','listbox', 'Position', [20 340 500 40], 'String', message);
    151151end
    152 %----------------------------------------------------
    153 %  determine the list of index pairs of processing file
    154 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    155 function [num_i1,num_i2,num_j1,num_j2,num_i_out,num_j_out]=find_file_indices(num_i,num_j,ind_shift,NomType,mode)
    156 num_i1=num_i;% set of first image numbers by default
    157 num_i2=num_i;
    158 num_j1=num_j;
    159 num_j2=num_j;
    160 num_i_out=num_i;
    161 num_j_out=num_j;
    162 if isequal (NomType,'_i1-i2_j') |isequal (NomType,'_i1-i2')
    163     num_i1_line=num_i+ind_shift(3);% set of first image numbers
    164     num_i2_line=num_i+ind_shift(4);
    165     % adjust the first and last field number
    166         indsel=find(num_i1_line >= 1);
    167     num_i_out=num_i(indsel);
    168     num_i1_line=num_i1_line(indsel);
    169     num_i2_line=num_i2_line(indsel);
    170     num_j1=meshgrid(num_j,ones(size(num_i1_line)));
    171     num_j2=meshgrid(num_j,ones(size(num_i1_line)));
    172     [xx,num_i1]=meshgrid(num_j,num_i1_line);
    173     [xx,num_i2]=meshgrid(num_j,num_i2_line);
    174 elseif isequal (NomType,'_i_j1-j2') || isequal (NomType,'#_ab')
    175     if isequal(mode,'bursts') %case of bursts (png_old or png_2D)
    176         num_j1=ind_shift(1)*ones(size(num_i));
    177         num_j2=ind_shift(2)*ones(size(num_i));
    178     else
    179         num_j1_col=num_j+ind_shift(1);% set of first image numbers
    180         num_j2_col=num_j+ind_shift(2);
    181         % adjust the first field number
    182         indsel=find((num_j1_col >= 1));   
    183         num_j_out=num_j(indsel);
    184         num_j1_col=num_j1_col(indsel);
    185         num_j2_col=num_j2_col(indsel);
    186         [num_i1,num_j1]=meshgrid(num_i,num_j1_col);
    187         [num_i2,num_j2]=meshgrid(num_i,num_j2_col);
    188     end   
    189 end
Note: See TracChangeset for help on using the changeset viewer.