Changeset 1099 for trunk/src/series


Ignore:
Timestamp:
Apr 14, 2021, 4:02:22 PM (3 years ago)
Author:
sommeria
Message:

extrtact_multitif cleaned

Location:
trunk/src/series
Files:
2 deleted
1 edited

Legend:

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

    r1093 r1099  
    1 %'ima2netcdf': read image series and transform to netcdf
     1%'extract_multitif': read image series from PCO cameras (tiff image series) and write .png images
    22%------------------------------------------------------------------------
    3 
    43   
    5 % function ParamOut=ima2netcdf(Param)
     4% function ParamOut=extract_multitif(Param)
    65%
    76%%%%%%%%%%% GENERAL TO ALL SERIES ACTION FCTS %%%%%%%%%%%%%%%%%%%%%%%%%%%
     
    5958%=======================================================================
    6059
    61 function ParamOut=extract_multitif_parallel(Param)
     60function ParamOut=extract_multitif(Param)
    6261
    6362%%%%%%%%%%%%%%%%%    INPUT PREPARATION MODE (no RUN)    %%%%%%%%%%%%%%%%%
     
    7473    ParamOut.OutputFileMode='NbSlice';% '=NbInput': 1 output file per input file index, '=NbInput_i': 1 file per input file index i, '=NbSlice': 1 file per slice
    7574    ParamOut.CheckOverwriteVisible='on'; % manage the overwrite of existing files (default=1)
    76     ParamOut.CPUTime=7;% expected time for writting one image ( in minute)
     75    ParamOut.CPUTime=10;% expected time for writting the output of one source image ( in minute)
    7776    %% root input file(s) and type
    7877    % check the existence of the first file in the series
     
    126125ImagesPerLevel=size(XmlInput.Time,2)-1;%100;%use the xmlinformation to get the nbre of j indices
    127126
    128 %% create the xml file of PCO camera if it does not exist
    129 Newxml=fullfile(Param.InputTable{1,1},[Param.InputTable{1,2} '.xml']);
    130 if ~exist(Newxml,'file')
    131 XmlInput.Camera.CameraName='PCO';
    132 XmlInput=rmfield(XmlInput,'Time');
    133 XmlInput=rmfield(XmlInput,'TimeUnit');
    134 t=struct2xml(XmlInput);
    135 t=set(t,1,'name','ImaDoc');
    136 save(t,Newxml);
     127%% create the xml file for timing if it does not exist : example to adapt
     128TEST=0;
     129if TEST
     130    count0=14;
     131    Dtj=0.05;% time interval between frames
     132    ImagesPerLevel=455;% total number of images per position, ImagesPerLevel-Nbj images skiiped during motion between two positions
     133    Nbj=390; %Nbre of images kept at a given position
     134    Dti=Dtj*ImagesPerLevel;
     135    NbLevel=11;
     136    NbScan=3;
     137    TimeReturn=268.5; %time needed to return back to the first position (in sec)
     138    NbReturn=round(TimeReturn/Dtj);
     139    NbSkipReturn=NbReturn+1-NbLevel*ImagesPerLevel;
     140   
     141    Newxml=fullfile(Param.InputTable{1,1},[Param.InputTable{1,2} '.xml']);
     142    if ~exist(Newxml,'file')
     143        XmlInput.Camera.CameraName='PCO';
     144        XmlInput.Camera.TimeUnit='s';
     145        XmlInput.Camera.BurstTiming.FrameFrequency=1;
     146        XmlInput.Camera.BurstTiming.Time=0;% for 200
     147        XmlInput.Camera.BurstTiming.Dtj=Dtj;
     148        XmlInput.Camera.BurstTiming.NbDtj=Nbj-1;
     149        XmlInput.Camera.BurstTiming.Dti=Dti;
     150        XmlInput.Camera.BurstTiming.NbDti=NbLevel-1;
     151        XmlInput.Camera.BurstTiming.Dtk=TimeReturn;
     152        XmlInput.Camera.BurstTiming.NbDtk=NbScan-1;
     153        t=struct2xml(XmlInput);
     154        t=set(t,1,'name','ImaDoc');
     155        save(t,Newxml);
     156    end
    137157end
    138158
Note: See TracChangeset for help on using the changeset viewer.