source: trunk/src/find_field_cells.m @ 1041

Last change on this file since 1041 was 1041, checked in by sommeria, 6 years ago

psmn parameters updated + miscenaleous

File size: 17.1 KB
[581]1%'find_field_cells': analyse the field structure for input in uvmat functions, grouping the variables  into 'fields' with common coordinates
[530]3% function  [CellInfo,NbDim,errormsg]=find_field_cells(Data)
[530]6% CellInfo: cell of structures describing field cells
7%     .CoordType:  type of coordinates for each field cell = 'scattered','grid','tps';
8%     .CoordIndex: array of the indices of the variables representing the coordinates (in the order z,y,x)
9%     .CoordSize: array of the nbre of values for each  coordinate in a grid, nbre of points in the unstructured case
[581]10%     .NbCentres_tps:
[530]11%     .SubRange_tps
12%     .VarIndex: arrays of the variable indices in the field cell
13%     .VarIndex_ancillary: indices of ancillary variables
14%              _color : color image, the last index, which is not a coordinate variable, represent the 3 color components rgb
15%              _discrete: like scalar, but set of data points without continuity, represented as dots in a usual plot, instead of continuous lines otherwise
16%              _errorflag: index of error flag
17%              _image   : B/W image, (behaves like scalar)
18%              _vector_x,_y,_z: indices of variables giving the vector components x, y, z
19%              _warnflag: index of warnflag   
[581]20%      .ProjModeRequest= 'interp_lin', 'interp_tps' indicate whether lin interpolation  or derivatives (tps) is needed to calculate the requested field
[576]21%      .FieldName = operation to be performed to finalise the field cell after projection
[530]22%      .SubCheck=0 /1 indicate that the field must be substracted (second  entry in uvmat)
23% NbDim: array with the length of CellVarIndex, giving the space dimension of each field cell
[514]24% errormsg: error message
26% INPUT:
27% Data: structure representing fields, output of check_field_structure
[575]28%            .ListGlobalAttributes
[514]29%            .ListVarName: cell array listing the names (cahr strings) of the variables
30%            .VarDimName: cell array of cells containing the set of dimension names for each variable of .ListVarName
31%            .VarAttribute: cell array of structures containing the variable attributes:
32%                     .VarAttribute{ilist}.key=value, where ilist is the variable
33%                     index, key is the name of the attribute, value its value (char string or number)
[575]34%            .Attr1, .Attr2....
35%        case of actual data:
36%            .Var1, .Var2...
37%        case of data structure, without the  data themselves
38%            .LisDimName: list of dimension names
39%            .DimValue: list of corresponding dimension values
[514]40% HELP:
41% to get the dimensions of arrays common to the field #icell
42%         VarIndex=CellVarIndex{icell}; % list of variable indices
43%         DimIndex=Data.VarDimIndex{VarIndex(1)} % list of dimensions for each variable in the cell #icell
[1027]46% Copyright 2008-2018, LEGI UMR 5519 / CNRS UGA G-INP, Grenoble, France
48%   Joel.Sommeria - Joel.Sommeria (A)
50%     This file is part of the toolbox UVMAT.
[514]52%     UVMAT is free software; you can redistribute it and/or modify
[809]53%     it under the terms of the GNU General Public License as published
54%     by the Free Software Foundation; either version 2 of the license,
55%     or (at your option) any later version.
[514]57%     UVMAT is distributed in the hope that it will be useful,
58%     but WITHOUT ANY WARRANTY; without even the implied warranty of
[809]60%     GNU General Public License (see LICENSE.txt) for more details.
[514]63%    group the variables  into 'fields' with common dimensions
[530]65function [CellInfo,NbDim,errormsg]=find_field_cells(Data)
[530]69if ~isfield(Data,'ListVarName'), errormsg='the list of variables .ListVarName is missing';return;end
70if ~isfield(Data,'VarDimName'), errormsg='the list of dimensions .VarDimName is missing';return;end
[581]71nbvar=numel(Data.ListVarName);%number of variables in the field structure
[530]72if ~isequal(numel(Data.VarDimName),nbvar), errormsg='.ListVarName and .VarDimName have unequal length';return;end
[535]73% check the existence of variable data
75for ilist=1:numel(Data.ListVarName)
76    if ~isfield(Data,Data.ListVarName{ilist})
[1009]77        check_var=0;% dimensions of array defined, but the corresponding array is not given
[535]78        break
[530]79    end
[515]82%% role of variables and list of requested operations
84%   'ancillary','image','color','discrete','scalar','coord_tps'};% rmq vector_x_tps and vector_y_tps to be replaced by vector_x and vector_y
[514]85Role=num2cell(blanks(nbvar));%initialize a cell array of nbvar blanks
[581]86ProjModeRequest=regexprep(Role,' ',''); % fieldRequest set to '' by default
[576]87FieldName=cell(size(Role)); % fieldRequest set to {} by default
[515]88CheckSub=zeros(size(Role));% =1 for fields to substract
[514]89Role=regexprep(Role,' ','scalar'); % Role set to 'scalar' by default
90if isfield(Data,'VarAttribute')
91    for ivar=1:numel(Data.VarAttribute)
92        if isfield(Data.VarAttribute{ivar},'Role')
93            Role{ivar}=Data.VarAttribute{ivar}.Role;
94        end
[581]95        if isfield(Data.VarAttribute{ivar},'ProjModeRequest')
96            ProjModeRequest{ivar}=Data.VarAttribute{ivar}.ProjModeRequest;
[515]97        end
[576]98        if isfield(Data.VarAttribute{ivar},'FieldName')
99            FieldName{ivar}=Data.VarAttribute{ivar}.FieldName;
[515]100        end
101        if isfield(Data.VarAttribute{ivar},'CheckSub')
102            CheckSub(ivar)=Data.VarAttribute{ivar}.CheckSub;
103        end
[514]104    end
[530]107%% find scattered (unstructured) coordinates
[674]108ivar_coord_x=find(strcmp('coord_x',Role));%find variables with Role='coord_x'
[674]113% loop on unstructured coordinate x -> different field cells
[530]114for icell=1:numel(ivar_coord_x)
[674]115    DimCell=Data.VarDimName{ivar_coord_x(icell)};% cell of dimension names for ivar_coord_x(icell)
116    if ischar(DimCell),DimCell={DimCell};end % transform char to cell for a single dimension
117    % look for variables sharing dimension(s) with ivar_coord_x(icell)
[530]118    check_cell=zeros(numel(DimCell),nbvar);
[1024]119    for idim=1:numel(DimCell);% for each variable with role coord_x, look at which other variables contain the same dimension
[530]120        for ivar=1:nbvar
121            check_cell(idim,ivar)=max(strcmp(DimCell{idim},Data.VarDimName{ivar}));
122        end
[514]123    end
[530]124    check_cell=sum(check_cell,1)==numel(DimCell);%logical array=1 for variables belonging to the current cell
[674]125    VarIndex=find(check_cell);% list of detected variable indices
[530]126    if ~(numel(VarIndex)==1 && numel(DimCell)==1)% exclude case of isolated coord_x variable (treated later)
127        if ~(numel(VarIndex)==1 && numel(DimCell)>1)% a variable is associated to coordinate
128            CellInfo{icell}.CoordIndex=ivar_coord_x(icell);
129            % size of coordinate var
[535]130            if check_var
131                CellInfo{icell}.CoordSize=numel(Data.(Data.ListVarName{ivar_coord_x(icell)}));
132            else
[530]133                for idim=1:numel(DimCell)
[674]134                    check_index= strcmp(DimCell{idim},Data.ListDimName);
135                    CellInfo{icell}.CoordSize(idim)=Data.DimValue(check_index);
[530]136                end
137                CellInfo{icell}.CoordSize=prod(CellInfo{icell}.CoordSize);
138            end
[1041]139            ind_scalar=find(strcmp('scalar',Role(VarIndex)));
140            ind_vector_x=find(strcmp('vector_x',Role(VarIndex)));
141            ind_vector_y=find(strcmp('vector_y',Role(VarIndex)));
[530]142            ind_y=find(strcmp('coord_y',Role(VarIndex)));
[1041]143            if numel([ind_scalar ind_vector_x ind_vector_y])==0           
144%             if numel(VarIndex)==2||isempty(ind_y)% no variable, except possibly y
[530]145                NbDim(icell)=1;
146            else
147                CellInfo{icell}.CoordType='scattered';
148                ind_z=find(strcmp('coord_z',Role(VarIndex)));
149                if numel(VarIndex)==3||isempty(ind_z)% no z variable, except possibly as a fct z(x,y)
150                    CellInfo{icell}.CoordIndex=[VarIndex(ind_y) CellInfo{icell}.CoordIndex];
151                    NbDim(icell)=2;
152                else
153                    CellInfo{icell}.CoordIndex=[VarIndex(ind_z) CellInfo{icell}.CoordIndex];
154                    NbDim(icell)=3;
155                end
156            end
[514]157        end
[530]158        CellInfo{icell}.VarIndex=VarIndex;
159        check_select=check_select|check_cell;
[514]160    end
163%% look for tps coordinates
[535]164ivar_remain=find(~check_select);% indices of remaining variables (not already selected)
[530]165check_coord_tps= strcmp('coord_tps',Role(~check_select));
[535]166ivar_tps=ivar_remain(check_coord_tps);% variable indices corresponding to tps coordinates
168% loop on the tps coordinate sets
[530]169for icell_tps=1:numel(ivar_tps)
[581]170    check_cell=zeros(1,nbvar);% =1 for the variables selected in the current cell
[582]171    check_cell(ivar_tps(icell_tps))=1;% mark the coordinate variable as selected
[535]172    DimCell=Data.VarDimName{ivar_tps(icell_tps)};% dimension names for the current tps coordinate variable
173    icell=numel(CellInfo)+icell_tps; % new field cell index
174    CellInfo{icell}.CoordIndex=ivar_tps(icell_tps);% index of the  tps coordinate variable
[530]175    if numel(DimCell)==3
176        VarDimName=Data.VarDimName(~check_select);
177        for ivardim=1:numel(VarDimName)
178            if strcmp(VarDimName{ivardim},DimCell{3})
[581]179                CellInfo{icell}.NbCentres_tps= ivar_remain(ivardim);% nbre of sites for each tps subdomain
180                check_cell(ivar_remain(ivardim))=1;% mark the variable as selected
[675]181            elseif strcmp(VarDimName{ivardim}{1},DimCell{2}) && numel(VarDimName{ivardim})>=3 && strcmp(VarDimName{ivardim}{3},DimCell{3})
[581]182                CellInfo{icell}.SubRange_tps=ivar_remain(ivardim);% subrange definiton for tps
183                check_cell(ivar_remain(ivardim))=1;% mark the variable as selected
184            elseif strcmp(VarDimName{ivardim}{1},DimCell{1}) && strcmp(VarDimName{ivardim}{2},DimCell{3})% variable
185                check_cell(ivar_remain(ivardim))=1;% mark the variable as selected
[530]186            end
187        end
[514]188    end
[530]189    CellInfo{icell}.CoordType='tps';
190    CellInfo{icell}.VarIndex=find(check_cell);
[535]191    if check_var
192        NbDim(icell)=size(Data.(Data.ListVarName{CellInfo{icell}.CoordIndex}),2);
193        CellInfo{icell}.CoordSize=size(Data.(Data.ListVarName{CellInfo{icell}.CoordIndex}),1)*size(Data.(Data.ListVarName{CellInfo{icell}.CoordIndex}),3);
194    else
[530]195        check_index_1= strcmp(DimCell{1},Data.ListDimName);
196        check_index_2= strcmp(DimCell{2},Data.ListDimName);
197        NbDim(icell)=Data.DimValue(check_index_2);
[648]198        if numel(DimCell)>=3
199        check_index_3= strcmp(DimCell{3},Data.ListDimName); 
[530]200        CellInfo{icell}.CoordSize=Data.DimValue(check_index_1)*Data.DimValue(check_index_3);
[648]201        end
[515]202    end
[530]203    check_select=check_select|check_cell;
[535]206%% look for coordinate variables and corresponding gridded data:
207% coordinate variables are variables associated with a single dimension, defining the coordinate values
208% two cases: 1)the coordiante variable represents the set of coordiante values
209%            2)the coordinate variable contains only two elements, representing the coordinate bounds for the dimension with the same name as the cordinate
210ivar_remain=find(~check_select);% indices of remaining variables, not already taken into account
[890]211ListVarName=Data.ListVarName(~check_select);%list of names of remaining variables
[530]212VarDimName=Data.VarDimName(~check_select);%dimensions of remaining variables
[575]213check_coord_select= cellfun(@numel,VarDimName)==1|cellfun(@ischar,VarDimName)==1;% find remaining variables with a single dimension
[690]214check_coord_select=check_coord_select & ~strcmp('ancillary',Role(~check_select));% do not select ancillary variables as coordinates
[580]215ListCoordIndex=ivar_remain(check_coord_select);% indices of remaining variables with a single dimension
216ListCoordName=ListVarName(check_coord_select);% corresponding names of remaining variables with a single dimension
217ListDimName=VarDimName(check_coord_select);% dimension names of remaining variables with a single dimension
219%remove redondant variables -> keep only one variable per dimension
[530]221for idim=1:numel(ListDimName)
[890]222    prev_ind=find(strcmp(ListDimName{idim},ListDimName(1:idim-1)));% check whether the dimension is already taken into account
223    if ~isempty(prev_ind)% in case of multiple coord variable
[535]224        if strcmp(ListCoordName{idim},ListDimName{idim}) %variable with the same name as the coordinate taken in priority
[890]225            check_keep(prev_ind)=0;% choose a variable with the same name as coordinate in priority
[530]226        else
[890]227           check_keep(idim)=0; %keep the first coordiante variable found
[514]228        end
229    end
[581]231ListCoordIndex=ListCoordIndex(check_keep);% list of coordinate variable indices
[535]232ListCoordName=ListCoordName(check_keep);% list of coordinate variable names
233ListDimName=ListDimName(check_keep);% list of coordinate dimension names
[535]235% determine dimension sizes
[530]237for ilist=1:numel(ListCoordIndex)
238    if iscell(ListDimName{ilist})
239        ListDimName(ilist)=ListDimName{ilist};%transform cell to string
[514]240    end
[535]241    if check_var% if the list of dimensions has been directly defined, no variable data available
242        CoordSize(ilist)=numel(Data.(ListCoordName{ilist}));% number of elements in the variable corresponding to the dimension #ilist
[530]243    else
[535]244        check_index= strcmp(ListDimName{ilist},Data.ListDimName);% find the  index in the list of dimensions
245        CoordSize(ilist)=Data.DimValue(check_index);% find the  corresponding dimension value
[530]246    end
[535]247    if CoordSize(ilist)==2% case of uniform grid coordinate defined by lower and upper bounds only
248        ListDimName{ilist}=ListCoordName{ilist};% replace the dimension name by the coordinate variable name
[514]249    end
[675]252%% group the remaining variables in cells sharing the same coordinate variables
[890]256for ivardim=1:numel(VarDimName) % loop at the list of dimensions for the remaining variables
[535]257    DimCell=VarDimName{ivardim};% dimension names of the current variable
[530]258    if ischar(DimCell), DimCell={DimCell}; end %transform char to cell if needed
259    DimIndices=[];
260    for idim=1:numel(DimCell)
[747]261        ind_dim=find(strcmp(DimCell{idim},ListDimName));%find the dim index in the list of dimensions ListDimName
[530]262        if ~isempty(ind_dim)
[747]263            DimIndices=[DimIndices ind_dim]; %update the list of dim indices included in DimCell
[944]264            if check_var & CoordSize(ind_dim)==2 % determine the size of the coordinate in case of coordinate definition limited to lower and upper bounds
[543]265                if isvector(Data.(ListVarName{ivardim}))
[535]266                    if numel(Data.(ListVarName{ivardim}))>2
267                        CoordSize(ind_dim)=numel(Data.(ListVarName{ivardim}));
268                    end
269                else
270                    CoordSize(ind_dim)=size(Data.(ListVarName{ivardim}),idim);
271                end
272            end
[514]273        end
[530]274    end
[535]275    % look for cells of variables with the same coordinate variables
[530]276    check_previous=0;
277    for iprev=1:numel(NewCellInfo)
278        if isequal(DimIndices,NewCellDimIndex{iprev})
279            check_previous=1;
280            NewCellInfo{iprev}.VarIndex=[NewCellInfo{iprev}.VarIndex ivar_remain(ivardim)];%append the current variable index to the found field cell
281            break
[514]282        end
[530]283    end
[535]284    % create a new cell if no previous one contains the coordinate variables
[530]285    if ~check_previous
286        nbcell=numel(NewCellInfo)+1;
287        NewCellDimIndex{nbcell}=DimIndices;
288        NewCellInfo{nbcell}.VarIndex=ivar_remain(ivardim);% create a new field cell with the current variable index
289        NewNbDim(nbcell)=numel(DimIndices);
290        NewCellInfo{nbcell}.CoordType='grid';
[747]291        NewCellInfo{nbcell}.CoordSize=CoordSize(DimIndices);
[530]292        NewCellInfo{nbcell}.CoordIndex=ListCoordIndex(DimIndices);
293    end
295NbDim=[NbDim NewNbDim];
296CellInfo=[CellInfo NewCellInfo];
[569]298%% suppress empty cells or cells with a single coordinate variable
300for icell=1:numel(check_remove)
[1023]301    if isempty(CellInfo{icell})||(numel(CellInfo{icell}.VarIndex)==1 && numel(check_coord)>=icell && check_coord(icell))
[569]302        check_remove(icell)=1;
303    end
308%% document roles of non-coordinate variables
309for icell=1:numel(CellInfo)
310    VarIndex=CellInfo{icell}.VarIndex;
311    for ivar=VarIndex
312        if isfield(CellInfo{icell},['VarIndex_' Role{ivar}])
313            CellInfo{icell}.(['VarIndex_' Role{ivar}])=[CellInfo{icell}.(['VarIndex_' Role{ivar}]) ivar];
314        else
315            CellInfo{icell}.(['VarIndex_' Role{ivar}])= ivar;
[514]316        end
[581]317        if ~isempty(ProjModeRequest{ivar})
318            CellInfo{icell}.ProjModeRequest=ProjModeRequest{ivar};
[514]319        end
[576]320        if ~isempty(FieldName{ivar})
321            CellInfo{icell}.FieldName=FieldName{ivar};
[514]322        end
[530]323        if CheckSub(ivar)==1
324            CellInfo{icell}.CheckSub=1;
[514]325        end
326    end
Note: See TracBrowser for help on using the repository browser.