Context Navigation

source: trunk/src/uvmat_doc/FUNCTIONS_DOC/proj_field.html @ 37

Last change on this file since 37 was 37, checked in by sommeria, 14 years ago
create_grid.fig , uvmat_doc and all the included files added
File size: 142.4 KB

Rev	Line
[37]	1	<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN"
	2	"http://www.w3.org/TR/REC-html40/loose.dtd">
	3	<html>
	4	<head>
	5	<title>Description of proj_field</title>
	6	<meta name="keywords" content="proj_field">
	7	<meta name="description" content="'proj_field': projects the field on a projection object">
	8	<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
	9	<meta name="generator" content="m2html © 2003 Guillaume Flandin">
	10	<meta name="robots" content="index, follow">
	11	<link type="text/css" rel="stylesheet" href="../m2html.css">
	12	</head>
	13	<body>
	14	<a name="_top"></a>
	15	<div><a href="../index.html">Home</a> > <a href="index.html">.</a> > proj_field.m</div>
	16
	17	<!--<table width="100%"><tr><td align="left"><a href="../index.html"><img alt="<" border="0" src="../left.png"> Master index</a></td>
	18	<td align="right"><a href="index.html">Index for . <img alt=">" border="0" src="../right.png"></a></td></tr></table>-->
	19
	20	<h1>proj_field
	21	</h1>
	22
	23	<h2><a name="_name"></a>PURPOSE <a href="#_top"><img alt="^" border="0" src="../up.png"></a></h2>
	24	<div class="box"><strong>'proj_field': projects the field on a projection object</strong></div>
	25
	26	<h2><a name="_synopsis"></a>SYNOPSIS <a href="#_top"><img alt="^" border="0" src="../up.png"></a></h2>
	27	<div class="box"><strong>function [ProjData,errormsg]=proj_field(FieldData,ObjectData,IndexObj) </strong></div>
	28
	29	<h2><a name="_description"></a>DESCRIPTION <a href="#_top"><img alt="^" border="0" src="../up.png"></a></h2>
	30	<div class="fragment"><pre class="comment">'proj_field': projects the field on a projection object
	31	--------------------------------------------------------------------------
	32	function [ProjData,errormsg]=proj_field(FieldData,ObjectData,IndexObj)
	33
	34	OUTPUT:
	35	ProjData structure containing the fields of the input field FieldData,
	36	transmitted or projected on the object, plus the additional fields
	37	.UMax, .UMin, .VMax, .VMin: min and max of velocity components in a domain
	38	.UMean,VMean: mean of the velocity components in a domain
	39	.AMin, AMax: min and max of a scalar
	40	.AMean: mean of a scalar in a domain
	41	.NbPix;
	42	.DimName= names of the matrix dimensions (matlab cell)
	43	.DimValue= values of the matricx dimensions (matlab vector, same length as .DimName);
	44	.VarName= names of the variables [ProjData.VarName {'A','AMean','AMin','AMax'}];
	45	.VarDimNameIndex= dimensions of the variables, indicated by indices in the list .DimName;
	46
	47	INPUT
	48	ObjectData: structure characterizing the projection object
	49	.Style : style of projection object
	50	.ProjMode=type of projection ;
	51	.CoordType: 'px' or 'phys' type of coordinates defining the object position
	52	.Phi angle of rotation (=0 by default)
	53	.ProjAngle=angle of projection;
	54	.DX,.DY,.DZ=increments along each coordinate
	55	.Coord(nbpoints,3): set of coordinates defining the object position;</pre></div>
	56
	57	<!-- crossreference -->
	58	<h2><a name="_cross"></a>CROSS-REFERENCE INFORMATION <a href="#_top"><img alt="^" border="0" src="../up.png"></a></h2>
	59	This function calls:
	60	<ul style="list-style-image:url(../matlabicon.gif)">
	61	<li><a href="find_field_indices.html" class="code" title="function [CellVarIndex,NbDim,VarType,errormsg]=find_field_indices(Data)">find_field_indices</a> 'find_file_indices': group the variables of a nc-formated Matlab structure into 'fields' with common dimensions</li><li><a href="griddata_uvmat.html" class="code" title="function ZI = griddata_uvmat(X,Y,Z,XI,YI)">griddata_uvmat</a> 'griddata_uvmat': function griddata_uvmat(vec2_X,vec2_Y,vec2_U,vec_X,vec_Y,'linear')</li><li><a href="warndlg_uvmat.html" class="code" title="function hwarn=warndlg_uvmat(warntext,title)">warndlg_uvmat</a> 'warndlg_uvmat': display warning message (error, warning, confirmation) in a given figure</li></ul>
	62	This function is called by:
	63	<ul style="list-style-image:url(../matlabicon.gif)">
	64	<li><a href="series.html" class="code" title="function varargout = series(varargin)">series</a> 'series': master function associated to the GUI series.m for analysis field series</li><li><a href="update_obj.html" class="code" title="function Object_out=update_obj(UvData,IndexObj,ObjectData,PlotHandles)">update_obj</a> 'update_obj': update the object graph representation and its projection field, record it in the uvmat interface</li><li><a href="uvmat.html" class="code" title="function varargout = uvmat(varargin)">uvmat</a> 'uvmat': function associated with the GUI 'uvmat.fig' for images and data field visualization</li></ul>
	65	<!-- crossreference -->
	66
	67	<h2><a name="_subfunctions"></a>SUBFUNCTIONS <a href="#_top"><img alt="^" border="0" src="../up.png"></a></h2>
	68	<ul style="list-style-image:url(../matlabicon.gif)">
	69	<li><a href="#_sub1" class="code">function [ProjData,errormsg]=proj_points(FieldData,ObjectData)</a></li><li><a href="#_sub2" class="code">function [ProjData,errormsg]=proj_patch(FieldData,ObjectData)</a></li><li><a href="#_sub3" class="code">function [ProjData,errormsg] = proj_line(FieldData, ObjectData)</a></li><li><a href="#_sub4" class="code">function [ProjData,errormsg] = proj_plane(FieldData, ObjectData)</a></li><li><a href="#_sub5" class="code">function [ProjData,errormsg]=proj_heading(FieldData,ObjectData)</a></li></ul>
	70	<h2><a name="_source"></a>SOURCE CODE <a href="#_top"><img alt="^" border="0" src="../up.png"></a></h2>
	71	<div class="fragment"><pre>0001 <span class="comment">%'proj_field': projects the field on a projection object</span>
	72	0002 <span class="comment">%--------------------------------------------------------------------------</span>
	73	0003 <span class="comment">% function [ProjData,errormsg]=proj_field(FieldData,ObjectData,IndexObj)</span>
	74	0004 <span class="comment">%</span>
	75	0005 <span class="comment">% OUTPUT:</span>
	76	0006 <span class="comment">% ProjData structure containing the fields of the input field FieldData,</span>
	77	0007 <span class="comment">% transmitted or projected on the object, plus the additional fields</span>
	78	0008 <span class="comment">% .UMax, .UMin, .VMax, .VMin: min and max of velocity components in a domain</span>
	79	0009 <span class="comment">% .UMean,VMean: mean of the velocity components in a domain</span>
	80	0010 <span class="comment">% .AMin, AMax: min and max of a scalar</span>
	81	0011 <span class="comment">% .AMean: mean of a scalar in a domain</span>
	82	0012 <span class="comment">% .NbPix;</span>
	83	0013 <span class="comment">% .DimName= names of the matrix dimensions (matlab cell)</span>
	84	0014 <span class="comment">% .DimValue= values of the matricx dimensions (matlab vector, same length as .DimName);</span>
	85	0015 <span class="comment">% .VarName= names of the variables [ProjData.VarName {'A','AMean','AMin','AMax'}];</span>
	86	0016 <span class="comment">% .VarDimNameIndex= dimensions of the variables, indicated by indices in the list .DimName;</span>
	87	0017 <span class="comment">%</span>
	88	0018 <span class="comment">%INPUT</span>
	89	0019 <span class="comment">% ObjectData: structure characterizing the projection object</span>
	90	0020 <span class="comment">% .Style : style of projection object</span>
	91	0021 <span class="comment">% .ProjMode=type of projection ;</span>
	92	0022 <span class="comment">% .CoordType: 'px' or 'phys' type of coordinates defining the object position</span>
	93	0023 <span class="comment">% .Phi angle of rotation (=0 by default)</span>
	94	0024 <span class="comment">% .ProjAngle=angle of projection;</span>
	95	0025 <span class="comment">% .DX,.DY,.DZ=increments along each coordinate</span>
	96	0026 <span class="comment">% .Coord(nbpoints,3): set of coordinates defining the object position;</span>
	97	0027
	98	0028 <span class="comment">%FieldData: data of the field to be projected on the projection object, with optional fields</span>
	99	0029 <span class="comment">% .Txt: error message, transmitted to the projection</span>
	100	0030 <span class="comment">% .CoordType: 'px' or 'phys' type of coordinates of the field, must be the same as for the projection object, transmitted</span>
	101	0031 <span class="comment">% .Mesh: typical distance between data points (used for mouse action or display), transmitted</span>
	102	0032 <span class="comment">% .CoordUnit, .TimeUnit, .dt: transmitted</span>
	103	0033 <span class="comment">% standardised description of fields, nc-formated Matlab structure with fields:</span>
	104	0034 <span class="comment">% .ListGlobalAttribute: cell listing the names of the global attributes</span>
	105	0035 <span class="comment">% .Att_1,Att_2... : values of the global attributes</span>
	106	0036 <span class="comment">% .ListDimName: cell listing the names of the array dimensions</span>
	107	0037 <span class="comment">% .DimValue: array dimension values (Matlab vector with the same length as .ListDimName</span>
	108	0038 <span class="comment">% .ListVarName: cell listing the names of the variables</span>
	109	0039 <span class="comment">% .VarDimIndex: cell containing the set of dimension indices (in list .ListDimName) for each variable of .ListVarName</span>
	110	0040 <span class="comment">% .VarAttribute: cell of structures s containing names and values of variable attributes (s.name=value) for each variable of .ListVarName</span>
	111	0041 <span class="comment">% .Var1, .Var2....: variables (Matlab arrays) with names listed in .ListVarName</span>
	112	0042 <span class="comment">% The variables are grouped in 'fields', made of a set of variables with common dimensions (using the function find_field_indices)</span>
	113	0043 <span class="comment">% The variable attribute 'Role' is used to define the role for plotting:</span>
	114	0044 <span class="comment">% Role = 'scalar': (default) represents a scalar field</span>
	115	0045 <span class="comment">% = 'coord': represents a set of unstructured coordinates, whose</span>
	116	0046 <span class="comment">% space dimension is given by the last array dimension (called 'nb_dim').</span>
	117	0047 <span class="comment">% = 'coord_x', 'coord_y', 'coord_z': represents a separate set of</span>
	118	0048 <span class="comment">% unstructured coordinate x, y or z</span>
	119	0049 <span class="comment">% = 'vector': represents a vector field whose number of components</span>
	120	0050 <span class="comment">% is given by the last dimension (called 'nb_dim')</span>
	121	0051 <span class="comment">% = 'vector_x', 'vector_y', 'vector_z' :represents the x, y or z component of a vector</span>
	122	0052 <span class="comment">% = 'warnflag' : provides a warning flag about the quality of data in a 'Field', default=0, no warning</span>
	123	0053 <span class="comment">% = 'errorflag': provides an error flag marking false data,</span>
	124	0054 <span class="comment">% default=0, no error. Different non zero values can represent different criteria of elimination.</span>
	125	0055 <span class="comment">%</span>
	126	0056 <span class="comment">% Default role of variables (by name)</span>
	127	0057 <span class="comment">% vector field:</span>
	128	0058 <span class="comment">% .X,.Y: position of the velocity vectors, projected on the object</span>
	129	0059 <span class="comment">% .U, .V, .W: velocity components, projected on the object</span>
	130	0060 <span class="comment">% .C, .CName: scalar associated to the vector</span>
	131	0061 <span class="comment">% .F : equivalent to 'warnflag'</span>
	132	0062 <span class="comment">% .FF: equivalent to 'errorflag'</span>
	133	0063 <span class="comment">% scalar field or image:</span>
	134	0064 <span class="comment">% .AName: name of a scalar (to be calculated from velocity fields after projection), transmitted</span>
	135	0065 <span class="comment">% .A: scalar, projected on the object</span>
	136	0066 <span class="comment">% .AX, .AY: positions for the scalar</span>
	137	0067 <span class="comment">% case of a structured grid: A is a dim 2 matrix and .AX=[first last] (length 2 vector) represents the first and last abscissa of the grid</span>
	138	0068 <span class="comment">% case of an unstructured scalar: A is a vector, AX and AY the corresponding coordinates</span>
	139	0069 <span class="comment">%</span>
	140	0070 <span class="comment">%AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA</span>
	141	0071 <span class="comment">% Copyright Joel Sommeria, 2008, LEGI / CNRS-UJF-INPG, sommeria@coriolis-legi.org.</span>
	142	0072 <span class="comment">%AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA</span>
	143	0073 <span class="comment">% This file is part of the toolbox UVMAT.</span>
	144	0074 <span class="comment">%</span>
	145	0075 <span class="comment">% UVMAT is free software; you can redistribute it and/or modify</span>
	146	0076 <span class="comment">% it under the terms of the GNU General Public License as published by</span>
	147	0077 <span class="comment">% the Free Software Foundation; either version 2 of the License, or</span>
	148	0078 <span class="comment">% (at your option) any later version.</span>
	149	0079 <span class="comment">%</span>
	150	0080 <span class="comment">% UVMAT is distributed in the hope that it will be useful,</span>
	151	0081 <span class="comment">% but WITHOUT ANY WARRANTY; without even the implied warranty of</span>
	152	0082 <span class="comment">% MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the</span>
	153	0083 <span class="comment">% GNU General Public License (file UVMAT/COPYING.txt) for more details.</span>
	154	0084 <span class="comment">%AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA</span>
	155	0085
	156	0086 <a name="_sub0" href="#_subfunctions" class="code">function [ProjData,errormsg]=proj_field(FieldData,ObjectData,IndexObj)</a>
	157	0087
	158	0088 <span class="keyword">if</span> isfield(ObjectData,<span class="string">'ProjMode'</span>) && (isequal(ObjectData.ProjMode,<span class="string">'none'</span>)\|\|isequal(ObjectData.ProjMode,<span class="string">'mask_inside'</span>)\|\|isequal(ObjectData.ProjMode,<span class="string">'mask_outside'</span>))
	159	0089 ProjData=[];
	160	0090 <span class="keyword">return</span>
	161	0091 <span class="keyword">end</span>
	162	0092 <span class="comment">%introduce default field properties (reading old standards)</span>
	163	0093 <span class="keyword">if</span> ~isfield(ObjectData,<span class="string">'Style'</span>)\|\|~isfield(ObjectData,<span class="string">'Coord'</span>)\|\|~isfield(ObjectData,<span class="string">'ProjMode'</span>)
	164	0094 ProjData=FieldData;
	165	0095 <span class="keyword">return</span>
	166	0096 <span class="keyword">end</span>
	167	0097
	168	0098 <span class="comment">% OBSOLETE</span>
	169	0099 <span class="keyword">if</span> isfield(ObjectData,<span class="string">'XMax'</span>) && ~isempty(ObjectData.XMax)
	170	0100 ObjectData.RangeX(1)=ObjectData.XMax;
	171	0101 <span class="keyword">end</span>
	172	0102 <span class="keyword">if</span> isfield(ObjectData,<span class="string">'XMin'</span>) && ~isempty(ObjectData.XMin)
	173	0103 ObjectData.RangeX(2)=ObjectData.XMin;
	174	0104 <span class="keyword">end</span>
	175	0105 <span class="keyword">if</span> isfield(ObjectData,<span class="string">'YMax'</span>) && ~isempty(ObjectData.YMax)
	176	0106 ObjectData.RangeY(1)=ObjectData.YMax;
	177	0107 <span class="keyword">end</span>
	178	0108 <span class="keyword">if</span> isfield(ObjectData,<span class="string">'YMin'</span>) && ~isempty(ObjectData.YMin)
	179	0109 ObjectData.RangeY(2)=ObjectData.YMin;
	180	0110 <span class="keyword">end</span>
	181	0111 <span class="keyword">if</span> isfield(ObjectData,<span class="string">'ZMax'</span>) && ~isempty(ObjectData.ZMax)
	182	0112 ObjectData.RangeZ(1)=ObjectData.ZMax;
	183	0113 <span class="keyword">end</span>
	184	0114 <span class="keyword">if</span> isfield(ObjectData,<span class="string">'ZMin'</span>) && ~isempty(ObjectData.ZMin)
	185	0115 ObjectData.RangeZ(2)=ObjectData.ZMin;
	186	0116 <span class="keyword">end</span>
	187	0117 <span class="comment">%%%%%%%%%%</span>
	188	0118
	189	0119 <span class="comment">% FieldData=document_field(FieldData);%transform FieldData to the standard format</span>
	190	0120 <span class="comment">% if ~isfield(FieldData,'VarAttribute')</span>
	191	0121 <span class="comment">% FieldData.VarAttribute={};</span>
	192	0122 <span class="comment">% end</span>
	193	0123
	194	0124 <span class="keyword">if</span> isequal(ObjectData.Style,<span class="string">'points'</span>)
	195	0125 [ProjData,errormsg]=<a href="#_sub1" class="code" title="subfunction [ProjData,errormsg]=proj_points(FieldData,ObjectData)">proj_points</a>(FieldData,ObjectData);
	196	0126 <span class="keyword">elseif</span> isequal(ObjectData.Style,<span class="string">'line'</span>)\|isequal(ObjectData.Style,<span class="string">'polyline'</span>)
	197	0127 [ProjData,errormsg] = <a href="#_sub3" class="code" title="subfunction [ProjData,errormsg] = proj_line(FieldData, ObjectData)">proj_line</a>(FieldData,ObjectData);
	198	0128 <span class="keyword">elseif</span> isequal(ObjectData.Style,<span class="string">'polygon'</span>)\|isequal(ObjectData.Style,<span class="string">'rectangle'</span>)\|isequal(ObjectData.Style,<span class="string">'ellipse'</span>)
	199	0129 <span class="keyword">if</span> isequal(ObjectData.ProjMode,<span class="string">'inside'</span>)\|isequal(ObjectData.ProjMode,<span class="string">'outside'</span>)
	200	0130 [ProjData,errormsg] = <a href="#_sub2" class="code" title="subfunction [ProjData,errormsg]=proj_patch(FieldData,ObjectData)">proj_patch</a>(FieldData,ObjectData);
	201	0131 <span class="keyword">else</span>
	202	0132 [ProjData,errormsg] = <a href="#_sub3" class="code" title="subfunction [ProjData,errormsg] = proj_line(FieldData, ObjectData)">proj_line</a>(FieldData,ObjectData);
	203	0133 <span class="keyword">end</span>
	204	0134 <span class="comment">%A FAIRE : GERER MASK</span>
	205	0135 <span class="keyword">elseif</span> isequal(ObjectData.Style,<span class="string">'plane'</span>)
	206	0136 <span class="comment">% if isfield(FieldData,'NbDim') & isequal(FieldData.NbDim,3)</span>
	207	0137 <span class="comment">% ProjData= proj_plane3D(FieldData,ObjectData);%</span>
	208	0138 <span class="comment">% else</span>
	209	0139 [ProjData,errormsg] = <a href="#_sub4" class="code" title="subfunction [ProjData,errormsg] = proj_plane(FieldData, ObjectData)">proj_plane</a>(FieldData,ObjectData);
	210	0140 <span class="comment">% end</span>
	211	0141 <span class="keyword">end</span>
	212	0142 <span class="keyword">if</span> exist(<span class="string">'IndexObj'</span>,<span class="string">'var'</span>)
	213	0143 ProjData.IndexObj=IndexObj;<span class="comment">%transfer object index</span>
	214	0144 <span class="keyword">end</span>
	215	0145
	216	0146 <span class="comment">%-----------------------------------------------------------------</span>
	217	0147 <span class="comment">%project on a set of points</span>
	218	0148 <a name="_sub1" href="#_subfunctions" class="code">function [ProjData,errormsg]=proj_points(FieldData,ObjectData)</a><span class="comment">%%</span>
	219	0149 <span class="comment">%-------------------------------------------------------------------</span>
	220	0150
	221	0151 siz=size(ObjectData.Coord);
	222	0152 width=0;
	223	0153 <span class="keyword">if</span> isfield(ObjectData,<span class="string">'Range'</span>)
	224	0154 width=ObjectData.Range(1,2);
	225	0155 <span class="keyword">end</span>
	226	0156 <span class="keyword">if</span> isfield(ObjectData,<span class="string">'RangeX'</span>)&~isempty(ObjectData.RangeX)
	227	0157 width=max(ObjectData.RangeX);
	228	0158 <span class="keyword">end</span>
	229	0159 <span class="keyword">if</span> isfield(ObjectData,<span class="string">'RangeY'</span>)&~isempty(ObjectData.RangeY)
	230	0160 width=max(width,max(ObjectData.RangeY));
	231	0161 <span class="keyword">end</span>
	232	0162 <span class="keyword">if</span> isfield(ObjectData,<span class="string">'RangeZ'</span>)&~isempty(ObjectData.RangeZ)
	233	0163 width=max(width,max(ObjectData.RangeZ));
	234	0164 <span class="keyword">end</span>
	235	0165 <span class="keyword">if</span> isequal(ObjectData.ProjMode,<span class="string">'projection'</span>)
	236	0166 <span class="keyword">if</span> width==0
	237	0167 errormsg=<span class="string">'projection range around points needed'</span>;
	238	0168 <span class="keyword">return</span>
	239	0169 <span class="keyword">end</span>
	240	0170 <span class="keyword">elseif</span> ~isequal(ObjectData.ProjMode,<span class="string">'interp'</span>)
	241	0171 errormsg=([<span class="string">'ProjMode option '</span> ObjectData.ProjMode <span class="string">' not available in proj_field'</span>]);
	242	0172 <span class="keyword">return</span>
	243	0173 <span class="keyword">end</span>
	244	0174 ProjData=<a href="#_sub5" class="code" title="subfunction [ProjData,errormsg]=proj_heading(FieldData,ObjectData)">proj_heading</a>(FieldData,ObjectData);
	245	0175 ProjData.NbDim=0;
	246	0176 ProjData.ListDimName= {<span class="string">'nb_points'</span>};
	247	0177 ProjData.DimValue=siz(1); <span class="comment">%nbre of projection points</span>
	248	0178
	249	0179
	250	0180 <span class="comment">% idimvar=0;</span>
	251	0181 [CellVarIndex,NbDim,VarTypeCell,errormsg]=<a href="find_field_indices.html" class="code" title="function [CellVarIndex,NbDim,VarType,errormsg]=find_field_indices(Data)">find_field_indices</a>(FieldData);
	252	0182 <span class="keyword">if</span> ~isempty(errormsg)
	253	0183 errormsg=[<span class="string">'error in proj_field/proj_points:'</span> errormsg];
	254	0184 <span class="keyword">return</span>
	255	0185 <span class="keyword">end</span>
	256	0186 <span class="comment">%LOOP ON GROUPS OF VARIABLES SHARING THE SAME DIMENSIONS</span>
	257	0187 <span class="comment">% CellVarIndex=cells of variable index arrays</span>
	258	0188 <span class="comment">% ivar_new=0; % index of the current variable in the projected field</span>
	259	0189 <span class="comment">% icoord=0;</span>
	260	0190 <span class="keyword">for</span> icell=1:length(CellVarIndex)
	261	0191 <span class="keyword">if</span> NbDim(icell)==1
	262	0192 <span class="keyword">continue</span>
	263	0193 <span class="keyword">end</span>
	264	0194 VarIndex=CellVarIndex{icell};<span class="comment">% indices of the selected variables in the list FieldData.ListVarName</span>
	265	0195 VarType=VarTypeCell{icell};
	266	0196 ivar_X=VarType.coord_x;
	267	0197 ivar_Y=VarType.coord_y;
	268	0198 ivar_Z=VarType.coord_z;
	269	0199 <span class="comment">% ivar_U=VarType.vector_x;</span>
	270	0200 <span class="comment">% ivar_V=VarType.vector_y;</span>
	271	0201 <span class="comment">% ivar_W=VarType.vector_z;</span>
	272	0202 <span class="comment">% ivar_C=VarType.scalar ;</span>
	273	0203 ivar_Anc=VarType.ancillary;
	274	0204 <span class="comment">% test_anc=zeros(size(VarIndex));</span>
	275	0205 test_anc(ivar_Anc)=ones(size(ivar_Anc));
	276	0206 ivar_F=VarType.warnflag;
	277	0207 ivar_FF=VarType.errorflag;
	278	0208 VarIndex([ivar_X ivar_Y ivar_Z ivar_Anc ivar_F ivar_FF])=[];<span class="comment">% not projected variables removed frlom list</span>
	279	0209 <span class="keyword">if</span> isempty(ivar_X)
	280	0210 test_grid=1;<span class="comment">%test for input data on regular grid (e.g. image)coordinates</span>
	281	0211
	282	0212 <span class="keyword">else</span>
	283	0213 <span class="keyword">if</span> length(ivar_X)>1 \| length(ivar_Y)>1 \| length(ivar_Z)>1
	284	0214 <a href="warndlg_uvmat.html" class="code" title="function hwarn=warndlg_uvmat(warntext,title)">warndlg_uvmat</a>(<span class="string">'multiple coordinate input in proj_field.m'</span>,<span class="string">'ERROR'</span>)
	285	0215 <span class="keyword">return</span>
	286	0216 <span class="keyword">end</span>
	287	0217 <span class="keyword">if</span> length(ivar_Y)~=1
	288	0218 <a href="warndlg_uvmat.html" class="code" title="function hwarn=warndlg_uvmat(warntext,title)">warndlg_uvmat</a>(<span class="string">'y coordinate not defined in proj_field.m'</span>,<span class="string">'ERROR'</span>)
	289	0219 <span class="keyword">return</span>
	290	0220 <span class="keyword">end</span>
	291	0221 test_grid=0;
	292	0222 <span class="keyword">end</span>
	293	0223 ProjData.ListVarName={<span class="string">'Y'</span>,<span class="string">'X'</span>,<span class="string">'NbVal'</span>};
	294	0224 ProjData.VarDimIndex={[1],[1],[1]};
	295	0225 ProjData.VarAttribute{1}.Role=<span class="string">'ancillary'</span>;
	296	0226 ProjData.VarAttribute{2}.Role=<span class="string">'ancillary'</span>;
	297	0227 ProjData.VarAttribute{3}.Role=<span class="string">'ancillary'</span>;
	298	0228 <span class="keyword">for</span> ivar=VarIndex
	299	0229 VarName=FieldData.ListVarName{ivar};
	300	0230 ProjData.ListVarName=[ProjData.ListVarName {VarName}];
	301	0231 ProjData.VarDimIndex=[ProjData.VarDimIndex {[1]}];
	302	0232 <span class="keyword">end</span>
	303	0233 <span class="keyword">if</span> ~test_grid
	304	0234 eval([<span class="string">'coord_x=FieldData.'</span> FieldData.ListVarName{ivar_X} <span class="string">';'</span>])
	305	0235 eval([<span class="string">'coord_y=FieldData.'</span> FieldData.ListVarName{ivar_Y} <span class="string">';'</span>])
	306	0236 test3D=0;<span class="comment">% TEST 3D CASE : NOT COMPLETED , 3D CASE : NOT COMPLETED</span>
	307	0237 <span class="keyword">if</span> length(ivar_Z)==1
	308	0238 eval([<span class="string">'coord_z=FieldData.'</span> FieldData.ListVarName{ivar_Z} <span class="string">';'</span>])
	309	0239 test3D=1;
	310	0240 <span class="keyword">end</span>
	311	0241 <span class="comment">% if length(ivar_U)>1 \| length(ivar_V)>1 \| length(ivar_W)>1</span>
	312	0242 <span class="comment">% warndlg_uvmat('multiple vector input in proj_field.m','ERROR')</span>
	313	0243 <span class="comment">% return</span>
	314	0244 <span class="comment">% end</span>
	315	0245 <span class="keyword">if</span> length(ivar_F)>1 \| length(ivar_FF)>1
	316	0246 <a href="warndlg_uvmat.html" class="code" title="function hwarn=warndlg_uvmat(warntext,title)">warndlg_uvmat</a>(<span class="string">'multiple flag input in proj_field.m'</span>,<span class="string">'ERROR'</span>)
	317	0247 <span class="keyword">return</span>
	318	0248 <span class="keyword">end</span>
	319	0249
	320	0250 <span class="keyword">for</span> ipoint=1:siz(1)
	321	0251 Xpoint=ObjectData.Coord(ipoint,:);
	322	0252 distX=coord_x-Xpoint(1);
	323	0253 distY=coord_y-Xpoint(2);
	324	0254 dist=distX.distX+distY.distY;
	325	0255 indsel=find(dist<width*width);
	326	0256 ProjData.X(ipoint,1)=Xpoint(1);
	327	0257 ProjData.Y(ipoint,1)=Xpoint(2);
	328	0258 <span class="keyword">if</span> isequal(length(ivar_FF),1)
	329	0259 FFName=FieldData.ListVarName{ivar_FF};
	330	0260 eval([<span class="string">'FF=FieldData.'</span> FFName <span class="string">'(indsel);'</span>])
	331	0261 ind_indsel=find(~FF);
	332	0262 indsel=indsel(ind_indsel);
	333	0263 <span class="keyword">end</span>
	334	0264 ProjData.NbVal(ipoint,1)=length(indsel);
	335	0265 <span class="keyword">for</span> ivar=VarIndex
	336	0266 VarName=FieldData.ListVarName{ivar};
	337	0267 <span class="keyword">if</span> isempty(indsel)
	338	0268 eval([<span class="string">'ProjData.'</span> VarName <span class="string">'(ipoint,1)=NaN;'</span>])
	339	0269 <span class="keyword">else</span>
	340	0270 eval([<span class="string">'Var=FieldData.'</span> VarName <span class="string">'(indsel);'</span>])
	341	0271 eval([<span class="string">'ProjData.'</span> VarName <span class="string">'(ipoint,1)=mean(Var);'</span>])
	342	0272 <span class="keyword">if</span> isequal(ObjectData.ProjMode,<span class="string">'interp'</span>)
	343	0273 eval([<span class="string">'ProjData.'</span> VarName <span class="string">'(ipoint,1)=griddata_uvmat(coord_x(indsel),coord_y(indsel),Var,Xpoint(1),Xpoint(2)))'</span>;])
	344	0274 <span class="keyword">end</span>
	345	0275 <span class="keyword">end</span>
	346	0276 <span class="keyword">end</span>
	347	0277 <span class="keyword">end</span>
	348	0278 <span class="keyword">else</span>
	349	0279 DimIndices=FieldData.VarDimIndex{VarIndex(1)};<span class="comment">%indices of the dimensions of the first variable (common to all variables in the cell)</span>
	350	0280 <span class="comment">%input field on a structured grid A REVOIR A REVOIR A REVOIR</span>
	351	0281 <span class="keyword">if</span> ~isempty(find(DimIndices))
	352	0282 DimValue=FieldData.DimValue(DimIndices);<span class="comment">%set of dimension values</span>
	353	0283 ListDimName=FieldData.ListDimName(DimIndices);
	354	0284 <span class="comment">% nbcolor=1; %default</span>
	355	0285 <span class="keyword">for</span> idim=1:length(ListDimName)
	356	0286 DimName=ListDimName{idim};
	357	0287 <span class="keyword">if</span> isequal(DimName,<span class="string">'rgb'</span>)\|isequal(DimName,<span class="string">'nb_coord'</span>)\|isequal(DimName,<span class="string">'nb_coord_i'</span>)
	358	0288 nbcolor=DimValue(idim);
	359	0289 DimIndices(idim)=[];
	360	0290 DimValue(idim)=[];
	361	0291 <span class="keyword">end</span>
	362	0292 <span class="keyword">if</span> isequal(DimName,<span class="string">'nb_coord_j'</span>)<span class="comment">% NOTE: CASE OF TENSOR NOT TREATED</span>
	363	0293 DimIndices(idim)=[];
	364	0294 DimValue(idim)=[];
	365	0295 <span class="keyword">end</span>
	366	0296 <span class="keyword">end</span>
	367	0297 ind_1=find(DimValue==1);
	368	0298 DimIndices(ind_1)=[]; <span class="comment">%suppress singleton dimensions</span>
	369	0299 indxy=find(DimVarIndex(DimIndices));<span class="comment">%select dimension variables (DimIndices non zero)</span>
	370	0300 nb_dim=length(DimIndices);<span class="comment">%number of space dimensions</span>
	371	0301 Coord_z=[];
	372	0302 Coord_y=[];
	373	0303 Coord_x=[];
	374	0304 <span class="keyword">for</span> idim=1:nb_dim <span class="comment">%loop on space dimensions</span>
	375	0305 test_interp(idim)=0;<span class="comment">%test for coordiate interpolation (non regular grid), =0 by default</span>
	376	0306 test_coord(idim)=0;<span class="comment">%test for defined coordinates, =0 by default</span>
	377	0307 ivar=DimVarIndex(DimIndices(idim));<span class="comment">% index of the variable corresponding to the current dimension</span>
	378	0308 <span class="keyword">if</span> ~isequal(ivar,0)<span class="comment">% a variable corresponds to the current dimension</span>
	379	0309 eval([<span class="string">'Coord{idim}=FieldData.'</span> FieldData.ListVarName{ivar} <span class="string">';'</span>]) ;<span class="comment">% position for the first index</span>
	380	0310 DCoord=diff(Coord{idim});
	381	0311 DCoord_min(idim)=min(DCoord);
	382	0312 DCoord_max=max(DCoord);
	383	0313 test_direct(idim)=DCoord_max>0;<span class="comment">% =1 for increasing values, 0 otherwise</span>
	384	0314 test_direct_min=DCoord_min(idim)>0;<span class="comment">% =1 for increasing values, 0 otherwise</span>
	385	0315 <span class="keyword">if</span> ~isequal(test_direct(idim),test_direct_min)
	386	0316 <a href="warndlg_uvmat.html" class="code" title="function hwarn=warndlg_uvmat(warntext,title)">warndlg_uvmat</a>([<span class="string">'non monotonic dimension variable # '</span> num2str(idim) <span class="string">' in proj_field.m'</span>],<span class="string">'ERROR'</span>)
	387	0317 <span class="keyword">return</span>
	388	0318 <span class="keyword">end</span>
	389	0319 test_interp(idim)=(DCoord_max-DCoord_min(idim))> 0.0001*abs(DCoord_max);<span class="comment">% test grid regularity</span>
	390	0320 test_coord(idim)=1;
	391	0321
	392	0322 <span class="keyword">else</span> <span class="comment">% no variable associated with the first dimension, look fo variable attributes Coord_1, _2 or _3</span>
	393	0323 Coord_i_str=[<span class="string">'Coord_'</span> num2str(idim)];
	394	0324 DCoord_min(idim)=1;<span class="comment">%default</span>
	395	0325 Coord{idim}=[0.5 DimValue(idim)];
	396	0326 test_direct(idim)=1;
	397	0327 <span class="comment">% for ivar=VarIndex</span>
	398	0328 <span class="comment">% if length(FieldData.VarAttribute)>=ivar && isfield(FieldData.VarAttribute{ivar},Coord_i_str)% if there is a variable attribute named Coord_1, _2 or _3</span>
	399	0329 <span class="comment">% eval(['Coord_i=FieldData.VarAttribute{ivar}.' Coord_i_str ';']);%'range x</span>
	400	0330 <span class="comment">% if isnumeric(Coord_i)</span>
	401	0331 <span class="comment">% if length(Coord_i)>=2</span>
	402	0332 <span class="comment">% Coord{idim}=[Coord_i(1) Coord_i(end)];</span>
	403	0333 <span class="comment">% %test_direct(idim)=(Coord{idim}(2)>Coord{idim}(1));</span>
	404	0334 <span class="comment">% else</span>
	405	0335 <span class="comment">% warndlg_uvmat(['two values needed for ' Coord_i_str 'in proj_field.m'],'ERROR')</span>
	406	0336 <span class="comment">% return</span>
	407	0337 <span class="comment">% end</span>
	408	0338 <span class="comment">% else</span>
	409	0339 <span class="comment">% warndlg_uvmat(['non numerical coordinate attributes' Coord_i_str 'in proj_field.m'],'ERROR')</span>
	410	0340 <span class="comment">% return</span>
	411	0341 <span class="comment">% end</span>
	412	0342 <span class="comment">% %test_coord(idim)=1;</span>
	413	0343 <span class="comment">% DCoord_min(idim)=(Coord{idim}(end)-Coord{idim}(1))/(DimValue(idim)-1);</span>
	414	0344 <span class="comment">% end</span>
	415	0345 <span class="comment">% end</span>
	416	0346 <span class="keyword">end</span>
	417	0347 <span class="keyword">end</span>
	418	0348 DX=DCoord_min(2);
	419	0349 DY=DCoord_min(1);
	420	0350 <span class="keyword">for</span> ipoint=1:siz(1)
	421	0351 xwidth=width/(abs(DX));
	422	0352 ywidth=width/(abs(DY));
	423	0353 i_min=round((ObjectData.Coord(ipoint,1)-Coord{2}(1))/DX+0.5-xwidth); <span class="comment">%minimum index of the selected region</span>
	424	0354 i_min=max(1,i_min);<span class="comment">%restrict to field limit</span>
	425	0355 i_plus=round((ObjectData.Coord(ipoint,1)-Coord{2}(1))/DX+0.5+xwidth);
	426	0356 i_plus=min(DimValue(2),i_plus); <span class="comment">%restrict to field limit</span>
	427	0357 j_min=round((ObjectData.Coord(ipoint,2)-Coord{1}(1))/DY-ywidth+0.5);
	428	0358 j_min=max(1,j_min);
	429	0359 j_plus=round((ObjectData.Coord(ipoint,2)-Coord{1}(1))/DY+ywidth+0.5);
	430	0360 j_plus=min(DimValue(1),j_plus);
	431	0361 ProjData.X(ipoint,1)=ObjectData.Coord(ipoint,1);
	432	0362 ProjData.Y(ipoint,1)=ObjectData.Coord(ipoint,2);
	433	0363 i_int=[i_min:i_plus];
	434	0364 j_int=[j_min:j_plus];
	435	0365 ProjData.NbVal(ipoint,1)=length(j_int)*length(i_int);
	436	0366 <span class="keyword">if</span> isempty(i_int) \| isempty(j_int)
	437	0367 <span class="keyword">for</span> ivar=VarIndex
	438	0368 eval([<span class="string">'ProjData.'</span> FieldData.ListVarName{ivar} <span class="string">'(ipoint,:)=NaN;'</span>]);
	439	0369 <span class="keyword">end</span>
	440	0370 errormsg=[<span class="string">'no data points in the selected projection range '</span> num2str(width) ];
	441	0371 <span class="keyword">else</span>
	442	0372 <span class="comment">%TODO: introduce circle in the selected subregion</span>
	443	0373 <span class="comment">%[I,J]=meshgrid([1:j_int],[1:i_int]);</span>
	444	0374 <span class="keyword">for</span> ivar=VarIndex
	445	0375 eval([<span class="string">'Avalue=FieldData.'</span> FieldData.ListVarName{ivar} <span class="string">'(j_int,i_int,:);'</span>]);
	446	0376 eval([<span class="string">'ProjData.'</span> FieldData.ListVarName{ivar} <span class="string">'(ipoint,:)=mean(mean(Avalue));'</span>]);
	447	0377 <span class="keyword">end</span>
	448	0378 <span class="keyword">end</span>
	449	0379 <span class="keyword">end</span>
	450	0380 <span class="keyword">end</span>
	451	0381 <span class="keyword">end</span>
	452	0382 <span class="keyword">end</span>
	453	0383 <span class="comment">% if isfield(FieldData,'U')& isfield(FieldData,'V')& ~isempty(FieldData.U)& ~isempty(FieldData.V)</span>
	454	0384 <span class="comment">% ProjData.ListVarName={'Y','X','U','V','NbVec'};</span>
	455	0385 <span class="comment">% ProjData.VarDimIndex={[1],[1],[1],[1],[1]};</span>
	456	0386 <span class="comment">% Role={'coord_y','coord_x','vector_x','vector_y','ancillary'};</span>
	457	0387 <span class="comment">% for ivar=1:5</span>
	458	0388 <span class="comment">% ProjData.VarAttribute{ivar}.Role=Role{ivar};</span>
	459	0389 <span class="comment">% end</span>
	460	0390 <span class="comment">% if isfield(FieldData,'FF')</span>
	461	0391 <span class="comment">% indsel=find(FieldData.FF==0);</span>
	462	0392 <span class="comment">% if isempty(indsel)</span>
	463	0393 <span class="comment">% FieldData.U=[];</span>
	464	0394 <span class="comment">% FieldData.V=[];</span>
	465	0395 <span class="comment">% FieldData.X=[];</span>
	466	0396 <span class="comment">% FieldData.Y=[];</span>
	467	0397 <span class="comment">% else</span>
	468	0398 <span class="comment">% FieldData.U=FieldData.U(indsel);</span>
	469	0399 <span class="comment">% FieldData.V=FieldData.V(indsel);</span>
	470	0400 <span class="comment">% FieldData.X=FieldData.X(indsel);</span>
	471	0401 <span class="comment">% FieldData.Y=FieldData.Y(indsel);</span>
	472	0402 <span class="comment">% end</span>
	473	0403 <span class="comment">% end</span>
	474	0404 <span class="comment">% if isequal(ObjectData.ProjMode,'projection')</span>
	475	0405 <span class="comment">% if ~isfield(ObjectData,'Range') & ~isfield(ObjectData,'RangeX')& ~isfield(ObjectData,'RangeY')</span>
	476	0406 <span class="comment">% errormsg='projection range needed';</span>
	477	0407 <span class="comment">% return</span>
	478	0408 <span class="comment">% end</span>
	479	0409 <span class="comment">%</span>
	480	0410 <span class="comment">% for ipoint=1:siz(1)</span>
	481	0411 <span class="comment">% Xpoint=ObjectData.Coord(ipoint,:);</span>
	482	0412 <span class="comment">% % Ypoint=ObjectData.YObject(ipoint);</span>
	483	0413 <span class="comment">% distX=FieldData.X-Xpoint(1);</span>
	484	0414 <span class="comment">% distY=FieldData.Y-Xpoint(2);</span>
	485	0415 <span class="comment">% dist=distX.distX+distY.distY;</span>
	486	0416 <span class="comment">% % indsel=find(dist<ObjectData.YMax*ObjectData.YMax);</span>
	487	0417 <span class="comment">% indsel=find(dist<width*width);</span>
	488	0418 <span class="comment">% ProjData.X(ipoint,1)=Xpoint(1);</span>
	489	0419 <span class="comment">% ProjData.Y(ipoint,1)=Xpoint(2);</span>
	490	0420 <span class="comment">% ProjData.U(ipoint,1)=NaN;</span>
	491	0421 <span class="comment">% ProjData.V(ipoint,1)=NaN;</span>
	492	0422 <span class="comment">% ProjData.NbVec(ipoint,1)=length(indsel);</span>
	493	0423 <span class="comment">% if ~isempty(indsel)</span>
	494	0424 <span class="comment">% vector_x=FieldData.U(indsel);</span>
	495	0425 <span class="comment">% vector_y=FieldData.V(indsel);</span>
	496	0426 <span class="comment">% ProjData.U(ipoint,1)=sum(vector_x)/ProjData.NbVec(ipoint,1);</span>
	497	0427 <span class="comment">% ProjData.V(ipoint,1)=sum(vector_y)/ProjData.NbVec(ipoint,1);</span>
	498	0428 <span class="comment">% % ProjData.UMax(ipoint)=max(vector_x);</span>
	499	0429 <span class="comment">% % ProjData.UMin(ipoint)=min(vector_x);</span>
	500	0430 <span class="comment">% % ProjData.VMax(ipoint)=max(vector_y);</span>
	501	0431 <span class="comment">% % ProjData.VMin(ipoint)=min(vector_y);</span>
	502	0432 <span class="comment">% end</span>
	503	0433 <span class="comment">% end</span>
	504	0434 <span class="comment">% elseif isequal(ObjectData.ProjMode,'interp')</span>
	505	0435 <span class="comment">% ProjData.U=(griddata_uvmat(FieldData.X,FieldData.Y,FieldData.U,ProjData.Coord(:,1),ProjData.Coord(:,2)))';</span>
	506	0436 <span class="comment">% ProjData.V=(griddata_uvmat(FieldData.X,FieldData.Y,FieldData.V,ProjData.Coord(:,1),ProjData.Coord(:,2)))';</span>
	507	0437 <span class="comment">% else</span>
	508	0438 <span class="comment">% errormsg=(['ProjMode option ' ObjectData.ProjMode ' not available in proj_field'])</span>
	509	0439 <span class="comment">% return</span>
	510	0440 <span class="comment">% end</span>
	511	0441 <span class="comment">% end</span>
	512	0442 <span class="comment">% if isfield(FieldData,'AName')</span>
	513	0443 <span class="comment">% ProjData.ListVarName={'AY','AX','A'};</span>
	514	0444 <span class="comment">% ProjData.VarDimIndex={[1],[1],[1]};</span>
	515	0445 <span class="comment">% % ProjData.ListVarAttribute={'Role'};</span>
	516	0446 <span class="comment">% % ProjData.Role={'coord_y','coord_x',[]};</span>
	517	0447 <span class="comment">% % ProjData.VarAttribute{1}.Role='coord_y';</span>
	518	0448 <span class="comment">% % ProjData.VarAttribute{2}.Role='coord_x';</span>
	519	0449 <span class="comment">%</span>
	520	0450 <span class="comment">% % if ~isfield(FieldData,'A')</span>
	521	0451 <span class="comment">% % ProjData.A=calc_scal(FieldData.AName,ProjData);</span>
	522	0452 <span class="comment">% % ProjData.AName=FieldData.AName;</span>
	523	0453 <span class="comment">% % else</span>
	524	0454 <span class="comment">% siz=size(ObjectData.Coord);</span>
	525	0455 <span class="comment">% if isequal(size(FieldData.A),size(FieldData.AX))%case of unstructured field at points (AX,AY)</span>
	526	0456 <span class="comment">% for ipoint=1:siz(1)</span>
	527	0457 <span class="comment">% Xpoint=ObjectData.Coord(ipoint,1);</span>
	528	0458 <span class="comment">% Ypoint=ObjectData.Coord(ipoint,2);</span>
	529	0459 <span class="comment">% distX=FieldData.AX-Xpoint;</span>
	530	0460 <span class="comment">% distY=FieldData.AY-Ypoint;</span>
	531	0461 <span class="comment">% dist=distX.distX+distY.distY;</span>
	532	0462 <span class="comment">% % indsel=find(dist<ObjectData.YMax*ObjectData.YMax);</span>
	533	0463 <span class="comment">% indsel=find(dist<width*width);</span>
	534	0464 <span class="comment">% ProjData.NbScal=length(indsel);</span>
	535	0465 <span class="comment">% A=FieldData.A(indsel);</span>
	536	0466 <span class="comment">% ProjData.A(ipoint)=sum(A)/ProjData.NbScal;</span>
	537	0467 <span class="comment">% end</span>
	538	0468 <span class="comment">% elseif length(size(FieldData.A))>=2% A is already a matrix</span>
	539	0469 <span class="comment">% siz_A=size(FieldData.A);</span>
	540	0470 <span class="comment">% npx=siz_A(2);</span>
	541	0471 <span class="comment">% npy=siz_A(1);</span>
	542	0472 <span class="comment">% if length(siz_A)==3</span>
	543	0473 <span class="comment">% npcolor=siz_A(3);</span>
	544	0474 <span class="comment">% ProjData.ListDimName=[ProjData.ListDimName {'rgb'} ];</span>
	545	0475 <span class="comment">% ProjData.VarDimIndex{3}=[1 2];</span>
	546	0476 <span class="comment">% ProjData.DimValue=[siz(1) 3];</span>
	547	0477 <span class="comment">% else</span>
	548	0478 <span class="comment">% npcolor=1;</span>
	549	0479 <span class="comment">% end</span>
	550	0480 <span class="comment">% DX=abs(FieldData.AX(end)-FieldData.AX(1))/(npx-1);</span>
	551	0481 <span class="comment">% DY=abs(FieldData.AY(end)-FieldData.AY(1))/(npy-1);</span>
	552	0482 <span class="comment">% for ipoint=1:siz(1)</span>
	553	0483 <span class="comment">% i_min=round((ObjectData.Coord(ipoint,1)-FieldData.AX(1)-width)/DX+0.5);</span>
	554	0484 <span class="comment">% i_min=max(1,i_min);</span>
	555	0485 <span class="comment">% i_plus=round((ObjectData.Coord(ipoint,1)-FieldData.AX(1)+width)/DX+0.5);</span>
	556	0486 <span class="comment">% i_plus=min(npx,i_plus);</span>
	557	0487 <span class="comment">% j_min=-round((ObjectData.Coord(ipoint,2)-FieldData.AY(1)+width)/DY+0.5);</span>
	558	0488 <span class="comment">% j_min=max(1,j_min);</span>
	559	0489 <span class="comment">% j_plus=-round((ObjectData.Coord(ipoint,2)-FieldData.AY(1)-width)/DY+0.5);</span>
	560	0490 <span class="comment">% j_plus=min(npy,j_plus);</span>
	561	0491 <span class="comment">% ProjData.AX(ipoint,1)=ObjectData.Coord(ipoint,1);</span>
	562	0492 <span class="comment">% ProjData.AY(ipoint,1)=ObjectData.Coord(ipoint,2);</span>
	563	0493 <span class="comment">% i_int=[i_min:i_plus];</span>
	564	0494 <span class="comment">% j_int=[j_min:j_plus];</span>
	565	0495 <span class="comment">% if ~isempty(i_int) & ~isempty(j_int)</span>
	566	0496 <span class="comment">% Avalue=FieldData.A(j_int,i_int,:);</span>
	567	0497 <span class="comment">% ProjData.A(ipoint,:)=mean(mean(Avalue));</span>
	568	0498 <span class="comment">% else</span>
	569	0499 <span class="comment">% ProjData.A(ipoint,:)=NaN;</span>
	570	0500 <span class="comment">% end</span>
	571	0501 <span class="comment">% end</span>
	572	0502 <span class="comment">% end</span>
	573	0503 <span class="comment">% end</span>
	574	0504
	575	0505 <span class="comment">%-----------------------------------------------------------------</span>
	576	0506 <span class="comment">%project in a patch</span>
	577	0507 <a name="_sub2" href="#_subfunctions" class="code">function [ProjData,errormsg]=proj_patch(FieldData,ObjectData)</a><span class="comment">%%</span>
	578	0508 <span class="comment">%-------------------------------------------------------------------</span>
	579	0509 ProjData=<a href="#_sub5" class="code" title="subfunction [ProjData,errormsg]=proj_heading(FieldData,ObjectData)">proj_heading</a>(FieldData,ObjectData);
	580	0510
	581	0511 objectfield=fieldnames(ObjectData);
	582	0512 widthx=0;
	583	0513 widthy=0;
	584	0514 <span class="keyword">if</span> isfield(ObjectData,<span class="string">'RangeX'</span>)&~isempty(ObjectData.RangeX)
	585	0515 widthx=max(ObjectData.RangeX);
	586	0516 <span class="keyword">end</span>
	587	0517 <span class="keyword">if</span> isfield(ObjectData,<span class="string">'RangeY'</span>)&~isempty(ObjectData.RangeY)
	588	0518 widthy=max(ObjectData.RangeY);
	589	0519 <span class="keyword">end</span>
	590	0520
	591	0521 <span class="comment">%A REVOIR, GENERALISER: UTILISER proj_line</span>
	592	0522 ProjData.NbDim=0;
	593	0523 ProjData.ListDimName={};<span class="comment">%name of dimension</span>
	594	0524 ProjData.DimValue=[];<span class="comment">%values of dimension (nbre of vectors)</span>
	595	0525 ProjData.ListVarName={};
	596	0526 <span class="comment">%ProjData.VarDimIndex={};</span>
	597	0527 ProjData.VarDimName={};
	598	0528 <span class="keyword">if</span> isfield (FieldData,<span class="string">'ListVarAttribute'</span>)
	599	0529 ProjData.ListVarAttribute=FieldData.ListVarAttribute;<span class="comment">%list of variable attribute names</span>
	600	0530 <span class="keyword">for</span> iattr=1:length(ProjData.ListVarAttribute)<span class="comment">%initialization of variable attribute values</span>
	601	0531 AttrName=ProjData.ListVarAttribute{iattr};
	602	0532 eval([<span class="string">'ProjData.'</span> AttrName <span class="string">'={};'</span>])
	603	0533 <span class="keyword">end</span>;
	604	0534 <span class="keyword">end</span>
	605	0535
	606	0536 <span class="comment">%group the variables (fields of 'FieldData') in cells of variables with the same dimensions</span>
	607	0537 testfalse=0;
	608	0538 ListIndex={};
	609	0539 <span class="comment">% DimVarIndex=0;%initilise list of indices for dimension variables</span>
	610	0540 idimvar=0;
	611	0541 [CellVarIndex,NbDim,VarTypeCell,errormsg]=<a href="find_field_indices.html" class="code" title="function [CellVarIndex,NbDim,VarType,errormsg]=find_field_indices(Data)">find_field_indices</a>(FieldData);
	612	0542 <span class="keyword">if</span> ~isempty(errormsg)
	613	0543 errormsg=[<span class="string">'error in proj_field/proj_patch:'</span> errormsg];
	614	0544 <span class="keyword">return</span>
	615	0545 <span class="keyword">end</span>
	616	0546
	617	0547 <span class="comment">%LOOP ON GROUPS OF VARIABLES SHARING THE SAME DIMENSIONS</span>
	618	0548 dimcounter=0;
	619	0549 <span class="keyword">for</span> icell=1:length(CellVarIndex)
	620	0550 testX=0;
	621	0551 testY=0;
	622	0552 test_Amat=0;
	623	0553 testfalse=0;
	624	0554 VarIndex=CellVarIndex{icell};<span class="comment">% indices of the selected variables in the list FieldData.ListVarName</span>
	625	0555 VarType=VarTypeCell{icell};
	626	0556 DimIndices=FieldData.VarDimIndex{VarIndex(1)};<span class="comment">%indices of the dimensions of the first variable (common to all variables in the cell)</span>
	627	0557 <span class="keyword">if</span> NbDim(icell)~=2<span class="comment">% proj_patch acts only on fields of space dimension 2</span>
	628	0558 <span class="keyword">continue</span>
	629	0559 <span class="keyword">end</span>
	630	0560 testX=~isempty(VarType.coord_x) && ~isempty(VarType.coord_y);
	631	0561 testfalse=~isempty(VarType.errorflag);
	632	0562 testproj(VarIndex)=zeros(size(VarIndex));<span class="comment">%default</span>
	633	0563 testproj(VarType.scalar)=1;
	634	0564 testproj(VarType.vector_x)=1;
	635	0565 testproj(VarType.vector_y)=1;
	636	0566 testproj(VarType.vector_z)=1;
	637	0567 testproj(VarType.image)=1;
	638	0568 testproj(VarType.color)=1;
	639	0569 VarIndex=VarIndex(find(testproj(VarIndex)));<span class="comment">%select only the projected variables</span>
	640	0570 <span class="keyword">if</span> testX <span class="comment">%case of unstructured coordinates</span>
	641	0571 eval([<span class="string">'nbpoint=numel(FieldData.'</span> FieldData.ListVarName{VarIndex(1)} <span class="string">');'</span>])
	642	0572 <span class="keyword">for</span> ivar=[VarIndex VarType.coord_x VarType.coord_y VarType.errorflag]
	643	0573 VarName=FieldData.ListVarName{ivar};
	644	0574 eval([<span class="string">'FieldData.'</span> VarName <span class="string">'=reshape(FieldData.'</span> VarName <span class="string">',nbpoint,1);'</span>])
	645	0575 <span class="keyword">end</span>
	646	0576 XName=FieldData.ListVarName{VarType.coord_x};
	647	0577 YName=FieldData.ListVarName{VarType.coord_y};
	648	0578 eval([<span class="string">'coord_x=FieldData.'</span> XName <span class="string">';'</span>])
	649	0579 eval([<span class="string">'coord_y=FieldData.'</span> YName <span class="string">';'</span>])
	650	0580 <span class="keyword">end</span>
	651	0581 <span class="keyword">if</span> testfalse
	652	0582 FFName=FieldData.ListVarName{VarType.errorflag};
	653	0583 eval([<span class="string">'errorflag=FieldData.'</span> FFName <span class="string">';'</span>])
	654	0584 <span class="keyword">end</span>
	655	0585 <span class="comment">% image or 2D matrix</span>
	656	0586 <span class="keyword">if</span> numel(VarType.coord)>=2 & VarType.coord(1:2) > 0;
	657	0587 test_Amat=1;<span class="comment">%image or 2D matrix</span>
	658	0588 AYName=FieldData.ListVarName{VarType.coord(1)};
	659	0589 AXName=FieldData.ListVarName{VarType.coord(2)};
	660	0590 eval([<span class="string">'AX=FieldData.'</span> AXName <span class="string">';'</span>])
	661	0591 eval([<span class="string">'AY=FieldData.'</span> AYName <span class="string">';'</span>])
	662	0592 VarName=FieldData.ListVarName{VarIndex(1)};
	663	0593 eval([<span class="string">'DimValue=size(FieldData.'</span> VarName <span class="string">');'</span>])
	664	0594 testcolor=find(numel(DimValue)==3);
	665	0595 <span class="comment">% errormsg='multicomponent field not projected';</span>
	666	0596
	667	0597 <span class="keyword">for</span> idim=1:length(DimIndices)
	668	0598 Coord_i_str=[<span class="string">'Coord_'</span> num2str(idim)];
	669	0599 DCoord_min(idim)=1;<span class="comment">%default</span>
	670	0600 Coord{idim}=[0.5 DimValue(idim)];
	671	0601 test_direct(idim)=1;
	672	0602 <span class="keyword">if</span> isfield(FieldData,<span class="string">'VarAttribute'</span>)
	673	0603 <span class="keyword">for</span> ivar=VarIndex
	674	0604 <span class="keyword">if</span> length(FieldData.VarAttribute)>=ivar & isfield(FieldData.VarAttribute{ivar},Coord_i_str)<span class="comment">% if there is a variable attribute named Coord_1, _2 or _3</span>
	675	0605 eval([<span class="string">'Coord_i=FieldData.VarAttribute{ivar}.'</span> Coord_i_str <span class="string">';'</span>]);<span class="comment">%'range x</span>
	676	0606 <span class="keyword">if</span> isnumeric(Coord_i)
	677	0607 <span class="keyword">if</span> length(Coord_i)>=2
	678	0608 Coord{idim}=[Coord_i(1) Coord_i(end)];
	679	0609 <span class="comment">%test_direct(idim)=(Coord{idim}(2)>Coord{idim}(1));</span>
	680	0610 <span class="keyword">else</span>
	681	0611 <a href="warndlg_uvmat.html" class="code" title="function hwarn=warndlg_uvmat(warntext,title)">warndlg_uvmat</a>([<span class="string">'two values needed for '</span> Coord_i_str <span class="string">'in proj_field.m'</span>],<span class="string">'ERROR'</span>)
	682	0612 <span class="keyword">return</span>
	683	0613 <span class="keyword">end</span>
	684	0614 <span class="keyword">else</span>
	685	0615 <a href="warndlg_uvmat.html" class="code" title="function hwarn=warndlg_uvmat(warntext,title)">warndlg_uvmat</a>([<span class="string">'non numerical coordinate attributes'</span> Coord_i_str <span class="string">'in proj_field.m'</span>],<span class="string">'ERROR'</span>)
	686	0616 <span class="keyword">return</span>
	687	0617 <span class="keyword">end</span>
	688	0618 <span class="comment">%test_coord(idim)=1;</span>
	689	0619 DCoord_min(idim)=(Coord{idim}(end)-Coord{idim}(1))/(DimValue(idim)-1);
	690	0620 <span class="keyword">end</span>
	691	0621 <span class="keyword">end</span>
	692	0622 <span class="keyword">end</span>
	693	0623 <span class="keyword">end</span>
	694	0624 AX=linspace(Coord{2}(1),Coord{2}(2),DimValue(2));
	695	0625 AY=linspace(Coord{1}(1),Coord{1}(2),DimValue(1)); <span class="comment">%TODO : 3D case</span>
	696	0626 <span class="keyword">if</span> length(DimValue)==3
	697	0627 testcolor=1;
	698	0628 npxy(3)=3;
	699	0629 <span class="keyword">else</span>
	700	0630 testcolor=0;
	701	0631 <span class="keyword">end</span>
	702	0632 [Xi,Yi]=meshgrid(AX,AY);
	703	0633 npxy(1)=length(AY);
	704	0634 npxy(2)=length(AX);
	705	0635 Xi=reshape(Xi,npxy(1)*npxy(2),1);
	706	0636 Yi=reshape(Yi,npxy(1)*npxy(2),1);
	707	0637 <span class="keyword">for</span> ivar=1:length(VarIndex)
	708	0638 VarName=FieldData.ListVarName{VarIndex(ivar)};
	709	0639 eval([<span class="string">'FieldData.'</span> VarName <span class="string">'=reshape(FieldData.'</span> VarName <span class="string">',npxy(1)*npxy(2),npxy(3));'</span>]); <span class="comment">% keep only non false vectors</span>
	710	0640 <span class="keyword">end</span>
	711	0641 <span class="keyword">end</span>
	712	0642 <span class="comment">%select the indices in the range of action</span>
	713	0643 testin=[];<span class="comment">%default</span>
	714	0644 <span class="keyword">if</span> isequal(ObjectData.Style,<span class="string">'rectangle'</span>)
	715	0645 <span class="comment">% if ~isfield(ObjectData,'RangeX')\|~isfield(ObjectData,'RangeY')</span>
	716	0646 <span class="comment">% errormsg='rectangle half sides RangeX and RangeY needed'</span>
	717	0647 <span class="comment">% return</span>
	718	0648 <span class="comment">% end</span>
	719	0649 <span class="keyword">if</span> testX
	720	0650 distX=abs(coord_x-ObjectData.Coord(1,1));
	721	0651 distY=abs(coord_y-ObjectData.Coord(1,2));
	722	0652 testin=distX<widthx & distY<widthy;
	723	0653 <span class="keyword">elseif</span> test_Amat
	724	0654 distX=abs(Xi-ObjectData.Coord(1,1));
	725	0655 distY=abs(Yi-ObjectData.Coord(1,2));
	726	0656 testin=distX<widthx & distY<widthy;
	727	0657 <span class="keyword">end</span>
	728	0658 <span class="keyword">elseif</span> isequal(ObjectData.Style,<span class="string">'polygon'</span>)
	729	0659 <span class="keyword">if</span> testX
	730	0660 testin=inpolygon(coord_x,coord_y,ObjectData.Coord(:,1),ObjectData.Coord(:,2));
	731	0661 <span class="keyword">elseif</span> test_Amat
	732	0662 testin=inpolygon(Xi,Yi,ObjectData.Coord(:,1),ObjectData.Coord(:,2));
	733	0663 <span class="keyword">else</span><span class="comment">%calculate the scalar</span>
	734	0664 testin=[]; <span class="comment">%A REVOIR</span>
	735	0665 <span class="keyword">end</span>
	736	0666 <span class="keyword">elseif</span> isequal(ObjectData.Style,<span class="string">'ellipse'</span>)
	737	0667 X2Max=widthx*widthx;
	738	0668 Y2Max=(widthy)*(widthy);
	739	0669 <span class="keyword">if</span> testX
	740	0670 distX=(coord_x-ObjectData.Coord(1,1));
	741	0671 distY=(coord_y-ObjectData.Coord(1,2));
	742	0672 testin=(distX.distX/X2Max+distY.distY/Y2Max)<1;
	743	0673 <span class="keyword">elseif</span> test_Amat <span class="comment">%case of usual 2x2 matrix</span>
	744	0674 distX=(Xi-ObjectData.Coord(1,1));
	745	0675 distY=(Yi-ObjectData.Coord(1,2));
	746	0676 testin=(distX.distX/X2Max+distY.distY/Y2Max)<1;
	747	0677 <span class="keyword">end</span>
	748	0678 <span class="keyword">end</span>
	749	0679 <span class="comment">%selected indices</span>
	750	0680 <span class="keyword">if</span> isequal(ObjectData.ProjMode,<span class="string">'outside'</span>)
	751	0681 testin=~testin;
	752	0682 <span class="keyword">end</span>
	753	0683 <span class="keyword">if</span> testfalse
	754	0684 testin=testin & (errorflag==0); <span class="comment">% keep only non false vectors</span>
	755	0685 <span class="keyword">end</span>
	756	0686 indsel=find(testin);
	757	0687 <span class="keyword">if</span> length(DimIndices)==1
	758	0688 ProjDimName=FieldData.ListDimName(DimIndices);<span class="comment">% reproduce the previous dim name</span>
	759	0689 <span class="keyword">elseif</span> dimcounter==0
	760	0690 ProjDimName={<span class="string">'nb_point'</span>};
	761	0691 dimcounter=dimcounter+1;
	762	0692 <span class="keyword">else</span>
	763	0693 ProjDimName={[<span class="string">'nb_point_'</span> num2str(dimcounter)]};
	764	0694 <span class="keyword">end</span>
	765	0695 <span class="comment">% ProjData.ListDimName=[ProjData.ListDimName ProjDimName];</span>
	766	0696 <span class="comment">%ProjData.DimValue=[ProjData.DimValue length(indsel)];</span>
	767	0697 <span class="keyword">for</span> ivar=VarIndex
	768	0698 <span class="keyword">if</span> testproj(ivar)
	769	0699 VarName=FieldData.ListVarName{ivar};
	770	0700 eval([<span class="string">'ProjData.'</span> VarName <span class="string">'=FieldData.'</span> VarName <span class="string">'(indsel,:);'</span>]); <span class="comment">% keep only non false vectors</span>
	771	0701 ProjData.ListVarName=[ProjData.ListVarName {VarName}];
	772	0702 <span class="comment">% ProjData.VarDimIndex=[ProjData.VarDimIndex {length(ProjData.DimValue)}];</span>
	773	0703 ProjData.VarDimName=[ProjData.VarDimName ProjDimName];
	774	0704 <span class="keyword">end</span>
	775	0705 <span class="keyword">end</span>
	776	0706 <span class="keyword">if</span> test_Amat & testcolor
	777	0707 <span class="comment">%ProjData.ListDimName=[ProjData.ListDimName {'rgb'}];</span>
	778	0708 <span class="comment">% ProjData.DimValue=[ProjData.DimValue 3];</span>
	779	0709 <span class="comment">% ProjData.VarDimIndex={[1 2]};</span>
	780	0710 ProjData.VarDimName={{<span class="string">'nb_point'</span>,<span class="string">'rgb'</span>}};
	781	0711 <span class="keyword">end</span>
	782	0712 <span class="keyword">end</span>
	783	0713
	784	0714
	785	0715
	786	0716 <span class="comment">%-----------------------------------------------------------------</span>
	787	0717 <span class="comment">%project on a line</span>
	788	0718 <span class="comment">% AJOUTER flux,circul,error</span>
	789	0719 <a name="_sub3" href="#_subfunctions" class="code">function [ProjData,errormsg] = proj_line(FieldData, ObjectData)</a>
	790	0720 <span class="comment">%-----------------------------------------------------------------</span>
	791	0721 ProjData=<a href="#_sub5" class="code" title="subfunction [ProjData,errormsg]=proj_heading(FieldData,ObjectData)">proj_heading</a>(FieldData,ObjectData);<span class="comment">%transfer global attributes</span>
	792	0722 ProjData.NbDim=1;
	793	0723
	794	0724 <span class="comment">%initialisation of the input parameters</span>
	795	0725 ProjMode=<span class="string">'projection'</span>;<span class="comment">%direct projection on the line by default</span>
	796	0726 <span class="keyword">if</span> isfield(ObjectData,<span class="string">'ProjMode'</span>),ProjMode=ObjectData.ProjMode; <span class="keyword">end</span>;
	797	0727 ProjAngle=90; <span class="comment">%90 degrees projection by default</span>
	798	0728 <span class="keyword">if</span> isfield(FieldData,<span class="string">'ProjAngle'</span>),ProjAngle=ObjectData.ProjAngle; <span class="keyword">end</span>;
	799	0729 width=0;<span class="comment">%default width of the projection band</span>
	800	0730
	801	0731 <span class="keyword">if</span> isfield(ObjectData,<span class="string">'Range'</span>)&size(ObjectData.Range,2)>=2
	802	0732 width=abs(ObjectData.Range(1,2));
	803	0733 <span class="keyword">end</span>
	804	0734 <span class="keyword">if</span> isfield(ObjectData,<span class="string">'RangeY'</span>)
	805	0735 width=max(ObjectData.RangeY);
	806	0736 <span class="keyword">end</span>
	807	0737 error=0;<span class="comment">%default</span>
	808	0738 Xline=[];
	809	0739
	810	0740 flux=0;
	811	0741 circul=0;
	812	0742 liny=ObjectData.Coord(:,2);
	813	0743 siz_line=size(ObjectData.Coord);
	814	0744 <span class="keyword">if</span> siz_line(1)<2
	815	0745 <span class="keyword">return</span><span class="comment">% line needs at least 2 points to be defined</span>
	816	0746 <span class="keyword">end</span>
	817	0747 testfalse=0;
	818	0748 ListIndex={};
	819	0749
	820	0750 <span class="comment">%angles of the polyline and boundaries of action</span>
	821	0751 dlinx=diff(ObjectData.Coord(:,1));
	822	0752 dliny=diff(ObjectData.Coord(:,2));
	823	0753 theta=angle(dlinx+i*dliny);<span class="comment">%angle of each segment</span>
	824	0754 theta(siz_line(1))=theta(siz_line(1)-1);
	825	0755 <span class="comment">% determine a rectangles at +-width from the line (only used for the ProjMode='projection or 'filter')</span>
	826	0756 <span class="keyword">if</span> isequal(ProjMode,<span class="string">'projection'</span>) \|\| isequal(ProjMode,<span class="string">'filter'</span>)
	827	0757 xsup(1)=ObjectData.Coord(1,1)-width*sin(theta(1));
	828	0758 xinf(1)=ObjectData.Coord(1,1)+width*sin(theta(1));
	829	0759 ysup(1)=ObjectData.Coord(1,2)+width*cos(theta(1));
	830	0760 yinf(1)=ObjectData.Coord(1,2)-width*cos(theta(1));
	831	0761 <span class="keyword">for</span> ip=2:siz_line(1)
	832	0762 xsup(ip)=ObjectData.Coord(ip,1)-width*sin((theta(ip)+theta(ip-1))/2)/cos((theta(ip-1)-theta(ip))/2);
	833	0763 xinf(ip)=ObjectData.Coord(ip,1)+width*sin((theta(ip)+theta(ip-1))/2)/cos((theta(ip-1)-theta(ip))/2);
	834	0764 ysup(ip)=ObjectData.Coord(ip,2)+width*cos((theta(ip)+theta(ip-1))/2)/cos((theta(ip-1)-theta(ip))/2);
	835	0765 yinf(ip)=ObjectData.Coord(ip,2)-width*cos((theta(ip)+theta(ip-1))/2)/cos((theta(ip-1)-theta(ip))/2);
	836	0766 <span class="keyword">end</span>
	837	0767 <span class="keyword">end</span>
	838	0768
	839	0769 <span class="comment">%group the variables (fields of 'FieldData') in cells of variables with the same dimensions</span>
	840	0770 [CellVarIndex,NbDim,VarTypeCell,errormsg]=<a href="find_field_indices.html" class="code" title="function [CellVarIndex,NbDim,VarType,errormsg]=find_field_indices(Data)">find_field_indices</a>(FieldData);
	841	0771 <span class="keyword">if</span> ~isempty(errormsg)
	842	0772 errormsg=[<span class="string">'error in proj_field/proj_line:'</span> errormsg];
	843	0773 <span class="keyword">return</span>
	844	0774 <span class="keyword">end</span>
	845	0775
	846	0776 ProjData.ListVarName={};
	847	0777 ProjData.VarDimName={};
	848	0778 <span class="keyword">for</span> icell=1:length(CellVarIndex)
	849	0779 VarIndex=CellVarIndex{icell};<span class="comment">% indices of the selected variables in the list FieldData.ListVarName</span>
	850	0780 VarType=VarTypeCell{icell};
	851	0781 <span class="comment">%DimCell=FieldData.VarDimName{VarIndex(1)}; % list of dimension names fore the cell</span>
	852	0782 <span class="keyword">if</span> NbDim(icell)~=2<span class="comment">% proj_line acts only on fields of space dimension 2, TODO: check 3D case</span>
	853	0783 <span class="keyword">continue</span>
	854	0784 <span class="keyword">end</span>
	855	0785 testX=~isempty(VarType.coord_x) && ~isempty(VarType.coord_y);
	856	0786 testU=~isempty(VarType.vector_x) && ~isempty(VarType.vector_y);
	857	0787 testfalse=~isempty(VarType.errorflag);
	858	0788 testproj(VarIndex)=zeros(size(VarIndex));<span class="comment">%default</span>
	859	0789 testproj(VarType.scalar)=1;
	860	0790 <span class="comment">% testproj(VarType.vector_x)=1;</span>
	861	0791 <span class="comment">% testproj(VarType.vector_y)=1;</span>
	862	0792 testproj(VarType.vector_z)=1;
	863	0793 testproj(VarType.image)=1;
	864	0794 testproj(VarType.color)=1;
	865	0795 VarIndex=VarIndex(find(testproj(VarIndex)));<span class="comment">%select only the projected variables</span>
	866	0796 <span class="keyword">if</span> testU
	867	0797 VarIndex=[VarIndex VarType.vector_x VarType.vector_y];<span class="comment">%put u and v at the end of the list of variables</span>
	868	0798 <span class="keyword">end</span>
	869	0799 <span class="comment">%identify vector components</span>
	870	0800 <span class="keyword">if</span> testU
	871	0801 UName=FieldData.ListVarName{VarType.vector_x};
	872	0802 VName=FieldData.ListVarName{VarType.vector_y};
	873	0803 eval([<span class="string">'vector_x=FieldData.'</span> UName <span class="string">';'</span>])
	874	0804 eval([<span class="string">'vector_y=FieldData.'</span> VName <span class="string">';'</span>])
	875	0805 <span class="keyword">end</span>
	876	0806
	877	0807 <span class="comment">%identify error flag</span>
	878	0808 <span class="keyword">if</span> testfalse
	879	0809 FFName=FieldData.ListVarName{VarType.errorflag};
	880	0810 eval([<span class="string">'errorflag=FieldData.'</span> FFName <span class="string">';'</span>])
	881	0811 <span class="keyword">end</span>
	882	0812
	883	0813 <span class="comment">%identify coordinates</span>
	884	0814 <span class="comment">% ivar_x=VarType.coord_x % .ccord_x and .coord_y exist if NbDim(icell)>=2</span>
	885	0815 <span class="comment">% DimCell_x=FieldData.VarDimName{ivar_x};</span>
	886	0816 <span class="comment">% if ischar(DimCell_x)</span>
	887	0817 <span class="comment">% DimCell_x={DimCell_x};</span>
	888	0818 <span class="comment">% end</span>
	889	0819
	890	0820 <span class="comment">% case of unstructured positions (position given by the variables with role coord_x, coord_y</span>
	891	0821 <span class="keyword">if</span> testX
	892	0822 <span class="keyword">if</span> ~isequal(ProjMode,<span class="string">'interp'</span>)
	893	0823 <span class="keyword">if</span> width==0
	894	0824 errormsg=<span class="string">'range of the projection object is missing'</span>;
	895	0825 <span class="keyword">return</span>
	896	0826 <span class="keyword">else</span>
	897	0827 lambda=2/(width*width); <span class="comment">%smoothing factor used for filter: weight exp(-2) at distance width from the line</span>
	898	0828 <span class="keyword">end</span>
	899	0829 <span class="keyword">end</span>
	900	0830 <span class="keyword">if</span> ~isequal(ProjMode,<span class="string">'projection'</span>)
	901	0831 <span class="keyword">if</span> isfield(ObjectData,<span class="string">'DX'</span>)&~isempty(ObjectData.DX)
	902	0832 DX=abs(ObjectData.DX);<span class="comment">%mesh of interpolation points along the line</span>
	903	0833 <span class="keyword">else</span>
	904	0834 errormsg=<span class="string">'DX missing'</span>;
	905	0835 <span class="keyword">return</span>
	906	0836 <span class="keyword">end</span>
	907	0837 <span class="keyword">end</span>
	908	0838 XName= FieldData.ListVarName{VarType.coord_x};
	909	0839 YName= FieldData.ListVarName{VarType.coord_y};
	910	0840 eval([<span class="string">'coord_x=FieldData.'</span> XName <span class="string">';'</span>])
	911	0841 eval([<span class="string">'coord_y=FieldData.'</span> YName <span class="string">';'</span>])
	912	0842 <span class="keyword">end</span>
	913	0843
	914	0844 <span class="comment">%initiate projection</span>
	915	0845 <span class="keyword">for</span> ivar=1:length(VarIndex)
	916	0846 ProjLine{ivar}=[];
	917	0847 <span class="keyword">end</span>
	918	0848 XLine=[];
	919	0849 linelengthtot=0;
	920	0850
	921	0851 <span class="comment">% circul=0;</span>
	922	0852 <span class="comment">% flux=0;</span>
	923	0853 <span class="comment">%%%%%%% % A FAIRE CALCULER MEAN DES QUANTITES %%%%%%</span>
	924	0854 <span class="comment">%case of unstructured coordinates</span>
	925	0855 <span class="keyword">if</span> testX
	926	0856 <span class="keyword">for</span> ip=1:siz_line(1)-1 <span class="comment">%Loop on the segments of the polyline</span>
	927	0857 linelength=sqrt(dlinx(ip)dlinx(ip)+dliny(ip)dliny(ip));
	928	0858 <span class="comment">%select the vector indices in the range of action</span>
	929	0859 <span class="keyword">if</span> testfalse
	930	0860 flagsel=(errorflag==0); <span class="comment">% keep only non false vectors</span>
	931	0861 <span class="keyword">else</span>
	932	0862 flagsel=ones(size(coord_x));
	933	0863 <span class="keyword">end</span>
	934	0864 <span class="keyword">if</span> isequal(ProjMode,<span class="string">'projection'</span>) \| isequal(ProjMode,<span class="string">'filter'</span>)
	935	0865 flagsel=flagsel & ((coord_y -yinf(ip))(xinf(ip+1)-xinf(ip))>(coord_x-xinf(ip))(yinf(ip+1)-yinf(ip))) <span class="keyword">...</span>
	936	0866 & ((coord_y -ysup(ip))(xsup(ip+1)-xsup(ip))<(coord_x-xsup(ip))(ysup(ip+1)-ysup(ip))) <span class="keyword">...</span>
	937	0867 & ((coord_y -yinf(ip+1))(xsup(ip+1)-xinf(ip+1))>(coord_x-xinf(ip+1))(ysup(ip+1)-yinf(ip+1))) <span class="keyword">...</span>
	938	0868 & ((coord_y -yinf(ip))(xsup(ip)-xinf(ip))<(coord_x-xinf(ip))(ysup(ip)-yinf(ip)));
	939	0869 <span class="keyword">end</span>
	940	0870 indsel=find(flagsel);<span class="comment">%indsel =indices of good vectors</span>
	941	0871 X_sel=coord_x(indsel);
	942	0872 Y_sel=coord_y(indsel);
	943	0873 nbvar=0;
	944	0874 <span class="keyword">for</span> iselect=1:numel(VarIndex)-2*testU
	945	0875 VarName=FieldData.ListVarName{VarIndex(iselect)};
	946	0876 eval([<span class="string">'ProjVar{iselect}=FieldData.'</span> VarName <span class="string">'(indsel);'</span>]);<span class="comment">%scalar value</span>
	947	0877 <span class="keyword">end</span>
	948	0878 <span class="keyword">if</span> testU
	949	0879 ProjVar{numel(VarIndex)-1}=cos(theta(ip))vector_x(indsel)+sin(theta(ip))vector_y(indsel);<span class="comment">% longitudinal component</span>
	950	0880 ProjVar{numel(VarIndex)}=-sin(theta(ip))vector_x(indsel)+cos(theta(ip))vector_y(indsel);<span class="comment">%transverse component</span>
	951	0881 <span class="keyword">end</span>
	952	0882 <span class="keyword">if</span> isequal(ProjMode,<span class="string">'projection'</span>)
	953	0883 sintheta=sin(theta(ip));
	954	0884 costheta=cos(theta(ip));
	955	0885 Xproj=(X_sel-ObjectData.Coord(ip,1))costheta + (Y_sel-ObjectData.Coord(ip,2))sintheta; <span class="comment">%projection on the line</span>
	956	0886 [Xproj,indsort]=sort(Xproj);
	957	0887 <span class="keyword">for</span> ivar=1:numel(ProjVar)
	958	0888 <span class="keyword">if</span> ~isempty(ProjVar{ivar})
	959	0889 ProjVar{ivar}=ProjVar{ivar}(indsort);
	960	0890 <span class="keyword">end</span>
	961	0891 <span class="keyword">end</span>
	962	0892 <span class="keyword">elseif</span> isequal(ProjMode,<span class="string">'interp'</span>) <span class="comment">%linear interpolation:</span>
	963	0893 npoint=floor(linelength/DX)+1;<span class="comment">% nbre of points in the profile (interval DX)</span>
	964	0894 Xproj=[linelength/(2npoint):linelength/npoint:linelength-linelength/(2npoint)];
	965	0895 xreg=cos(theta(ip))*Xproj+ObjectData.Coord(ip,1);
	966	0896 yreg=sin(theta(ip))*Xproj+ObjectData.Coord(ip,2);
	967	0897 <span class="keyword">for</span> ivar=1:numel(ProjVar)
	968	0898 <span class="keyword">if</span> ~isempty(ProjVar{ivar})
	969	0899 ProjVar{ivar}=<a href="griddata_uvmat.html" class="code" title="function ZI = griddata_uvmat(X,Y,Z,XI,YI)">griddata_uvmat</a>(X_sel,Y_sel,ProjVar{ivar},xreg,yreg);
	970	0900 <span class="keyword">end</span>
	971	0901 <span class="keyword">end</span>
	972	0902 <span class="keyword">elseif</span> isequal(ProjMode,<span class="string">'filter'</span>) <span class="comment">%filtering</span>
	973	0903 npoint=floor(linelength/DX)+1;<span class="comment">% nbre of points in the profile (interval DX)</span>
	974	0904 Xproj=[linelength/(2npoint):linelength/npoint:linelength-linelength/(2npoint)];
	975	0905 siz=size(X_sel);
	976	0906 xregij=cos(theta(ip))Xproj'ones(1,siz(2))+ObjectData.Coord(ip,1);
	977	0907 yregij=sin(theta(ip))Xproj'ones(1,siz(2))+ObjectData.Coord(ip,2);
	978	0908 xij=ones(npoint,1)*X_sel;
	979	0909 yij=ones(npoint,1)*Y_sel;
	980	0910 Aij=exp(-lambda((xij-xregij).(xij-xregij)+(yij-yregij).*(yij-yregij)));
	981	0911 norm=ones(1,siz(2))*Aij';
	982	0912 <span class="keyword">for</span> ivar=1:numel(ProjVar)
	983	0913 <span class="keyword">if</span> ~isempty(ProjVar{ivar})
	984	0914 ProjVar{ivar}=ProjVar{ivar}*Aij'./norm;
	985	0915 <span class="keyword">end</span>
	986	0916 <span class="keyword">end</span>
	987	0917 <span class="keyword">end</span>
	988	0918 <span class="comment">%prolongate the total record</span>
	989	0919 <span class="keyword">for</span> ivar=1:numel(ProjVar)
	990	0920 <span class="keyword">if</span> ~isempty(ProjVar{ivar})
	991	0921 ProjLine{ivar}=[ProjLine{ivar}; ProjVar{ivar}];
	992	0922 <span class="keyword">end</span>
	993	0923 <span class="keyword">end</span>
	994	0924 XLine=[XLine ;(Xproj+linelengthtot)];<span class="comment">%along line abscissa</span>
	995	0925 linelengthtot=linelengthtot+linelength;
	996	0926 <span class="comment">% circul=circul+(sum(U_sel))*linelength/npoint;</span>
	997	0927 <span class="comment">% flux=flux+(sum(V_sel))*linelength/npoint;</span>
	998	0928 <span class="keyword">end</span>
	999	0929 ProjData.X=XLine';
	1000	0930 <span class="comment">% eval(['ProjData.' XName '=XLine'';']); % x coordinate</span>
	1001	0931 <span class="comment">%ProjData.ListDimName={XName};</span>
	1002	0932 <span class="comment">%ProjData.DimValue=[length(XLine)];%values of dimension (nbre of vectors)</span>
	1003	0933 cur_index=1;
	1004	0934 ProjData.ListVarName=[ProjData.ListVarName {XName}];
	1005	0935 ProjData.VarDimName=[ProjData.VarDimName {XName}];
	1006	0936 <span class="comment">%ProjData.VarDimIndex{1}=[1];</span>
	1007	0937 ProjData.VarAttribute{1}.long_name=<span class="string">'abscissa along line'</span>;
	1008	0938
	1009	0939 <span class="comment">%ProjData.VarAttribute{1}.Role='coord_x';</span>
	1010	0940 <span class="keyword">for</span> iselect=1:numel(VarIndex)
	1011	0941 VarName=FieldData.ListVarName{VarIndex(iselect)};
	1012	0942 eval([<span class="string">'ProjData.'</span> VarName <span class="string">'=ProjLine{iselect};'</span>])
	1013	0943 ProjData.ListVarName=[ProjData.ListVarName {VarName}];
	1014	0944 ProjData.VarDimName=[ProjData.VarDimName {XName}];
	1015	0945 <span class="comment">%ProjData.VarDimIndex=[ProjData.VarDimIndex {[1]}];</span>
	1016	0946 ProjData.VarAttribute{iselect}=FieldData.VarAttribute{VarIndex(iselect)};
	1017	0947 <span class="keyword">end</span>
	1018	0948
	1019	0949 <span class="comment">%case of structured coordinates</span>
	1020	0950 <span class="keyword">elseif</span> numel(VarType.coord)>=2 & VarType.coord(1:2) > 0;
	1021	0951 <span class="keyword">if</span> ~isequal(ObjectData.Style,<span class="string">'line'</span>)<span class="comment">% exclude polyline</span>
	1022	0952 errormsg=[<span class="string">'no projection available on '</span> ObjectData.Style <span class="string">'for structured coordinates'</span>]; <span class="comment">%</span>
	1023	0953 <span class="keyword">else</span>
	1024	0954 test_Amat=1;<span class="comment">%image or 2D matrix</span>
	1025	0955 test_interp2=0;<span class="comment">%default</span>
	1026	0956 <span class="comment">% if ~isempty(VarType.coord_y)</span>
	1027	0957 AYName=FieldData.ListVarName{VarType.coord(1)};
	1028	0958 AXName=FieldData.ListVarName{VarType.coord(2)};
	1029	0959 eval([<span class="string">'AX=FieldData.'</span> AXName <span class="string">';'</span>]);<span class="comment">% set of x positions</span>
	1030	0960 eval([<span class="string">'AY=FieldData.'</span> AYName <span class="string">';'</span>]);<span class="comment">% set of y positions</span>
	1031	0961 AName=FieldData.ListVarName{VarType.scalar};
	1032	0962 eval([<span class="string">'A=FieldData.'</span> AName <span class="string">';'</span>]);<span class="comment">% scalar</span>
	1033	0963 npxy=size(A);
	1034	0964 npx=npxy(2);
	1035	0965 npy=npxy(1);
	1036	0966 <span class="keyword">if</span> numel(AX)==2
	1037	0967 DX=(AX(2)-AX(1))/(npx-1);
	1038	0968 <span class="keyword">else</span>
	1039	0969 DX_vec=diff(AX);
	1040	0970 DX=max(DX_vec);
	1041	0971 DX_min=min(DX_vec);
	1042	0972 <span class="keyword">if</span> (DX-DX_min)>0.0001*abs(DX)
	1043	0973 test_interp2=1;
	1044	0974 DX=DX_min;
	1045	0975 <span class="keyword">end</span>
	1046	0976 <span class="keyword">end</span>
	1047	0977 <span class="keyword">if</span> numel(AY)==2
	1048	0978 DY=(AY(2)-AY(1))/(npy-1);
	1049	0979 <span class="keyword">else</span>
	1050	0980 DY_vec=diff(AY);
	1051	0981 DY=max(DY_vec);
	1052	0982 DY_min=min(DY_vec);
	1053	0983 <span class="keyword">if</span> (DY-DY_min)>0.0001*abs(DY)
	1054	0984 test_interp2=1;
	1055	0985 DY=DY_min;
	1056	0986 <span class="keyword">end</span>
	1057	0987 <span class="keyword">end</span>
	1058	0988 AXI=linspace(AX(1),AX(end), npx);<span class="comment">%set of x positions for the interpolated input data</span>
	1059	0989 AYI=linspace(AY(1),AY(end), npy);<span class="comment">%set of x positions for the interpolated input data</span>
	1060	0990 <span class="comment">% else % coordinates defined by first and last values (regular grids)</span>
	1061	0991 <span class="comment">% if isempty(Coord_x)</span>
	1062	0992 <span class="comment">% Coord_x=[0.5 npx-0.5];% default coordinate equal to index</span>
	1063	0993 <span class="comment">% end</span>
	1064	0994 <span class="comment">% if isempty(Coord_y)</span>
	1065	0995 <span class="comment">% Coord_y=[npy-0.5 0.5];% default coordinate equal to reversed index (image convention)</span>
	1066	0996 <span class="comment">% end</span>
	1067	0997 <span class="comment">% AXI=linspace(Coord_x(1),Coord_x(2), npx);%set of y positions in the input data</span>
	1068	0998 <span class="comment">% DX=(Coord_x(2)-Coord_x(1))/(npx-1);%grid mesh</span>
	1069	0999 <span class="comment">% AYI=linspace(Coord_y(1),Coord_y(2), npy);%set of y positions in the input data</span>
	1070	1000 <span class="comment">% DY=(Coord_y(2)-Coord_y(1))/(npy-1);%grid mesh</span>
	1071	1001 <span class="comment">% % indcolor([1 2])=[]; %indices of non-coordinate dimensions (e.g; color component)</span>
	1072	1002 <span class="comment">% end</span>
	1073	1003 <span class="keyword">if</span> isfield(ObjectData,<span class="string">'DX'</span>)
	1074	1004 DXY_line=ObjectData.DX;<span class="comment">%mesh on the projection line</span>
	1075	1005 <span class="keyword">else</span>
	1076	1006 DXY_line=sqrt(abs(DX*DY));<span class="comment">% mesh on the projection line</span>
	1077	1007 <span class="keyword">end</span>
	1078	1008 dlinx=ObjectData.Coord(2,1)-ObjectData.Coord(1,1);
	1079	1009 dliny=ObjectData.Coord(2,2)-ObjectData.Coord(1,2);
	1080	1010 linelength=sqrt(dlinxdlinx+dlinydliny);
	1081	1011 theta=angle(dlinx+i*dliny);<span class="comment">%angle of the line</span>
	1082	1012 <span class="keyword">if</span> isfield(FieldData,<span class="string">'RangeX'</span>)
	1083	1013 XMin=min(FieldData.RangeX);<span class="comment">%shift of the origin on the line</span>
	1084	1014 <span class="keyword">else</span>
	1085	1015 XMin=0;
	1086	1016 <span class="keyword">end</span>
	1087	1017 ProjData.AX=linspace(XMin,XMin+linelength,linelength/DXY_line+1);<span class="comment">%abscissa of the new pixels along the line</span>
	1088	1018 y=linspace(-width,width,2*width/DXY_line+1);<span class="comment">%ordintes of the new pixels (coordinate across the line)</span>
	1089	1019 npX=length(ProjData.AX);
	1090	1020 npY=length(y); <span class="comment">%TODO: utiliser proj_grid</span>
	1091	1021 [X,Y]=meshgrid(ProjData.AX,y);<span class="comment">%grid in the line coordinates</span>
	1092	1022 XIMA=ObjectData.Coord(1,1)+(X-XMin)cos(theta)-Ysin(theta);
	1093	1023 YIMA=ObjectData.Coord(1,2)+(X-XMin)sin(theta)+Ycos(theta);
	1094	1024 XIMA=(XIMA-AX(1))/DX+1;<span class="comment">% index of the original image along x</span>
	1095	1025 YIMA=(YIMA-AY(1))/DY+1;<span class="comment">% index of the original image along y</span>
	1096	1026 XIMA=reshape(round(XIMA),1,npX*npY);<span class="comment">%indices reorganized in 'line'</span>
	1097	1027 YIMA=reshape(round(YIMA),1,npX*npY);
	1098	1028 flagin=XIMA>=1 & XIMA<=npx & YIMA >=1 & YIMA<=npy;<span class="comment">%flagin=1 inside the original image</span>
	1099	1029 ind_in=find(flagin);
	1100	1030 ind_out=find(~flagin);
	1101	1031 ICOMB=(XIMA-1)*npy+YIMA;
	1102	1032 ICOMB=ICOMB(flagin);<span class="comment">%index corresponding to XIMA and YIMA in the aligned original image vec_A</span>
	1103	1033 nbcolor=1; <span class="comment">%color images</span>
	1104	1034 <span class="keyword">if</span> numel(npxy)==2
	1105	1035 nbcolor=1;
	1106	1036 <span class="keyword">elseif</span> length(npxy)==3
	1107	1037 nbcolor=npxy(3);
	1108	1038 <span class="keyword">else</span>
	1109	1039 errormsg=<span class="string">'multicomponent field not projected'</span>;
	1110	1040 display(errormsg)
	1111	1041 <span class="keyword">return</span>
	1112	1042 <span class="keyword">end</span>
	1113	1043 nbvar=length(ProjData.ListVarName);<span class="comment">% number of var from previous cells</span>
	1114	1044 <span class="comment">% nbdim=length(ProjData.ListDimName);%number of dim from previous cells</span>
	1115	1045 ProjData.ListVarName=[ProjData.ListVarName {AXName}];
	1116	1046 ProjData.VarDimName=[ProjData.VarDimName {AXName}];
	1117	1047 <span class="keyword">for</span> ivar=VarIndex
	1118	1048 VarName{ivar}=FieldData.ListVarName{ivar};
	1119	1049 <span class="keyword">if</span> test_interp2<span class="comment">% interpolate on new grid</span>
	1120	1050 eval([<span class="string">'FieldData.'</span> VarName{ivar} <span class="string">'=interp2(FieldData.'</span> AXName <span class="string">',FieldData.'</span> AYName <span class="string">',FieldData.'</span> VarName{ivar} <span class="string">',AXI,AYI'');'</span>]) <span class="comment">%TO TEST</span>
	1121	1051 <span class="keyword">end</span>
	1122	1052 eval([<span class="string">'vec_A=reshape(squeeze(FieldData.'</span> VarName{ivar} <span class="string">'),npx*npy,nbcolor);'</span>]) <span class="comment">%put the original image in colum</span>
	1123	1053 <span class="keyword">if</span> nbcolor==1
	1124	1054 vec_B(ind_in)=vec_A(ICOMB);
	1125	1055 vec_B(ind_out)=zeros(size(ind_out));
	1126	1056 A_out=reshape(vec_B,npY,npX);
	1127	1057 eval([<span class="string">'ProjData.'</span> VarName{ivar} <span class="string">'=((sum(A_out,1)/npY))'';'</span>]);
	1128	1058 <span class="keyword">elseif</span> nbcolor==3
	1129	1059 vec_B(ind_in,[1:3])=vec_A(ICOMB,:);
	1130	1060 vec_B(ind_out,1)=zeros(size(ind_out));
	1131	1061 vec_B(ind_out,2)=zeros(size(ind_out));
	1132	1062 vec_B(ind_out,3)=zeros(size(ind_out));
	1133	1063 A_out=reshape(vec_B,npY,npX,nbcolor);
	1134	1064 eval([<span class="string">'ProjData.'</span> VarName{ivar} <span class="string">'=squeeze(sum(A_out,1)/npY);'</span>]);
	1135	1065 <span class="keyword">end</span>
	1136	1066 ProjData.ListVarName=[ProjData.ListVarName VarName{ivar} ];
	1137	1067 ProjData.VarDimName=[ProjData.VarDimName {AXName}];<span class="comment">%to generalize with the initial name of the x coordinate</span>
	1138	1068 <span class="comment">%ProjData.VarDimIndex=[ProjData.VarDimIndex nbdim+1];</span>
	1139	1069 <span class="keyword">end</span>
	1140	1070 <span class="keyword">if</span> testU
	1141	1071 eval([<span class="string">'vector_x =ProjData.'</span> VarName{ivar_u} <span class="string">';'</span>])
	1142	1072 eval([<span class="string">'vector_y =ProjData.'</span> VarName{ivar_v} <span class="string">';'</span>])
	1143	1073 eval([<span class="string">'ProjData.'</span> VarName{ivar_u} <span class="string">'=cos(theta)vector_x+sin(theta)vector_y;'</span>])
	1144	1074 eval([<span class="string">'ProjData.'</span> VarName{ivar_v} <span class="string">'=-sin(theta)vector_x+cos(theta)vector_y;'</span>])
	1145	1075 <span class="keyword">end</span>
	1146	1076 <span class="comment">% ProjData.AX=(AXI)';</span>
	1147	1077 <span class="comment">%ProjData.ListDimName=[ProjData.ListDimName 'AX'];</span>
	1148	1078 <span class="comment">%ProjData.DimValue=[ProjData.DimValue length(ProjData.AX)];%values of dimension (nbre of vectors)</span>
	1149	1079 <span class="comment">%ProjData.VarDimIndex=[ProjData.VarDimIndex length(ProjData.ListDimName)];</span>
	1150	1080 ProjData.VarAttribute{nbvar+1}.long_name=<span class="string">'abscissa along line'</span>;
	1151	1081 <span class="comment">%ProjData.VarAttribute{nbvar+1}.Role='coord_x';</span>
	1152	1082 <span class="keyword">if</span> nbcolor==3
	1153	1083 <span class="comment">%ProjData.ListDimName=[ProjData.ListDimName 'rgb'];%GENERALIZE TO OTHER COMPONENT DIMENSIONS</span>
	1154	1084 <span class="comment">%ProjData.DimValue=[ProjData.DimValue 3];%values of dimension</span>
	1155	1085 <span class="comment">%ProjData.VarDimIndex{end}=[length(ProjData.ListDimName)-1 length(ProjData.ListDimName)];</span>
	1156	1086 ProjData.VarDimName{end}={XName,<span class="string">'rgb'</span>};
	1157	1087 <span class="keyword">end</span>
	1158	1088 <span class="keyword">end</span>
	1159	1089 <span class="keyword">end</span>
	1160	1090 <span class="keyword">end</span>
	1161	1091
	1162	1092 <span class="comment">% %shorter case for horizontal or vertical line (A FAIRE</span>
	1163	1093 <span class="comment">% % Rangx=[0.5 npx-0.5];%image coordiantes of corners</span>
	1164	1094 <span class="comment">% % Rangy=[npy-0.5 0.5];</span>
	1165	1095 <span class="comment">% % if isfield(Calib,'Pxcmx')&isfield(Calib,'Pxcmy')%old calib</span>
	1166	1096 <span class="comment">% % Rangx=Rangx/Calib.Pxcmx;</span>
	1167	1097 <span class="comment">% % Rangy=Rangy/Calib.Pxcmy;</span>
	1168	1098 <span class="comment">% % else</span>
	1169	1099 <span class="comment">% % [Rangx]=phys_XYZ(Calib,Rangx,[0.5 0.5],[0 0]);%case of translations without rotation and quadratic deformation</span>
	1170	1100 <span class="comment">% % [xx,Rangy]=phys_XYZ(Calib,[0.5 0.5],Rangy,[0 0]);</span>
	1171	1101 <span class="comment">% % end</span>
	1172	1102 <span class="comment">%</span>
	1173	1103 <span class="comment">% % test_scal=0;%default% 3- 'UserData':(get(handles.Tag,'UserData')</span>
	1174	1104
	1175	1105
	1176	1106 <span class="comment">%-----------------------------------------------------------------</span>
	1177	1107 <span class="comment">%project on a plane</span>
	1178	1108 <span class="comment">% AJOUTER flux,circul,error</span>
	1179	1109 <span class="keyword">function</span> [ProjData,errormsg] = <a href="#_sub4" class="code" title="subfunction [ProjData,errormsg] = proj_plane(FieldData, ObjectData)">proj_plane</a>(FieldData, ObjectData)
	1180	1110 <span class="comment">%-----------------------------------------------------------------</span>
	1181	1111
	1182	1112 <span class="comment">%initialisation of the input parameters of the projection plane</span>
	1183	1113 <span class="comment">%-----------------------------------------------------------------</span>
	1184	1114 ProjMode=<span class="string">'projection'</span>;<span class="comment">%direct projection by default</span>
	1185	1115 <span class="keyword">if</span> isfield(ObjectData,<span class="string">'ProjMode'</span>),ProjMode=ObjectData.ProjMode; <span class="keyword">end</span>;
	1186	1116
	1187	1117 <span class="comment">%axis origin</span>
	1188	1118 <span class="keyword">if</span> isempty(ObjectData.Coord)
	1189	1119 ObjectData.Coord(1,1)=0;<span class="comment">%origin of the plane set to [0 0] by default</span>
	1190	1120 ObjectData.Coord(1,2)=0;
	1191	1121 ObjectData.Coord(1,3)=0;
	1192	1122 <span class="keyword">end</span>
	1193	1123
	1194	1124 <span class="comment">%rotation angles</span>
	1195	1125 Phi=0;<span class="comment">%default</span>
	1196	1126 Theta=0;
	1197	1127 Psi=0;
	1198	1128 <span class="keyword">if</span> isfield(ObjectData,<span class="string">'Phi'</span>)&& ~isempty(ObjectData.Phi)
	1199	1129 Phi=(pi/180)*ObjectData.Phi;<span class="comment">%first Euler angle in radian</span>
	1200	1130 <span class="keyword">end</span>
	1201	1131 <span class="keyword">if</span> isfield(ObjectData,<span class="string">'Theta'</span>)&& ~isempty(ObjectData.Theta)
	1202	1132 Theta=(pi/180)*ObjectData.Theta;<span class="comment">%second Euler angle in radian</span>
	1203	1133 <span class="keyword">end</span>
	1204	1134 <span class="keyword">if</span> isfield(ObjectData,<span class="string">'Psi'</span>)&& ~isempty(ObjectData.Psi)
	1205	1135 Psi=(pi/180)*ObjectData.Psi;<span class="comment">%third Euler angle in radian</span>
	1206	1136 <span class="keyword">end</span>
	1207	1137
	1208	1138 <span class="comment">%components of the unity vector normal to the projection plane</span>
	1209	1139 NormVec_X=-sin(Phi)*sin(Theta);
	1210	1140 NormVec_Y=cos(Phi)*sin(Theta);
	1211	1141 NormVec_Z=cos(Theta);
	1212	1142
	1213	1143 <span class="comment">% test for 3D fields</span>
	1214	1144 test3D=0;
	1215	1145 <span class="keyword">if</span> isfield(FieldData,<span class="string">'NbDim'</span>)
	1216	1146 test3D=isequal(FieldData.NbDim,3);
	1217	1147 <span class="keyword">end</span>
	1218	1148 test3C=test3D; <span class="comment">%default 3 vel components</span>
	1219	1149
	1220	1150 <span class="comment">%mesh sizes DX and DY</span>
	1221	1151 DX=0;
	1222	1152 DY=0; <span class="comment">%default</span>
	1223	1153 <span class="keyword">if</span> isfield(ObjectData,<span class="string">'DX'</span>)&~isempty(ObjectData.DX)
	1224	1154 DX=abs(ObjectData.DX);<span class="comment">%mesh of interpolation points</span>
	1225	1155 <span class="keyword">end</span>
	1226	1156 <span class="keyword">if</span> isfield(ObjectData,<span class="string">'DY'</span>)&~isempty(ObjectData.DY)
	1227	1157 DY=abs(ObjectData.DY);<span class="comment">%mesh of interpolation points</span>
	1228	1158 <span class="keyword">end</span>
	1229	1159 <span class="keyword">if</span> ~isequal(ProjMode,<span class="string">'projection'</span>) & (DX==0\|DY==0)
	1230	1160 errormsg=<span class="string">'DX or DY missing'</span>;
	1231	1161 display(errormsg)
	1232	1162 <span class="keyword">return</span>
	1233	1163 <span class="keyword">end</span>
	1234	1164
	1235	1165 <span class="comment">%extrema along each axis</span>
	1236	1166 testXMin=0;
	1237	1167 testXMax=0;
	1238	1168 testYMin=0;
	1239	1169 testYMax=0;
	1240	1170 <span class="keyword">if</span> isfield(ObjectData,<span class="string">'RangeX'</span>)
	1241	1171 XMin=min(ObjectData.RangeX);
	1242	1172 XMax=max(ObjectData.RangeX);
	1243	1173 testXMin=XMax>XMin;
	1244	1174 testXMax=1;
	1245	1175 <span class="keyword">end</span>
	1246	1176 <span class="keyword">if</span> isfield(ObjectData,<span class="string">'RangeY'</span>)
	1247	1177 YMin=min(ObjectData.RangeY);
	1248	1178 YMax=max(ObjectData.RangeY);
	1249	1179 testYMin=YMax>YMin;
	1250	1180 testYMax=1;
	1251	1181 <span class="keyword">end</span>
	1252	1182 width=0;<span class="comment">%default width of the projection band</span>
	1253	1183 <span class="keyword">if</span> isfield(ObjectData,<span class="string">'RangeZ'</span>)
	1254	1184 width=max(ObjectData.RangeZ);
	1255	1185 <span class="keyword">end</span>
	1256	1186
	1257	1187 <span class="comment">% initiate Matlab structure for physical field</span>
	1258	1188 ProjData=<a href="#_sub5" class="code" title="subfunction [ProjData,errormsg]=proj_heading(FieldData,ObjectData)">proj_heading</a>(FieldData,ObjectData);
	1259	1189 ProjData.NbDim=2;
	1260	1190 ProjData.ListDimName={};<span class="comment">%name of dimension</span>
	1261	1191 ProjData.DimValue=[];<span class="comment">%values of dimension (nbre of vectors)</span>
	1262	1192 ProjData.ListVarName={};
	1263	1193 ProjData.VarDimIndex={};
	1264	1194
	1265	1195 error=0;<span class="comment">%default</span>
	1266	1196 flux=0;
	1267	1197 testfalse=0;
	1268	1198 ListIndex={};
	1269	1199
	1270	1200 <span class="comment">%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%</span>
	1271	1201 <span class="comment">%group the variables (fields of 'FieldData') in cells of variables with the same dimensions</span>
	1272	1202 <span class="comment">%-----------------------------------------------------------------</span>
	1273	1203 idimvar=0;
	1274	1204 [CellVarIndex,NbDim,VarTypeCell,errormsg]=<a href="find_field_indices.html" class="code" title="function [CellVarIndex,NbDim,VarType,errormsg]=find_field_indices(Data)">find_field_indices</a>(FieldData);
	1275	1205 <span class="keyword">if</span> ~isempty(errormsg)
	1276	1206 errormsg=[<span class="string">'error in proj_field/proj_plane:'</span> errormsg];
	1277	1207 <span class="keyword">return</span>
	1278	1208 <span class="keyword">end</span>
	1279	1209 <span class="comment">%LOOP ON GROUPS OF VARIABLES SHARING THE SAME DIMENSIONS</span>
	1280	1210 <span class="comment">% CellVarIndex=cells of variable index arrays</span>
	1281	1211 ivar_new=0; <span class="comment">% index of the current variable in the projected field</span>
	1282	1212 icoord=0;
	1283	1213 <span class="keyword">for</span> icell=1:length(CellVarIndex)
	1284	1214 <span class="keyword">if</span> NbDim(icell)==1
	1285	1215 <span class="keyword">continue</span>
	1286	1216 <span class="keyword">end</span>
	1287	1217 VarIndex=CellVarIndex{icell};<span class="comment">% indices of the selected variables in the list FieldData.ListVarName</span>
	1288	1218 VarType=VarTypeCell{icell};
	1289	1219 ivar_X=VarType.coord_x;
	1290	1220 ivar_Y=VarType.coord_y;
	1291	1221 ivar_Z=VarType.coord_z;
	1292	1222 ivar_U=VarType.vector_x;
	1293	1223 ivar_V=VarType.vector_y;
	1294	1224 ivar_W=VarType.vector_z;
	1295	1225 ivar_C=VarType.scalar ;
	1296	1226 ivar_Anc=VarType.ancillary;
	1297	1227 test_anc=zeros(size(VarIndex));
	1298	1228 test_anc(ivar_Anc)=ones(size(ivar_Anc));
	1299	1229 ivar_F=VarType.warnflag;
	1300	1230 ivar_FF=VarType.errorflag;
	1301	1231 <span class="keyword">if</span> isempty(ivar_X)
	1302	1232 test_grid=1;<span class="comment">%test for input data on regular grid (e.g. image)coordinates</span>
	1303	1233 <span class="keyword">else</span>
	1304	1234 <span class="keyword">if</span> length(ivar_Y)~=1
	1305	1235 <a href="warndlg_uvmat.html" class="code" title="function hwarn=warndlg_uvmat(warntext,title)">warndlg_uvmat</a>(<span class="string">'y coordinate missing in proj_field.m'</span>,<span class="string">'ERROR'</span>)
	1306	1236 <span class="keyword">return</span>
	1307	1237 <span class="keyword">end</span>
	1308	1238 test_grid=0;
	1309	1239 <span class="keyword">end</span>
	1310	1240 DimIndices=FieldData.VarDimIndex{VarIndex(1)};<span class="comment">%indices of the dimensions of the first variable (common to all variables in the cell)</span>
	1311	1241
	1312	1242 <span class="comment">%case of input fields with unstructured coordinates</span>
	1313	1243 <span class="keyword">if</span> ~test_grid
	1314	1244 <span class="keyword">if</span> length(ivar_X)==1
	1315	1245 eval([<span class="string">'coord_x=FieldData.'</span> FieldData.ListVarName{ivar_X} <span class="string">';'</span>])
	1316	1246 <span class="keyword">end</span>
	1317	1247 <span class="keyword">if</span> length(ivar_Y)==1
	1318	1248 eval([<span class="string">'coord_y=FieldData.'</span> FieldData.ListVarName{ivar_Y} <span class="string">';'</span>])
	1319	1249 <span class="keyword">end</span>
	1320	1250 <span class="keyword">if</span> length(ivar_Z)==1
	1321	1251 eval([<span class="string">'coord_z=FieldData.'</span> FieldData.ListVarName{ivar_Z} <span class="string">';'</span>])
	1322	1252 <span class="keyword">end</span>
	1323	1253 <span class="keyword">if</span> length(ivar_U)>1 \| length(ivar_V)>1 \| length(ivar_W)>1
	1324	1254 <a href="warndlg_uvmat.html" class="code" title="function hwarn=warndlg_uvmat(warntext,title)">warndlg_uvmat</a>(<span class="string">'multiple vector input in proj_field.m'</span>,<span class="string">'ERROR'</span>)
	1325	1255 <span class="keyword">return</span>
	1326	1256 <span class="keyword">end</span>
	1327	1257 <span class="keyword">if</span> length(ivar_X)>1 \| length(ivar_Y)>1 \| length(ivar_Z)>1
	1328	1258 <a href="warndlg_uvmat.html" class="code" title="function hwarn=warndlg_uvmat(warntext,title)">warndlg_uvmat</a>(<span class="string">'multiple coordinate input in proj_field.m'</span>,<span class="string">'ERROR'</span>)
	1329	1259 <span class="keyword">return</span>
	1330	1260 <span class="keyword">end</span>
	1331	1261 <span class="keyword">if</span> length(ivar_F)>1 \| length(ivar_FF)>1
	1332	1262 <a href="warndlg_uvmat.html" class="code" title="function hwarn=warndlg_uvmat(warntext,title)">warndlg_uvmat</a>(<span class="string">'multiple flag input in proj_field.m'</span>,<span class="string">'ERROR'</span>)
	1333	1263 <span class="keyword">return</span>
	1334	1264 <span class="keyword">end</span>
	1335	1265 <span class="keyword">if</span> length(ivar_Y)~=1
	1336	1266 <a href="warndlg_uvmat.html" class="code" title="function hwarn=warndlg_uvmat(warntext,title)">warndlg_uvmat</a>(<span class="string">'y coordinate missing in proj_field.m'</span>,<span class="string">'ERROR'</span>)
	1337	1267 <span class="keyword">return</span>
	1338	1268 <span class="keyword">end</span>
	1339	1269 <span class="comment">% translate initial coordinates</span>
	1340	1270 coord_x=coord_x-ObjectData.Coord(1,1);
	1341	1271 coord_y=coord_y-ObjectData.Coord(1,2);
	1342	1272 <span class="keyword">if</span> ~isempty(ivar_Z)
	1343	1273 coord_z=coord_z-ObjectData.Coord(1,3);
	1344	1274 <span class="keyword">end</span>
	1345	1275 <span class="comment">% selection of the vectors in the projection range (3D case)</span>
	1346	1276 <span class="keyword">if</span> length(ivar_Z)==1 & width > 0
	1347	1277 <span class="comment">%components of the unitiy vector normal to the projection plane</span>
	1348	1278 fieldZ=NormVec_Xcoord_x + NormVec_Ycoord_y+ NormVec_Z*coord_z;<span class="comment">% distance to the plane</span>
	1349	1279 indcut=find(abs(fieldZ) <= width);
	1350	1280 <span class="keyword">for</span> ivar=VarIndex
	1351	1281 VarName=FieldData.ListVarName{ivar};
	1352	1282 eval([<span class="string">'FieldData.'</span> VarName <span class="string">'=FieldData.'</span> VarName <span class="string">'(indcut);'</span>])
	1353	1283 <span class="comment">% end</span>
	1354	1284 <span class="comment">% A VOIR : CAS DE VAR STRUCTUREE MAIS PAS GRILLE REGULIERE : INTERPOLER SUR GRILLE REGULIERE</span>
	1355	1285 <span class="keyword">end</span>
	1356	1286 coord_x=coord_x(indcut);
	1357	1287 coord_y=coord_y(indcut);
	1358	1288 coord_z=coord_z(indcut);
	1359	1289 <span class="keyword">end</span>
	1360	1290
	1361	1291 <span class="comment">%rotate coordinates if needed</span>
	1362	1292 <span class="keyword">if</span> isequal(Phi,0)
	1363	1293 coord_X=coord_x;
	1364	1294 coord_Y=coord_y;
	1365	1295 <span class="keyword">if</span> ~isequal(Theta,0)
	1366	1296 coord_Y=coord_Y *cos(Theta);
	1367	1297 <span class="keyword">end</span>
	1368	1298 <span class="keyword">else</span>
	1369	1299 coord_X=(coord_x cos(Phi) + coord_y sin(Phi));
	1370	1300 coord_Y=(-coord_x sin(Phi) + coord_y cos(Phi))*cos(Theta);
	1371	1301 <span class="keyword">end</span>
	1372	1302 <span class="keyword">if</span> ~isempty(ivar_Z)
	1373	1303 coord_Y=coord_Y+coord_z *sin(Theta);
	1374	1304 <span class="keyword">end</span>
	1375	1305 <span class="keyword">if</span> ~isequal(Psi,0)
	1376	1306 coord_X=(coord_X cos(Psi) - coord_Y sin(Psi));<span class="comment">%A VERIFIER</span>
	1377	1307 coord_Y=(coord_X sin(Psi) + coord_Y cos(Psi));
	1378	1308 <span class="keyword">end</span>
	1379	1309
	1380	1310 <span class="comment">%restriction to the range of x and y if imposed</span>
	1381	1311 testin=ones(size(coord_X)); <span class="comment">%default</span>
	1382	1312 testbound=0;
	1383	1313 <span class="keyword">if</span> testXMin
	1384	1314 testin=testin & (coord_X >= XMin);
	1385	1315 testbound=1;
	1386	1316 <span class="keyword">end</span>
	1387	1317 <span class="keyword">if</span> testXMax
	1388	1318 testin=testin & (coord_X <= XMax);
	1389	1319 testbound=1;
	1390	1320 <span class="keyword">end</span>
	1391	1321 <span class="keyword">if</span> testYMin
	1392	1322 testin=testin & (coord_Y >= YMin);
	1393	1323 testbound=1;
	1394	1324 <span class="keyword">end</span>
	1395	1325 <span class="keyword">if</span> testYMin
	1396	1326 testin=testin & (coord_Y <= YMax);
	1397	1327 testbound=1;
	1398	1328 <span class="keyword">end</span>
	1399	1329 <span class="keyword">if</span> testbound
	1400	1330 indcut=find(testin);
	1401	1331 <span class="keyword">for</span> ivar=VarIndex
	1402	1332 VarName=FieldData.ListVarName{ivar};
	1403	1333 eval([<span class="string">'FieldData.'</span> VarName <span class="string">'=FieldData.'</span> VarName <span class="string">'(indcut);'</span>])
	1404	1334 <span class="keyword">end</span>
	1405	1335 coord_X=coord_X(indcut);
	1406	1336 coord_Y=coord_Y(indcut);
	1407	1337 <span class="keyword">if</span> length(ivar_Z)==1
	1408	1338 coord_Z=coord_Z(indcut);
	1409	1339 <span class="keyword">end</span>
	1410	1340 <span class="keyword">end</span>
	1411	1341 <span class="comment">% different cases of projection</span>
	1412	1342 <span class="keyword">if</span> isequal(ObjectData.ProjMode,<span class="string">'projection'</span>)
	1413	1343 ProjData.ListDimName=[ProjData.ListDimName FieldData.ListDimName(DimIndices(1))];<span class="comment">%add the point index to the list of dimensions</span>
	1414	1344 ProjData.DimValue=[ProjData.DimValue length(coord_X)];
	1415	1345 nbvar=0;
	1416	1346 <span class="keyword">for</span> ivar=VarIndex <span class="comment">%transfer variables to the projection plane</span>
	1417	1347 VarName=FieldData.ListVarName{ivar};
	1418	1348 <span class="keyword">if</span> ivar==ivar_X <span class="comment">%x coordinate</span>
	1419	1349 eval([<span class="string">'ProjData.'</span> VarName <span class="string">'=coord_X;'</span>])
	1420	1350 <span class="keyword">elseif</span> ivar==ivar_Y <span class="comment">% y coordinate</span>
	1421	1351 eval([<span class="string">'ProjData.'</span> VarName <span class="string">'=coord_Y;'</span>])
	1422	1352 <span class="keyword">elseif</span> isempty(ivar_Z) \|\| ivar~=ivar_Z <span class="comment">% other variables (except Z coordinate wyhich is not reproduced)</span>
	1423	1353 eval([<span class="string">'ProjData.'</span> VarName <span class="string">'=FieldData.'</span> VarName <span class="string">';'</span>])
	1424	1354 <span class="keyword">end</span>
	1425	1355 <span class="keyword">if</span> isempty(ivar_Z) \|\| ivar~=ivar_Z
	1426	1356 ProjData.ListVarName=[ProjData.ListVarName VarName];
	1427	1357 ProjData.VarDimIndex=[ProjData.VarDimIndex DimIndices(1)];
	1428	1358 nbvar=nbvar+1;
	1429	1359 <span class="keyword">if</span> isfield(FieldData,<span class="string">'VarAttribute'</span>) & length(FieldData.VarAttribute) >=ivar
	1430	1360 ProjData.VarAttribute{nbvar}=FieldData.VarAttribute{ivar};
	1431	1361 <span class="keyword">end</span>
	1432	1362 <span class="keyword">end</span>
	1433	1363 <span class="keyword">end</span>
	1434	1364 <span class="keyword">elseif</span> isequal(ObjectData.ProjMode,<span class="string">'interp'</span>)\|\|isequal(ObjectData.ProjMode,<span class="string">'filter'</span>)<span class="comment">%interpolate data on a regular grid</span>
	1435	1365 coord_x_proj=[XMin:DX:XMax];
	1436	1366 coord_y_proj=[YMin:DY:YMax];
	1437	1367 ProjData.ListDimName={<span class="string">'coord_y'</span>,<span class="string">'coord_x'</span>};
	1438	1368 ProjData.DimValue=[length(coord_y_proj) length(coord_x_proj)];
	1439	1369 ProjData.ListVarName={};
	1440	1370 ProjData.VarDimIndex={};
	1441	1371 <span class="keyword">if</span> isempty(ivar_X), ivar_X=0; <span class="keyword">end</span>;
	1442	1372 <span class="keyword">if</span> isempty(ivar_Y), ivar_Y=0; <span class="keyword">end</span>;
	1443	1373 <span class="keyword">if</span> isempty(ivar_Z), ivar_Z=0; <span class="keyword">end</span>;
	1444	1374 <span class="keyword">if</span> isempty(ivar_U), ivar_U=0; <span class="keyword">end</span>;
	1445	1375 <span class="keyword">if</span> isempty(ivar_V), ivar_V=0; <span class="keyword">end</span>;
	1446	1376 <span class="keyword">if</span> isempty(ivar_W), ivar_W=0; <span class="keyword">end</span>;
	1447	1377 <span class="keyword">if</span> isempty(ivar_F), ivar_F=0; <span class="keyword">end</span>;
	1448	1378 <span class="keyword">if</span> isempty(ivar_FF), ivar_FF=0; <span class="keyword">end</span>;
	1449	1379 <span class="keyword">if</span> ~isequal(ivar_FF,0)
	1450	1380 VarName_FF=FieldData.ListVarName{ivar_FF};
	1451	1381 eval([<span class="string">'indsel=find(FieldData.'</span> VarName_FF <span class="string">'==0);'</span>])
	1452	1382 coord_X=coord_X(indsel);
	1453	1383 coord_Y=coord_Y(indsel);
	1454	1384 <span class="keyword">end</span>
	1455	1385 FF=zeros(1,length(coord_y_proj)*length(coord_x_proj));
	1456	1386 testFF=0;
	1457	1387 <span class="keyword">for</span> ivar=VarIndex
	1458	1388 VarName=FieldData.ListVarName{ivar};
	1459	1389 <span class="keyword">if</span> ~( ivar==ivar_X \| ivar==ivar_Y \| ivar==ivar_Z \| ivar==ivar_F \| ivar==ivar_FF \| test_anc(ivar)==1)
	1460	1390 ivar_new=ivar_new+1;
	1461	1391 ProjData.ListVarName=[ProjData.ListVarName {VarName}];
	1462	1392 ProjData.VarDimIndex=[ProjData.VarDimIndex {[1 2]}];
	1463	1393 <span class="keyword">if</span> isfield(FieldData,<span class="string">'VarAttribute'</span>) & length(FieldData.VarAttribute) >=ivar
	1464	1394 ProjData.VarAttribute{ivar_new}=FieldData.VarAttribute{ivar};
	1465	1395 <span class="keyword">end</span>
	1466	1396 ProjData.VarAttribute{ivar_new}.Coord_2=[XMin XMax];
	1467	1397 ProjData.VarAttribute{ivar_new}.Coord_1=[YMin YMax];
	1468	1398 <span class="keyword">if</span> ~isequal(ivar_FF,0)
	1469	1399 eval([<span class="string">'FieldData.'</span> VarName <span class="string">'=FieldData.'</span> VarName <span class="string">'(indsel);'</span>])
	1470	1400 <span class="keyword">end</span>
	1471	1401 eval([<span class="string">'ProjData.'</span> VarName <span class="string">'=griddata_uvmat(coord_X,coord_Y,FieldData.'</span> VarName <span class="string">',coord_x_proj,coord_y_proj'');'</span>])
	1472	1402 eval([<span class="string">'varline=reshape(ProjData.'</span> VarName <span class="string">',1,length(coord_y_proj)*length(coord_x_proj));'</span>])
	1473	1403 FFlag= isnan(varline); <span class="comment">%detect undefined values NaN</span>
	1474	1404 indnan=find(FFlag);
	1475	1405 <span class="keyword">if</span>~isempty(indnan)
	1476	1406 varline(indnan)=zeros(size(indnan));
	1477	1407 eval([<span class="string">'ProjData.'</span> VarName <span class="string">'=reshape(varline,length(coord_y_proj),length(coord_x_proj));'</span>])
	1478	1408 FF(indnan)=ones(size(indnan));
	1479	1409 <span class="comment">% eval(['ProjData.' VarName '(indnan)=zeros(size(indnan));'])%put NaN to 0</span>
	1480	1410 <span class="comment">% FF=FF\|FFlag;</span>
	1481	1411 testFF=1;
	1482	1412 <span class="keyword">end</span>
	1483	1413 <span class="keyword">if</span> ivar==ivar_U
	1484	1414 ivar_U=ivar_new;
	1485	1415 <span class="keyword">end</span>
	1486	1416 <span class="keyword">if</span> ivar==ivar_V
	1487	1417 ivar_V=ivar_new;
	1488	1418 <span class="keyword">end</span>
	1489	1419 <span class="keyword">if</span> ivar==ivar_W
	1490	1420 ivar_W=ivar_new;
	1491	1421 <span class="keyword">end</span>
	1492	1422 <span class="keyword">end</span>
	1493	1423 <span class="keyword">end</span>
	1494	1424 <span class="keyword">if</span> testFF
	1495	1425 ProjData.FF=reshape(FF,length(coord_y_proj),length(coord_x_proj));
	1496	1426 ProjData.ListVarName=[ProjData.ListVarName {<span class="string">'FF'</span>}];
	1497	1427 ProjData.VarDimIndex=[ProjData.VarDimIndex {[1 2]}];
	1498	1428 ProjData.VarAttribute{ivar_new+1}.Role=<span class="string">'errorflag'</span>;
	1499	1429 <span class="keyword">end</span>
	1500	1430 <span class="keyword">end</span>
	1501	1431 <span class="comment">%case of fields defined on a structured grid</span>
	1502	1432 <span class="keyword">else</span>
	1503	1433 DimValue=FieldData.DimValue(DimIndices);<span class="comment">%set of dimension values</span>
	1504	1434 ListDimName=FieldData.ListDimName(DimIndices);
	1505	1435 nbcolor=1; <span class="comment">%default</span>
	1506	1436 <span class="keyword">for</span> idim=1:length(ListDimName)
	1507	1437 DimName=ListDimName{idim};
	1508	1438 <span class="keyword">if</span> isequal(DimName,<span class="string">'rgb'</span>)\|isequal(DimName,<span class="string">'nb_coord'</span>)\|isequal(DimName,<span class="string">'nb_coord_i'</span>)
	1509	1439 nbcolor=DimValue(idim);
	1510	1440 DimIndices(idim)=[];
	1511	1441 DimValue(idim)=[];
	1512	1442 <span class="keyword">end</span>
	1513	1443 <span class="keyword">if</span> isequal(DimName,<span class="string">'nb_coord_j'</span>)<span class="comment">% NOTE: CASE OF TENSOR NOT TREATED</span>
	1514	1444 DimIndices(idim)=[];
	1515	1445 DimValue(idim)=[];
	1516	1446 <span class="keyword">end</span>
	1517	1447 <span class="keyword">end</span>
	1518	1448 ind_1=find(DimValue==1);
	1519	1449 DimIndices(ind_1)=[]; <span class="comment">%suppress singleton dimensions</span>
	1520	1450 indxy=find(DimVarIndex(DimIndices));<span class="comment">%select dimension variables (DimIndices non zero)</span>
	1521	1451 nb_dim=length(DimIndices);<span class="comment">%number of space dimensions</span>
	1522	1452 Coord_z=[];
	1523	1453 Coord_y=[];
	1524	1454 Coord_x=[];
	1525	1455
	1526	1456 <span class="keyword">for</span> idim=1:nb_dim <span class="comment">%loop on space dimensions</span>
	1527	1457 test_interp(idim)=0;<span class="comment">%test for coordiate interpolation (non regular grid), =0 by default</span>
	1528	1458 test_coord(idim)=0;<span class="comment">%test for defined coordinates, =0 by default</span>
	1529	1459 ivar=DimVarIndex(DimIndices(idim));<span class="comment">% index of the variable corresponding to the current dimension</span>
	1530	1460 <span class="keyword">if</span> ~isequal(ivar,0)<span class="comment">% a variable corresponds to the current dimension</span>
	1531	1461 eval([<span class="string">'Coord{idim}=FieldData.'</span> FieldData.ListVarName{ivar} <span class="string">';'</span>]) ;<span class="comment">% position for the first index</span>
	1532	1462 DCoord=diff(Coord{idim});
	1533	1463 DCoord_min(idim)=min(DCoord);
	1534	1464 DCoord_max=max(DCoord);
	1535	1465 test_direct(idim)=DCoord_max>0;<span class="comment">% =1 for increasing values, 0 otherwise</span>
	1536	1466 test_direct_min=DCoord_min(idim)>0;<span class="comment">% =1 for increasing values, 0 otherwise</span>
	1537	1467 <span class="keyword">if</span> ~isequal(test_direct(idim),test_direct_min)
	1538	1468 <a href="warndlg_uvmat.html" class="code" title="function hwarn=warndlg_uvmat(warntext,title)">warndlg_uvmat</a>([<span class="string">'non monotonic dimension variable # '</span> num2str(idim) <span class="string">' in proj_field.m'</span>],<span class="string">'ERROR'</span>)
	1539	1469 <span class="keyword">return</span>
	1540	1470 <span class="keyword">end</span>
	1541	1471 test_interp(idim)=(DCoord_max-DCoord_min(idim))> 0.0001*abs(DCoord_max);<span class="comment">% test grid regularity</span>
	1542	1472 test_coord(idim)=1;
	1543	1473
	1544	1474 <span class="keyword">else</span> <span class="comment">% no variable associated with the first dimension, look for variable attributes Coord_1, _2 or _3</span>
	1545	1475 Coord_i_str=[<span class="string">'Coord_'</span> num2str(idim)];
	1546	1476 DCoord_min(idim)=1;<span class="comment">%default</span>
	1547	1477 Coord{idim}=[0.5 DimValue(idim)-0.5];
	1548	1478 test_direct(idim)=1;
	1549	1479 <span class="keyword">for</span> ivar=VarIndex
	1550	1480 <span class="keyword">if</span> isfield(FieldData.VarAttribute{ivar},Coord_i_str)<span class="comment">% if there is a variable attribute named Coord_1, _2 or _3</span>
	1551	1481 eval([<span class="string">'Coord{idim}=FieldData.VarAttribute{ivar}.'</span> Coord_i_str <span class="string">';'</span>]);<span class="comment">%'range x</span>
	1552	1482 <span class="keyword">if</span> isnumeric(Coord{idim})
	1553	1483 <span class="keyword">if</span> length(Coord{idim})>=2
	1554	1484 test_direct(idim)=(Coord{idim}(2)>Coord{idim}(1));
	1555	1485 <span class="keyword">else</span>
	1556	1486 <a href="warndlg_uvmat.html" class="code" title="function hwarn=warndlg_uvmat(warntext,title)">warndlg_uvmat</a>([<span class="string">'two values needed for '</span> Coord_i_str <span class="string">'in proj_field.m'</span>],<span class="string">'ERROR'</span>)
	1557	1487 <span class="keyword">return</span>
	1558	1488 <span class="keyword">end</span>
	1559	1489 <span class="keyword">else</span>
	1560	1490 <a href="warndlg_uvmat.html" class="code" title="function hwarn=warndlg_uvmat(warntext,title)">warndlg_uvmat</a>([<span class="string">'non numerical coordinate attributes'</span> Coord_i_str <span class="string">'in proj_field.m'</span>],<span class="string">'ERROR'</span>)
	1561	1491 <span class="keyword">return</span>
	1562	1492 <span class="keyword">end</span>
	1563	1493 DCoord_min(idim)=(Coord{idim}(end)-Coord{idim}(1))/(DimValue(idim)-1);
	1564	1494 <span class="keyword">end</span>
	1565	1495 <span class="keyword">end</span>
	1566	1496 <span class="keyword">end</span>
	1567	1497 <span class="keyword">end</span>
	1568	1498 <span class="keyword">if</span> nb_dim==2
	1569	1499 <span class="keyword">if</span> DY==0
	1570	1500 DY=abs(DCoord_min(1));
	1571	1501 <span class="keyword">end</span>
	1572	1502 npY=1+round(abs(Coord{1}(end)-Coord{1}(1))/DY);<span class="comment">%nbre of points after interpolation</span>
	1573	1503 npy=1+round(abs(Coord{1}(end)-Coord{1}(1))/abs(DCoord_min(1)));<span class="comment">%nbre of points after possible interpolation on a regular grid</span>
	1574	1504 <span class="keyword">if</span> DX==0
	1575	1505 DX=abs(DCoord_min(2));
	1576	1506 <span class="keyword">end</span>
	1577	1507 npX=1+round(abs(Coord{2}(end)-Coord{2}(1))/DX);<span class="comment">%nbre of points after interpol</span>
	1578	1508 npx=1+round(abs(Coord{2}(end)-Coord{2}(1))/abs(DCoord_min(2)));<span class="comment">%nbre of points after possible interpolation on a regular grid</span>
	1579	1509 Coord_y=linspace(Coord{1}(1),Coord{1}(end),npY);
	1580	1510 test_direct_y=test_direct(1);
	1581	1511 Coord_x=linspace(Coord{2}(1),Coord{2}(end),npX);
	1582	1512 test_direct_x=test_direct(2);
	1583	1513 DAX=DCoord_min(2);
	1584	1514 DAY=DCoord_min(1);
	1585	1515 <span class="keyword">elseif</span> nb_dim==3
	1586	1516 DZ=abs(DCoord_min(1));
	1587	1517 npz=1+round(abs(Coord{1}(end)-Coord{1}(1))/DZ);<span class="comment">%nbre of points after interpolation</span>
	1588	1518 <span class="keyword">if</span> DY==0
	1589	1519 DY=abs(DCoord_min(2));
	1590	1520 <span class="keyword">end</span>
	1591	1521 npY=1+round(abs(Coord{2}(end)-Coord{2}(1))/DY);<span class="comment">%nbre of points after interpol</span>
	1592	1522 npy=1+round(abs(Coord{2}(end)-Coord{2}(1))/abs(DCoord_min(2)));<span class="comment">%nbre of points before interpol</span>
	1593	1523 <span class="keyword">if</span> DX==0
	1594	1524 DX=abs(DCoord_min(3));
	1595	1525 <span class="keyword">end</span>
	1596	1526 npX=1+round(abs(Coord{3}(end)-Coord{3}(1))/DX);<span class="comment">%nbre of points after interpol</span>
	1597	1527 npx=1+round(abs(Coord{3}(end)-Coord{3}(1))/abs(DCoord_min(3)));<span class="comment">%nbre of points before interpol</span>
	1598	1528 Coord_z=linspace(Coord{1}(1),Coord{1}(end),npz);
	1599	1529 test_direct_z=test_direct(1);
	1600	1530 Coord_y=linspace(Coord{2}(1),Coord{2}(end),npY);
	1601	1531 test_direct_y=test_direct(2);
	1602	1532 Coord_x=linspace(Coord{3}(1),Coord{3}(end),npX);
	1603	1533 test_direct_x=test_direct(3);
	1604	1534 <span class="keyword">end</span>
	1605	1535 minAX=min(Coord_x);
	1606	1536 maxAX=max(Coord_x);
	1607	1537 minAY=min(Coord_y);
	1608	1538 maxAY=max(Coord_y);
	1609	1539 xcorner=[minAX maxAX minAX maxAX]-ObjectData.Coord(1,1);
	1610	1540 ycorner=[maxAY maxAY minAY minAY]-ObjectData.Coord(1,2);
	1611	1541 xcor_new=xcornercos(Phi)+ycornersin(Phi);<span class="comment">%coord new frame</span>
	1612	1542 ycor_new=-xcornersin(Phi)+ycornercos(Phi);
	1613	1543 <span class="keyword">if</span> ~testXMax
	1614	1544 XMax=max(xcor_new);
	1615	1545 <span class="keyword">end</span>
	1616	1546 <span class="keyword">if</span> ~testXMin
	1617	1547 XMin=min(xcor_new);
	1618	1548 <span class="keyword">end</span>
	1619	1549 <span class="keyword">if</span> ~testYMax
	1620	1550 YMax=max(ycor_new);
	1621	1551 <span class="keyword">end</span>
	1622	1552 <span class="keyword">if</span> ~testYMin
	1623	1553 YMin=min(ycor_new);
	1624	1554 <span class="keyword">end</span>
	1625	1555 DXinit=(maxAX-minAX)/(npx-1);
	1626	1556 DYinit=(maxAY-minAY)/(npy-1);
	1627	1557 <span class="keyword">if</span> DX==0
	1628	1558 DX=DXinit;
	1629	1559 <span class="keyword">end</span>
	1630	1560 <span class="keyword">if</span> DY==0
	1631	1561 DY=DYinit;
	1632	1562 <span class="keyword">end</span>
	1633	1563 npX=floor((XMax-XMin)/DX+1);
	1634	1564 npY=floor((YMax-YMin)/DY+1);
	1635	1565 <span class="keyword">if</span> test_direct_y
	1636	1566 coord_y_proj=linspace(YMin,YMax,npY);<span class="comment">%abscissa of the new pixels along the line</span>
	1637	1567 <span class="keyword">else</span>
	1638	1568 coord_y_proj=linspace(YMax,YMin,npY);<span class="comment">%abscissa of the new pixels along the line</span>
	1639	1569 <span class="keyword">end</span>
	1640	1570 <span class="keyword">if</span> test_direct_x
	1641	1571 coord_x_proj=linspace(XMin,XMax,npX);<span class="comment">%abscissa of the new pixels along the line</span>
	1642	1572 <span class="keyword">else</span>
	1643	1573 coord_x_proj=linspace(XMax,XMin,npX);<span class="comment">%abscissa of the new pixels along the line</span>
	1644	1574 <span class="keyword">end</span>
	1645	1575
	1646	1576 <span class="comment">% case with no rotation and interpolation</span>
	1647	1577 <span class="keyword">if</span> isequal(ProjMode,<span class="string">'projection'</span>) && isequal(Phi,0) && isequal(Theta,0) && isequal(Psi,0)
	1648	1578 <span class="keyword">if</span> test_direct(1)
	1649	1579 min_ind1=ceil((YMin-Coord{1}(1))/DYinit)+1;
	1650	1580 max_ind1=floor((YMax-Coord{1}(1))/DYinit)+1;
	1651	1581 Ybound(1)=Coord{1}(1)+DYinit*(min_ind1-1);
	1652	1582 Ybound(2)=Coord{1}(1)+DYinit*(max_ind1-1);
	1653	1583 <span class="keyword">else</span>
	1654	1584 min_ind1=ceil((Coord{1}(1)-YMax)/DYinit)+1;
	1655	1585 max_ind1=floor((Coord{1}(1)-YMin)/DYinit)+1;
	1656	1586 Ybound(2)=Coord{1}(1)-DYinit*(max_ind1-1);
	1657	1587 Ybound(1)=Coord{1}(1)-DYinit*(min_ind1-1);
	1658	1588 <span class="keyword">end</span>
	1659	1589 <span class="keyword">if</span> test_direct(2)==1
	1660	1590 min_ind2=ceil((XMin-Coord{2}(1))/DXinit)+1;
	1661	1591 max_ind2=floor((XMax-Coord{2}(1))/DXinit)+1;
	1662	1592 Xbound(1)=Coord{2}(1)+DXinit*(min_ind2-1);
	1663	1593 Xbound(2)=Coord{2}(1)+DXinit*(max_ind2-1);
	1664	1594 <span class="keyword">else</span>
	1665	1595 min_ind2=ceil((Coord{2}(1)-XMax)/DXinit)+1;
	1666	1596 max_ind2=floor((Coord{2}(1)-XMin)/DXinit)+1;
	1667	1597 Xbound(2)=Coord{2}(1)+DXinit*(max_ind2-1);
	1668	1598 Xbound(1)=Coord{2}(1)+DXinit*(min_ind2-1);
	1669	1599 <span class="keyword">end</span>
	1670	1600 min_ind1=max(min_ind1,1);<span class="comment">% deals with margin (bound lower than the first index)</span>
	1671	1601 min_ind2=max(min_ind2,1);
	1672	1602 max_ind1=min(max_ind1,npy);
	1673	1603 max_ind2=min(max_ind2,npx);
	1674	1604 <span class="keyword">for</span> ivar=VarIndex
	1675	1605 VarName=FieldData.ListVarName{ivar};
	1676	1606 <span class="keyword">if</span> isequal(ObjectData.ProjMode,<span class="string">'interp'</span>)\|\|isequal(ObjectData.ProjMode,<span class="string">'filter'</span>)<span class="comment">% coordinates common to all fields</span>
	1677	1607 ProjData.ListDimName={<span class="string">'coord_y'</span>,<span class="string">'coord_x'</span>};
	1678	1608 ProjData.DimValue=[length(coord_y_proj) length(coord_x_proj)];
	1679	1609 <span class="keyword">else</span>
	1680	1610 icoord=icoord+1;
	1681	1611 ProjData.ListDimName=[ProjData.ListDimName {[<span class="string">'coord_y_'</span> num2str(icoord)],[<span class="string">'coord_x_'</span> num2str(icoord)]}];
	1682	1612 ProjData.DimValue=[ProjData.DimValue length(coord_y_proj) length(coord_x_proj)];
	1683	1613 <span class="keyword">end</span>
	1684	1614 ProjData.ListVarName=[ProjData.ListVarName VarName];
	1685	1615 ProjData.VarDimIndex=[ProjData.VarDimIndex [nb_dim-1 nb_dim]];
	1686	1616 <span class="keyword">if</span> length(FieldData.VarAttribute)>=ivar
	1687	1617 ProjData.VarAttribute{length(ProjData.ListVarName)}=FieldData.VarAttribute{ivar};
	1688	1618 <span class="keyword">end</span>
	1689	1619 ProjData.VarAttribute{length(ProjData.ListVarName)}.Coord_2=Xbound;
	1690	1620 ProjData.VarAttribute{length(ProjData.ListVarName)}.Coord_1=Ybound;
	1691	1621 eval([<span class="string">'ProjData.'</span> VarName <span class="string">'=FieldData.'</span> VarName <span class="string">'(min_ind1:max_ind1,min_ind2:max_ind2) ;'</span>]);
	1692	1622 <span class="keyword">end</span>
	1693	1623 <span class="keyword">else</span>
	1694	1624 <span class="comment">% case with rotation and/or interpolation</span>
	1695	1625 <span class="keyword">if</span> isempty(Coord_z) <span class="comment">%2D case</span>
	1696	1626 [X,Y]=meshgrid(coord_x_proj,coord_y_proj);<span class="comment">%grid in the new coordinates</span>
	1697	1627 XIMA=ObjectData.Coord(1,1)+(X)cos(Phi)-Ysin(Phi);<span class="comment">%corresponding coordinates in the original image</span>
	1698	1628 YIMA=ObjectData.Coord(1,2)+(X)sin(Phi)+Ycos(Phi);
	1699	1629 XIMA=(XIMA-minAX)/DXinit+1;<span class="comment">% image index along x</span>
	1700	1630 YIMA=(-YIMA+maxAY)/DYinit+1;<span class="comment">% image index along y</span>
	1701	1631 XIMA=reshape(round(XIMA),1,npX*npY);<span class="comment">%indices reorganized in 'line'</span>
	1702	1632 YIMA=reshape(round(YIMA),1,npX*npY);
	1703	1633 flagin=XIMA>=1 & XIMA<=npx & YIMA >=1 & YIMA<=npy;<span class="comment">%flagin=1 inside the original image</span>
	1704	1634 <span class="keyword">if</span> isequal(ObjectData.ProjMode,<span class="string">'filter'</span>)
	1705	1635 npx_filter=ceil(abs(DX/DAX));
	1706	1636 npy_filter=ceil(abs(DY/DAY));
	1707	1637 Mfilter=ones(npy_filter,npx_filter)/(npx_filter*npy_filter);
	1708	1638 test_filter=1;
	1709	1639 <span class="keyword">else</span>
	1710	1640 test_filter=0;
	1711	1641 <span class="keyword">end</span>
	1712	1642 <span class="keyword">for</span> ivar=VarIndex
	1713	1643 VarName=FieldData.ListVarName{ivar} ;
	1714	1644 <span class="keyword">if</span> test_interp(1) \| test_interp(2)<span class="comment">%interpolate on a regular grid</span>
	1715	1645 eval([<span class="string">'FieldData.'</span> VarName <span class="string">'=interp2(Coord{2},Coord{1},FieldData.'</span> VarName <span class="string">',Coord_x,Coord_y'');'</span>]) <span class="comment">%TO TEST</span>
	1716	1646 <span class="keyword">end</span>
	1717	1647 <span class="comment">%filter the field (image) if option 'filter' is used</span>
	1718	1648 <span class="keyword">if</span> test_filter
	1719	1649 Aclass=class(FieldData.A);
	1720	1650 eval([<span class="string">'FieldData.'</span> VarName <span class="string">'=filter2(Mfilter,FieldData.'</span> VarName <span class="string">',''valid'');'</span>])
	1721	1651 <span class="keyword">if</span> ~isequal(Aclass,<span class="string">'double'</span>)
	1722	1652 eval([<span class="string">'FieldData.'</span> VarName <span class="string">'='</span> Aclass <span class="string">'(FieldData.'</span> VarName <span class="string">');'</span>])<span class="comment">%revert to integer values</span>
	1723	1653 <span class="keyword">end</span>
	1724	1654 <span class="keyword">end</span>
	1725	1655 eval([<span class="string">'vec_A=reshape(FieldData.'</span> VarName <span class="string">',npx*npy,nbcolor);'</span>])<span class="comment">%put the original image in line</span>
	1726	1656 ind_in=find(flagin);
	1727	1657 ind_out=find(~flagin);
	1728	1658 ICOMB=(XIMA-1)*npy+YIMA;
	1729	1659 ICOMB=ICOMB(flagin);<span class="comment">%index corresponding to XIMA and YIMA in the aligned original image vec_A</span>
	1730	1660 vec_B(ind_in,[1:nbcolor])=vec_A(ICOMB,:);
	1731	1661 <span class="keyword">for</span> icolor=1:nbcolor
	1732	1662 vec_B(ind_out,icolor)=zeros(size(ind_out));
	1733	1663 <span class="keyword">end</span>
	1734	1664 <span class="keyword">if</span> isequal(ObjectData.ProjMode,<span class="string">'interp'</span>) \|\| isequal(ObjectData.ProjMode,<span class="string">'filter'</span>)<span class="comment">% coordinates common to all fields</span>
	1735	1665 ProjData.ListDimName={<span class="string">'coord_y'</span>,<span class="string">'coord_x'</span>};
	1736	1666 ProjData.DimValue=[length(coord_y_proj) length(coord_x_proj)];
	1737	1667 <span class="keyword">else</span>
	1738	1668 icoord=icoord+1;
	1739	1669 ProjData.ListDimName=[ProjData.ListDimName {[<span class="string">'coord_y_'</span> num2str(icoord)],[<span class="string">'coord_x_'</span> num2str(icoord)]}];
	1740	1670 ProjData.DimValue=[ProjData.DimValue length(coord_y_proj) length(coord_x_proj)];
	1741	1671 <span class="keyword">end</span>
	1742	1672 ProjData.ListVarName=[ProjData.ListVarName VarName];
	1743	1673 ProjData.VarDimIndex=[ProjData.VarDimIndex [nb_dim-1 nb_dim]];
	1744	1674 <span class="keyword">if</span> length(FieldData.VarAttribute)>=ivar
	1745	1675 ProjData.VarAttribute{length(ProjData.ListVarName)}=FieldData.VarAttribute{ivar};
	1746	1676 <span class="keyword">end</span>
	1747	1677 <span class="keyword">if</span> test_direct(2)==1
	1748	1678 ProjData.VarAttribute{length(ProjData.ListVarName)}.Coord_2=[XMin XMax];
	1749	1679 <span class="keyword">else</span>
	1750	1680 ProjData.VarAttribute{length(ProjData.ListVarName)}.Coord_2=[XMax XMin];
	1751	1681 <span class="keyword">end</span>
	1752	1682 <span class="keyword">if</span> test_direct(1)==1
	1753	1683 ProjData.VarAttribute{length(ProjData.ListVarName)}.Coord_1=[YMin YMax];
	1754	1684 <span class="keyword">else</span>
	1755	1685 ProjData.VarAttribute{length(ProjData.ListVarName)}.Coord_1=[YMax YMin];
	1756	1686 <span class="keyword">end</span>
	1757	1687 eval([<span class="string">'ProjData.'</span> VarName <span class="string">'=reshape(vec_B,npY,npX,nbcolor);'</span>]);
	1758	1688 <span class="keyword">end</span>
	1759	1689 ProjData.FF=reshape(~flagin,npY,npX);<span class="comment">%false flag A FAIRE: tenir compte d'un flga antérieur</span>
	1760	1690 ProjData.ListVarName=[ProjData.ListVarName <span class="string">'FF'</span>];
	1761	1691 ProjData.VarDimIndex=[ProjData.VarDimIndex [nb_dim-1 nb_dim]];
	1762	1692 ProjData.VarAttribute{length(ProjData.ListVarName)}.Role=<span class="string">'errorflag'</span>;
	1763	1693 <span class="keyword">else</span> <span class="comment">%3D case</span>
	1764	1694 <span class="keyword">if</span> isequal(Theta,0) & isequal(Phi,0)
	1765	1695 test_sup=(Coord{1}>=ObjectData.Coord(1,3));
	1766	1696 iz_sup=find(test_sup);
	1767	1697 iz=iz_sup(1);
	1768	1698 <span class="keyword">if</span> iz>=1 & iz<=npz
	1769	1699 ProjData.ListDimName=[ProjData.ListDimName ListDimName(2:end)];
	1770	1700 ProjData.DimValue=[ProjData.DimValue npY npX];
	1771	1701 <span class="keyword">for</span> ivar=VarIndex
	1772	1702 VarName=FieldData.ListVarName{ivar};
	1773	1703 ProjData.ListVarName=[ProjData.ListVarName VarName];
	1774	1704 ProjData.VarDimIndex=[ProjData.VarDimIndex [nb_dim-2 nb_dim-1]];
	1775	1705 ProjData.VarAttribute{length(ProjData.ListVarName)}=FieldData.VarAttribute{ivar}; <span class="comment">%reproduce the variable attributes</span>
	1776	1706 <span class="keyword">if</span> test_direct_x==1
	1777	1707 ProjData.VarAttribute{length(ProjData.ListVarName)}.Coord_2=[XMin XMax];
	1778	1708 <span class="keyword">else</span>
	1779	1709 ProjData.VarAttribute{length(ProjData.ListVarName)}.Coord_2=[XMax XMin];
	1780	1710 <span class="keyword">end</span>
	1781	1711 <span class="keyword">if</span> test_direct_y==1
	1782	1712 ProjData.VarAttribute{length(ProjData.ListVarName)}.Coord_1=[YMin YMax];
	1783	1713 <span class="keyword">else</span>
	1784	1714 ProjData.VarAttribute{length(ProjData.ListVarName)}.Coord_1=[YMax YMin];
	1785	1715 <span class="keyword">end</span>
	1786	1716 eval([<span class="string">'ProjData.'</span> VarName <span class="string">'=squeeze(FieldData.'</span> VarName <span class="string">'(iz,:,:));'</span>])<span class="comment">% select the z index iz</span>
	1787	1717 <span class="comment">%TODO : do a vertical average for a thick plane</span>
	1788	1718 <span class="keyword">if</span> test_interp(2) \| test_interp(3)
	1789	1719 eval([<span class="string">'ProjData.'</span> VarName <span class="string">'=interp2(Coord{3},Coord{2},ProjData.'</span> VarName <span class="string">',Coord_x,Coord_y'');'</span>])
	1790	1720 <span class="keyword">end</span>
	1791	1721 <span class="keyword">end</span>
	1792	1722 <span class="keyword">end</span>
	1793	1723 <span class="keyword">else</span>
	1794	1724 <a href="warndlg_uvmat.html" class="code" title="function hwarn=warndlg_uvmat(warntext,title)">warndlg_uvmat</a>(<span class="string">'projection of structured coordinates on oblique plane not yet implemented'</span>,<span class="string">'ERROR'</span>)
	1795	1725 <span class="comment">%TODO: use interp3</span>
	1796	1726 <span class="keyword">return</span>
	1797	1727 <span class="keyword">end</span>
	1798	1728 <span class="keyword">end</span>
	1799	1729 <span class="keyword">end</span>
	1800	1730 <span class="keyword">end</span>
	1801	1731 <span class="comment">%projection of velocity components in the rotated coordinates</span>
	1802	1732 <span class="keyword">if</span> ~isequal(Phi,0) & length(ivar_U)==1
	1803	1733 <span class="keyword">if</span> isempty(ivar_V)
	1804	1734 <a href="warndlg_uvmat.html" class="code" title="function hwarn=warndlg_uvmat(warntext,title)">warndlg_uvmat</a>(<span class="string">'v velocity component missing in proj_field.m'</span>,<span class="string">'ERROR'</span>)
	1805	1735 <span class="keyword">return</span>
	1806	1736 <span class="keyword">end</span>
	1807	1737 UName=FieldData.ListVarName{ivar_U};
	1808	1738 VName=FieldData.ListVarName{ivar_V};
	1809	1739 eval([<span class="string">'ProjData.'</span> UName <span class="string">'=cos(Phi)ProjData.'</span> UName <span class="string">'+ sin(Phi)ProjData.'</span> VName <span class="string">';'</span>])
	1810	1740 eval([<span class="string">'ProjData.'</span> VName <span class="string">'=cos(Theta)(-sin(Phi)ProjData.'</span> UName <span class="string">'+ cos(Phi)*ProjData.'</span> VName <span class="string">');'</span>])
	1811	1741 <span class="keyword">if</span> ~isempty(ivar_W)
	1812	1742 WName=FieldData.ListVarName{ivar_W};
	1813	1743 eval([<span class="string">'ProjData.'</span> VName <span class="string">'=ProjData.'</span> VName <span class="string">'+ ProjData.'</span> WName <span class="string">'*sin(Theta);'</span>])<span class="comment">%</span>
	1814	1744 eval([<span class="string">'ProjData.'</span> WName <span class="string">'=NormVec_XProjData.'</span> UName <span class="string">'+ NormVec_YProjData.'</span> VName <span class="string">'+ NormVec_Z* ProjData.'</span> WName <span class="string">';'</span>]);
	1815	1745 <span class="keyword">end</span>
	1816	1746 <span class="keyword">if</span> ~isequal(Psi,0)
	1817	1747 eval([<span class="string">'ProjData.'</span> UName <span class="string">'=cos(Psi)* ProjData.'</span> UName <span class="string">'- sin(Psi)*ProjData.'</span> VName <span class="string">';'</span>]);
	1818	1748 eval([<span class="string">'ProjData.'</span> VName <span class="string">'=sin(Psi)* ProjData.'</span> UName <span class="string">'+ cos(Psi)*ProjData.'</span> VName <span class="string">';'</span>]);
	1819	1749 <span class="keyword">end</span>
	1820	1750 <span class="keyword">end</span>
	1821	1751 <span class="keyword">end</span>
	1822	1752
	1823	1753 <span class="comment">%-----------------------------------------------------------------</span>
	1824	1754 <span class="comment">%transmit the global attributes</span>
	1825	1755 <a name="_sub4" href="#_subfunctions" class="code">function [ProjData,errormsg]=proj_heading(FieldData,ObjectData)</a>
	1826	1756 <span class="comment">%-----------------------------------------------------------------</span>
	1827	1757 <span class="comment">% ProjData=FieldData;</span>
	1828	1758 ProjData=[];<span class="comment">%default</span>
	1829	1759 <span class="keyword">if</span> ~isfield(FieldData,<span class="string">'ListGlobalAttribute'</span>)
	1830	1760 ProjData.ListGlobalAttribute={};
	1831	1761 <span class="keyword">else</span>
	1832	1762 ProjData.ListGlobalAttribute=FieldData.ListGlobalAttribute;
	1833	1763 <span class="keyword">end</span>
	1834	1764 <span class="keyword">if</span> isfield(FieldData,<span class="string">'Txt'</span>)
	1835	1765 errormsg=FieldData.Txt; <span class="comment">%transmit erreur message</span>
	1836	1766 <span class="keyword">return</span>;
	1837	1767 <span class="keyword">end</span>
	1838	1768 <span class="keyword">for</span> iattr=1:length(ProjData.ListGlobalAttribute)
	1839	1769 AttrName=ProjData.ListGlobalAttribute{iattr};
	1840	1770 <span class="keyword">if</span> isfield(FieldData,AttrName)
	1841	1771 eval([<span class="string">'ProjData.'</span> AttrName <span class="string">'=FieldData.'</span> AttrName <span class="string">';'</span>]);
	1842	1772 <span class="keyword">end</span>
	1843	1773 <span class="keyword">end</span>
	1844	1774 <span class="keyword">if</span> isfield(FieldData,<span class="string">'CoordType'</span>)
	1845	1775 <span class="keyword">if</span> isfield(ObjectData,<span class="string">'CoordType'</span>)&~isequal(FieldData.CoordType,ObjectData.CoordType)
	1846	1776 errormsg=[ObjectData.Style <span class="string">' in '</span> ObjectData.CoordType <span class="string">' coordinates, while field in '</span> FieldData.CoordType <span class="string">' coordinates'</span>];
	1847	1777 <span class="keyword">return</span>
	1848	1778 <span class="keyword">else</span>
	1849	1779 ProjData.CoordType=FieldData.CoordType;
	1850	1780 <span class="keyword">end</span>
	1851	1781 <span class="keyword">end</span>
	1852	1782
	1853	1783 ListObject={<span class="string">'Style'</span>,<span class="string">'ProjMode'</span>,<span class="string">'RangeX'</span>,<span class="string">'RangeY'</span>,<span class="string">'RangeZ'</span>,<span class="string">'Phi'</span>,<span class="string">'Theta'</span>,<span class="string">'Psi'</span>,<span class="string">'Coord'</span>};
	1854	1784 <span class="keyword">for</span> ilist=1:length(ListObject)
	1855	1785 <span class="keyword">if</span> isfield(ObjectData,ListObject{ilist})
	1856	1786 eval([<span class="string">'val=ObjectData.'</span> ListObject{ilist} <span class="string">';'</span>])
	1857	1787 <span class="keyword">if</span> ~isempty(val)
	1858	1788 eval([<span class="string">'ProjData.Object'</span> ListObject{ilist} <span class="string">'=val;'</span>]);
	1859	1789 ProjData.ListGlobalAttribute=[ProjData.ListGlobalAttribute {[<span class="string">'Object'</span> ListObject{ilist}]}];
	1860	1790 <span class="keyword">end</span>
	1861	1791 <span class="keyword">end</span>
	1862	1792 <span class="keyword">end</span></pre></div>
	1863	<hr><address>Generated on Fri 13-Nov-2009 11:17:03 by <strong><a href="http://www.artefact.tk/software/matlab/m2html/">m2html</a></strong> © 2003</address>
	1864	</body>
	1865	</html>

Note: See TracBrowser for help on using the repository browser.

Download in other formats: