source: trunk/src/read_field.m @ 1071

Last change on this file since 1071 was 1071, checked in by g7moreau, 5 years ago
  • Update COPYRIGHT
File size: 14.3 KB
RevLine 
[497]1%'read_field': read the fields from files in different formats (netcdf files, images, video)
[181]2%--------------------------------------------------------------------------
[450]3%  function [Field,ParamOut,errormsg] = read_field(FileName,FileType,ParamIn,num)
[181]4%
5% OUTPUT:
6% Field: matlab structure representing the field
7% ParamOut: structure representing parameters:
8%        .FieldName; field name
9%        .VelType
[1045]10%        .CivStage: stage of civx processing (=0, not Civx, =1 (civ1), =2  (fix1)....
[181]11%        .Npx,.Npy: for images, nbre of pixels in x and y
12% errormsg: error message, ='' by default
13%
14%INPUT
[466]15% FileName: name of the input file
[784]16% FileType: type of file, as determined by the function get_file_info.m
[450]17% ParamIn: movie object or Matlab structure of input parameters
[445]18%     .FieldName: name (char string) of the input field (for Civx data)
19%     .VelType: char string giving the type of velocity data ('civ1', 'filter1', 'civ2'...)
[181]20%     .ColorVar: variable used for vector color
21%     .Npx, .Npy: nbre of pixels along x and y (used for .vol input files)
[1045]22%     .TimeDimName: name of the dimension considered as 'time', selected index value then set by input 'num'
[466]23% num: frame number for movies
[497]24%
25% see also read_image.m,read_civxdata.m,read_civdata.m,
[466]26
[809]27%=======================================================================
[1071]28% Copyright 2008-2020, LEGI UMR 5519 / CNRS UGA G-INP, Grenoble, France
[809]29%   http://www.legi.grenoble-inp.fr
30%   Joel.Sommeria - Joel.Sommeria (A) legi.cnrs.fr
31%
32%     This file is part of the toolbox UVMAT.
33%
34%     UVMAT is free software; you can redistribute it and/or modify
35%     it under the terms of the GNU General Public License as published
36%     by the Free Software Foundation; either version 2 of the license,
37%     or (at your option) any later version.
38%
39%     UVMAT is distributed in the hope that it will be useful,
40%     but WITHOUT ANY WARRANTY; without even the implied warranty of
41%     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
42%     GNU General Public License (see LICENSE.txt) for more details.
43%=======================================================================
44
[450]45function [Field,ParamOut,errormsg] = read_field(FileName,FileType,ParamIn,num)
[748]46%% default output and check input
[181]47Field=[];
[354]48if ~exist('num','var')
49    num=1;
50end
[1033]51if isempty(num)
52    num=1;
53end
[354]54if ~exist('ParamIn','var')
55    ParamIn=[];
56end
57ParamOut=ParamIn;%default
[181]58errormsg='';
[1022]59if isempty(regexp(FileName,'^http://'))&& ~exist(FileName,'file')
[747]60    errormsg=['input file ' FileName ' does not exist'];
[635]61    return
62end
[334]63A=[];
[575]64InputField={};
65check_colorvar=0;
[526]66if isstruct(ParamIn)
[575]67    if isfield(ParamIn,'FieldName')
68        if ischar(ParamIn.FieldName)
[576]69            InputField={ParamIn.FieldName};
[575]70        else
71            InputField= ParamIn.FieldName;
72        end
[526]73    end
[654]74    check_colorvar=zeros(size(InputField));
[684]75    if isfield(ParamIn,'ColorVar')&&~isempty(ParamIn.ColorVar)
[526]76        InputField=[ParamIn.FieldName {ParamIn.ColorVar}];
[654]77        check_colorvar(numel(InputField))=1;
[526]78    end
[521]79end
[575]80
[334]81%% distingush different input file types
[527]82switch FileType
[748]83    case 'civdata'% new format for civ results
[527]84        [Field,ParamOut.VelType,errormsg]=read_civdata(FileName,InputField,ParamIn.VelType);
[907]85        if ~isempty(errormsg),errormsg=['read_civdata / ' errormsg];return,end
[527]86        ParamOut.CivStage=Field.CivStage;
[958]87    case 'pivdata_fluidimage'
88        [Field,ParamOut.VelType,errormsg]=read_pivdata_fluidimage(FileName,InputField,ParamIn.VelType);
89        ParamOut.CivStage=Field.CivStage;
[748]90    case 'civx'% old (obsolete) format for civ results
[527]91        ParamOut.FieldName='velocity';%Civx data found, set .FieldName='velocity' by default
92        [Field,ParamOut.VelType,errormsg]=read_civxdata(FileName,InputField,ParamIn.VelType);
93        if ~isempty(errormsg),errormsg=['read_civxdata / ' errormsg];return,end
94        ParamOut.CivStage=Field.CivStage;
[748]95    case 'netcdf'% general netcdf file (not recognized as civ)
[1045]96        ListVarName={};
[675]97        Role={};
98        ProjModeRequest={};
[1045]99        %ListInputField={};
100        %checkU=0;
101        %checkV=0;
[1041]102        % scan the list InputField
[1045]103        Operator=cell(1,numel(InputField));
104        InputVar=cell(1,numel(InputField));
[527]105        for ilist=1:numel(InputField)
[748]106            % look for input variables to read
[527]107            r=regexp(InputField{ilist},'(?<Operator>(^vec|^norm))\((?<UName>.+),(?<VName>.+)\)$','names');
[654]108            if isempty(r)%  no operator used
[1045]109                ListVarName=[ListVarName InputField(ilist)];%append the variable name
110                %InputVar{ilist}=InputField(ilist);
111                %ListInputField=[ListInputField InputField(ilist)];
112                if check_colorvar(ilist)% case of field used for vector color
113                    Role{numel(ListVarName)}='ancillary';% not projected with interpolation
114                    ProjModeRequest{numel(ListVarName)}='';
115                else
116                    Role{numel(ListVarName)}='scalar';
117                    ProjModeRequest{numel(ListVarName)}='interp_lin';%scalar field (requires interpolation for plot)
[681]118                end
[1045]119                Operator{numel(ListVarName)}='';
120            else  % an operator 'vec' or 'norm' is used
121                ListVarName=[ListVarName {r.UName}]; % append the variable in the list if not previously listed
122                if  strcmp(r.Operator,'norm')
123                    if check_colorvar(ilist)
124                    Role=[Role {'ancillary'}];
125                    else
126                       Role=[Role {'scalar'}];
[907]127                    end
[654]128                else
[1045]129                     Role=[Role {'vector_x'}];
[654]130                end
[1045]131                %ListInputField=[ListInputField InputField(ilist)];
132                ListVarName=[ListVarName {r.VName}];% append the variable in the list if not previously listed
133                Role=[Role {'vector_y'}];
134                %ListInputField=[ListInputField InputField(ilist)];
135                Operator{numel(ListVarName)-1}=r.Operator;
136                Operator{numel(ListVarName)}='';           
[667]137                if ~check_colorvar(ilist) && strcmp(r.Operator,'norm')
[1045]138                    ProjModeRequest{numel(ListVarName)}='interp_lin';%scalar field (requires interpolation for plot)
139                    ProjModeRequest{numel(ListVarName)-1}='interp_lin';%scalar field (requires interpolation for plot)
[667]140                else
[1045]141                    ProjModeRequest{numel(ListVarName)}='';
142                    ProjModeRequest{numel(ListVarName)-1}='';
[667]143                end
[493]144            end
[527]145        end
[748]146        if ~isfield(ParamIn,'Coord_z')
147            ParamIn.Coord_z=[];
148        end
149        NbCoord=~isempty(ParamIn.Coord_x)+~isempty(ParamIn.Coord_y)+~isempty(ParamIn.Coord_z);
[648]150        if isfield(ParamIn,'TimeDimName')% case of reading of a single time index in a multidimensional array
[1045]151            [Field,var_detect,ichoice,errormsg]=nc2struct(FileName,'TimeDimName',ParamIn.TimeDimName,num,[ParamIn.Coord_x ParamIn.Coord_y ParamIn.Coord_z ListVarName]);
[747]152        elseif isfield(ParamIn,'TimeVarName')% case of reading of a single time  in a multidimensional array
[1045]153            [Field,var_detect,ichoice,errormsg]=nc2struct(FileName,'TimeVarName',ParamIn.TimeVarName,num,[ParamIn.Coord_x ParamIn.Coord_y ParamIn.Coord_z ListVarName]);
[907]154            if numel(num)~=1
[1001]155                NbCoord=NbCoord+1;% adds time coordinate, except if a single time has been selected
[907]156            end
[648]157        else
[1045]158            [Field,var_detect,ichoice,errormsg]=nc2struct(FileName,[ParamIn.Coord_x ParamIn.Coord_y ParamIn.Coord_z ListVarName]);
[648]159        end
[747]160        if ~isempty(errormsg)
[534]161            return
162        end
[1001]163        CheckStructured=1;
[1045]164        %scan all the variables
165        NbCoord=0;
166        if ~isempty(ParamIn.Coord_x)
167            index_Coord_x=find(strcmp(ParamIn.Coord_x,Field.ListVarName));
168            Field.VarAttribute{index_Coord_x}.Role='coord_x';%
169            NbCoord=NbCoord+1;
[675]170        end
[1045]171        if ~isempty(ParamIn.Coord_y)
172            if ischar(ParamIn.Coord_y)
173                index_Coord_y=find(strcmp(ParamIn.Coord_y,Field.ListVarName));
174                Field.VarAttribute{index_Coord_y}.Role='coord_y';%
175                NbCoord=NbCoord+1;
176            else
177                for icoord_y=1:numel(ParamIn.Coord_y)
178                    index_Coord_y=find(strcmp(ParamIn.Coord_y{icoord_y},Field.ListVarName));
179                    Field.VarAttribute{index_Coord_y}.Role='coord_y';%
180                    NbCoord=NbCoord+1;
[1001]181                end
182            end
[1002]183        end
[1045]184        NbDim=1;
185        if ~isempty(ParamIn.Coord_z)
186            index_Coord_z=find(strcmp(ParamIn.Coord_z,Field.ListVarName));
187            Field.VarAttribute{index_Coord_z}.Role='coord_z';%
188            NbCoord=NbCoord+1;
189            NbDim=3;
190        elseif ~isempty(ParamIn.FieldName)
191            NbDim=2;
192        end
[684]193        NormName='';
194        UName='';
195        VName='';
[1045]196        if numel(Field.ListVarName)>NbCoord % if there are variables beyond coord (exclude 1 D plots)
197            VarAttribute=cell(1,numel(ListVarName));
198            for ilist=1:numel(ListVarName)
199                index_var=find(strcmp(ListVarName{ilist},Field.ListVarName),1);
200                VarDimName{ilist}=Field.VarDimName{index_var};
201                DimOrder=[];
202                if NbDim ==2
203                    DimOrder=[find(strcmp(ParamIn.Coord_y,VarDimName{ilist})) find(strcmp(ParamIn.Coord_x,VarDimName{ilist}))];
204                elseif NbDim ==3
205                    DimOrder=[find(strcmp(ParamIn.Coord_z,VarDimName{ilist}))...
206                        find(strcmp(ParamIn.Coord_y,VarDimName{ilist}))...
207                        find(strcmp(ParamIn.Coord_x,VarDimName{ilist}))];
[907]208                end
[1045]209                if ~isempty(DimOrder)
210                    Field.(ListVarName{ilist})=permute(Field.(ListVarName{ilist}),DimOrder);
211                    VarDimName{ilist}=VarDimName{ilist}(DimOrder);
[907]212                end
[1045]213                if numel(Field.VarAttribute)>=index_var
214                VarAttribute{ilist}=Field.VarAttribute{index_var};% read var attributes from input if exist
[684]215                end
[681]216            end
[1045]217            check_remove=false(1,numel(Field.ListVarName));
218            for ilist=1:numel(ListVarName)
219                VarAttribute{ilist}.Role=Role{ilist};
220                VarAttribute{ilist}.ProjModeRequest=ProjModeRequest{ilist};
221                if isfield(ParamIn,'FieldName')
222                    VarAttribute{ilist}.Operator=Operator{ilist};
[907]223                end
[1045]224                if strcmp(Operator{ilist},'norm')
225                    UName=ListVarName{ilist};
226                    VName=ListVarName{ilist+1};
227                    ListVarName{ilist}='norm';
228                    Field.norm=Field.(UName).*Field.(UName)+Field.(VName).*Field.(VName);
229                    Field.norm=sqrt(Field.norm);
230                    check_remove(ilist+1)=true;
231                    VarAttribute{ilist}.Operator='';
232                end
[769]233            end
[1045]234            ListVarName(check_remove)=[];
235            VarDimName(check_remove)=[];
236            VarAttribute(check_remove)=[];
237            Field.ListVarName=[Field.ListVarName(1:NbCoord) ListVarName];% complement the list of vqriables, which may be listed twice
238            Field.VarDimName=[Field.VarDimName(1:NbCoord) VarDimName];
239            Field.VarAttribute=[Field.VarAttribute(1:NbCoord) VarAttribute];
[684]240        end
[1045]241       
[527]242    case 'video'
243        if strcmp(class(ParamIn),'VideoReader')
244            A=read(ParamIn,num);
245        else
246            ParamOut=VideoReader(FileName);
247            A=read(ParamOut,num);
248        end
249    case 'mmreader'
250        if strcmp(class(ParamIn),'mmreader')
251            A=read(ParamIn,num);
252        else
253            ParamOut=mmreader(FileName);
254            A=read(ParamOut,num);
255        end
256    case 'vol'
257        A=imread(FileName);
258        Npz=size(A,1)/ParamIn.Npy;
259        A=reshape(A',ParamIn.Npx,ParamIn.Npy,Npz);
260        A=permute(A,[3 2 1]);
261    case 'multimage'
[907]262        %  warning 'off'
[527]263        A=imread(FileName,num);
264    case 'image'
265        A=imread(FileName);
[784]266    case 'rdvision'
[1058]267        [A,FileInfo,timestamps,errormsg]=read_rdvision(FileName,num);
[1007]268    case 'image_DaVis'
269        Input=readimx(FileName);
[1033]270        if numel(Input.Frames)==1
271            num=1;
272        end
[1007]273        A=Input.Frames{num}.Components{1}.Planes{1}';
[1033]274        for ilist=1:numel(Input.Frames{1}.Attributes)
275            if strcmp(Input.Frames{1}.Attributes{ilist}.Name,'AcqTimeSeries')
[1045]276                timestamps=str2num(Input.Frames{1}.Attributes{ilist}.Value(1:end-3))/1000000;
277                break
[1033]278            end
279        end
[980]280    case 'cine_phantom'
281        [A,FileInfo] = read_cine_phantom(FileName,num );
[784]282    otherwise
283        errormsg=[ FileType ': invalid input file type for uvmat'];
[1045]284       
[527]285end
[445]286
[334]287%% case of image
288if ~isempty(A)
[1032]289    if strcmp(FileType,'rdvision')||strcmp(FileType,'image_DaVis')
[784]290        Field.Time=timestamps;
291    end
[452]292    if isstruct(ParamOut)
[580]293        ParamOut.FieldName='image';
[452]294    end
[182]295    Npz=1;%default
[181]296    npxy=size(A);
297    Field.NbDim=2;%default
298    Field.AName='image';
[782]299    Field.ListVarName={'Coord_y','Coord_x','A'}; %
[1024]300    Field.VarAttribute{1}.Unit='pixel';
301    Field.VarAttribute{2}.Unit='pixel';
[1045]302    Field.VarAttribute{1}.Role='coord_y';
303    Field.VarAttribute{2}.Role='coord_x';
[1024]304    Field.VarAttribute{3}.Role='scalar';
[181]305    if ndims(A)==3
306        if Npz==1;%color
[782]307            Field.VarDimName={'Coord_y','Coord_x',{'Coord_y','Coord_x','rgb'}}; %
308            Field.Coord_y=[npxy(1)-0.5 0.5];
309            Field.Coord_x=[0.5 npxy(2)-0.5]; % coordinates of the first and last pixel centers
[452]310            if isstruct(ParamOut)
[580]311                ParamOut.Npx=npxy(2);% display image size on the interface
312                ParamOut.Npy=npxy(1);
[452]313            end
[181]314        else
315            Field.NbDim=3;
316            Field.ListVarName=['AZ' Field.ListVarName];
[782]317            Field.VarDimName={'AZ','Coord_y','Coord_x',{'AZ','Coord_y','Coord_x'}};
[181]318            Field.AZ=[npxy(1)-0.5 0.5];
[782]319            Field.Coord_y=[npxy(2)-0.5 0.5];
320            Field.Coord_x=[0.5 npxy(3)-0.5]; % coordinates of the first and last pixel centers
[452]321            if isstruct(ParamOut)
[580]322                ParamOut.Npx=npxy(3);% display image size on the interface
323                ParamOut.Npy=npxy(2);
324            end
[181]325        end
326    else
[782]327        Field.VarDimName={'Coord_y','Coord_x',{'Coord_y','Coord_x'}}; %
328        Field.Coord_y=[npxy(1)-0.5 0.5];
329        Field.Coord_x=[0.5 npxy(2)-0.5]; % coordinates of the first and last pixel centers
[181]330    end
331    Field.A=A;
332    Field.CoordUnit='pixel'; %used for mouse_motion
333end
334
335
[334]336
Note: See TracBrowser for help on using the repository browser.