Mercurial > octave-nkf

comparison scripts/plot/draw/isocolors.m @ 17572:7bb76a22cde1

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
maint: Split scripts/plot directory into 4 pieces. scripts/gui : user-interface functions scripts/plot/appearance : functions controlling plot appearance scripts/plot/draw : plotting functions which produce graphs scripts/plot/util : low-level plotting functions and utilities. * scripts/gui/guidata.m, scripts/gui/guihandles.m, scripts/gui/module.mk, scripts/gui/private/__file_filter__.m, scripts/gui/private/__fltk_file_filter__.m, scripts/gui/private/__is_function__.m, scripts/gui/private/__uigetdir_fltk__.m, scripts/gui/private/__uigetfile_fltk__.m, scripts/gui/private/__uiobject_split_args__.m, scripts/gui/private/__uiputfile_fltk__.m, scripts/gui/uicontextmenu.m, scripts/gui/uicontrol.m, scripts/gui/uigetdir.m, scripts/gui/uigetfile.m, scripts/gui/uimenu.m, scripts/gui/uipanel.m, scripts/gui/uipushtool.m, scripts/gui/uiputfile.m, scripts/gui/uiresume.m, scripts/gui/uitoggletool.m, scripts/gui/uitoolbar.m, scripts/gui/uiwait.m, scripts/gui/waitbar.m, scripts/gui/waitforbuttonpress.m: Moved from scripts/plot to scripts/gui * scripts/plot/appearance/__clabel__.m, scripts/plot/appearance/__getlegenddata__.m, scripts/plot/appearance/axis.m, scripts/plot/appearance/box.m, scripts/plot/appearance/caxis.m, scripts/plot/appearance/clabel.m, scripts/plot/appearance/daspect.m, scripts/plot/appearance/diffuse.m, scripts/plot/appearance/grid.m, scripts/plot/appearance/gtext.m, scripts/plot/appearance/hidden.m, scripts/plot/appearance/legend.m, scripts/plot/appearance/orient.m, scripts/plot/appearance/pbaspect.m, scripts/plot/appearance/private/__axis_label__.m, scripts/plot/appearance/private/__axis_limits__.m, scripts/plot/appearance/shading.m, scripts/plot/appearance/specular.m, scripts/plot/appearance/text.m, scripts/plot/appearance/title.m, scripts/plot/appearance/view.m, scripts/plot/appearance/whitebg.m, scripts/plot/appearance/xlabel.m, scripts/plot/appearance/xlim.m, scripts/plot/appearance/ylabel.m, scripts/plot/appearance/ylim.m, scripts/plot/appearance/zlabel.m, scripts/plot/appearance/zlim.m: Moved from scripts/plot to subdir appearance. * scripts/plot/draw/area.m, scripts/plot/draw/bar.m, scripts/plot/draw/barh.m, scripts/plot/draw/colorbar.m, scripts/plot/draw/comet.m, scripts/plot/draw/comet3.m, scripts/plot/draw/compass.m, scripts/plot/draw/contour.m, scripts/plot/draw/contour3.m, scripts/plot/draw/contourc.m, scripts/plot/draw/contourf.m, scripts/plot/draw/cylinder.m, scripts/plot/draw/ellipsoid.m, scripts/plot/draw/errorbar.m, scripts/plot/draw/ezcontour.m, scripts/plot/draw/ezcontourf.m, scripts/plot/draw/ezmesh.m, scripts/plot/draw/ezmeshc.m, scripts/plot/draw/ezplot.m, scripts/plot/draw/ezplot3.m, scripts/plot/draw/ezpolar.m, scripts/plot/draw/ezsurf.m, scripts/plot/draw/ezsurfc.m, scripts/plot/draw/feather.m, scripts/plot/draw/fill.m, scripts/plot/draw/fplot.m, scripts/plot/draw/hist.m, scripts/plot/draw/isocolors.m, scripts/plot/draw/isonormals.m, scripts/plot/draw/isosurface.m, scripts/plot/draw/line.m, scripts/plot/draw/loglog.m, scripts/plot/draw/loglogerr.m, scripts/plot/draw/mesh.m, scripts/plot/draw/meshc.m, scripts/plot/draw/meshz.m, scripts/plot/draw/pareto.m, scripts/plot/draw/patch.m, scripts/plot/draw/pcolor.m, scripts/plot/draw/peaks.m, scripts/plot/draw/pie.m, scripts/plot/draw/pie3.m, scripts/plot/draw/plot.m, scripts/plot/draw/plot3.m, scripts/plot/draw/plotmatrix.m, scripts/plot/draw/plotyy.m, scripts/plot/draw/polar.m, scripts/plot/draw/private/__add_datasource__.m, scripts/plot/draw/private/__bar__.m, scripts/plot/draw/private/__contour__.m, scripts/plot/draw/private/__errcomm__.m, scripts/plot/draw/private/__errplot__.m, scripts/plot/draw/private/__ezplot__.m, scripts/plot/draw/private/__interp_cube__.m, scripts/plot/draw/private/__line__.m, scripts/plot/draw/private/__marching_cube__.m, scripts/plot/draw/private/__patch__.m, scripts/plot/draw/private/__pie__.m, scripts/plot/draw/private/__plt__.m, scripts/plot/draw/private/__quiver__.m, scripts/plot/draw/private/__scatter__.m, scripts/plot/draw/private/__stem__.m, scripts/plot/draw/quiver.m, scripts/plot/draw/quiver3.m, scripts/plot/draw/rectangle.m, scripts/plot/draw/ribbon.m, scripts/plot/draw/rose.m, scripts/plot/draw/scatter.m, scripts/plot/draw/scatter3.m, scripts/plot/draw/semilogx.m, scripts/plot/draw/semilogxerr.m, scripts/plot/draw/semilogy.m, scripts/plot/draw/semilogyerr.m, scripts/plot/draw/shrinkfaces.m, scripts/plot/draw/slice.m, scripts/plot/draw/sombrero.m, scripts/plot/draw/sphere.m, scripts/plot/draw/stairs.m, scripts/plot/draw/stem.m, scripts/plot/draw/stem3.m, scripts/plot/draw/stemleaf.m, scripts/plot/draw/surf.m, scripts/plot/draw/surface.m, scripts/plot/draw/surfc.m, scripts/plot/draw/surfl.m, scripts/plot/draw/surfnorm.m, scripts/plot/draw/tetramesh.m, scripts/plot/draw/trimesh.m, scripts/plot/draw/triplot.m, scripts/plot/draw/trisurf.m, scripts/plot/draw/waterfall.m: Moved from plot/ to subdir draw. * scripts/plot/util/__actual_axis_position__.m, scripts/plot/util/__default_plot_options__.m, scripts/plot/util/__gnuplot_drawnow__.m, scripts/plot/util/__next_line_color__.m, scripts/plot/util/__next_line_style__.m, scripts/plot/util/__plt_get_axis_arg__.m, scripts/plot/util/__pltopt__.m, scripts/plot/util/allchild.m, scripts/plot/util/ancestor.m, scripts/plot/util/axes.m, scripts/plot/util/cla.m, scripts/plot/util/clf.m, scripts/plot/util/close.m, scripts/plot/util/closereq.m, scripts/plot/util/colstyle.m, scripts/plot/util/copyobj.m, scripts/plot/util/figure.m, scripts/plot/util/findall.m, scripts/plot/util/findfigs.m, scripts/plot/util/findobj.m, scripts/plot/util/gca.m, scripts/plot/util/gcbf.m, scripts/plot/util/gcbo.m, scripts/plot/util/gcf.m, scripts/plot/util/gco.m, scripts/plot/util/ginput.m, scripts/plot/util/gnuplot_binary.in, scripts/plot/util/graphics_toolkit.m, scripts/plot/util/hdl2struct.m, scripts/plot/util/hggroup.m, scripts/plot/util/hold.m, scripts/plot/util/isaxes.m, scripts/plot/util/isfigure.m, scripts/plot/util/ishghandle.m, scripts/plot/util/ishold.m, scripts/plot/util/isprop.m, scripts/plot/util/linkprop.m, scripts/plot/util/meshgrid.m, scripts/plot/util/ndgrid.m, scripts/plot/util/newplot.m, scripts/plot/util/print.m, scripts/plot/util/printd.m, scripts/plot/util/private/__add_default_menu__.m, scripts/plot/util/private/__fltk_ginput__.m, scripts/plot/util/private/__fltk_print__.m, scripts/plot/util/private/__ghostscript__.m, scripts/plot/util/private/__gnuplot_get_var__.m, scripts/plot/util/private/__gnuplot_ginput__.m, scripts/plot/util/private/__gnuplot_has_feature__.m, scripts/plot/util/private/__gnuplot_has_terminal__.m, scripts/plot/util/private/__gnuplot_open_stream__.m, scripts/plot/util/private/__gnuplot_print__.m, scripts/plot/util/private/__gnuplot_version__.m, scripts/plot/util/private/__go_draw_axes__.m, scripts/plot/util/private/__go_draw_figure__.m, scripts/plot/util/private/__print_parse_opts__.m, scripts/plot/util/private/__tight_eps_bbox__.m, scripts/plot/util/refresh.m, scripts/plot/util/refreshdata.m, scripts/plot/util/saveas.m, scripts/plot/util/shg.m, scripts/plot/util/struct2hdl.m, scripts/plot/util/subplot.m: Moved from plot to subdir util. * etc/HACKING: Updated directory structure info. * scripts/Makefile.am, scripts/plot/appearance/module.mk, scripts/plot/draw/module.mk, scripts/plot/util/module.mk: Added new directories to build system.
author Rik <rik@octave.org>
date Fri, 04 Oct 2013 17:09:08 -0700
parents scripts/plot/isocolors.m@eaab03308c0b
children d63878346099
comparison
equal deleted inserted replaced
17571:6e4ea5c8a4bb 17572:7bb76a22cde1
1 ## Copyright (C) 2009-2012 Martin Helm
2 ##
3 ## This file is part of Octave.
4 ##
5 ## Octave is free software; you can redistribute it and/or modify it
6 ## under the terms of the GNU General Public License as published by
7 ## the Free Software Foundation; either version 3 of the License, or (at
8 ## your option) any later version.
9 ##
10 ## Octave is distributed in the hope that it will be useful, but
11 ## WITHOUT ANY WARRANTY; without even the implied warranty of
12 ## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
13 ## General Public License for more details.
14 ##
15 ## You should have received a copy of the GNU General Public License
16 ## along with Octave; see the file COPYING. If not, see
17 ## <http://www.gnu.org/licenses/>.
18
19 ## -*- texinfo -*-
20 ## @deftypefn {Function File} {[@var{cd}] =} isocolors (@var{c}, @var{v})
21 ## @deftypefnx {Function File} {[@var{cd}] =} isocolors (@var{x}, @var{y}, @var{z}, @var{c}, @var{v})
22 ## @deftypefnx {Function File} {[@var{cd}] =} isocolors (@var{x}, @var{y}, @var{z}, @var{r}, @var{g}, @var{b}, @var{v})
23 ## @deftypefnx {Function File} {[@var{cd}] =} isocolors (@var{r}, @var{g}, @var{b}, @var{v})
24 ## @deftypefnx {Function File} {[@var{cd}] =} isocolors (@dots{}, @var{p})
25 ## @deftypefnx {Function File} {} isocolors (@dots{})
26 ##
27 ## If called with one output argument and the first input argument
28 ## @var{c} is a three-dimensional array that contains color values and
29 ## the second input argument @var{v} keeps the vertices of a geometry
30 ## then return a matrix @var{cd} with color data information for the
31 ## geometry at computed points
32 ## @command{[x, y, z] = meshgrid (1:l, 1:m, 1:n)}. The output argument
33 ## @var{cd} can be taken to manually set FaceVertexCData of a patch.
34 ##
35 ## If called with further input arguments @var{x}, @var{y} and @var{z}
36 ## which are three--dimensional arrays of the same size than @var{c}
37 ## then the color data is taken at those given points. Instead of the
38 ## color data @var{c} this function can also be called with RGB values
39 ## @var{r}, @var{g}, @var{b}. If input argumnets @var{x}, @var{y},
40 ## @var{z} are not given then again @command{meshgrid} computed values
41 ## are taken.
42 ##
43 ## Optionally, the patch handle @var{p} can be given as the last input
44 ## argument to all variations of function calls instead of the vertices
45 ## data @var{v}. Finally, if no output argument is given then directly
46 ## change the colors of a patch that is given by the patch handle
47 ## @var{p}.
48 ##
49 ## For example:
50 ##
51 ## @example
52 ## function [] = isofinish (p)
53 ## set (gca, "PlotBoxAspectRatioMode", "manual", ...
54 ## "PlotBoxAspectRatio", [1 1 1]);
55 ## set (p, "FaceColor", "interp");
56 ## ## set (p, "FaceLighting", "flat");
57 ## ## light ("Position", [1 1 5]); ## Available with JHandles
58 ## endfunction
59 ##
60 ## N = 15; # Increase number of vertices in each direction
61 ## iso = .4; # Change isovalue to .1 to display a sphere
62 ## lin = linspace (0, 2, N);
63 ## [x, y, z] = meshgrid (lin, lin, lin);
64 ## c = abs ((x-.5).^2 + (y-.5).^2 + (z-.5).^2);
65 ## figure (); # Open another figure window
66 ##
67 ## subplot (2,2,1); view (-38, 20);
68 ## [f, v] = isosurface (x, y, z, c, iso);
69 ## p = patch ("Faces", f, "Vertices", v, "EdgeColor", "none");
70 ## cdat = rand (size (c)); # Compute random patch color data
71 ## isocolors (x, y, z, cdat, p); # Directly set colors of patch
72 ## isofinish (p); # Call user function isofinish
73 ##
74 ## subplot (2,2,2); view (-38, 20);
75 ## p = patch ("Faces", f, "Vertices", v, "EdgeColor", "none");
76 ## [r, g, b] = meshgrid (lin, 2-lin, 2-lin);
77 ## cdat = isocolors (x, y, z, c, v); # Compute color data vertices
78 ## set (p, "FaceVertexCData", cdat); # Set color data manually
79 ## isofinish (p);
80 ##
81 ## subplot (2,2,3); view (-38, 20);
82 ## p = patch ("Faces", f, "Vertices", v, "EdgeColor", "none");
83 ## cdat = isocolors (r, g, b, c, p); # Compute color data patch
84 ## set (p, "FaceVertexCData", cdat); # Set color data manually
85 ## isofinish (p);
86 ##
87 ## subplot (2,2,4); view (-38, 20);
88 ## p = patch ("Faces", f, "Vertices", v, "EdgeColor", "none");
89 ## r = g = b = repmat ([1:N] / N, [N, 1, N]); # Black to white
90 ## cdat = isocolors (x, y, z, r, g, b, v);
91 ## set (p, "FaceVertexCData", cdat);
92 ## isofinish (p);
93 ## @end example
94 ##
95 ## @seealso{isosurface, isonormals}
96 ## @end deftypefn
97
98 ## Author: Martin Helm <martin@mhelm.de>
99
100 function varargout = isocolors (varargin)
101 calc_rgb = false;
102 switch (nargin)
103 case 2
104 c = varargin{1};
105 vp = varargin{2};
106 x = 1:size (c, 2);
107 y = 1:size (c, 1);
108 z = 1:size (c, 3);
109 case 4
110 calc_rgb = true;
111 R = varargin{1};
112 G = varargin{2};
113 B = varargin{3};
114 vp = varargin{4};
115 x = 1:size (R, 1);
116 y = 1:size (R, 2);
117 z = 1:size (R, 3);
118 case 5
119 x = varargin{1};
120 y = varargin{2};
121 z = varargin{3};
122 c = varargin{4};
123 vp = varargin{5};
124 case 7
125 calc_rgb = true;
126 x = varargin{1};
127 y = varargin{2};
128 z = varargin{3};
129 R = varargin{4};
130 G = varargin{5};
131 B = varargin{6};
132 vp = varargin{7};
133 otherwise
134 print_usage ();
135 endswitch
136 if (ismatrix (vp) && columns (vp) == 3)
137 pa = [];
138 v = vp;
139 elseif ( ishandle (vp) )
140 pa = vp;
141 v = get (pa, "Vertices");
142 else
143 error ("isocolors: last argument is not a vertex list or patch handle");
144 endif
145 if (calc_rgb)
146 new_col = zeros (rows (v), 3);
147 new_col(:,1) = __interp_cube__ (x, y, z, R, v, "values" );
148 new_col(:,2) = __interp_cube__ (x, y, z, G, v, "values" );
149 new_col(:,3) = __interp_cube__ (x, y, z, B, v, "values" );
150 else
151 new_col = __interp_cube__ (x, y, z, c, v, "values" );
152 endif
153 switch (nargout)
154 case 0
155 if (!isempty (pa))
156 set (pa, "FaceVertexCData", new_col);
157 endif
158 case 1
159 varargout = {new_col};
160 otherwise
161 print_usage ();
162 endswitch
163 endfunction
164
165
166 %!test
167 %! [x, y, z] = meshgrid (0:.5:2, 0:.5:2, 0:.5:2);
168 %! c = (x-.5).^2 + (y-.5).^2 + (z-.5).^2;
169 %! [f, v] = isosurface (x, y, z, c, .4);
170 %! cdat = isocolors (x, y, z, c, v);
171 %! assert (rows (cdat) == rows (v));
172 ## Can't create a patch handle for tests without a figure
173