Context Navigation

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

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

Line
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: