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

file relabeling implemented for civ_series

File:
1 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);
Note: See TracChangeset for help on using the changeset viewer.