Changeset 1188 for trunk/src/series


Ignore:
Timestamp:
Dec 11, 2025, 8:45:31 AM (2 months ago)
Author:
sommeria
Message:

file relabeling implemented for civ_series

Location:
trunk/src/series
Files:
2 edited

Legend:

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

    r1187 r1188  
    163163        MaxIndex_i=find(SeriesData.i1_series{iview_image}(1,2,:),1,'last' )-1;%max ref index i detected in the series (corresponding to the last non-zero value of i1_series)
    164164    else
    165         ref_i=squeeze(max(SeriesData.i1_series{iview_image}(1,:,:),[],2));% select ref_j index for each ref_i
    166         ref_j=squeeze(max(SeriesData.j1_series{iview_image}(1,:,:),[],3));% select ref_i index for each ref_j
    167         MinIndex_i=min(find(ref_i))-1;
    168         MaxIndex_i=max(find(ref_i))-1;
    169         MaxIndex_j=max(find(ref_j))-1;
    170         MinIndex_j=min(find(ref_j))-1;
     165       
     166        if ndims(SeriesData.j1_series{iview_image})==3% usuual file series input
     167            ref_i=squeeze(max(SeriesData.i1_series{iview_image}(1,:,:),[],2));% select ref_j index for each ref_i
     168            ref_j=squeeze(max(SeriesData.j1_series{iview_image}(1,:,:),[],3));% select ref_i index for each ref_j
     169            MinIndex_i=min(find(ref_i))-1;
     170            MaxIndex_i=max(find(ref_i))-1;
     171            MaxIndex_j=max(find(ref_j))-1;
     172            MinIndex_j=min(find(ref_j))-1;
     173        else %case with relabeling
     174            MinIndex_i=1;
     175            MaxIndex_i=numel(SeriesData.i1_series{iview_image});% case relabel
     176            MinIndex_j=1;
     177            MaxIndex_j=max(1,numel(SeriesData.j1_series{iview_image}));% =1 if j1_series empty
     178        end
    171179    end
    172180end
     
    176184time=[];
    177185TimeUnit='frame'; %default
    178 CoordUnit='';%default
    179 pxcm_search=1;
     186% CoordUnit='';%default
     187% pxcm_search=1;
    180188if isfield(SeriesData,'Time') &&numel(SeriesData.Time')>=1 && ~isempty(SeriesData.Time{1})
    181189    time=SeriesData.Time{1};
     
    184192    TimeUnit=Param.IndexRange.TimeUnit;
    185193end
    186 if isfield(SeriesData,'GeometryCalib')
    187     tsai=SeriesData.GeometryCalib;
    188     if isfield(tsai,'fx_fy')
    189         pxcm_search=max(tsai.fx_fy(1),tsai.fx_fy(2));%pixels:cm estimated for the search range
    190     end
    191     if isfield(tsai,'CoordUnit')
    192         CoordUnit=tsai.CoordUnit;
    193     end
    194 end
     194% if isfield(SeriesData,'GeometryCalib')
     195%     tsai=SeriesData.GeometryCalib;
     196%     if isfield(tsai,'fx_fy')
     197%         pxcm_search=max(tsai.fx_fy(1),tsai.fx_fy(2));%pixels:cm estimated for the search range
     198%     end
     199%     if isfield(tsai,'CoordUnit')
     200%         CoordUnit=tsai.CoordUnit;
     201%     end
     202% end
    195203
    196204%% timing display
     
    824832%reproduce by default the chosen pair in the checkciv2 menu
    825833set(handles.ListPairCiv2,'Value',get(handles.ListPairCiv1,'Value'))%civ2 selection the same as civ1 by default
    826 %ListPairCiv2_Callback(hObject, eventdata, handles)
    827834
    828835% %------------------------------------------------------------------------
     
    978985        end
    979986        if index_pair ~=0
    980         [tild,indsort]=sort(dt);
     987        [~,indsort]=sort(dt);
    981988        displ_pair=displ_pair(indsort);
    982989        displ_pair_dt=displ_pair_dt(indsort);
  • trunk/src/series/civ_series.m

    r1187 r1188  
    256256
    257257%% File relabeling documented by the xml file (e.g. PCO)
    258 CheckRelabel=isfield(Param,'FileSeries' );%=true for index relabeling (PCO)
    259 
    260 % %% introduce input image transform
    261 % transform_fct=[];%default, no transform
    262 % if isfield(Param,'FieldTransform')&&~isempty(Param.FieldTransform.TransformName)
    263 %         currentdir=pwd;
    264 %     cd(Param.FieldTransform.TransformPath)
    265 %     transform_fct=str2func(Param.FieldTransform.TransformName);
    266 %     cd (currentdir)
    267 % end
     258CheckRelabel=isfield(Param.IndexRange,'Relabel' )&& Param.IndexRange.Relabel;%=true for index relabeling (PCO)
    268259
    269260%%%%% MAIN LOOP %%%%%%
     
    291282    OutputPath=fullfile(Param.OutputPath,Param.Experiment,Param.Device);
    292283    if CheckRelabel
    293          RootFileOut=index2filename(Param.FileSeries,1,1,MaxIndex_j);
     284         RootFileOut=index2filename(XmlData.FileSeries,1,1,MaxIndex_j);
    294285    else
    295286        RootFileOut=RootFile_A;
     
    343334                FrameIndex_A=1;
    344335            elseif CheckRelabel
    345             [RootFile,FrameIndex_A]=index2filename(Param.FileSeries,i1_series_Civ1(ifield),j1_series_Civ1(ifield),MaxIndex_j);
    346             ImageName_A=fullfile(RootPath_A,SubDir_A,RootFile);
     336                [RootFile,FrameIndex_A]=index2filename(XmlData.FileSeries,i1_series_Civ1(ifield),j1_series_Civ1(ifield),MaxIndex_j);
     337                ImageName_A=fullfile(RootPath_A,SubDir_A,RootFile);
    347338            else
    348339                ImageName_A=fullfile_uvmat(RootPath_A,SubDir_A,RootFile_A,FileExt_A,NomType_A,i1_series_Civ1(ifield),[],j1_series_Civ1(ifield));
     
    382373            end
    383374            if CheckRelabel
    384                 [RootFile,FrameIndex_B]=index2filename(Param.FileSeries,i2_series_Civ1(ifield),j2_series_Civ1(ifield),MaxIndex_j);
     375                [RootFile,FrameIndex_B]=index2filename(XmlData.FileSeries,i2_series_Civ1(ifield),j2_series_Civ1(ifield),MaxIndex_j);
    385376                ImageName_B=fullfile(RootPath_B,SubDir_B,RootFile);
    386377            else
     
    632623            ImageName_A_Civ2=Param.ActionInput.RefFile;
    633624        elseif CheckRelabel
    634             [RootFile,FrameIndex_A_2]=index2filename(Param.FileSeries,i1_series_Civ2(ifield),j1_series_Civ2(ifield),MaxIndex_j);
     625            [RootFile,FrameIndex_A_2]=index2filename(XmlData.FileSeries,i1_series_Civ2(ifield),j1_series_Civ2(ifield),MaxIndex_j);
    635626            ImageName_A_Civ2=fullfile(RootPath_A,SubDir_A,RootFile);
    636627        else
     
    644635        end
    645636        if CheckRelabel
    646             [RootFile,FrameIndex_B_2]=index2filename(Param.FileSeries,i2_civ2,j2_civ2,MaxIndex_j);
     637            [RootFile,FrameIndex_B_2]=index2filename(XmlData.FileSeries,i2_civ2,j2_civ2,MaxIndex_j);
    647638            ImageName_B_Civ2=fullfile(RootPath_B,SubDir_B,RootFile);
    648639        else
Note: See TracChangeset for help on using the changeset viewer.