Changeset 1180 for trunk/src/series/sub_background.m

Timestamp:

Mar 27, 2025, 5:59:21 PM (3 weeks ago)

Author:

sommeria

Message:

various bugs repaired,Relabeling frames installed for multitif

File:

• trunk/src/series/sub_background.m (modified) (5 diffs)

trunk/src/series/sub_background.m

@@ -78 +78 @@
     ParamOut.VelType='off';% menu for selecting the velocity type (options 'off'/'one'/'two',  'off' by default)
     ParamOut.FieldName='off';% menu for selecting the field (s) in the input file(options 'off'/'one'/'two', 'off' by default)
-    ParamOut.FieldTransform = 'off';%can use a transform function
-    ParamOut.ProjObject='off';%cannot use projection object(option 'off'/'on',
-    ParamOut.Mask='on';%can use mask option   (option 'off'/'on', 'off' by default)
+    ParamOut.FieldTransform = 'off';%can use a transform function (option 'off'/'on','off' by default)
+    ParamOut.ProjObject='off';%cannot use projection object(option 'off'/'on','off' by default)
+    ParamOut.Mask='off';%cannot use mask option   (option 'off'/'on', 'off' by default)
     ParamOut.OutputDirExt='.sback';%set the output dir extension
     ParamOut.OutputFileMode='NbInput';% '=NbInput': 1 output file per input file index, '=NbInput_i': 1 file per input file index i, '=NbSlice': 1 file per slice
@@ -250 +250 @@
 %% Input file info
 if CheckRelabel
-      [RootFileOut,FileIndexString]=index2filename(Param.FileSeries,Param.IndexRange.first_i,j_indices(1),NbField_j);
-       FirstFileName=fullfile(RootPath,SubDir,[RootFileOut FileIndexString FileExt]);
-else
-FirstFileName=fullfile_uvmat(RootPath,SubDir,RootFile,FileExt,NomType,Param.IndexRange.first_i,[],j_indices(1));%get first file name
-RootFileOut=RootFile;
+    [RootFileOut,frame_index]=index2filename(Param.FileSeries,Param.IndexRange.first_i,j_indices(1),NbField_j);
+    FirstFileName=fullfile(RootPath,SubDir,RootFileOut);
+else
+    FirstFileName=fullfile_uvmat(RootPath,SubDir,RootFile,FileExt,NomType,Param.IndexRange.first_i,[],j_indices(1));%get first file name
+    RootFileOut=RootFile;
 end
 [FileInfo,MovieObject]=get_file_info(FirstFileName);
 FileType=FileInfo.FileType;
-if isfield(FileInfo,'NumberOfFrames') && FileInfo.NumberOfFrames >1
-    if isempty(regexp(NomType,'1$', 'once'))% no file indexing
-        frame_index=i_indices;% the index i denotes the frame number in a movie, no index j
+if ~CheckRelabel
+    if isfield(FileInfo,'NumberOfFrames') && FileInfo.NumberOfFrames >1
+        if isempty(regexp(NomType,'1$', 'once'))% no file indexing
+            frame_index=i_indices;% the index i denotes the frame number in a movie, no index j
+        else
+            frame_index=j_indices;% the index j denotes the frame number in a movie
+            MovieObject=[]; %not a single video object
+        end
     else
-        frame_index=j_indices;% the index j denotes the frame number in a movie
-        MovieObject=[]; %not a single video object
-    end
-else
-    frame_index=ones(1,nbfield);
+        frame_index=ones(1,nbfield);
+    end
 end
 
@@ -288 +290 @@
 
 %% prealocate memory for the sliding background
-Ak=zeros(FileInfo.Height,FileInfo.Width,nbaver_ima,['uint' num2str(FileInfo.BitDepth)]); %prealocate memory    
+Ak=zeros(FileInfo.Height,FileInfo.Width,nbaver_ima,['uint' num2str(FileInfo.BitDepth)]); %prealocate memory
 
 %% selection of frame indices
-if Param.ActionInput.CheckVolume 
+if Param.ActionInput.CheckVolume
     nbfield=floor(nbfield/NbSlice_j)*NbSlice_j;% truncate the total number of frames in case of incomplete series
     indselect=1:nbfield;
-     indselect=reshape(indselect,NbSlice_j,[]);
-      NbSlice=NbSlice_j;
-else
-       NbSlice=NbSlice_i;
+    indselect=reshape(indselect,NbSlice_j,[]);
+    NbSlice=NbSlice_j;
+else
+    NbSlice=NbSlice_i;
     nbfield=floor(nbfield/NbSlice)*NbSlice;% truncate the total number of frames in case of incomplete series
     indselect=reshape(1:nbfield,NbSlice,[]);
     for j_slice=1:NbSlice
-    indselect(j_slice,:)=j_slice:NbSlice:nbfield;% select file indices of the slice
+        indselect(j_slice,:)=j_slice:NbSlice:nbfield;% select file indices of the slice
     end
 end
@@ -313 +315 @@
         %filename=filecell{1,ifile};
         if CheckRelabel
-            [RootFile,FileIndexString,FrameIndex]=index2filename(Param.FileSeries,i_indices(ifile),j_indices(ifile),NbField_j);
-            filename=fullfile(RootPath,SubDir,[RootFile FileIndexString FileExt]);
+            [filename,FrameIndex]=index2filename(Param.FileSeries,i_indices(ifile),j_indices(ifile),NbField_j);
+             filename=fullfile(RootPath,SubDir,filename);
         else
             filename=fullfile_uvmat(RootPath,SubDir,RootFile,FileExt,NomType,i_indices(ifile),[],j_indices(ifile));
@@ -368 +370 @@
                 ifile=indselect(j_slice,ifield+iburst+step*halfnbaver);
                 if CheckRelabel
-                    [RootFile,FileIndexString,FrameIndex]=index2filename(Param.FileSeries,i_indices(ifile),j_indices(ifile),NbField_j);
-                    filename=fullfile(RootPath,SubDir,[RootFile FileIndexString FileExt]);
+                    [filename,FrameIndex]=index2filename(Param.FileSeries,i_indices(ifile),j_indices(ifile),NbField_j);
+                    filename=fullfile(RootPath,SubDir,filename);
                 else
                     filename=fullfile_uvmat(RootPath,SubDir,RootFile,FileExt,NomType,i_indices(ifile),[],j_indices(ifile));