[37] | 1 | <!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN" |
---|
| 2 | "http://www.w3.org/TR/REC-html40/loose.dtd"> |
---|
| 3 | <html> |
---|
| 4 | <head> |
---|
| 5 | <title>Description of plot_field</title> |
---|
| 6 | <meta name="keywords" content="plot_field"> |
---|
| 7 | <meta name="description" content="'plot_field': plot any field with the structure defined in the uvmat package"> |
---|
| 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> > plot_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>plot_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>'plot_field': plot any field with the structure defined in the uvmat package</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 [PlotType,PlotParamOut,haxes]= plot_field(Data,haxes,PlotParam,KeepLim,PosColorbar) </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">'plot_field': plot any field with the structure defined in the uvmat package |
---|
| 31 | ------------------------------------------------------------------------ |
---|
| 32 | |
---|
| 33 | This function is used by uvmat to plot fields. It automatically chooses the representation |
---|
| 34 | appropriate to the input field structure: |
---|
| 35 | 2D vector fields are represented by arrows, 2D scalar fiedlds by grey scale images or contour plots, 1D fields are represented by usual plot with (abscissa, ordinate). |
---|
| 36 | The input field structure is first tested by check_field_structure.m, |
---|
| 37 | then split into blocks of related variables by find_field_indices.m. |
---|
| 38 | The dimensionality of each block is obtained by this fuction |
---|
| 39 | considering the presence of variables with the attribute .Role='coord_x' |
---|
| 40 | and/or coord_y and/or coord_z (case of unstructured coordinates), or |
---|
| 41 | dimension variables (case of matrices). |
---|
| 42 | |
---|
| 43 | function [PlotType,PlotParamOut,haxes]= plot_field(Data,haxes,PlotParam,KeepLim,PosColorbar) |
---|
| 44 | |
---|
| 45 | OUPUT: |
---|
| 46 | PlotType: type of plot: 'text','line'(curve plot),'plane':2D view,'volume' |
---|
| 47 | PlotParamOut: structure, representing the updated plotting parameters, in case of automatic scaling |
---|
| 48 | haxes: handle of the plotting axis, when a new figure is created. |
---|
| 49 | |
---|
| 50 | INPUT |
---|
| 51 | Data: structure describing the field to plot |
---|
| 52 | (optional) .ListGlobalAttribute: cell listing the names of the global attributes |
---|
| 53 | .Att_1,Att_2... : values of the global attributes |
---|
| 54 | (requested) .ListVarName: list of variable names to select (cell array of char strings {'VarName1', 'VarName2',...} ) |
---|
| 55 | (requested) .VarDimName: list of dimension names for each element of .ListVarName (cell array of string cells) |
---|
| 56 | .VarAttribute: cell of attributes for each element of .ListVarName (cell array of structures of the form VarAtt.key=value) |
---|
| 57 | (requested) .Var1, .Var2....: variables (Matlab arrays) with names listed in .ListVarName |
---|
| 58 | |
---|
| 59 | Variable attribute .Role : |
---|
| 60 | The only variable attribute used for plotting purpose is .Role which can take |
---|
| 61 | the values |
---|
| 62 | Role = 'scalar': (default) represents a scalar field |
---|
| 63 | = 'coord_x', 'coord_y', 'coord_z': represents a separate set of |
---|
| 64 | unstructured coordinate x, y or z |
---|
| 65 | = 'vector': represents a vector field whose number of components |
---|
| 66 | is given by the last dimension (called 'nb_dim') |
---|
| 67 | = 'vector_x', 'vector_y', 'vector_z' :represents the x, y or z component of a vector |
---|
| 68 | = 'warnflag' : provides a warning flag about the quality of data in a 'Field', default=0, no warning |
---|
| 69 | = 'errorflag': provides an error flag marking false data, |
---|
| 70 | default=0, no error. Different non zero values can represent different criteria of elimination.</pre></div> |
---|
| 71 | |
---|
| 72 | <!-- crossreference --> |
---|
| 73 | <h2><a name="_cross"></a>CROSS-REFERENCE INFORMATION <a href="#_top"><img alt="^" border="0" src="../up.png"></a></h2> |
---|
| 74 | This function calls: |
---|
| 75 | <ul style="list-style-image:url(../matlabicon.gif)"> |
---|
| 76 | <li><a href="check_field_structure.html" class="code" title="function [DataOut,errormsg]=check_field_structure(Data)">check_field_structure</a> 'check_field_structure': check the representation of fields by a Matlab structure</li><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="msgbox_uvmat.html" class="code" title="function varargout = msgbox_uvmat(varargin)">msgbox_uvmat</a> 'msgbox_uvmat': associated with GUI msgbox_uvmat.fig to display message boxes, for error, warning or input calls</li><li><a href="plot_text.html" class="code" title="function hdisplay=plot_text(FieldData,hdisplay_in)">plot_text</a> 'plot_text': function for displaying the content of a Matlab structure in a figure</li><li><a href="proj_grid.html" class="code" title="function [A,rangx,rangy]=proj_grid(vec_X,vec_Y,vec_A,rgx_in,rgy_in,npxy_in)">proj_grid</a> 'proj_grid': project fields with unstructured coordinantes on a regular grid</li><li><a href="read_plot_param.html" class="code" title="function PlotParam=read_plot_param(handles)">read_plot_param</a> 'read_plot_param': read plotting parameters from the interface uvmat</li><li><a href="set_col_vec.html" class="code" title="function [colorlist,col_vec,colcode_out]=set_col_vec(colcode,vec_C)">set_col_vec</a> 'set_col_vec': sets the color code for vectors depending on a scalar vec_C and parameters given by the struct colcode</li><li><a href="set_title.html" class="code" title="function TITLE=set_title(Style,ProjMode)">set_title</a> 'set_title': defines the 'TITLE' of a projection object</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> |
---|
| 77 | This function is called by: |
---|
| 78 | <ul style="list-style-image:url(../matlabicon.gif)"> |
---|
| 79 | <li><a href="geometry_calib.html" class="code" title="function varargout = geometry_calib(varargin)">geometry_calib</a> 'geometry_calib': performs geometric calibration from a set of reference points</li><li><a href="get_field.html" class="code" title="function varargout = get_field(varargin)">get_field</a> 'get_field': display variables and attributes from a Netcdf file, and plot selected fields</li><li><a href="mouse_down.html" class="code" title="function xy=mouse_down(hObject,eventdata)">mouse_down</a> 'mouse_down': function activated when the mouse button is pressed on a figure (callback for 'WindowButtonDownFcn'</li><li><a href="probe_calib.html" class="code" title="function varargout = probe_calib(varargin)">probe_calib</a> 'probe_calib': performs geometric calibration from a set of reference points</li><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> |
---|
| 80 | <!-- crossreference --> |
---|
| 81 | |
---|
| 82 | <h2><a name="_subfunctions"></a>SUBFUNCTIONS <a href="#_top"><img alt="^" border="0" src="../up.png"></a></h2> |
---|
| 83 | <ul style="list-style-image:url(../matlabicon.gif)"> |
---|
| 84 | <li><a href="#_sub1" class="code">function [AxeData,haxes]=plot_hist(Data,haxes,PlotParam)</a></li><li><a href="#_sub2" class="code">function [AxeData,haxes]=plot_profile(data,CellVarIndex,VarType,haxes,PlotParam)</a></li><li><a href="#_sub3" class="code">function [AxeData,haxes,PlotParamOut,PlotType]=plot_plane(Data,CellVarIndex,VarTypeCell,haxes,PlotParam,KeepLim,PosColorbar)</a></li><li><a href="#_sub4" class="code">function quiresetn(haxes,x,y,u,v,scale,colorlist,col_vec)</a></li><li><a href="#_sub5" class="code">function YTick=colbartick(MinA,MaxA)</a></li></ul> |
---|
| 85 | <h2><a name="_source"></a>SOURCE CODE <a href="#_top"><img alt="^" border="0" src="../up.png"></a></h2> |
---|
| 86 | <div class="fragment"><pre>0001 <span class="comment">%'plot_field': plot any field with the structure defined in the uvmat package</span> |
---|
| 87 | 0002 <span class="comment">%------------------------------------------------------------------------</span> |
---|
| 88 | 0003 <span class="comment">%</span> |
---|
| 89 | 0004 <span class="comment">% This function is used by uvmat to plot fields. It automatically chooses the representation</span> |
---|
| 90 | 0005 <span class="comment">% appropriate to the input field structure:</span> |
---|
| 91 | 0006 <span class="comment">% 2D vector fields are represented by arrows, 2D scalar fiedlds by grey scale images or contour plots, 1D fields are represented by usual plot with (abscissa, ordinate).</span> |
---|
| 92 | 0007 <span class="comment">% The input field structure is first tested by check_field_structure.m,</span> |
---|
| 93 | 0008 <span class="comment">% then split into blocks of related variables by find_field_indices.m.</span> |
---|
| 94 | 0009 <span class="comment">% The dimensionality of each block is obtained by this fuction</span> |
---|
| 95 | 0010 <span class="comment">% considering the presence of variables with the attribute .Role='coord_x'</span> |
---|
| 96 | 0011 <span class="comment">% and/or coord_y and/or coord_z (case of unstructured coordinates), or</span> |
---|
| 97 | 0012 <span class="comment">% dimension variables (case of matrices).</span> |
---|
| 98 | 0013 <span class="comment">%</span> |
---|
| 99 | 0014 <span class="comment">% function [PlotType,PlotParamOut,haxes]= plot_field(Data,haxes,PlotParam,KeepLim,PosColorbar)</span> |
---|
| 100 | 0015 <span class="comment">%</span> |
---|
| 101 | 0016 <span class="comment">% OUPUT:</span> |
---|
| 102 | 0017 <span class="comment">% PlotType: type of plot: 'text','line'(curve plot),'plane':2D view,'volume'</span> |
---|
| 103 | 0018 <span class="comment">% PlotParamOut: structure, representing the updated plotting parameters, in case of automatic scaling</span> |
---|
| 104 | 0019 <span class="comment">% haxes: handle of the plotting axis, when a new figure is created.</span> |
---|
| 105 | 0020 <span class="comment">%</span> |
---|
| 106 | 0021 <span class="comment">%INPUT</span> |
---|
| 107 | 0022 <span class="comment">% Data: structure describing the field to plot</span> |
---|
| 108 | 0023 <span class="comment">% (optional) .ListGlobalAttribute: cell listing the names of the global attributes</span> |
---|
| 109 | 0024 <span class="comment">% .Att_1,Att_2... : values of the global attributes</span> |
---|
| 110 | 0025 <span class="comment">% (requested) .ListVarName: list of variable names to select (cell array of char strings {'VarName1', 'VarName2',...} )</span> |
---|
| 111 | 0026 <span class="comment">% (requested) .VarDimName: list of dimension names for each element of .ListVarName (cell array of string cells)</span> |
---|
| 112 | 0027 <span class="comment">% .VarAttribute: cell of attributes for each element of .ListVarName (cell array of structures of the form VarAtt.key=value)</span> |
---|
| 113 | 0028 <span class="comment">% (requested) .Var1, .Var2....: variables (Matlab arrays) with names listed in .ListVarName</span> |
---|
| 114 | 0029 <span class="comment">%</span> |
---|
| 115 | 0030 <span class="comment">% Variable attribute .Role :</span> |
---|
| 116 | 0031 <span class="comment">% The only variable attribute used for plotting purpose is .Role which can take</span> |
---|
| 117 | 0032 <span class="comment">% the values</span> |
---|
| 118 | 0033 <span class="comment">% Role = 'scalar': (default) represents a scalar field</span> |
---|
| 119 | 0034 <span class="comment">% = 'coord_x', 'coord_y', 'coord_z': represents a separate set of</span> |
---|
| 120 | 0035 <span class="comment">% unstructured coordinate x, y or z</span> |
---|
| 121 | 0036 <span class="comment">% = 'vector': represents a vector field whose number of components</span> |
---|
| 122 | 0037 <span class="comment">% is given by the last dimension (called 'nb_dim')</span> |
---|
| 123 | 0038 <span class="comment">% = 'vector_x', 'vector_y', 'vector_z' :represents the x, y or z component of a vector</span> |
---|
| 124 | 0039 <span class="comment">% = 'warnflag' : provides a warning flag about the quality of data in a 'Field', default=0, no warning</span> |
---|
| 125 | 0040 <span class="comment">% = 'errorflag': provides an error flag marking false data,</span> |
---|
| 126 | 0041 <span class="comment">% default=0, no error. Different non zero values can represent different criteria of elimination.</span> |
---|
| 127 | 0042 |
---|
| 128 | 0043 |
---|
| 129 | 0044 <span class="comment">% additional elements characterizing the projection object (should not be necessary)--</span> |
---|
| 130 | 0045 <span class="comment">% Data.Style : style of projection object</span> |
---|
| 131 | 0046 <span class="comment">% Data.XObject,.YObject: set of coordinates defining the object position;</span> |
---|
| 132 | 0047 <span class="comment">% Data.ProjMode=type of projection ;</span> |
---|
| 133 | 0048 <span class="comment">% Data.ProjAngle=angle of projection;</span> |
---|
| 134 | 0049 <span class="comment">% Data.DX,.DY,.DZ=increments;</span> |
---|
| 135 | 0050 <span class="comment">% Data.MaxY,MinY: min and max Y</span> |
---|
| 136 | 0051 |
---|
| 137 | 0052 <span class="comment">% haxes: handle of the plotting axes to update with the new plot. If this input is absent or not a valid axes handle, a new figure is created.</span> |
---|
| 138 | 0053 <span class="comment">%</span> |
---|
| 139 | 0054 <span class="comment">% PlotParam: parameters for plotting, as read on the uvmat interface (by function 'read_plot_param.m')</span> |
---|
| 140 | 0055 <span class="comment">% --scalars--</span> |
---|
| 141 | 0056 <span class="comment">% .Scalar.MaxA: upper bound (saturation color) for the scalar representation, max(field) by default</span> |
---|
| 142 | 0057 <span class="comment">% .Scalar.MinA: lower bound (saturation) for the scalar representation, min(field) by default</span> |
---|
| 143 | 0058 <span class="comment">% .Scalar.AutoScal: =1 (default) lower and upper bounds of the scalar representation set to the min and max of the field</span> |
---|
| 144 | 0059 <span class="comment">% =0 lower and upper bound imposed by .AMax and .MinA</span> |
---|
| 145 | 0060 <span class="comment">% .Scalar.BW= 1 black and white representation imposed, =0 by default.</span> |
---|
| 146 | 0061 <span class="comment">% .Scalar.Contours= 1: represent scalars by contour plots (Matlab function 'contour'); =0 by default</span> |
---|
| 147 | 0062 <span class="comment">% .IncrA : contour interval</span> |
---|
| 148 | 0063 <span class="comment">% -- vectors--</span> |
---|
| 149 | 0064 <span class="comment">% .Vectors.VecScale: scale for the vector representation</span> |
---|
| 150 | 0065 <span class="comment">% .Vectors.AutoVec: =0 (default) automatic length for vector representation, =1: length set by .VecScale</span> |
---|
| 151 | 0066 <span class="comment">% .Vectors.HideFalse= 0 (default) false vectors represented in magenta, =1: false vectors not represented;</span> |
---|
| 152 | 0067 <span class="comment">% .Vectors.HideWarning= 0 (default) vectors marked by warnflag~=0 marked in black, 1: no warning representation;</span> |
---|
| 153 | 0068 <span class="comment">% .Vectors.decimate4 = 0 (default) all vectors reprtesented, =1: half of the vectors represented along each coordinate</span> |
---|
| 154 | 0069 <span class="comment">% -- vector color--</span> |
---|
| 155 | 0070 <span class="comment">% .Vectors.ColorCode= 'black','white': imposed color (default ='blue')</span> |
---|
| 156 | 0071 <span class="comment">% 'rgb', : three colors red, blue, green depending</span> |
---|
| 157 | 0072 <span class="comment">% on thresholds .colcode1 and .colcode2 on the input scalar value (C)</span> |
---|
| 158 | 0073 <span class="comment">% 'brg': like rgb but reversed color order (blue, green, red)</span> |
---|
| 159 | 0074 <span class="comment">% '64 colors': continuous color from blue to red (multijet)</span> |
---|
| 160 | 0075 <span class="comment">% .Vectors.colcode1 : first threshold for rgb, first value for'continuous'</span> |
---|
| 161 | 0076 <span class="comment">% .Vectors.colcode2 : second threshold for rgb, last value (saturation) for 'continuous'</span> |
---|
| 162 | 0077 <span class="comment">% .Vectors.FixedCbounds; =0 (default): the bounds on C representation are min and max, =1: they are fixed by .Minc and .MaxC</span> |
---|
| 163 | 0078 <span class="comment">% .Vectors.MinC = imposed minimum of the scalar field used for vector color;</span> |
---|
| 164 | 0079 <span class="comment">% .Vectors.MaxC = imposed maximum of the scalar field used for vector color;</span> |
---|
| 165 | 0080 <span class="comment">%</span> |
---|
| 166 | 0081 <span class="comment">% KeepLim:=0 (default) adjust axes limit to the X,Y data, =1: preserves the previous axes limits</span> |
---|
| 167 | 0082 <span class="comment">% PosColorbar: if not empty, display a colorbar for B&W images</span> |
---|
| 168 | 0083 <span class="comment">% imposed position of the colorbar (ex [0.821 0.471 0.019 0.445])</span> |
---|
| 169 | 0084 <span class="comment">%</span> |
---|
| 170 | 0085 <span class="comment">%AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA</span> |
---|
| 171 | 0086 <span class="comment">% Copyright Joel Sommeria, 2008, LEGI / CNRS-UJF-INPG, sommeria@coriolis-legi.org.</span> |
---|
| 172 | 0087 <span class="comment">%AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA</span> |
---|
| 173 | 0088 <span class="comment">% This file is part of the toolbox UVMAT.</span> |
---|
| 174 | 0089 <span class="comment">%</span> |
---|
| 175 | 0090 <span class="comment">% UVMAT is free software; you can redistribute it and/or modify</span> |
---|
| 176 | 0091 <span class="comment">% it under the terms of the GNU General Public License as published by</span> |
---|
| 177 | 0092 <span class="comment">% the Free Software Foundation; either version 2 of the License, or</span> |
---|
| 178 | 0093 <span class="comment">% (at your option) any later version.</span> |
---|
| 179 | 0094 <span class="comment">%</span> |
---|
| 180 | 0095 <span class="comment">% UVMAT is distributed in the hope that it will be useful,</span> |
---|
| 181 | 0096 <span class="comment">% but WITHOUT ANY WARRANTY; without even the implied warranty of</span> |
---|
| 182 | 0097 <span class="comment">% MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the</span> |
---|
| 183 | 0098 <span class="comment">% GNU General Public License (file UVMAT/COPYING.txt) for more details.</span> |
---|
| 184 | 0099 <span class="comment">%AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA</span> |
---|
| 185 | 0100 |
---|
| 186 | 0101 <a name="_sub0" href="#_subfunctions" class="code">function [PlotType,PlotParamOut,haxes]= plot_field(Data,haxes,PlotParam,KeepLim,PosColorbar)</a> |
---|
| 187 | 0102 <span class="comment">%default output</span> |
---|
| 188 | 0103 <span class="keyword">if</span> ~exist(<span class="string">'PlotParam'</span>,<span class="string">'var'</span>),PlotParam=[];<span class="keyword">end</span>; |
---|
| 189 | 0104 <span class="keyword">if</span> ~exist(<span class="string">'KeepLim'</span>,<span class="string">'var'</span>),KeepLim=0;<span class="keyword">end</span>; |
---|
| 190 | 0105 <span class="keyword">if</span> ~exist(<span class="string">'PosColorbar'</span>,<span class="string">'var'</span>),PosColorbar=[];<span class="keyword">end</span>; |
---|
| 191 | 0106 PlotType=<span class="string">'text'</span>; <span class="comment">%default</span> |
---|
| 192 | 0107 PlotParamOut=PlotParam;<span class="comment">%default</span> |
---|
| 193 | 0108 |
---|
| 194 | 0109 <span class="comment">% check input structure</span> |
---|
| 195 | 0110 [Data,errormsg]=<a href="check_field_structure.html" class="code" title="function [DataOut,errormsg]=check_field_structure(Data)">check_field_structure</a>(Data); |
---|
| 196 | 0111 |
---|
| 197 | 0112 <span class="keyword">if</span> ~isempty(errormsg) |
---|
| 198 | 0113 <a href="msgbox_uvmat.html" class="code" title="function varargout = msgbox_uvmat(varargin)">msgbox_uvmat</a>(<span class="string">'ERROR'</span>,[<span class="string">'input of plot_field/check_field_structure: '</span> errormsg]) |
---|
| 199 | 0114 display([<span class="string">'input of plot_field/check_field_structure:: '</span> errormsg]) |
---|
| 200 | 0115 <span class="keyword">return</span> |
---|
| 201 | 0116 <span class="keyword">end</span> |
---|
| 202 | 0117 |
---|
| 203 | 0118 testnewfig=1;<span class="comment">%test to create a new figure (default)</span> |
---|
| 204 | 0119 testzoomaxes=0;<span class="comment">%test for the existence of a zoom secondary figure attached to the plotting axes</span> |
---|
| 205 | 0120 <span class="keyword">if</span> exist(<span class="string">'haxes'</span>,<span class="string">'var'</span>) |
---|
| 206 | 0121 <span class="keyword">if</span> ishandle(haxes) |
---|
| 207 | 0122 <span class="keyword">if</span> isequal(get(haxes,<span class="string">'Type'</span>),<span class="string">'axes'</span>) |
---|
| 208 | 0123 <span class="comment">% hfig=get(haxes,'Parent');</span> |
---|
| 209 | 0124 axes(haxes) |
---|
| 210 | 0125 testnewfig=0; |
---|
| 211 | 0126 AxeData=get(haxes,<span class="string">'UserData'</span>); |
---|
| 212 | 0127 <span class="keyword">if</span> isfield(AxeData,<span class="string">'ZoomAxes'</span>)&& ishandle(AxeData.ZoomAxes) |
---|
| 213 | 0128 <span class="keyword">if</span> isequal(get(AxeData.ZoomAxes,<span class="string">'Type'</span>),<span class="string">'axes'</span>) |
---|
| 214 | 0129 testzoomaxes=1; |
---|
| 215 | 0130 zoomaxes=AxeData.ZoomAxes; |
---|
| 216 | 0131 <span class="keyword">end</span> |
---|
| 217 | 0132 <span class="keyword">end</span> |
---|
| 218 | 0133 <span class="keyword">end</span> |
---|
| 219 | 0134 <span class="keyword">end</span> |
---|
| 220 | 0135 <span class="keyword">end</span> |
---|
| 221 | 0136 <span class="keyword">if</span> testnewfig<span class="comment">% create a new figure and axes if the plotting axes does not exist</span> |
---|
| 222 | 0137 hfig=figure; |
---|
| 223 | 0138 <span class="keyword">if</span> isfield(Data,<span class="string">'IndexObj'</span>)&isfield(Data,<span class="string">'Style'</span>)&isfield(Data,<span class="string">'ProjMode'</span>) |
---|
| 224 | 0139 figname=[num2str(Data.IndexObj) <span class="string">'-'</span> <a href="set_title.html" class="code" title="function TITLE=set_title(Style,ProjMode)">set_title</a>(Data.Style,Data.ProjMode)]; |
---|
| 225 | 0140 set(hfig,<span class="string">'Name'</span>,figname) |
---|
| 226 | 0141 <span class="keyword">end</span> |
---|
| 227 | 0142 testhandle=0; |
---|
| 228 | 0143 <span class="keyword">if</span> isfield(PlotParam,<span class="string">'text_display_1'</span>)& ishandle(PlotParam.text_display_1) |
---|
| 229 | 0144 set(hfig,<span class="string">'UserData'</span>,PlotParam) |
---|
| 230 | 0145 testhandle=1; |
---|
| 231 | 0146 <span class="keyword">end</span> |
---|
| 232 | 0147 set(hfig,<span class="string">'Units'</span>,<span class="string">'normalized'</span>) |
---|
| 233 | 0148 set(hfig,<span class="string">'WindowButtonDownFcn'</span>,<span class="string">'mouse_down'</span>) |
---|
| 234 | 0149 <span class="comment">%set(hfig,'WindowButtonMotionFcn',{'mouse_motion',PlotParam})%set mouse action function</span> |
---|
| 235 | 0150 set(hfig,<span class="string">'WindowButtonMotionFcn'</span>,<span class="string">'mouse_motion'</span>)<span class="comment">%set mouse action function</span> |
---|
| 236 | 0151 set(hfig,<span class="string">'WindowButtonUpFcn'</span>,<span class="string">'mouse_up'</span>)<span class="comment">%set mouse action function</span> |
---|
| 237 | 0152 haxes=axes; |
---|
| 238 | 0153 set(haxes,<span class="string">'position'</span>,[0.13,0.2,0.775,0.73]) |
---|
| 239 | 0154 <span class="keyword">end</span> |
---|
| 240 | 0155 <span class="keyword">if</span> isfield(PlotParam,<span class="string">'text_display_1'</span>)& ishandle(PlotParam.text_display_1) |
---|
| 241 | 0156 PlotParam=<a href="read_plot_param.html" class="code" title="function PlotParam=read_plot_param(handles)">read_plot_param</a>(PlotParam); |
---|
| 242 | 0157 <span class="keyword">end</span> |
---|
| 243 | 0158 <span class="keyword">if</span> testnewfig |
---|
| 244 | 0159 PlotParam.NextPlot=<span class="string">'add'</span>; <span class="comment">%parameter for plot_profile and plot_hist</span> |
---|
| 245 | 0160 <span class="keyword">end</span> |
---|
| 246 | 0161 <span class="keyword">if</span> isfield(PlotParam,<span class="string">'Auto_xy'</span>) && isequal(PlotParam.Auto_xy,1) |
---|
| 247 | 0162 set(haxes,<span class="string">'DataAspectRatioMode'</span>,<span class="string">'auto'</span>)<span class="comment">%automatic aspect ratio</span> |
---|
| 248 | 0163 <span class="keyword">end</span> |
---|
| 249 | 0164 |
---|
| 250 | 0165 <span class="comment">% check the cells of fields :</span> |
---|
| 251 | 0166 testnbdim=1; |
---|
| 252 | 0167 [CellVarIndex,NbDim,VarType,errormsg]=<a href="find_field_indices.html" class="code" title="function [CellVarIndex,NbDim,VarType,errormsg]=find_field_indices(Data)">find_field_indices</a>(Data); |
---|
| 253 | 0168 |
---|
| 254 | 0169 <span class="keyword">if</span> ~isempty(errormsg) |
---|
| 255 | 0170 <a href="msgbox_uvmat.html" class="code" title="function varargout = msgbox_uvmat(varargin)">msgbox_uvmat</a>(<span class="string">'ERROR'</span>,[<span class="string">'input of plot_field/find_field_indices: '</span> errormsg]) |
---|
| 256 | 0171 display([<span class="string">'input of plot_field: '</span> errormsg]) |
---|
| 257 | 0172 <span class="keyword">return</span> |
---|
| 258 | 0173 <span class="keyword">end</span> |
---|
| 259 | 0174 <span class="keyword">if</span> ~isfield(Data,<span class="string">'NbDim'</span>) <span class="comment">%& ~isfield(Data,'Style')%determine the space dimensionb if not defined: choose the kind of plot</span> |
---|
| 260 | 0175 [Data.NbDim,imax]=max(NbDim); |
---|
| 261 | 0176 <span class="keyword">end</span> |
---|
| 262 | 0177 <span class="keyword">if</span> isequal(Data.NbDim,0) <span class="comment">% TODO: chech whether this function is still used, replace by plot_profile ?</span> |
---|
| 263 | 0178 <span class="keyword">if</span> isfield(Data,<span class="string">'Style'</span>) & isequal(Data.Style,<span class="string">'points'</span>) |
---|
| 264 | 0179 AxeData=<a href="plot_text.html" class="code" title="function hdisplay=plot_text(FieldData,hdisplay_in)">plot_text</a>(Data,haxes); |
---|
| 265 | 0180 PlotType=<span class="string">'text'</span>; |
---|
| 266 | 0181 <span class="keyword">else</span> |
---|
| 267 | 0182 [AxeData,haxes]=<a href="#_sub1" class="code" title="subfunction [AxeData,haxes]=plot_hist(Data,haxes,PlotParam)">plot_hist</a>(Data,haxes,PlotParam); |
---|
| 268 | 0183 <span class="keyword">end</span> |
---|
| 269 | 0184 <span class="keyword">elseif</span> isequal(Data.NbDim,1) |
---|
| 270 | 0185 [AxeData,haxes]=<a href="#_sub2" class="code" title="subfunction [AxeData,haxes]=plot_profile(data,CellVarIndex,VarType,haxes,PlotParam)">plot_profile</a>(Data,CellVarIndex,VarType,haxes,PlotParam);<span class="comment">%</span> |
---|
| 271 | 0186 <span class="keyword">if</span> testzoomaxes |
---|
| 272 | 0187 [AxeData,zoomaxes,PlotParamOut]=<a href="#_sub2" class="code" title="subfunction [AxeData,haxes]=plot_profile(data,CellVarIndex,VarType,haxes,PlotParam)">plot_profile</a>(Data,CellVarIndex,VarType,zoomaxes,PlotParam); |
---|
| 273 | 0188 AxeData.ZoomAxes=zoomaxes; |
---|
| 274 | 0189 <span class="keyword">end</span> |
---|
| 275 | 0190 PlotType=<span class="string">'line'</span>; |
---|
| 276 | 0191 <span class="keyword">elseif</span> isequal(Data.NbDim,2) |
---|
| 277 | 0192 ind_select=find(NbDim>=2); |
---|
| 278 | 0193 <span class="keyword">if</span> numel(ind_select)>2 |
---|
| 279 | 0194 <a href="msgbox_uvmat.html" class="code" title="function varargout = msgbox_uvmat(varargin)">msgbox_uvmat</a>(<span class="string">'ERROR'</span>,[<span class="string">'more than two fields to map'</span>]) |
---|
| 280 | 0195 display([<span class="string">'more than two fields to map'</span>]) |
---|
| 281 | 0196 <span class="keyword">return</span> |
---|
| 282 | 0197 <span class="keyword">end</span> |
---|
| 283 | 0198 [AxeData,haxes,PlotParamOut,PlotType]=<a href="#_sub3" class="code" title="subfunction [AxeData,haxes,PlotParamOut,PlotType]=plot_plane(Data,CellVarIndex,VarTypeCell,haxes,PlotParam,KeepLim,PosColorbar)">plot_plane</a>(Data,CellVarIndex(ind_select),VarType(ind_select),haxes,PlotParam,KeepLim,PosColorbar); |
---|
| 284 | 0199 <span class="keyword">if</span> testzoomaxes |
---|
| 285 | 0200 [AxeData,zoomaxes,PlotParamOut]=<a href="#_sub3" class="code" title="subfunction [AxeData,haxes,PlotParamOut,PlotType]=plot_plane(Data,CellVarIndex,VarTypeCell,haxes,PlotParam,KeepLim,PosColorbar)">plot_plane</a>(Data,CellVarIndex(ind_select),VarType(ind_select),zoomaxes,PlotParam,1,PosColorbar); |
---|
| 286 | 0201 AxeData.ZoomAxes=zoomaxes; |
---|
| 287 | 0202 <span class="keyword">end</span> |
---|
| 288 | 0203 <span class="keyword">elseif</span> isequal(Data.NbDim,3) |
---|
| 289 | 0204 <a href="msgbox_uvmat.html" class="code" title="function varargout = msgbox_uvmat(varargin)">msgbox_uvmat</a>(<span class="string">'ERROR'</span>,<span class="string">'volume plot not implemented yet'</span>) |
---|
| 290 | 0205 <span class="keyword">return</span> |
---|
| 291 | 0206 <span class="comment">%plot_volume(haxes,Data,PlotParam)% A FAIRE</span> |
---|
| 292 | 0207 <span class="comment">%PlotType='volume';</span> |
---|
| 293 | 0208 <span class="keyword">else</span> |
---|
| 294 | 0209 testnbdim=0; |
---|
| 295 | 0210 <span class="keyword">end</span> |
---|
| 296 | 0211 |
---|
| 297 | 0212 <span class="comment">%display (or delete) error message</span> |
---|
| 298 | 0213 htext=findobj(haxes,<span class="string">'Tag'</span>,<span class="string">'hTxt'</span>); |
---|
| 299 | 0214 <span class="keyword">if</span> isfield(Data,<span class="string">'Txt'</span>) |
---|
| 300 | 0215 <span class="keyword">if</span> isempty(htext) |
---|
| 301 | 0216 Xlim=get(haxes,<span class="string">'XLim'</span>); |
---|
| 302 | 0217 Ylim=get(haxes,<span class="string">'YLim'</span>); |
---|
| 303 | 0218 htext=text(Xlim(1),(Ylim(1)+Ylim(2))/2,Data.Txt,<span class="string">'Tag'</span>,<span class="string">'hTxt'</span>,<span class="string">'Color'</span>,<span class="string">'r'</span>); |
---|
| 304 | 0219 set(htext,<span class="string">'Interpreter'</span>,<span class="string">'none'</span>) |
---|
| 305 | 0220 <span class="keyword">else</span> |
---|
| 306 | 0221 set(htext,<span class="string">'String'</span>,Data.Txt) |
---|
| 307 | 0222 <span class="keyword">end</span> |
---|
| 308 | 0223 <span class="keyword">elseif</span> ~isempty(htext) |
---|
| 309 | 0224 delete(htext) |
---|
| 310 | 0225 <span class="keyword">end</span> |
---|
| 311 | 0226 |
---|
| 312 | 0227 <span class="comment">% set graph aspect ratio</span> |
---|
| 313 | 0228 |
---|
| 314 | 0229 |
---|
| 315 | 0230 set(haxes,<span class="string">'UserData'</span>,AxeData) |
---|
| 316 | 0231 <span class="comment">%set(haxes,'Tag','uvmat');</span> |
---|
| 317 | 0232 |
---|
| 318 | 0233 |
---|
| 319 | 0234 <span class="comment">%-------------------------------------------</span> |
---|
| 320 | 0235 <a name="_sub1" href="#_subfunctions" class="code">function [AxeData,haxes]=plot_hist(Data,haxes,PlotParam)</a><span class="comment">% TODO: chech whether this function is still used, replace by plot_profile ?</span> |
---|
| 321 | 0236 <span class="comment">%------------------------------------------</span> |
---|
| 322 | 0237 AxeData=get(haxes,<span class="string">'UserData'</span>); <span class="comment">%defau</span> |
---|
| 323 | 0238 hfig=get(haxes,<span class="string">'parent'</span>); |
---|
| 324 | 0239 <span class="keyword">if</span> ~isfield(Data,<span class="string">'ListVarName'</span>) |
---|
| 325 | 0240 <span class="keyword">return</span> |
---|
| 326 | 0241 <span class="keyword">end</span> |
---|
| 327 | 0242 ColorOrder=[1 0 0;0 0.5 0;0 0 1;0 0.75 0.75;0.75 0 0.75;0.75 0.75 0;0.25 0.25 0.25]; |
---|
| 328 | 0243 set(haxes,<span class="string">'ColorOrder'</span>,ColorOrder) |
---|
| 329 | 0244 <span class="keyword">if</span> isfield(PlotParam,<span class="string">'NextPlot'</span>) |
---|
| 330 | 0245 set(haxes,<span class="string">'NextPlot'</span>,PlotParam.NextPlot) |
---|
| 331 | 0246 <span class="keyword">end</span> |
---|
| 332 | 0247 charplot=<span class="string">'''-'''</span>; |
---|
| 333 | 0248 iplot=0; |
---|
| 334 | 0249 legend_str={}; |
---|
| 335 | 0250 label_str=<span class="string">''</span>; |
---|
| 336 | 0251 textmean={}; |
---|
| 337 | 0252 plotstr=<span class="string">'plot('</span>; |
---|
| 338 | 0253 <span class="keyword">for</span> ilist=1:length(Data.ListVarName) |
---|
| 339 | 0254 VarName=Data.ListVarName{ilist}; |
---|
| 340 | 0255 eval([<span class="string">'['</span> VarName <span class="string">'hist,'</span> VarName <span class="string">'val]=hist(double(Data.'</span> VarName <span class="string">'),100);'</span>]);<span class="comment">%coordinate variable set as c</span> |
---|
| 341 | 0256 plotstr=[plotstr VarName <span class="string">'val,'</span> VarName <span class="string">'hist,'</span> charplot <span class="string">','</span>]; |
---|
| 342 | 0257 eval([<span class="string">'nbcomponent2=size(Data.'</span> VarName <span class="string">',2);'</span>]); |
---|
| 343 | 0258 eval([<span class="string">'nbcomponent1=size(Data.'</span> VarName <span class="string">',1);'</span>]); |
---|
| 344 | 0259 eval([<span class="string">'varmean=mean(double(Data.'</span> VarName <span class="string">'));'</span>]);<span class="comment">%mean value</span> |
---|
| 345 | 0260 textmean=[textmean; {[VarName <span class="string">'mean= '</span> num2str(varmean,4)]}]; |
---|
| 346 | 0261 <span class="keyword">if</span> nbcomponent1==1| nbcomponent2==1 |
---|
| 347 | 0262 legend_str=[legend_str {VarName}]; <span class="comment">%variable with one component</span> |
---|
| 348 | 0263 <span class="keyword">else</span> |
---|
| 349 | 0264 <span class="keyword">for</span> ic=1:min(nbcomponent1,nbcomponent2) |
---|
| 350 | 0265 legend_str=[legend_str [VarName <span class="string">'_'</span> num2str(ic)]]; <span class="comment">%variable with severals components</span> |
---|
| 351 | 0266 <span class="comment">% labeled by their index (e.g. color component)</span> |
---|
| 352 | 0267 <span class="keyword">end</span> |
---|
| 353 | 0268 <span class="keyword">end</span> |
---|
| 354 | 0269 label_str=[label_str <span class="string">' '</span> VarName]; |
---|
| 355 | 0270 <span class="keyword">end</span> |
---|
| 356 | 0271 <span class="keyword">if</span> ~isequal(plotstr,<span class="string">'plot('</span>) |
---|
| 357 | 0272 plotstr(end)=<span class="string">')'</span>; |
---|
| 358 | 0273 eval(plotstr) |
---|
| 359 | 0274 hlegend=findobj(hfig,<span class="string">'Tag'</span>,<span class="string">'legend'</span>);<span class="comment">%find existing legend on the plot</span> |
---|
| 360 | 0275 <span class="keyword">if</span> ~isempty(hlegend) |
---|
| 361 | 0276 legend_old=get(hlegend,<span class="string">'String'</span>); |
---|
| 362 | 0277 <span class="keyword">if</span> isequal(size(legend_old,1),size(legend_str,1)) |
---|
| 363 | 0278 legend_str=[legend_old legend_str]; |
---|
| 364 | 0279 <span class="keyword">end</span> |
---|
| 365 | 0280 <span class="keyword">end</span> |
---|
| 366 | 0281 legend(legend_str) |
---|
| 367 | 0282 xlabel(label_str) |
---|
| 368 | 0283 ylabel(<span class="string">'nb values'</span>) |
---|
| 369 | 0284 grid on |
---|
| 370 | 0285 title_str=<span class="string">''</span>; |
---|
| 371 | 0286 <span class="keyword">if</span> isfield(Data,<span class="string">'filename'</span>) |
---|
| 372 | 0287 [Path, title_str, ext]=fileparts(Data.filename); |
---|
| 373 | 0288 title_str=[title_str ext]; |
---|
| 374 | 0289 <span class="keyword">end</span> |
---|
| 375 | 0290 <span class="keyword">if</span> isfield(Data,<span class="string">'Action'</span>) |
---|
| 376 | 0291 <span class="keyword">if</span> ~isequal(title_str,<span class="string">''</span>) |
---|
| 377 | 0292 title_str=[title_str <span class="string">', '</span>]; |
---|
| 378 | 0293 <span class="keyword">end</span> |
---|
| 379 | 0294 title_str=[title_str Data.Action]; |
---|
| 380 | 0295 <span class="keyword">end</span> |
---|
| 381 | 0296 htitle=title(title_str); |
---|
| 382 | 0297 set(htitle,<span class="string">'Interpreter'</span>,<span class="string">'none'</span>)<span class="comment">% desable tex interpreter</span> |
---|
| 383 | 0298 hlist=findobj(gcf,<span class="string">'Style'</span>,<span class="string">'listbox'</span>); |
---|
| 384 | 0299 <span class="keyword">if</span> isempty(hlist) |
---|
| 385 | 0300 uicontrol(<span class="string">'Style'</span>,<span class="string">'popupmenu'</span>,<span class="string">'Position'</span>,[20 20 200 20],<span class="string">'String'</span>,textmean); |
---|
| 386 | 0301 <span class="keyword">else</span> |
---|
| 387 | 0302 set(hlist(1),<span class="string">'String'</span>,textmean) |
---|
| 388 | 0303 <span class="keyword">end</span> |
---|
| 389 | 0304 <span class="keyword">end</span> |
---|
| 390 | 0305 AxeData=Data; |
---|
| 391 | 0306 |
---|
| 392 | 0307 |
---|
| 393 | 0308 <span class="comment">%----------------------------------------------------------</span> |
---|
| 394 | 0309 <a name="_sub2" href="#_subfunctions" class="code">function [AxeData,haxes]=plot_profile(data,CellVarIndex,VarType,haxes,PlotParam)</a> |
---|
| 395 | 0310 <span class="comment">%-----------------------------------------------------------</span> |
---|
| 396 | 0311 axes(haxes) |
---|
| 397 | 0312 hfig=get(haxes,<span class="string">'parent'</span>); |
---|
| 398 | 0313 AxeData=data; |
---|
| 399 | 0314 ColorOrder=[1 0 0;0 0.5 0;0 0 1;0 0.75 0.75;0.75 0 0.75;0.75 0.75 0;0.25 0.25 0.25]; |
---|
| 400 | 0315 set(haxes,<span class="string">'ColorOrder'</span>,ColorOrder) |
---|
| 401 | 0316 <span class="keyword">if</span> isfield(PlotParam,<span class="string">'NextPlot'</span>) |
---|
| 402 | 0317 set(haxes,<span class="string">'NextPlot'</span>,PlotParam.NextPlot) |
---|
| 403 | 0318 <span class="keyword">end</span> |
---|
| 404 | 0319 legend_str={}; |
---|
| 405 | 0320 |
---|
| 406 | 0321 <span class="comment">%initiates string of the plot command</span> |
---|
| 407 | 0322 plotstr=<span class="string">'plot('</span>; |
---|
| 408 | 0323 textmean={}; |
---|
| 409 | 0324 abscissa_name=<span class="string">''</span>; |
---|
| 410 | 0325 coord_x_index=[]; |
---|
| 411 | 0326 <span class="keyword">for</span> icell=1:length(CellVarIndex) |
---|
| 412 | 0327 testfalse=0; |
---|
| 413 | 0328 VarIndex=CellVarIndex{icell};<span class="comment">% indices of the selected variables in the list data.ListVarName</span> |
---|
| 414 | 0329 DimCell=data.VarDimName{VarIndex(1)}; |
---|
| 415 | 0330 <span class="keyword">if</span> ischar(DimCell) |
---|
| 416 | 0331 DimCell={DimCell}; |
---|
| 417 | 0332 <span class="keyword">end</span> |
---|
| 418 | 0333 XName=DimCell{1}; <span class="comment">%first dimension considered as abscissa</span> |
---|
| 419 | 0334 coord_x_index_cell=VarType{icell}.coord(1); |
---|
| 420 | 0335 <span class="keyword">if</span> isequal(coord_x_index_cell,0) |
---|
| 421 | 0336 <span class="keyword">continue</span> <span class="comment">% the cell has no abscissa, skip it</span> |
---|
| 422 | 0337 <span class="keyword">end</span> |
---|
| 423 | 0338 <span class="keyword">if</span> ~isempty(coord_x_index)&&~isequal(coord_x_index_cell,coord_x_index) |
---|
| 424 | 0339 <span class="keyword">continue</span> <span class="comment">%all the selected variables must have the same first dimension</span> |
---|
| 425 | 0340 <span class="keyword">else</span> |
---|
| 426 | 0341 coord_x_index=coord_x_index_cell; |
---|
| 427 | 0342 <span class="keyword">end</span> |
---|
| 428 | 0343 testplot=ones(size(data.ListVarName));<span class="comment">%default test for plotted variables</span> |
---|
| 429 | 0344 testcoordvar=0; |
---|
| 430 | 0345 charplot_0=<span class="string">'''-'''</span>;<span class="comment">%default</span> |
---|
| 431 | 0346 <span class="keyword">if</span> isfield(data,<span class="string">'ObjectProjMode'</span>)& isequal(data.ObjectProjMode,<span class="string">'projection'</span>) |
---|
| 432 | 0347 charplot_0=<span class="string">'''+'''</span>; |
---|
| 433 | 0348 <span class="keyword">end</span> |
---|
| 434 | 0349 xtitle=<span class="string">''</span>; |
---|
| 435 | 0350 |
---|
| 436 | 0351 xtitle=data.ListVarName{coord_x_index}; |
---|
| 437 | 0352 eval([<span class="string">'coord_x{icell}=data.'</span> data.ListVarName{coord_x_index} <span class="string">';'</span>]);<span class="comment">%coordinate variable set as coord_x</span> |
---|
| 438 | 0353 <span class="keyword">if</span> isfield(data,<span class="string">'VarAttribute'</span>)&& numel(data.VarAttribute)>=coord_x_index && isfield(data.VarAttribute{coord_x_index},<span class="string">'units'</span>) |
---|
| 439 | 0354 xtitle=[xtitle <span class="string">'('</span> data.VarAttribute{coord_x_index}.units <span class="string">')'</span>]; |
---|
| 440 | 0355 <span class="keyword">end</span> |
---|
| 441 | 0356 eval([<span class="string">'coord_x{icell}=data.'</span> data.ListVarName{coord_x_index} <span class="string">';'</span>]);<span class="comment">%coordinate variable set as coord_x</span> |
---|
| 442 | 0357 testcoordvar=1; |
---|
| 443 | 0358 testplot(coord_x_index)=0; |
---|
| 444 | 0359 <span class="keyword">if</span> ~isempty(VarType{icell}.ancillary') |
---|
| 445 | 0360 testplot(VarType{icell}.ancillary)=0; |
---|
| 446 | 0361 <span class="keyword">end</span> |
---|
| 447 | 0362 <span class="keyword">if</span> ~isempty(VarType{icell}.warnflag') |
---|
| 448 | 0363 testplot(VarType{icell}.warnflag)=0; |
---|
| 449 | 0364 <span class="keyword">end</span> |
---|
| 450 | 0365 <span class="keyword">if</span> ~isempty(VarType{icell}.discrete') |
---|
| 451 | 0366 charplot_0=<span class="string">'''+'''</span>; |
---|
| 452 | 0367 <span class="keyword">else</span> |
---|
| 453 | 0368 charplot_0=<span class="string">'''-'''</span>; |
---|
| 454 | 0369 <span class="keyword">end</span> |
---|
| 455 | 0370 <span class="keyword">if</span> testcoordvar==0 |
---|
| 456 | 0371 coord_x{icell}=[1:data.DimValue(DimIndices(1))];<span class="comment">%abscissa by default if no coordinate variable</span> |
---|
| 457 | 0372 <span class="comment">% charplot_0='''-''';</span> |
---|
| 458 | 0373 <span class="keyword">end</span> |
---|
| 459 | 0374 <span class="keyword">if</span> isfield(data,<span class="string">'VarAttribute'</span>) |
---|
| 460 | 0375 VarAttribute=data.VarAttribute; |
---|
| 461 | 0376 <span class="keyword">for</span> ivar=1:length(VarIndex) |
---|
| 462 | 0377 <span class="keyword">if</span> length(VarAttribute)>=VarIndex(ivar) & isfield(VarAttribute{VarIndex(ivar)},<span class="string">'long_name'</span>) |
---|
| 463 | 0378 plotname{VarIndex(ivar)}=VarAttribute{VarIndex(ivar)}.long_name; |
---|
| 464 | 0379 <span class="keyword">else</span> |
---|
| 465 | 0380 plotname{VarIndex(ivar)}=data.ListVarName{VarIndex(ivar)};<span class="comment">%name for display in plot A METTRE</span> |
---|
| 466 | 0381 <span class="keyword">end</span> |
---|
| 467 | 0382 <span class="keyword">end</span> |
---|
| 468 | 0383 <span class="keyword">end</span> |
---|
| 469 | 0384 <span class="keyword">for</span> ivar=1:length(VarIndex) |
---|
| 470 | 0385 <span class="keyword">if</span> testplot(VarIndex(ivar)) |
---|
| 471 | 0386 VarName=data.ListVarName{VarIndex(ivar)}; |
---|
| 472 | 0387 eval([<span class="string">'data.'</span> VarName <span class="string">'=squeeze(data.'</span> VarName <span class="string">');'</span>]) |
---|
| 473 | 0388 <span class="keyword">if</span> isequal(VarName,<span class="string">'A'</span>) |
---|
| 474 | 0389 charplot=<span class="string">'''-'''</span>; |
---|
| 475 | 0390 <span class="keyword">else</span> |
---|
| 476 | 0391 charplot=charplot_0; |
---|
| 477 | 0392 <span class="keyword">end</span> |
---|
| 478 | 0393 plotstr=[plotstr <span class="string">'coord_x{'</span> num2str(icell) <span class="string">'},data.'</span> VarName <span class="string">','</span> charplot <span class="string">','</span>]; |
---|
| 479 | 0394 eval([<span class="string">'nbcomponent2=size(data.'</span> VarName <span class="string">',2);'</span>]); |
---|
| 480 | 0395 eval([<span class="string">'nbcomponent1=size(data.'</span> VarName <span class="string">',1);'</span>]); |
---|
| 481 | 0396 eval([<span class="string">'varmean=mean(double(data.'</span> VarName <span class="string">'));'</span>]);<span class="comment">%mean value</span> |
---|
| 482 | 0397 textmean=[textmean; {[VarName <span class="string">'mean= '</span> num2str(varmean,4)]}]; |
---|
| 483 | 0398 <span class="keyword">if</span> nbcomponent1==1| nbcomponent2==1 |
---|
| 484 | 0399 legend_str=[legend_str {VarName}]; <span class="comment">%variable with one component</span> |
---|
| 485 | 0400 <span class="keyword">else</span> <span class="comment">%variable with severals components</span> |
---|
| 486 | 0401 <span class="keyword">for</span> ic=1:min(nbcomponent1,nbcomponent2) |
---|
| 487 | 0402 legend_str=[legend_str [VarName <span class="string">'_'</span> num2str(ic)]]; <span class="comment">%variable with severals components</span> |
---|
| 488 | 0403 <span class="keyword">end</span> <span class="comment">% labeled by their index (e.g. color component)</span> |
---|
| 489 | 0404 <span class="keyword">end</span> |
---|
| 490 | 0405 <span class="keyword">end</span> |
---|
| 491 | 0406 <span class="keyword">end</span> |
---|
| 492 | 0407 <span class="keyword">end</span> |
---|
| 493 | 0408 <span class="keyword">if</span> ~isequal(plotstr,<span class="string">'plot('</span>) |
---|
| 494 | 0409 plotstr(end)=<span class="string">')'</span>; |
---|
| 495 | 0410 <span class="comment">%execute plot (instruction plotstr)</span> |
---|
| 496 | 0411 eval(plotstr) |
---|
| 497 | 0412 <span class="comment">%%%%%</span> |
---|
| 498 | 0413 grid on |
---|
| 499 | 0414 hxlabel=xlabel(xtitle); |
---|
| 500 | 0415 set(hxlabel,<span class="string">'Interpreter'</span>,<span class="string">'none'</span>)<span class="comment">% desable tex interpreter</span> |
---|
| 501 | 0416 <span class="keyword">if</span> length(legend_str)>=1 |
---|
| 502 | 0417 hylabel=ylabel(legend_str{end}); |
---|
| 503 | 0418 set(hylabel,<span class="string">'Interpreter'</span>,<span class="string">'none'</span>)<span class="comment">% desable tex interpreter</span> |
---|
| 504 | 0419 <span class="keyword">end</span> |
---|
| 505 | 0420 <span class="keyword">if</span> ~isempty(legend_str) |
---|
| 506 | 0421 hlegend=findobj(hfig,<span class="string">'Tag'</span>,<span class="string">'legend'</span>); |
---|
| 507 | 0422 <span class="keyword">if</span> isempty(hlegend) |
---|
| 508 | 0423 hlegend=legend(legend_str); |
---|
| 509 | 0424 txt=ver; |
---|
| 510 | 0425 Release=txt(1).Release; |
---|
| 511 | 0426 relnumb=str2num(Release(3:4)); |
---|
| 512 | 0427 <span class="keyword">if</span> relnumb >= 14 |
---|
| 513 | 0428 set(hlegend,<span class="string">'Interpreter'</span>,<span class="string">'none'</span>)<span class="comment">% desable tex interpreter</span> |
---|
| 514 | 0429 <span class="keyword">end</span> |
---|
| 515 | 0430 <span class="keyword">else</span> |
---|
| 516 | 0431 legend_old=get(hlegend,<span class="string">'String'</span>); |
---|
| 517 | 0432 <span class="keyword">if</span> isequal(size(legend_old,1),size(legend_str,1))&~isequal(legend_old,legend_str) |
---|
| 518 | 0433 set(hlegend,<span class="string">'String'</span>,[legend_old legend_str]); |
---|
| 519 | 0434 <span class="keyword">end</span> |
---|
| 520 | 0435 <span class="keyword">end</span> |
---|
| 521 | 0436 <span class="keyword">end</span> |
---|
| 522 | 0437 title_str=<span class="string">''</span>; |
---|
| 523 | 0438 <span class="keyword">if</span> isfield(data,<span class="string">'filename'</span>) |
---|
| 524 | 0439 [Path, title_str, ext]=fileparts(data.filename); |
---|
| 525 | 0440 title_str=[title_str ext]; |
---|
| 526 | 0441 <span class="keyword">end</span> |
---|
| 527 | 0442 <span class="keyword">if</span> isfield(data,<span class="string">'Action'</span>) |
---|
| 528 | 0443 <span class="keyword">if</span> ~isequal(title_str,<span class="string">''</span>) |
---|
| 529 | 0444 title_str=[title_str <span class="string">', '</span>]; |
---|
| 530 | 0445 <span class="keyword">end</span> |
---|
| 531 | 0446 title_str=[title_str data.Action]; |
---|
| 532 | 0447 <span class="keyword">end</span> |
---|
| 533 | 0448 htitle=title(title_str); |
---|
| 534 | 0449 txt=ver; |
---|
| 535 | 0450 Release=txt(1).Release; |
---|
| 536 | 0451 relnumb=str2num(Release(3:4)); |
---|
| 537 | 0452 <span class="keyword">if</span> relnumb >= 14 |
---|
| 538 | 0453 set(htitle,<span class="string">'Interpreter'</span>,<span class="string">'none'</span>)<span class="comment">% desable tex interpreter</span> |
---|
| 539 | 0454 <span class="keyword">end</span> |
---|
| 540 | 0455 <span class="comment">% A REPRENDRE Mean</span> |
---|
| 541 | 0456 <span class="comment">% hlist=findobj(gcf,'Style','listbox','Tag','liststat');</span> |
---|
| 542 | 0457 <span class="comment">% if isempty(hlist)</span> |
---|
| 543 | 0458 <span class="comment">% 'text'</span> |
---|
| 544 | 0459 <span class="comment">% textmean</span> |
---|
| 545 | 0460 <span class="comment">% set(gca,'position',[0.13,0.2,0.775,0.73])</span> |
---|
| 546 | 0461 <span class="comment">% uicontrol('Style','popupmenu','Position',[20 20 200 20],'String',textmean,'Tag','liststat');</span> |
---|
| 547 | 0462 <span class="comment">% else</span> |
---|
| 548 | 0463 <span class="comment">% set(hlist(1),'String',textmean)</span> |
---|
| 549 | 0464 <span class="comment">% end</span> |
---|
| 550 | 0465 <span class="keyword">end</span> |
---|
| 551 | 0466 |
---|
| 552 | 0467 |
---|
| 553 | 0468 <span class="comment">%---------------------------------------</span> |
---|
| 554 | 0469 <span class="comment">% plot_plane</span> |
---|
| 555 | 0470 <span class="comment">%----------------------------------------</span> |
---|
| 556 | 0471 <a name="_sub3" href="#_subfunctions" class="code">function [AxeData,haxes,PlotParamOut,PlotType]=plot_plane(Data,CellVarIndex,VarTypeCell,haxes,PlotParam,KeepLim,PosColorbar)</a> |
---|
| 557 | 0472 |
---|
| 558 | 0473 <span class="comment">%default plotting parameters</span> |
---|
| 559 | 0474 PlotType=<span class="string">'plane'</span>;<span class="comment">%default</span> |
---|
| 560 | 0475 <span class="keyword">if</span> ~exist(<span class="string">'PlotParam'</span>,<span class="string">'var'</span>) |
---|
| 561 | 0476 PlotParam=[]; |
---|
| 562 | 0477 <span class="keyword">end</span> |
---|
| 563 | 0478 <span class="keyword">if</span> ~isfield(PlotParam,<span class="string">'Scalar'</span>) |
---|
| 564 | 0479 PlotParam.Scalar=[]; |
---|
| 565 | 0480 <span class="keyword">end</span> |
---|
| 566 | 0481 <span class="keyword">if</span> ~isfield(PlotParam,<span class="string">'Vectors'</span>) |
---|
| 567 | 0482 PlotParam.Vectors=[]; |
---|
| 568 | 0483 <span class="keyword">end</span> |
---|
| 569 | 0484 PlotParamOut=PlotParam;<span class="comment">%default</span> |
---|
| 570 | 0485 |
---|
| 571 | 0486 <span class="comment">%plotting axes</span> |
---|
| 572 | 0487 hfig=get(haxes,<span class="string">'parent'</span>); |
---|
| 573 | 0488 hcol=findobj(hfig,<span class="string">'Tag'</span>,<span class="string">'Colorbar'</span>); <span class="comment">%look for colorbar axes</span> |
---|
| 574 | 0489 hima=findobj(haxes,<span class="string">'Tag'</span>,<span class="string">'ima'</span>);<span class="comment">% search existing image in the current axes</span> |
---|
| 575 | 0490 AxeData=get(haxes,<span class="string">'UserData'</span>); <span class="comment">%default</span> |
---|
| 576 | 0491 <span class="keyword">if</span> ~isstruct(AxeData)<span class="comment">% AxeData must be a structure</span> |
---|
| 577 | 0492 AxeData=[]; |
---|
| 578 | 0493 <span class="keyword">end</span> |
---|
| 579 | 0494 AxeData.NbDim=2; |
---|
| 580 | 0495 <span class="keyword">if</span> isfield(Data,<span class="string">'ObjectCoord'</span>) |
---|
| 581 | 0496 AxeData.ObjectCoord=Data.ObjectCoord; |
---|
| 582 | 0497 <span class="keyword">end</span> |
---|
| 583 | 0498 |
---|
| 584 | 0499 test_ima=0; <span class="comment">%default: test for image or map plot</span> |
---|
| 585 | 0500 test_vec=0; <span class="comment">%default: test for vector plots</span> |
---|
| 586 | 0501 test_black=0; |
---|
| 587 | 0502 test_false=0; |
---|
| 588 | 0503 test_C=0; |
---|
| 589 | 0504 <span class="keyword">for</span> icell=1:length(CellVarIndex) <span class="comment">% length(CellVarIndex) =1 or 2 (from the calling function)</span> |
---|
| 590 | 0505 <span class="comment">% VarIndex=CellVarIndex{icell};</span> |
---|
| 591 | 0506 VarType=VarTypeCell{icell}; |
---|
| 592 | 0507 ivar_X=VarType.coord_x; <span class="comment">% defines (unique) index for the variable representing unstructured x coordinate (default =[])</span> |
---|
| 593 | 0508 ivar_Y=VarType.coord_y; <span class="comment">% defines (unique)index for the variable representing unstructured y coordinate (default =[])</span> |
---|
| 594 | 0509 ivar_U=VarType.vector_x; <span class="comment">% defines (unique) index for the variable representing x vector component (default =[])</span> |
---|
| 595 | 0510 ivar_V=VarType.vector_y; <span class="comment">% defines (unique) index for the variable representing y vector component (default =[])</span> |
---|
| 596 | 0511 ivar_C=[VarType.scalar VarType.image VarType.color VarType.ancillary]; <span class="comment">%defines index (indices) for the scalar or ancillary fields</span> |
---|
| 597 | 0512 <span class="keyword">if</span> numel(ivar_C)>1 |
---|
| 598 | 0513 <a href="msgbox_uvmat.html" class="code" title="function varargout = msgbox_uvmat(varargin)">msgbox_uvmat</a>(<span class="string">'ERROR'</span>,<span class="string">'too many scalar inputs in plot_field.m'</span>) |
---|
| 599 | 0514 <span class="keyword">return</span> |
---|
| 600 | 0515 <span class="keyword">end</span> |
---|
| 601 | 0516 ivar_F=VarType.warnflag; <span class="comment">%defines index (unique) for warning flag variable</span> |
---|
| 602 | 0517 ivar_FF=VarType.errorflag; <span class="comment">%defines index (unique) for error flag variable</span> |
---|
| 603 | 0518 idim_Y=[]; |
---|
| 604 | 0519 test_grid=0; |
---|
| 605 | 0520 <span class="keyword">if</span> ~isempty(ivar_U) && ~isempty(ivar_V)<span class="comment">% vector components detected</span> |
---|
| 606 | 0521 <span class="keyword">if</span> test_vec |
---|
| 607 | 0522 <a href="msgbox_uvmat.html" class="code" title="function varargout = msgbox_uvmat(varargin)">msgbox_uvmat</a>(<span class="string">'ERROR'</span>,<span class="string">'attempt to plot two vector fields'</span>) |
---|
| 608 | 0523 <span class="keyword">return</span> |
---|
| 609 | 0524 <span class="keyword">else</span> |
---|
| 610 | 0525 test_vec=1; |
---|
| 611 | 0526 eval([<span class="string">'vec_U=Data.'</span> Data.ListVarName{ivar_U} <span class="string">';'</span>]) |
---|
| 612 | 0527 eval([<span class="string">'vec_V=Data.'</span> Data.ListVarName{ivar_V} <span class="string">';'</span>]) |
---|
| 613 | 0528 <span class="keyword">if</span> ~isempty(ivar_X) && ~isempty(ivar_Y)<span class="comment">% 2D field (with unstructured coordinates or structured ones (then ivar_X and ivar_Y empty)</span> |
---|
| 614 | 0529 eval([<span class="string">'vec_X=Data.'</span> Data.ListVarName{ivar_X} <span class="string">';'</span>]) |
---|
| 615 | 0530 eval([<span class="string">'vec_Y=Data.'</span> Data.ListVarName{ivar_Y} <span class="string">';'</span>]) |
---|
| 616 | 0531 <span class="keyword">elseif</span> numel(VarType.coord)==2 & VarType.coord~=[0 0];<span class="comment">%coordinates defines by dimension variables</span> |
---|
| 617 | 0532 eval([<span class="string">'y=Data.'</span> Data.ListVarName{VarType.coord(1)} <span class="string">';'</span>]) |
---|
| 618 | 0533 eval([<span class="string">'x=Data.'</span> Data.ListVarName{VarType.coord(2)} <span class="string">';'</span>]) |
---|
| 619 | 0534 <span class="keyword">if</span> numel(y)==2 <span class="comment">% y defined by first and last values on aregular mesh</span> |
---|
| 620 | 0535 y=linspace(y(1),y(2),size(vec_U,1)); |
---|
| 621 | 0536 <span class="keyword">end</span> |
---|
| 622 | 0537 <span class="keyword">if</span> numel(x)==2 <span class="comment">% y defined by first and last values on aregular mesh</span> |
---|
| 623 | 0538 x=linspace(x(1),x(2),size(vec_U,2)); |
---|
| 624 | 0539 <span class="keyword">end</span> |
---|
| 625 | 0540 [vec_X,vec_Y]=meshgrid(x,y); |
---|
| 626 | 0541 <span class="keyword">else</span> |
---|
| 627 | 0542 <a href="msgbox_uvmat.html" class="code" title="function varargout = msgbox_uvmat(varargin)">msgbox_uvmat</a>(<span class="string">'ERROR'</span>,<span class="string">'error in plot_field: bad coordinate definition for vector field'</span>) |
---|
| 628 | 0543 <span class="keyword">return</span> |
---|
| 629 | 0544 <span class="keyword">end</span> |
---|
| 630 | 0545 <span class="keyword">if</span> ~isempty(ivar_C) |
---|
| 631 | 0546 eval([<span class="string">'vec_C=Data.'</span> Data.ListVarName{ivar_C} <span class="string">';'</span>]) ; |
---|
| 632 | 0547 vec_C=reshape(vec_C,1,numel(vec_C)); |
---|
| 633 | 0548 test_C=1; |
---|
| 634 | 0549 <span class="keyword">end</span> |
---|
| 635 | 0550 <span class="keyword">if</span> ~(isfield(PlotParam.Vectors,<span class="string">'HideWarning'</span>)&& isequal(PlotParam.Vectors.HideWarning,1))&&~isempty(ivar_F) &&~test_black |
---|
| 636 | 0551 <span class="keyword">if</span> numel(ivar_F)>1 |
---|
| 637 | 0552 <a href="warndlg_uvmat.html" class="code" title="function hwarn=warndlg_uvmat(warntext,title)">warndlg_uvmat</a>(<span class="string">'multiple warning flags in plot_field.m, make a selection'</span>,<span class="string">'ERROR'</span>) |
---|
| 638 | 0553 <span class="keyword">return</span> |
---|
| 639 | 0554 <span class="keyword">elseif</span> test_vec |
---|
| 640 | 0555 eval([<span class="string">'AxeData.F=Data.'</span> Data.ListVarName{ivar_F} <span class="string">';'</span>]) <span class="comment">% flags for false vectors</span> |
---|
| 641 | 0556 test_black=1; |
---|
| 642 | 0557 <span class="keyword">end</span> |
---|
| 643 | 0558 <span class="keyword">end</span> |
---|
| 644 | 0559 <span class="keyword">if</span> ~isempty(ivar_FF) && ~test_false |
---|
| 645 | 0560 <span class="keyword">if</span> numel(ivar_FF)>1 |
---|
| 646 | 0561 <a href="warndlg_uvmat.html" class="code" title="function hwarn=warndlg_uvmat(warntext,title)">warndlg_uvmat</a>(<span class="string">'multiple error flags in plot_field.m, make a selection'</span>,<span class="string">'ERROR'</span>) |
---|
| 647 | 0562 <span class="keyword">return</span> |
---|
| 648 | 0563 <span class="keyword">else</span> <span class="comment">%A VOIR CAS IMAGES</span> |
---|
| 649 | 0564 eval([<span class="string">'AxeData.FF=Data.'</span> Data.ListVarName{ivar_FF} <span class="string">';'</span>]) <span class="comment">% flags for false vectors</span> |
---|
| 650 | 0565 test_false=1; |
---|
| 651 | 0566 <span class="keyword">end</span> |
---|
| 652 | 0567 <span class="keyword">end</span> |
---|
| 653 | 0568 <span class="keyword">end</span> |
---|
| 654 | 0569 <span class="keyword">elseif</span> ~isempty(ivar_C) <span class="comment">%scalar or image</span> |
---|
| 655 | 0570 <span class="keyword">if</span> test_ima |
---|
| 656 | 0571 <a href="msgbox_uvmat.html" class="code" title="function varargout = msgbox_uvmat(varargin)">msgbox_uvmat</a>(<span class="string">'ERROR'</span>,<span class="string">'attempt to plot two scalar fields or images'</span>) |
---|
| 657 | 0572 <span class="keyword">return</span> |
---|
| 658 | 0573 <span class="keyword">end</span> |
---|
| 659 | 0574 eval([<span class="string">'A=Data.'</span> Data.ListVarName{ivar_C} <span class="string">';'</span>]) ;<span class="comment">% scalar represented as color image</span> |
---|
| 660 | 0575 test_ima=1; |
---|
| 661 | 0576 <span class="keyword">if</span> ~isempty(ivar_X) && ~isempty(ivar_Y)<span class="comment">% 2D field (with unstructured coordinates or structured ones (then ivar_X and ivar_Y empty)</span> |
---|
| 662 | 0577 eval([<span class="string">'AX=Data.'</span> Data.ListVarName{ivar_X} <span class="string">';'</span>]) |
---|
| 663 | 0578 eval([<span class="string">'AY=Data.'</span> Data.ListVarName{ivar_Y} <span class="string">';'</span>]) |
---|
| 664 | 0579 [A,AX,AY]=<a href="proj_grid.html" class="code" title="function [A,rangx,rangy]=proj_grid(vec_X,vec_Y,vec_A,rgx_in,rgy_in,npxy_in)">proj_grid</a>(AX',AY',A',[],[],<span class="string">'np>256'</span>); <span class="comment">% interpolate on a grid</span> |
---|
| 665 | 0580 <span class="keyword">elseif</span> numel(VarType.coord)>=2 & VarType.coord(1:2)~=[0 0];<span class="comment">%structured coordinates</span> |
---|
| 666 | 0581 eval([<span class="string">'AY=Data.'</span> Data.ListVarName{VarType.coord(1)} <span class="string">';'</span>]) |
---|
| 667 | 0582 eval([<span class="string">'AX=Data.'</span> Data.ListVarName{VarType.coord(2)} <span class="string">';'</span>]) |
---|
| 668 | 0583 test_interp_X=0; <span class="comment">%default, regularly meshed X coordinate</span> |
---|
| 669 | 0584 test_interp_Y=0; <span class="comment">%default, regularly meshed Y coordinate</span> |
---|
| 670 | 0585 <span class="keyword">if</span> numel(AY)>2 |
---|
| 671 | 0586 DAY=diff(AY); |
---|
| 672 | 0587 DAY_min=min(DAY); |
---|
| 673 | 0588 DAY_max=max(DAY); |
---|
| 674 | 0589 <span class="keyword">if</span> sign(DAY_min)~=sign(DAY_max);<span class="comment">% =1 for increasing values, 0 otherwise</span> |
---|
| 675 | 0590 errormsg=[<span class="string">'errror in plot_field.m: non monotonic dimension variable # '</span> ListVarName{VarType.coord(1)} ]; |
---|
| 676 | 0591 <span class="keyword">return</span> |
---|
| 677 | 0592 <span class="keyword">end</span> |
---|
| 678 | 0593 test_interp_Y=(DAY_max-DAY_min)> 0.0001*abs(DAY_max); |
---|
| 679 | 0594 <span class="keyword">end</span> |
---|
| 680 | 0595 <span class="keyword">if</span> numel(AX)>2 |
---|
| 681 | 0596 DAX=diff(AX); |
---|
| 682 | 0597 DAX_min=min(DAX); |
---|
| 683 | 0598 DAX_max=max(DAX); |
---|
| 684 | 0599 <span class="keyword">if</span> sign(DAX_min)~=sign(DAX_max);<span class="comment">% =1 for increasing values, 0 otherwise</span> |
---|
| 685 | 0600 errormsg=[<span class="string">'errror in plot_field.m: non monotonic dimension variable # '</span> ListVarName{VarType.coord(2)} ]; |
---|
| 686 | 0601 <span class="keyword">return</span> |
---|
| 687 | 0602 <span class="keyword">end</span> |
---|
| 688 | 0603 test_interp_X=(DAX_max-DAX_min)> 0.0001*abs(DAX_max); |
---|
| 689 | 0604 <span class="keyword">end</span> |
---|
| 690 | 0605 <span class="keyword">if</span> test_interp_Y |
---|
| 691 | 0606 npxy(1)=max([256 floor((AY(end)-AY(1))/DAY_min) floor((AY(end)-AY(1))/DAY_max)]); |
---|
| 692 | 0607 yI=linspace(AY(1),AY(end),npxy(1)); |
---|
| 693 | 0608 <span class="keyword">if</span> ~test_interp_X |
---|
| 694 | 0609 xI=linspace(AX(1),AX(2),size(A,2));<span class="comment">%default</span> |
---|
| 695 | 0610 AX=xI; |
---|
| 696 | 0611 <span class="keyword">end</span> |
---|
| 697 | 0612 <span class="keyword">end</span> |
---|
| 698 | 0613 <span class="keyword">if</span> test_interp_X |
---|
| 699 | 0614 npxy(1)=max([256 floor((AX(end)-AX(1))/DAX_min) floor((AX(end)-AX(1))/DAX_max)]); |
---|
| 700 | 0615 xI=linspace(AX(1),AX(end),npxy(2)); |
---|
| 701 | 0616 <span class="keyword">if</span> ~test_interp_Y |
---|
| 702 | 0617 yI=linspace(AY(1),AY(end),size(A,1)); |
---|
| 703 | 0618 AY=yI; |
---|
| 704 | 0619 <span class="keyword">end</span> |
---|
| 705 | 0620 <span class="keyword">end</span> |
---|
| 706 | 0621 <span class="keyword">if</span> test_interp_X || test_interp_Y |
---|
| 707 | 0622 [AX2D,AY2D]=meshgrid(AX,AY); |
---|
| 708 | 0623 A=interp2(AX2D,AY2D,double(A),xI,yI'); |
---|
| 709 | 0624 <span class="keyword">end</span> |
---|
| 710 | 0625 AX=[AX(1) AX(end)];<span class="comment">% keep only the lower and upper bounds for image represnetation</span> |
---|
| 711 | 0626 AY=[AY(1) AY(end)]; |
---|
| 712 | 0627 <span class="keyword">else</span> |
---|
| 713 | 0628 <a href="msgbox_uvmat.html" class="code" title="function varargout = msgbox_uvmat(varargin)">msgbox_uvmat</a>(<span class="string">'ERROR'</span>,<span class="string">'bad coordinate definition for vector field'</span>) |
---|
| 714 | 0629 <span class="keyword">return</span> |
---|
| 715 | 0630 <span class="keyword">end</span> |
---|
| 716 | 0631 <span class="keyword">end</span> |
---|
| 717 | 0632 <span class="keyword">if</span> isfield(Data,<span class="string">'VarAttribute'</span>) |
---|
| 718 | 0633 VarAttribute=Data.VarAttribute; |
---|
| 719 | 0634 <span class="keyword">end</span> |
---|
| 720 | 0635 <span class="keyword">end</span> |
---|
| 721 | 0636 |
---|
| 722 | 0637 <span class="comment">%%%%%%%%%%%%%%%%%%%%% image or scalar plot %%%%%%%%%%%%%%%%%%%%%%%%%%</span> |
---|
| 723 | 0638 |
---|
| 724 | 0639 <span class="keyword">if</span> ~isfield(PlotParam.Scalar,<span class="string">'Contours'</span>) |
---|
| 725 | 0640 PlotParam.Scalar.Contours=0; <span class="comment">%default</span> |
---|
| 726 | 0641 <span class="keyword">end</span> |
---|
| 727 | 0642 PlotParamOut=PlotParam; <span class="comment">%default</span> |
---|
| 728 | 0643 <span class="keyword">if</span> test_ima |
---|
| 729 | 0644 <span class="comment">% distinguish B/W and color images</span> |
---|
| 730 | 0645 np=size(A);<span class="comment">%size of image</span> |
---|
| 731 | 0646 siz=size(np); |
---|
| 732 | 0647 <span class="comment">%set the color map</span> |
---|
| 733 | 0648 <span class="keyword">if</span> siz(2)==2 <span class="comment">%for black and white images</span> |
---|
| 734 | 0649 <span class="keyword">if</span> ~isfield(PlotParam.Scalar,<span class="string">'AutoScal'</span>) |
---|
| 735 | 0650 PlotParam.Scalar.AutoScal=0;<span class="comment">%default</span> |
---|
| 736 | 0651 <span class="keyword">end</span> |
---|
| 737 | 0652 <span class="keyword">if</span> ~isfield(PlotParam.Scalar,<span class="string">'MinA'</span>) |
---|
| 738 | 0653 PlotParam.Scalar.MinA=[];<span class="comment">%default</span> |
---|
| 739 | 0654 <span class="keyword">end</span> |
---|
| 740 | 0655 <span class="keyword">if</span> ~isfield(PlotParam.Scalar,<span class="string">'MaxA'</span>) |
---|
| 741 | 0656 PlotParam.Scalar.MaxA=[];<span class="comment">%default</span> |
---|
| 742 | 0657 <span class="keyword">end</span> |
---|
| 743 | 0658 <span class="keyword">if</span> isequal(PlotParam.Scalar.AutoScal,0)|isempty(PlotParam.Scalar.MinA)|~isa(PlotParam.Scalar.MinA,<span class="string">'double'</span>) <span class="comment">%correct if there is no numerical data in edit box</span> |
---|
| 744 | 0659 MinA=double(min(min(A))); |
---|
| 745 | 0660 <span class="keyword">else</span> |
---|
| 746 | 0661 MinA=PlotParam.Scalar.MinA; |
---|
| 747 | 0662 <span class="keyword">end</span>; |
---|
| 748 | 0663 <span class="keyword">if</span> isequal(PlotParam.Scalar.AutoScal,0)|isempty(PlotParam.Scalar.MaxA)|~isa(PlotParam.Scalar.MaxA,<span class="string">'double'</span>) <span class="comment">%correct if there is no numerical data in edit box</span> |
---|
| 749 | 0664 MaxA=double(max(max(A))); |
---|
| 750 | 0665 <span class="keyword">else</span> |
---|
| 751 | 0666 MaxA=PlotParam.Scalar.MaxA; |
---|
| 752 | 0667 <span class="keyword">end</span>; |
---|
| 753 | 0668 PlotParamOut.Scalar.MinA=MinA; |
---|
| 754 | 0669 PlotParamOut.Scalar.MaxA=MaxA; |
---|
| 755 | 0670 axes(haxes) |
---|
| 756 | 0671 <span class="keyword">if</span> isequal(PlotParam.Scalar.Contours,1) |
---|
| 757 | 0672 <span class="keyword">if</span> ~isempty(hima) & ishandle(hima) |
---|
| 758 | 0673 delete(hima) |
---|
| 759 | 0674 <span class="keyword">end</span> |
---|
| 760 | 0675 <span class="keyword">if</span> ~isfield(PlotParam.Scalar,<span class="string">'IncrA'</span>) |
---|
| 761 | 0676 PlotParam.Scalar.IncrA=[]; |
---|
| 762 | 0677 <span class="keyword">end</span> |
---|
| 763 | 0678 <span class="keyword">if</span> isempty(PlotParam.Scalar.IncrA)<span class="comment">% | PlotParam.Scalar.AutoScal==0</span> |
---|
| 764 | 0679 cont=<a href="#_sub5" class="code" title="subfunction YTick=colbartick(MinA,MaxA)">colbartick</a>(MinA,MaxA); |
---|
| 765 | 0680 intercont=cont(2)-cont(1);<span class="comment">%default</span> |
---|
| 766 | 0681 PlotParamOut.Scalar.IncrA=intercont; |
---|
| 767 | 0682 <span class="keyword">else</span> |
---|
| 768 | 0683 intercont=PlotParam.Scalar.IncrA; |
---|
| 769 | 0684 <span class="keyword">end</span> |
---|
| 770 | 0685 B=A; |
---|
| 771 | 0686 abscontmin=intercont*floor(MinA/intercont); |
---|
| 772 | 0687 abscontmax=intercont*ceil(MaxA/intercont); |
---|
| 773 | 0688 contmin=intercont*floor(min(min(B))/intercont); |
---|
| 774 | 0689 contmax=intercont*ceil(max(max(B))/intercont); |
---|
| 775 | 0690 cont_pos_plus=[0:intercont:contmax]; |
---|
| 776 | 0691 cont_pos_min=[double(contmin):intercont:-intercont]; |
---|
| 777 | 0692 cont_pos=[cont_pos_min cont_pos_plus]; |
---|
| 778 | 0693 sizpx=(AX(end)-AX(1))/(np(2)-1); |
---|
| 779 | 0694 sizpy=(AY(1)-AY(end))/(np(1)-1); |
---|
| 780 | 0695 x_cont=[AX(1):sizpx:AX(end)]; <span class="comment">% pixel x coordinates for image display</span> |
---|
| 781 | 0696 y_cont=[AY(1):-sizpy:AY(end)]; <span class="comment">% pixel x coordinates for image display</span> |
---|
| 782 | 0697 txt=ver;<span class="comment">%version of Matlab</span> |
---|
| 783 | 0698 Release=txt(1).Release; |
---|
| 784 | 0699 relnumb=str2num(Release(3:4)); |
---|
| 785 | 0700 <span class="keyword">if</span> relnumb >= 14 |
---|
| 786 | 0701 vec=linspace(0,1,(abscontmax-abscontmin)/intercont);<span class="comment">%define a greyscale colormap with steps intercont</span> |
---|
| 787 | 0702 map=[vec' vec' vec']; |
---|
| 788 | 0703 colormap(map); |
---|
| 789 | 0704 [var,hcontour]=contour(x_cont,y_cont,B,cont_pos); |
---|
| 790 | 0705 set(hcontour,<span class="string">'Fill'</span>,<span class="string">'on'</span>) |
---|
| 791 | 0706 set(hcontour,<span class="string">'LineStyle'</span>,<span class="string">'none'</span>) |
---|
| 792 | 0707 hold on |
---|
| 793 | 0708 <span class="keyword">end</span> |
---|
| 794 | 0709 [var_p,hcontour_p]=contour(x_cont,y_cont,B,cont_pos_plus,<span class="string">'k-'</span>); |
---|
| 795 | 0710 hold on |
---|
| 796 | 0711 [var_m,hcontour_m]=contour(x_cont,y_cont,B,cont_pos_min,<span class="string">':'</span>); |
---|
| 797 | 0712 set(hcontour_m,<span class="string">'LineColor'</span>,[1 1 1]) |
---|
| 798 | 0713 hold off |
---|
| 799 | 0714 caxis([abscontmin abscontmax]) |
---|
| 800 | 0715 colormap(map); |
---|
| 801 | 0716 <span class="keyword">end</span> |
---|
| 802 | 0717 <span class="keyword">if</span> ~isequal(PlotParam.Scalar.Contours,1) |
---|
| 803 | 0718 <span class="comment">% rescale the grey levels with min and max, put a grey scale colorbar</span> |
---|
| 804 | 0719 <span class="keyword">if</span> (isfield(PlotParam.Scalar,<span class="string">'BW'</span>)& (isequal(PlotParam.Scalar.BW,1))|isa(A,<span class="string">'uint8'</span>)| isa(A,<span class="string">'uint16'</span>))<span class="comment">%images</span> |
---|
| 805 | 0720 B=A; |
---|
| 806 | 0721 vec=linspace(0,1,255);<span class="comment">%define a linear greyscale colormap</span> |
---|
| 807 | 0722 map=[vec' vec' vec']; |
---|
| 808 | 0723 colormap(map); <span class="comment">%grey scale color map</span> |
---|
| 809 | 0724 <span class="keyword">else</span> |
---|
| 810 | 0725 B=A; |
---|
| 811 | 0726 colormap(<span class="string">'default'</span>); <span class="comment">% standard faulse colors for div, vort , scalar fields</span> |
---|
| 812 | 0727 <span class="keyword">end</span> |
---|
| 813 | 0728 <span class="keyword">end</span> |
---|
| 814 | 0729 <span class="keyword">elseif</span> siz(2)==3 <span class="comment">%color images</span> |
---|
| 815 | 0730 axes(haxes) |
---|
| 816 | 0731 B=uint8(A); |
---|
| 817 | 0732 MinA=0; |
---|
| 818 | 0733 MaxA=255; |
---|
| 819 | 0734 <span class="keyword">end</span> |
---|
| 820 | 0735 <span class="keyword">if</span> ~isequal(PlotParam.Scalar.Contours,1) |
---|
| 821 | 0736 <span class="comment">%interpolate to increase resolution</span> |
---|
| 822 | 0737 test_interp=1; |
---|
| 823 | 0738 <span class="keyword">if</span> max(np) <= 64 |
---|
| 824 | 0739 npxy=8*np;<span class="comment">% increase the resolution 8 times</span> |
---|
| 825 | 0740 <span class="keyword">elseif</span> max(np) <= 128 |
---|
| 826 | 0741 npxy=4*np;<span class="comment">% increase the resolution 4 times</span> |
---|
| 827 | 0742 <span class="keyword">elseif</span> max(np) <= 256 |
---|
| 828 | 0743 npxy=2*np;<span class="comment">% increase the resolution 2 times</span> |
---|
| 829 | 0744 <span class="keyword">else</span> |
---|
| 830 | 0745 npxy=np; |
---|
| 831 | 0746 test_interp=0; <span class="comment">% no interpolation done</span> |
---|
| 832 | 0747 <span class="keyword">end</span> |
---|
| 833 | 0748 <span class="keyword">if</span> test_interp==1<span class="comment">%if we interpolate</span> |
---|
| 834 | 0749 x=[AX(1):(AX(2)-AX(1))/(np(2)-1):AX(2)]; |
---|
| 835 | 0750 y=[AY(1):(AY(2)-AY(1))/(np(1)-1):AY(2)]; |
---|
| 836 | 0751 [X,Y]=meshgrid(x,y); |
---|
| 837 | 0752 <span class="comment">% xi=[AX(1):(AX(2)-AX(1))/(npxy(2)-1):AX(2)];</span> |
---|
| 838 | 0753 <span class="comment">% yi=[AY(1):(AY(2)-AY(1))/(npxy(1)-1):AY(2)];</span> |
---|
| 839 | 0754 |
---|
| 840 | 0755 xi=linspace(AX(1),AX(2),npxy(2)); |
---|
| 841 | 0756 yi=linspace(AY(1),AY(2),npxy(1)); |
---|
| 842 | 0757 B = interp2(X,Y,double(B),xi,yi'); |
---|
| 843 | 0758 <span class="keyword">end</span> |
---|
| 844 | 0759 <span class="keyword">if</span> isempty(hima) |
---|
| 845 | 0760 tag=get(haxes,<span class="string">'Tag'</span>); |
---|
| 846 | 0761 hima=imagesc(AX,AY,B,[MinA MaxA]); |
---|
| 847 | 0762 set(hima,<span class="string">'Tag'</span>,<span class="string">'ima'</span>,<span class="string">'HitTest'</span>,<span class="string">'off'</span>) |
---|
| 848 | 0763 set(haxes,<span class="string">'Tag'</span>,tag);<span class="comment">%preserve the axes tag (removed by image fct !!!)</span> |
---|
| 849 | 0764 <span class="comment">%axis(haxes,'image'); % preserve the image aspect ratio %TESTMOD</span> |
---|
| 850 | 0765 <span class="keyword">else</span> |
---|
| 851 | 0766 set(hima,<span class="string">'CData'</span>,B); |
---|
| 852 | 0767 <span class="keyword">if</span> MinA<MaxA |
---|
| 853 | 0768 caxis([MinA MaxA]) |
---|
| 854 | 0769 <span class="keyword">else</span> |
---|
| 855 | 0770 caxis([MinA MinA+1]) |
---|
| 856 | 0771 <span class="keyword">end</span> |
---|
| 857 | 0772 set(hima,<span class="string">'XData'</span>,AX); |
---|
| 858 | 0773 set(hima,<span class="string">'YData'</span>,AY); |
---|
| 859 | 0774 <span class="keyword">end</span> |
---|
| 860 | 0775 <span class="keyword">end</span> |
---|
| 861 | 0776 <span class="keyword">if</span> ~isstruct(AxeData) |
---|
| 862 | 0777 AxeData=[]; |
---|
| 863 | 0778 <span class="keyword">end</span> |
---|
| 864 | 0779 AxeData.A=A; |
---|
| 865 | 0780 AxeData.AX=[AX(1) AX(end)]; |
---|
| 866 | 0781 AxeData.AY=[AY(1) AY(end)]; |
---|
| 867 | 0782 test_ima=1; |
---|
| 868 | 0783 <span class="comment">%display the colorbar code for B/W images if Poscolorbar not empty</span> |
---|
| 869 | 0784 <span class="keyword">if</span> siz(2)==2 & exist(<span class="string">'PosColorbar'</span>,<span class="string">'var'</span>)& ~isempty(PosColorbar) |
---|
| 870 | 0785 <span class="keyword">if</span> isempty(hcol)|~ishandle(hcol) |
---|
| 871 | 0786 hcol=colorbar;<span class="comment">%create new colorbar</span> |
---|
| 872 | 0787 <span class="keyword">end</span> |
---|
| 873 | 0788 <span class="keyword">if</span> length(PosColorbar)==4 |
---|
| 874 | 0789 set(hcol,<span class="string">'Position'</span>,PosColorbar) |
---|
| 875 | 0790 <span class="keyword">end</span> |
---|
| 876 | 0791 YTick=0;<span class="comment">%default</span> |
---|
| 877 | 0792 <span class="keyword">if</span> MaxA>MinA |
---|
| 878 | 0793 <span class="keyword">if</span> isequal(PlotParam.Scalar.Contours,1) |
---|
| 879 | 0794 colbarlim=get(hcol,<span class="string">'YLim'</span>); |
---|
| 880 | 0795 scale_bar=(colbarlim(2)-colbarlim(1))/(abscontmax-abscontmin); |
---|
| 881 | 0796 YTick=cont_pos(2:end-1); |
---|
| 882 | 0797 YTick_scaled=colbarlim(1)+scale_bar*(YTick-abscontmin); |
---|
| 883 | 0798 set(hcol,<span class="string">'YTick'</span>,YTick_scaled); |
---|
| 884 | 0799 <span class="keyword">elseif</span> (isfield(PlotParam.Scalar,<span class="string">'BW'</span>) & isequal(PlotParam.Scalar.BW,1))|isa(A,<span class="string">'uint8'</span>)| isa(A,<span class="string">'uint16'</span>)<span class="comment">%images</span> |
---|
| 885 | 0800 hi=get(hcol,<span class="string">'children'</span>); |
---|
| 886 | 0801 <span class="keyword">if</span> iscell(hi)<span class="comment">%multiple images in colorbar</span> |
---|
| 887 | 0802 hi=hi{1}; |
---|
| 888 | 0803 <span class="keyword">end</span> |
---|
| 889 | 0804 set(hi,<span class="string">'YData'</span>,[MinA MaxA]) |
---|
| 890 | 0805 set(hi,<span class="string">'CData'</span>,[1:256]') |
---|
| 891 | 0806 set(hcol,<span class="string">'YLim'</span>,[MinA MaxA]) |
---|
| 892 | 0807 YTick=<a href="#_sub5" class="code" title="subfunction YTick=colbartick(MinA,MaxA)">colbartick</a>(MinA,MaxA); |
---|
| 893 | 0808 set(hcol,<span class="string">'YTick'</span>,YTick) |
---|
| 894 | 0809 <span class="keyword">else</span> |
---|
| 895 | 0810 hi=get(hcol,<span class="string">'children'</span>); |
---|
| 896 | 0811 <span class="keyword">if</span> iscell(hi)<span class="comment">%multiple images in colorbar</span> |
---|
| 897 | 0812 hi=hi{1}; |
---|
| 898 | 0813 <span class="keyword">end</span> |
---|
| 899 | 0814 set(hi,<span class="string">'YData'</span>,[MinA MaxA]) |
---|
| 900 | 0815 set(hi,<span class="string">'CData'</span>,[1:64]') |
---|
| 901 | 0816 YTick=<a href="#_sub5" class="code" title="subfunction YTick=colbartick(MinA,MaxA)">colbartick</a>(MinA,MaxA); |
---|
| 902 | 0817 set(hcol,<span class="string">'YLim'</span>,[MinA MaxA]) |
---|
| 903 | 0818 set(hcol,<span class="string">'YTick'</span>,YTick) |
---|
| 904 | 0819 <span class="keyword">end</span> |
---|
| 905 | 0820 set(hcol,<span class="string">'Yticklabel'</span>,num2str(YTick')); |
---|
| 906 | 0821 <span class="keyword">end</span> |
---|
| 907 | 0822 <span class="keyword">elseif</span> ishandle(hcol) |
---|
| 908 | 0823 delete(hcol); <span class="comment">%erase existing colorbar if not needed</span> |
---|
| 909 | 0824 <span class="keyword">end</span> |
---|
| 910 | 0825 <span class="keyword">else</span><span class="comment">%no scalar plot</span> |
---|
| 911 | 0826 <span class="keyword">if</span> ~isempty(hima) && ishandle(hima) |
---|
| 912 | 0827 delete(hima) |
---|
| 913 | 0828 <span class="keyword">end</span> |
---|
| 914 | 0829 <span class="keyword">if</span> ~isempty(hcol)&& ishandle(hcol) |
---|
| 915 | 0830 delete(hcol) |
---|
| 916 | 0831 <span class="keyword">end</span> |
---|
| 917 | 0832 AxeData.A=[]; |
---|
| 918 | 0833 AxeData.AX=[]; |
---|
| 919 | 0834 AxeData.AY=[]; |
---|
| 920 | 0835 PlotParamOut=rmfield(PlotParamOut,<span class="string">'Scalar'</span>); |
---|
| 921 | 0836 <span class="keyword">end</span> |
---|
| 922 | 0837 |
---|
| 923 | 0838 <span class="comment">%%%%%%%%%%%%%%%%%%%%% vector plot %%%%%%%%%%%%%%%%%%%%%%%%%%</span> |
---|
| 924 | 0839 <span class="keyword">if</span> test_vec |
---|
| 925 | 0840 <span class="comment">%vector scale representation</span> |
---|
| 926 | 0841 <span class="keyword">if</span> size(vec_U,1)==numel(vec_Y) && size(vec_U,2)==numel(vec_X); <span class="comment">% x, y coordinate variables</span> |
---|
| 927 | 0842 [vec_X,vec_Y]=meshgrid(vec_X,vec_Y); |
---|
| 928 | 0843 <span class="keyword">end</span> |
---|
| 929 | 0844 vec_X=reshape(vec_X,1,numel(vec_X));<span class="comment">%reshape in matlab vectors</span> |
---|
| 930 | 0845 vec_Y=reshape(vec_Y,1,numel(vec_Y)); |
---|
| 931 | 0846 vec_U=reshape(vec_U,1,numel(vec_U)); |
---|
| 932 | 0847 vec_V=reshape(vec_V,1,numel(vec_V)); |
---|
| 933 | 0848 MinMaxX=max(vec_X)-min(vec_X); |
---|
| 934 | 0849 MinMaxY=max(vec_Y)-min(vec_Y); |
---|
| 935 | 0850 AxeData.Mesh=sqrt((MinMaxX*MinMaxY)/length(vec_X)); |
---|
| 936 | 0851 <span class="keyword">if</span> ~isfield(PlotParam.Vectors,<span class="string">'AutoVec'</span>) || isequal(PlotParam.Vectors.AutoVec,0)|| ~isfield(PlotParam.Vectors,<span class="string">'VecScale'</span>)<span class="keyword">...</span> |
---|
| 937 | 0852 ||isempty(PlotParam.Vectors.VecScale)||~isa(PlotParam.Vectors.VecScale,<span class="string">'double'</span>) <span class="comment">%automatic vector scale</span> |
---|
| 938 | 0853 scale=[]; |
---|
| 939 | 0854 <span class="keyword">if</span> test_false <span class="comment">%remove false vectors</span> |
---|
| 940 | 0855 indsel=find(AxeData.FF==0);<span class="comment">%indsel =indices of good vectors</span> |
---|
| 941 | 0856 <span class="keyword">else</span> |
---|
| 942 | 0857 indsel=[1:numel(vec_X)];<span class="comment">%</span> |
---|
| 943 | 0858 <span class="keyword">end</span> |
---|
| 944 | 0859 <span class="keyword">if</span> isempty(vec_U) |
---|
| 945 | 0860 scale=1; |
---|
| 946 | 0861 <span class="keyword">else</span> |
---|
| 947 | 0862 <span class="keyword">if</span> isempty(indsel) |
---|
| 948 | 0863 MaxU=max(abs(vec_U)); |
---|
| 949 | 0864 MaxV=max(abs(vec_V)); |
---|
| 950 | 0865 <span class="keyword">else</span> |
---|
| 951 | 0866 MaxU=max(abs(vec_U(indsel))); |
---|
| 952 | 0867 MaxV=max(abs(vec_V(indsel))); |
---|
| 953 | 0868 <span class="keyword">end</span> |
---|
| 954 | 0869 scale=MinMaxX/(max(MaxU,MaxV)*50); |
---|
| 955 | 0870 PlotParam.Vectors.VecScale=scale;<span class="comment">%update the 'scale' display</span> |
---|
| 956 | 0871 <span class="keyword">end</span> |
---|
| 957 | 0872 <span class="keyword">else</span> |
---|
| 958 | 0873 scale=PlotParam.Vectors.VecScale; <span class="comment">%impose the length of vector representation</span> |
---|
| 959 | 0874 <span class="keyword">end</span>; |
---|
| 960 | 0875 |
---|
| 961 | 0876 <span class="comment">%record vectors on the plotting axes</span> |
---|
| 962 | 0877 <span class="keyword">if</span> test_C==0 |
---|
| 963 | 0878 vec_C=ones(1,numel(vec_X)); |
---|
| 964 | 0879 <span class="keyword">end</span> |
---|
| 965 | 0880 AxeData.X=vec_X'; |
---|
| 966 | 0881 AxeData.Y=vec_Y'; |
---|
| 967 | 0882 AxeData.U=vec_U'; |
---|
| 968 | 0883 AxeData.V=vec_V'; |
---|
| 969 | 0884 AxeData.C=vec_C'; |
---|
| 970 | 0885 <span class="keyword">if</span> isfield(Data,<span class="string">'W'</span>) |
---|
| 971 | 0886 AxeData.W=Data.W; |
---|
| 972 | 0887 <span class="keyword">end</span> |
---|
| 973 | 0888 |
---|
| 974 | 0889 <span class="comment">%decimate by a factor 2 in vector mesh(4 in nbre of vectors)</span> |
---|
| 975 | 0890 <span class="keyword">if</span> isfield(PlotParam.Vectors,<span class="string">'decimate4'</span>)&isequal(PlotParam.Vectors.decimate4,1) |
---|
| 976 | 0891 diffy=diff(vec_Y); <span class="comment">%difference dy=vec_Y(i+1)-vec_Y(i)</span> |
---|
| 977 | 0892 dy_thresh=max(abs(diffy))/2; |
---|
| 978 | 0893 ind_jump=find(abs(diffy) > dy_thresh); <span class="comment">%indices with diff(vec_Y)> max/2, detect change of line</span> |
---|
| 979 | 0894 ind_sel=[1:ind_jump(1)];<span class="comment">%select the first line</span> |
---|
| 980 | 0895 <span class="keyword">for</span> i=2:2:length(ind_jump)-1 |
---|
| 981 | 0896 ind_sel=[ind_sel [ind_jump(i)+1:ind_jump(i+1)]];<span class="comment">% select the odd lines</span> |
---|
| 982 | 0897 <span class="keyword">end</span> |
---|
| 983 | 0898 nb_sel=length(ind_sel); |
---|
| 984 | 0899 ind_sel=ind_sel([1:2:nb_sel]);<span class="comment">% take half the points on a line</span> |
---|
| 985 | 0900 vec_X=vec_X(ind_sel); |
---|
| 986 | 0901 vec_Y=vec_Y(ind_sel); |
---|
| 987 | 0902 vec_U=vec_U(ind_sel); |
---|
| 988 | 0903 vec_V=vec_V(ind_sel); |
---|
| 989 | 0904 vec_C=vec_C(ind_sel); |
---|
| 990 | 0905 <span class="keyword">if</span> isfield(Data,<span class="string">'F'</span>)&& isequal(size(Data.F),size(vec_X))sdfsdf |
---|
| 991 | 0906 Data.F=Data.F(ind_sel); |
---|
| 992 | 0907 <span class="keyword">end</span> |
---|
| 993 | 0908 <span class="keyword">if</span> isfield(Data,<span class="string">'FF'</span>)&& isequal(size(Data.FF),size(vec_X)) |
---|
| 994 | 0909 Data.FF=Data.FF(ind_sel); |
---|
| 995 | 0910 <span class="keyword">end</span> |
---|
| 996 | 0911 <span class="keyword">end</span> |
---|
| 997 | 0912 |
---|
| 998 | 0913 <span class="comment">%get main level color code</span> |
---|
| 999 | 0914 [colorlist,col_vec,PlotParamOut.Vectors]=<a href="set_col_vec.html" class="code" title="function [colorlist,col_vec,colcode_out]=set_col_vec(colcode,vec_C)">set_col_vec</a>(PlotParam.Vectors,vec_C); |
---|
| 1000 | 0915 <span class="comment">% take flags into account: add flag colors to the list of colors</span> |
---|
| 1001 | 0916 sizlist=size(colorlist); |
---|
| 1002 | 0917 nbcolor=sizlist(1); |
---|
| 1003 | 0918 <span class="keyword">if</span> test_black && isfield (Data,<span class="string">'F'</span>) |
---|
| 1004 | 0919 nbcolor=nbcolor+1; |
---|
| 1005 | 0920 colorlist(nbcolor,:)=[0 0 0]; <span class="comment">%add black to the list of colors</span> |
---|
| 1006 | 0921 <span class="keyword">if</span> isfield(Data,<span class="string">'FF'</span>) |
---|
| 1007 | 0922 ind_flag=find(Data.F~=1 & Data.FF==0); <span class="comment">%flag warning but not faulse</span> |
---|
| 1008 | 0923 <span class="keyword">else</span> |
---|
| 1009 | 0924 ind_flag=find(Data.F~=1); |
---|
| 1010 | 0925 <span class="keyword">end</span> |
---|
| 1011 | 0926 col_vec(ind_flag)=nbcolor; |
---|
| 1012 | 0927 <span class="keyword">end</span> |
---|
| 1013 | 0928 nbcolor=nbcolor+1; |
---|
| 1014 | 0929 <span class="keyword">if</span> test_false && isfield (Data,<span class="string">'FF'</span>) |
---|
| 1015 | 0930 ind_flag=find(Data.FF~=0); |
---|
| 1016 | 0931 <span class="keyword">if</span> isfield(PlotParam.Vectors,<span class="string">'HideFalse'</span>) & PlotParam.Vectors.HideFalse==1 |
---|
| 1017 | 0932 colorlist(nbcolor,:)=[NaN NaN NaN];<span class="comment">% no plot of false vectors</span> |
---|
| 1018 | 0933 <span class="keyword">else</span> |
---|
| 1019 | 0934 colorlist(nbcolor,:)=[1 0 1];<span class="comment">% magenta color</span> |
---|
| 1020 | 0935 <span class="keyword">end</span> |
---|
| 1021 | 0936 col_vec(ind_flag)=nbcolor; |
---|
| 1022 | 0937 <span class="keyword">end</span> |
---|
| 1023 | 0938 <span class="comment">%plot vectors:</span> |
---|
| 1024 | 0939 <a href="#_sub4" class="code" title="subfunction quiresetn(haxes,x,y,u,v,scale,colorlist,col_vec)">quiresetn</a>(haxes,vec_X,vec_Y,vec_U,vec_V,scale,colorlist,col_vec); |
---|
| 1025 | 0940 <span class="keyword">else</span> |
---|
| 1026 | 0941 hvec=findobj(haxes,<span class="string">'Tag'</span>,<span class="string">'vel'</span>); |
---|
| 1027 | 0942 <span class="keyword">if</span> ~isempty(hvec) |
---|
| 1028 | 0943 delete(hvec); |
---|
| 1029 | 0944 <span class="keyword">end</span> |
---|
| 1030 | 0945 AxeData.X=[]; |
---|
| 1031 | 0946 AxeData.Y=[]; |
---|
| 1032 | 0947 AxeData.U=[]; |
---|
| 1033 | 0948 AxeData.V=[]; |
---|
| 1034 | 0949 AxeData.C=[]; |
---|
| 1035 | 0950 AxeData.W=[]; |
---|
| 1036 | 0951 AxeData.Mesh=[]; |
---|
| 1037 | 0952 PlotParamOut=rmfield(PlotParamOut,<span class="string">'Vectors'</span>); |
---|
| 1038 | 0953 <span class="keyword">end</span> |
---|
| 1039 | 0954 <span class="keyword">if</span> isfield(Data,<span class="string">'Z'</span>) |
---|
| 1040 | 0955 AxeData.Z=Data.Z;<span class="comment">% A REVOIR</span> |
---|
| 1041 | 0956 <span class="keyword">end</span> |
---|
| 1042 | 0957 listfields={<span class="string">'AY'</span>,<span class="string">'AX'</span>,<span class="string">'A'</span>,<span class="string">'X'</span>,<span class="string">'Y'</span>,<span class="string">'U'</span>,<span class="string">'V'</span>,<span class="string">'C'</span>,<span class="string">'W'</span>,<span class="string">'F'</span>,<span class="string">'FF'</span>}; |
---|
| 1043 | 0958 listdim={<span class="string">'AY'</span>,<span class="string">'AX'</span>,{<span class="string">'AY'</span>,<span class="string">'AX'</span>},<span class="string">'nb_vectors'</span>,<span class="string">'nb_vectors'</span>,<span class="string">'nb_vectors'</span>,<span class="string">'nb_vectors'</span>,<span class="string">'nb_vectors'</span>,<span class="string">'nb_vectors'</span>,<span class="string">'nb_vectors'</span>,<span class="string">'nb_vectors'</span>}; |
---|
| 1044 | 0959 Role={<span class="string">'coord_y'</span>,<span class="string">'coord_x'</span>,<span class="string">'scalar'</span>,<span class="string">'coord_x'</span>,<span class="string">'coord_y'</span>,<span class="string">'vector_x'</span>,<span class="string">'vector_y'</span>,<span class="string">'scalar'</span>,<span class="string">'vector_z'</span>,<span class="string">'warnflag'</span>,<span class="string">'falseflag'</span>}; |
---|
| 1045 | 0960 ind_select=[]; |
---|
| 1046 | 0961 nbvar=0; |
---|
| 1047 | 0962 |
---|
| 1048 | 0963 <span class="keyword">for</span> ilist=1:numel(listfields) |
---|
| 1049 | 0964 eval([<span class="string">'testvar=isfield(AxeData,listfields{ilist}) && ~isempty(AxeData.'</span> listfields{ilist} <span class="string">');'</span>]) |
---|
| 1050 | 0965 <span class="keyword">if</span> testvar |
---|
| 1051 | 0966 nbvar=nbvar+1; |
---|
| 1052 | 0967 AxeData.ListVarName{nbvar}=listfields{ilist}; |
---|
| 1053 | 0968 AxeData.VarDimName{nbvar}=listdim{ilist}; |
---|
| 1054 | 0969 AxeData.VarAttribute{nbvar}.Role=Role{ilist}; |
---|
| 1055 | 0970 <span class="keyword">end</span> |
---|
| 1056 | 0971 <span class="keyword">end</span> |
---|
| 1057 | 0972 |
---|
| 1058 | 0973 <span class="comment">% adjust the size of the plot to include the whole field, except if KeepLim=1</span> |
---|
| 1059 | 0974 <span class="keyword">if</span> ~(exist(<span class="string">'KeepLim'</span>,<span class="string">'var'</span>) && isequal(KeepLim,1)) <span class="comment">%adjust the graph limits*</span> |
---|
| 1060 | 0975 test_lim=0; |
---|
| 1061 | 0976 <span class="keyword">if</span> test_vec |
---|
| 1062 | 0977 Xlim=[min(vec_X) max(vec_X)]; |
---|
| 1063 | 0978 Ylim=[min(vec_Y) max(vec_Y)]; |
---|
| 1064 | 0979 test_lim=1; |
---|
| 1065 | 0980 <span class="keyword">if</span> test_ima<span class="comment">%both background image and vectors coexist, take the wider bound</span> |
---|
| 1066 | 0981 Xlim(1)=min(AX(1),Xlim(1)); |
---|
| 1067 | 0982 Xlim(2)=max(AX(end),Xlim(2)); |
---|
| 1068 | 0983 Ylim(1)=min(AY(end),Ylim(1)); |
---|
| 1069 | 0984 Ylim(2)=max(AY(1),Ylim(2)); |
---|
| 1070 | 0985 <span class="keyword">end</span> |
---|
| 1071 | 0986 <span class="keyword">elseif</span> test_ima <span class="comment">%only image plot</span> |
---|
| 1072 | 0987 Xlim(1)=min(AX(1),AX(end)); |
---|
| 1073 | 0988 Xlim(2)=max(AX(1),AX(end)); |
---|
| 1074 | 0989 Ylim(1)=min(AY(1),AY(end)); |
---|
| 1075 | 0990 Ylim(2)=max(AY(1),AY(end)); |
---|
| 1076 | 0991 test_lim=1; |
---|
| 1077 | 0992 <span class="keyword">end</span> |
---|
| 1078 | 0993 <span class="keyword">if</span> test_lim |
---|
| 1079 | 0994 set(haxes,<span class="string">'XLim'</span>,Xlim);<span class="comment">% set x limits of frame in axes coordinates</span> |
---|
| 1080 | 0995 set(haxes,<span class="string">'YLim'</span>,Ylim);<span class="comment">% set y limits of frame in axes coordinate</span> |
---|
| 1081 | 0996 <span class="keyword">end</span> |
---|
| 1082 | 0997 <span class="keyword">end</span> |
---|
| 1083 | 0998 <span class="keyword">if</span> ~(isfield(PlotParam,<span class="string">'Auto_xy'</span>) && isequal(PlotParam.Auto_xy,1)) |
---|
| 1084 | 0999 set(haxes,<span class="string">'DataAspectRatio'</span>,[1 1 1]) |
---|
| 1085 | 1000 <span class="keyword">end</span> |
---|
| 1086 | 1001 set(haxes,<span class="string">'YDir'</span>,<span class="string">'normal'</span>) |
---|
| 1087 | 1002 |
---|
| 1088 | 1003 <span class="comment">%---------------------------------------------</span> |
---|
| 1089 | 1004 <span class="comment">%function for plotting vectors</span> |
---|
| 1090 | 1005 <span class="comment">%------------------------------------------------</span> |
---|
| 1091 | 1006 <span class="comment">%INPUT:</span> |
---|
| 1092 | 1007 <span class="comment">% haxes: handles of the plotting axes</span> |
---|
| 1093 | 1008 <span class="comment">%x,y,u,v: vectors coordinates and vector components to plot, arrays withb the same dimension</span> |
---|
| 1094 | 1009 <span class="comment">% scale: scaling factor for vector length representation</span> |
---|
| 1095 | 1010 <span class="comment">%colorlist(icolor,:): list of vector colors, dim (nbcolor,3), depending on color #i</span> |
---|
| 1096 | 1011 <span class="comment">%col_vec: matlab vector setting the color number #i for each velocity vector</span> |
---|
| 1097 | 1012 <a name="_sub4" href="#_subfunctions" class="code">function quiresetn(haxes,x,y,u,v,scale,colorlist,col_vec)</a> |
---|
| 1098 | 1013 |
---|
| 1099 | 1014 <span class="comment">%define arrows</span> |
---|
| 1100 | 1015 theta=0.5 ;<span class="comment">%angle arrow</span> |
---|
| 1101 | 1016 alpha=0.3 ;<span class="comment">%length arrow</span> |
---|
| 1102 | 1017 rot=alpha*[cos(theta) -sin(theta); sin(theta) cos(theta)]'; |
---|
| 1103 | 1018 <span class="comment">%find the existing lines</span> |
---|
| 1104 | 1019 <span class="comment">%h=findobj(gca,'Type','Line');% search existing lines in the current axes</span> |
---|
| 1105 | 1020 h=findobj(haxes,<span class="string">'Tag'</span>,<span class="string">'vel'</span>);<span class="comment">% search existing lines in the current axes</span> |
---|
| 1106 | 1021 sizh=size(h); |
---|
| 1107 | 1022 set(h,<span class="string">'EraseMode'</span>,<span class="string">'xor'</span>); |
---|
| 1108 | 1023 set(haxes,<span class="string">'NextPlot'</span>,<span class="string">'replacechildren'</span>); |
---|
| 1109 | 1024 <span class="comment">%htext=findobj('Tag','scalevec');</span> |
---|
| 1110 | 1025 |
---|
| 1111 | 1026 <span class="comment">%drawnow</span> |
---|
| 1112 | 1027 <span class="comment">%create lines (if no lines) or modify them</span> |
---|
| 1113 | 1028 <span class="keyword">if</span> ~isequal(size(col_vec),size(x)) |
---|
| 1114 | 1029 col_vec=ones(size(x));<span class="comment">% case of error in col_vec input</span> |
---|
| 1115 | 1030 <span class="keyword">end</span> |
---|
| 1116 | 1031 sizlist=size(colorlist); |
---|
| 1117 | 1032 ncolor=sizlist(1); |
---|
| 1118 | 1033 |
---|
| 1119 | 1034 <span class="keyword">for</span> icolor=1:ncolor |
---|
| 1120 | 1035 <span class="comment">%determine the line positions for each color icolor</span> |
---|
| 1121 | 1036 ind=find(col_vec==icolor); |
---|
| 1122 | 1037 xc=x(ind); |
---|
| 1123 | 1038 yc=y(ind); |
---|
| 1124 | 1039 uc=u(ind)*scale; |
---|
| 1125 | 1040 vc=v(ind)*scale; |
---|
| 1126 | 1041 n=size(xc); |
---|
| 1127 | 1042 xN=NaN*ones(size(xc)); |
---|
| 1128 | 1043 matx=[xc(:) xc(:)+uc(:) xN(:)]'; |
---|
| 1129 | 1044 matx=reshape(matx,1,3*n(2)); |
---|
| 1130 | 1045 maty=[yc(:) yc(:)+vc(:) xN(:)]'; |
---|
| 1131 | 1046 maty=reshape(maty,1,3*n(2)); |
---|
| 1132 | 1047 |
---|
| 1133 | 1048 <span class="comment">%determine arrow heads</span> |
---|
| 1134 | 1049 arrowplus=rot*[uc;vc]; |
---|
| 1135 | 1050 arrowmoins=rot'*[uc;vc]; |
---|
| 1136 | 1051 x1=xc+uc-arrowplus(1,:); |
---|
| 1137 | 1052 x2=xc+uc; |
---|
| 1138 | 1053 x3=xc+uc-arrowmoins(1,:); |
---|
| 1139 | 1054 y1=yc+vc-arrowplus(2,:); |
---|
| 1140 | 1055 y2=yc+vc; |
---|
| 1141 | 1056 y3=yc+vc-arrowmoins(2,:); |
---|
| 1142 | 1057 matxar=[x1(:) x2(:) x3(:) xN(:)]'; |
---|
| 1143 | 1058 matxar=reshape(matxar,1,4*n(2)); |
---|
| 1144 | 1059 matyar=[y1(:) y2(:) y3(:) xN(:)]'; |
---|
| 1145 | 1060 matyar=reshape(matyar,1,4*n(2)); |
---|
| 1146 | 1061 <span class="comment">%draw the line or modify the existing ones</span> |
---|
| 1147 | 1062 isn=isnan(colorlist(icolor,:));<span class="comment">%test if color NaN</span> |
---|
| 1148 | 1063 <span class="keyword">if</span> 2*icolor > sizh(1) <span class="comment">%if icolor exceeds the number of existing ones</span> |
---|
| 1149 | 1064 axes(haxes) |
---|
| 1150 | 1065 <span class="keyword">if</span> ~isn(1) <span class="comment">%if the vectors are visible color not nan</span> |
---|
| 1151 | 1066 <span class="keyword">if</span> n(2)>0 |
---|
| 1152 | 1067 hold on |
---|
| 1153 | 1068 line(matx,maty,<span class="string">'Color'</span>,colorlist(icolor,:),<span class="string">'Tag'</span>,<span class="string">'vel'</span>);<span class="comment">% plot new lines</span> |
---|
| 1154 | 1069 line(matxar,matyar,<span class="string">'Color'</span>,colorlist(icolor,:),<span class="string">'Tag'</span>,<span class="string">'vel'</span>);<span class="comment">% plot arrows</span> |
---|
| 1155 | 1070 <span class="keyword">end</span> |
---|
| 1156 | 1071 <span class="keyword">end</span> |
---|
| 1157 | 1072 <span class="keyword">else</span> |
---|
| 1158 | 1073 <span class="keyword">if</span> isn(1) |
---|
| 1159 | 1074 delete(h(2*icolor-1)) |
---|
| 1160 | 1075 delete(h(2*icolor)) |
---|
| 1161 | 1076 <span class="keyword">else</span> |
---|
| 1162 | 1077 set(h(2*icolor-1),<span class="string">'Xdata'</span>,matx,<span class="string">'Ydata'</span>,maty); |
---|
| 1163 | 1078 set(h(2*icolor-1),<span class="string">'Color'</span>,colorlist(icolor,:)); |
---|
| 1164 | 1079 set(h(2*icolor-1),<span class="string">'EraseMode'</span>,<span class="string">'xor'</span>); |
---|
| 1165 | 1080 set(h(2*icolor),<span class="string">'Xdata'</span>,matxar,<span class="string">'Ydata'</span>,matyar); |
---|
| 1166 | 1081 set(h(2*icolor),<span class="string">'Color'</span>,colorlist(icolor,:)); |
---|
| 1167 | 1082 set(h(2*icolor),<span class="string">'EraseMode'</span>,<span class="string">'xor'</span>); |
---|
| 1168 | 1083 <span class="keyword">end</span> |
---|
| 1169 | 1084 <span class="keyword">end</span> |
---|
| 1170 | 1085 <span class="keyword">end</span> |
---|
| 1171 | 1086 <span class="keyword">if</span> sizh(1) > 2*ncolor |
---|
| 1172 | 1087 <span class="keyword">for</span> icolor=ncolor+1 : sizh(1)/2<span class="comment">%delete additional objects</span> |
---|
| 1173 | 1088 delete(h(2*icolor-1)) |
---|
| 1174 | 1089 delete(h(2*icolor)) |
---|
| 1175 | 1090 <span class="keyword">end</span> |
---|
| 1176 | 1091 <span class="keyword">end</span> |
---|
| 1177 | 1092 |
---|
| 1178 | 1093 <span class="comment">%---------------------------------------</span> |
---|
| 1179 | 1094 <span class="comment">%determine tick positions for colorbar</span> |
---|
| 1180 | 1095 <span class="comment">%------------------------------------</span> |
---|
| 1181 | 1096 <a name="_sub5" href="#_subfunctions" class="code">function YTick=colbartick(MinA,MaxA)</a> |
---|
| 1182 | 1097 <span class="comment">%determine tick positions with "simple" values between MinA and MaxA</span> |
---|
| 1183 | 1098 YTick=0;<span class="comment">%default</span> |
---|
| 1184 | 1099 maxabs=max([abs(MinA) abs(MaxA)]); |
---|
| 1185 | 1100 <span class="keyword">if</span> maxabs>0 |
---|
| 1186 | 1101 ord=10^(floor(log10(maxabs)));<span class="comment">%order of magnitude</span> |
---|
| 1187 | 1102 div=1; |
---|
| 1188 | 1103 siz2=1; |
---|
| 1189 | 1104 <span class="keyword">while</span> siz2<2 |
---|
| 1190 | 1105 <span class="comment">% values=[-9:div:9];</span> |
---|
| 1191 | 1106 values=-10:div:10; |
---|
| 1192 | 1107 ind=find((ord*values-MaxA)<0 & (ord*values-MinA)>0);<span class="comment">%indices of 'values' such that MinA<ord*values<MaxA</span> |
---|
| 1193 | 1108 siz=size(ind); |
---|
| 1194 | 1109 <span class="keyword">if</span> siz(2)<4<span class="comment">%if there are less than 4 selected values (4 levels)</span> |
---|
| 1195 | 1110 values=[-9:0.5*div:9]; |
---|
| 1196 | 1111 ind=find((ord*values-MaxA)<0 & (ord*values-MinA)>0); |
---|
| 1197 | 1112 <span class="keyword">end</span> |
---|
| 1198 | 1113 siz2=size(ind,2); |
---|
| 1199 | 1114 <span class="comment">% siz2=siz(2)</span> |
---|
| 1200 | 1115 div=div/10; |
---|
| 1201 | 1116 <span class="keyword">end</span> |
---|
| 1202 | 1117 YTick=ord*values(ind); |
---|
| 1203 | 1118 <span class="keyword">end</span></pre></div> |
---|
| 1204 | <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> |
---|
| 1205 | </body> |
---|
| 1206 | </html> |
---|