|
6405
|
1 ## Copyright (C) 2005 John W. Eaton |
|
|
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 2, or (at your option) |
|
|
8 ## 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, write to the Free |
|
|
17 ## Software Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA |
|
|
18 ## 02110-1301, USA. |
|
|
19 |
|
|
20 ## -*- texinfo -*- |
|
|
21 ## @deftypefn {Function File} {} __go_draw_axes__ (@var{axis_obj}, @var{plot_stream}) |
|
|
22 ## Display the axes @var{axis_obj} on @var{plot_stream}. |
|
|
23 ## @end deftypefn |
|
|
24 |
|
|
25 ## Author: jwe |
|
|
26 |
|
|
27 function __go_draw_axes__ (h, plot_stream) |
|
|
28 |
|
|
29 if (nargin == 2) |
|
|
30 |
|
|
31 axis_obj = get (h); |
|
|
32 |
|
|
33 parent_figure_obj = get (axis_obj.parent); |
|
|
34 |
|
|
35 have_newer_gnuplot = compare_versions (__gnuplot_version__ (), "4.0", ">"); |
|
|
36 |
|
|
37 ## Set axis properties here? |
|
|
38 |
|
|
39 if (! isempty (axis_obj.outerposition)) |
|
|
40 pos = axis_obj.outerposition; |
|
|
41 fprintf (plot_stream, "set origin %g, %g;\n", pos(1), pos(2)); |
|
|
42 fprintf (plot_stream, "set size %g, %g;\n", pos(3), pos(4)); |
|
|
43 endif |
|
|
44 |
|
|
45 if (! isempty (axis_obj.position)) |
|
|
46 pos = axis_obj.position; |
|
|
47 fprintf (plot_stream, "set origin %g, %g;\n", pos(1), pos(2)); |
|
|
48 fprintf (plot_stream, "set size %g, %g;\n", pos(3), pos(4)); |
|
|
49 endif |
|
|
50 |
|
|
51 if (strcmp (axis_obj.dataaspectratiomode, "manual")) |
|
|
52 r = axis_obj.dataaspectratio; |
|
|
53 fprintf (plot_stream, "set size ratio %g;\n", -r(2)/r(1)); |
|
|
54 else |
|
|
55 fputs (plot_stream, "set size noratio;\n"); |
|
|
56 endif |
|
|
57 |
|
|
58 if (! isempty (axis_obj.title)) |
|
|
59 t = get (axis_obj.title); |
|
|
60 if (isempty (t.string)) |
|
|
61 fputs (plot_stream, "unset title;\n"); |
|
|
62 else |
|
|
63 fprintf (plot_stream, "set title \"%s\";\n", |
|
|
64 undo_string_escapes (t.string)); |
|
|
65 endif |
|
|
66 endif |
|
|
67 |
|
|
68 if (! isempty (axis_obj.xlabel)) |
|
|
69 t = get (axis_obj.xlabel); |
|
|
70 if (isempty (t.string)) |
|
|
71 fputs (plot_stream, "unset xlabel;\n"); |
|
|
72 else |
|
|
73 fprintf (plot_stream, "set xlabel \"%s\";\n", |
|
|
74 undo_string_escapes (t.string)); |
|
|
75 endif |
|
|
76 endif |
|
|
77 |
|
|
78 if (! isempty (axis_obj.ylabel)) |
|
|
79 t = get (axis_obj.ylabel); |
|
|
80 if (isempty (t.string)) |
|
|
81 fputs (plot_stream, "unset ylabel;\n"); |
|
|
82 else |
|
|
83 fprintf (plot_stream, "set ylabel \"%s\";\n", |
|
|
84 undo_string_escapes (t.string)); |
|
|
85 endif |
|
|
86 endif |
|
|
87 |
|
|
88 if (! isempty (axis_obj.zlabel)) |
|
|
89 t = get (axis_obj.zlabel); |
|
|
90 if (isempty (t.string)) |
|
|
91 fputs (plot_stream, "unset zlabel;\n"); |
|
|
92 else |
|
|
93 fprintf (plot_stream, "set zlabel \"%s\";\n", |
|
|
94 undo_string_escapes (t.string)); |
|
|
95 endif |
|
|
96 endif |
|
|
97 |
|
|
98 if (strcmp (axis_obj.xgrid, "on")) |
|
|
99 fputs (plot_stream, "set grid xtics;\n"); |
|
|
100 else |
|
|
101 fputs (plot_stream, "set grid noxtics;\n"); |
|
|
102 endif |
|
|
103 |
|
|
104 if (strcmp (axis_obj.ygrid, "on")) |
|
|
105 fputs (plot_stream, "set grid ytics;\n"); |
|
|
106 else |
|
|
107 fputs (plot_stream, "set grid noytics;\n"); |
|
|
108 endif |
|
|
109 |
|
|
110 if (strcmp (axis_obj.zgrid, "on")) |
|
|
111 fputs (plot_stream, "set grid ztics;\n"); |
|
|
112 else |
|
|
113 fputs (plot_stream, "set grid ztics;\n"); |
|
|
114 endif |
|
|
115 |
|
|
116 if (strcmp (axis_obj.xminorgrid, "on")) |
|
|
117 fputs (plot_stream, "set mxtics 5;\n"); |
|
|
118 fputs (plot_stream, "set grid mxtics;\n"); |
|
|
119 else |
|
|
120 fputs (plot_stream, "set grid nomxtics;\n"); |
|
|
121 endif |
|
|
122 |
|
|
123 if (strcmp (axis_obj.yminorgrid, "on")) |
|
|
124 fputs (plot_stream, "set mytics 5;\n"); |
|
|
125 fputs (plot_stream, "set grid mytics;\n"); |
|
|
126 else |
|
|
127 fputs (plot_stream, "set grid nomytics;\n"); |
|
|
128 endif |
|
|
129 |
|
|
130 if (strcmp (axis_obj.zminorgrid, "on")) |
|
|
131 fputs (plot_stream, "set mztics 5;\n"); |
|
|
132 fputs (plot_stream, "set grid mztics;\n"); |
|
|
133 else |
|
|
134 fputs (plot_stream, "set grid nomztics;\n"); |
|
|
135 endif |
|
|
136 |
|
|
137 if (strcmp (axis_obj.xtickmode, "manual")) |
|
|
138 xtic = axis_obj.xtick; |
|
|
139 if (isempty (xtic)) |
|
|
140 fputs (plot_stream, "unset xtics;\n"); |
|
|
141 else |
|
|
142 ## FIXME |
|
|
143 endif |
|
|
144 else |
|
|
145 fputs (plot_stream, "set xtics;\n"); |
|
|
146 endif |
|
|
147 |
|
|
148 if (strcmp (axis_obj.ytickmode, "manual")) |
|
|
149 ytic = axis_obj.ytick; |
|
|
150 if (isempty (ytic)) |
|
|
151 fputs (plot_stream, "unset ytics;\n"); |
|
|
152 else |
|
|
153 ## FIXME |
|
|
154 endif |
|
|
155 else |
|
|
156 fputs (plot_stream, "set ytics;\n"); |
|
|
157 endif |
|
|
158 |
|
|
159 if (strcmp (axis_obj.ztickmode, "manual")) |
|
|
160 ztic = axis_obj.ztick; |
|
|
161 if (isempty (ztic)) |
|
|
162 fputs (plot_stream, "unset ztics;\n"); |
|
|
163 else |
|
|
164 ## FIXME |
|
|
165 endif |
|
|
166 else |
|
|
167 fputs (plot_stream, "set ztics;\n"); |
|
|
168 endif |
|
|
169 |
|
|
170 if (strcmp (axis_obj.xticklabelmode, "manual")) |
|
|
171 ## FIXME -- we should be able to specify the actual tick labels, |
|
|
172 ## not just the format. |
|
|
173 xticlabel = axis_obj.xticklabel; |
|
|
174 fprintf (plot_stream, "set format x \"%s\";\n", xticlabel); |
|
|
175 else |
|
|
176 fputs (plot_stream, "set xtics;\n"); |
|
|
177 endif |
|
|
178 |
|
|
179 if (strcmp (axis_obj.yticklabelmode, "manual")) |
|
|
180 ## FIXME -- we should be able to specify the actual tick labels, |
|
|
181 ## not just the format. |
|
|
182 yticlabel = axis_obj.yticklabel; |
|
|
183 fprintf (plot_stream, "set format y \"%s\";\n", yticlabel); |
|
|
184 else |
|
|
185 fputs (plot_stream, "set ytics;\n"); |
|
|
186 endif |
|
|
187 |
|
|
188 if (strcmp (axis_obj.zticklabelmode, "manual")) |
|
|
189 ## FIXME -- we should be able to specify the actual tick labels, |
|
|
190 ## not just the format. |
|
|
191 zticlabel = axis_obj.zticklabel; |
|
|
192 fprintf (plot_stream, "set format z \"%s\";\n", zticlabel); |
|
|
193 else |
|
|
194 fputs (plot_stream, "set ztics;\n"); |
|
|
195 endif |
|
|
196 |
|
|
197 xlogscale = strcmp (axis_obj.xscale, "log"); |
|
|
198 if (xlogscale) |
|
|
199 fputs (plot_stream, "set logscale x;\n"); |
|
|
200 else |
|
|
201 fputs (plot_stream, "unset logscale x;\n"); |
|
|
202 endif |
|
|
203 |
|
|
204 ylogscale = strcmp (axis_obj.yscale, "log"); |
|
|
205 if (ylogscale) |
|
|
206 fputs (plot_stream, "set logscale y;\n"); |
|
|
207 else |
|
|
208 fputs (plot_stream, "unset logscale y;\n"); |
|
|
209 endif |
|
|
210 |
|
|
211 zlogscale = strcmp (axis_obj.zscale, "log"); |
|
|
212 if (zlogscale) |
|
|
213 fputs (plot_stream, "set logscale z;\n"); |
|
|
214 else |
|
|
215 fputs (plot_stream, "unset logscale z;\n"); |
|
|
216 endif |
|
|
217 |
|
|
218 xautoscale = strcmp (axis_obj.xlimmode, "auto"); |
|
|
219 yautoscale = strcmp (axis_obj.ylimmode, "auto"); |
|
|
220 zautoscale = strcmp (axis_obj.zlimmode, "auto"); |
|
|
221 |
|
|
222 kids = axis_obj.children; |
|
|
223 |
|
|
224 nd = 0; |
|
|
225 data_idx = 0; |
|
|
226 data = cell (); |
|
|
227 |
|
|
228 xminp = yminp = zminp = Inf; |
|
|
229 xmax = ymax = zmax = -Inf; |
|
|
230 xmin = ymin = zmin = Inf; |
|
|
231 |
|
|
232 palette_set = 0; |
|
|
233 |
|
|
234 [view_cmd, view_fcn, view_zoom] = image_viewer (); |
|
|
235 use_gnuplot_for_images = (ischar (view_fcn) |
|
|
236 && strcmp (view_fcn, "gnuplot_internal")); |
|
|
237 |
|
|
238 ximg_data = {}; |
|
|
239 ximg_data_idx = 0; |
|
|
240 |
|
|
241 for i = 1:length (kids) |
|
|
242 |
|
|
243 obj = get (kids(i)); |
|
|
244 |
|
|
245 switch (obj.type) |
|
|
246 case "image" |
|
|
247 ## FIXME - Is there a better way to determine if the plot |
|
|
248 ## command should be "plot" or "splot"????? Could have images |
|
|
249 ## projected into 3D so there is really no reason to limit |
|
|
250 ## this. |
|
|
251 if (nd == 0) |
|
|
252 nd = 2; |
|
|
253 endif |
|
|
254 |
|
|
255 img_data = obj.cdata; |
|
|
256 img_colormap = parent_figure_obj.colormap; |
|
|
257 img_xdata = obj.xdata; |
|
|
258 img_ydata = obj.ydata; |
|
|
259 |
|
|
260 if (use_gnuplot_for_images) |
|
|
261 |
|
|
262 data_idx++; |
|
|
263 |
|
|
264 [y_dim, x_dim] = size (img_data(:,:,1)); |
|
|
265 if (x_dim > 1) |
|
|
266 dx = abs (img_xdata(2)-img_xdata(1))/(x_dim-1); |
|
|
267 else |
|
|
268 dx = 1; |
|
|
269 endif |
|
|
270 if (y_dim > 1) |
|
|
271 dy = abs (img_ydata(2)-img_ydata(1))/(y_dim-1); |
|
|
272 else |
|
|
273 dy = 1; |
|
|
274 endif |
|
|
275 x_origin = min (img_xdata); |
|
|
276 y_origin = min (img_ydata); |
|
|
277 |
|
|
278 if (nd == 2) |
|
|
279 if (xautoscale) |
|
|
280 xmin = min (xmin, min (img_xdata) - dx/2); |
|
|
281 xmax = max (xmax, max (img_xdata) + dx/2); |
|
|
282 xminp = min (xminp, min (img_xdata((img_xdata - dx/2)>0)) - dx/2); |
|
|
283 endif |
|
|
284 if (yautoscale) |
|
|
285 ymin = min (ymin, min (img_ydata) - dy/2); |
|
|
286 ymax = max (ymax, max (img_ydata) + dy/2); |
|
|
287 yminp = min (yminp, min (img_ydata((img_ydata - dy/2)>0)) - dy/2); |
|
|
288 endif |
|
|
289 else |
|
|
290 ## Can have images in 3D, but the image routines don't seem |
|
|
291 ## to have a means of arbitrary projection. |
|
|
292 endif |
|
|
293 |
|
|
294 ## Let the file be deleted when Octave exits or `purge_tmp_files' |
|
|
295 ## is called. |
|
|
296 [img_fid, img_fname] = mkstemp (fullfile (P_tmpdir, "gpimageXXXXXX"), 1); |
|
|
297 if (ndims (img_data) == 3) |
|
|
298 fwrite (img_fid, permute (img_data, [3, 1, 2])(:), "float"); |
|
|
299 format = "1:2:3"; |
|
|
300 imagetype = "rgbimage"; |
|
|
301 else |
|
|
302 fwrite (img_fid, img_data(:), "float"); |
|
|
303 format = "1"; |
|
|
304 imagetype = "image"; |
|
|
305 ## Only need to set pallete once because it doesn't change |
|
|
306 ## on a figure. |
|
|
307 if (! palette_set) |
|
|
308 palette_set = 1; |
|
|
309 palette_size = rows (img_colormap); |
|
|
310 fprintf (plot_stream, |
|
|
311 "set palette positive color model RGB maxcolors %i;\n", |
|
|
312 palette_size); |
|
|
313 if (palette_size <= 128) |
|
|
314 ## Break up command to avoid buffer overflow. |
|
|
315 fprintf (plot_stream, "set palette file \"-\" using 1:2:3:4;\n"); |
|
|
316 for i = 1:palette_size |
|
|
317 fprintf (plot_stream, "%g %g %g %g;\n", |
|
|
318 1e-3*round (1e3*[(i-1)/(palette_size-1), img_colormap(i,:)])); |
|
|
319 endfor |
|
|
320 fprintf (plot_stream, "e;\n"); |
|
|
321 else |
|
|
322 ## Let the file be deleted when Octave exits or |
|
|
323 ## `purge_tmp_files' is called. |
|
|
324 [cmap_fid, cmap_fname, msg] = mkstemp (fullfile (P_tmpdir, "gpimageXXXXXX"), 1); |
|
|
325 fwrite (cmap_fid, img_colormap', "float32", 0, "ieee-le"); |
|
|
326 fclose (cmap_fid); |
|
|
327 fprintf (plot_stream, |
|
|
328 "set palette file \"%s\" binary record=%d using 1:2:3;\n", |
|
|
329 cmap_fname, palette_size); |
|
|
330 endif |
|
|
331 endif |
|
|
332 endif |
|
|
333 fclose (img_fid); |
|
|
334 |
|
|
335 filespec{data_idx} = img_fname; |
|
|
336 titlespec{data_idx} = ""; |
|
|
337 usingclause{data_idx} = sprintf ("binary array=%dx%d scan=yx flipy origin=(%g,%g) dx=%g dy=%g using %s", |
|
|
338 x_dim, y_dim, x_origin, y_origin, dx, dy, format); |
|
|
339 withclause{data_idx} = sprintf ("with %s", imagetype); |
|
|
340 |
|
|
341 ## Data in file, set to zero for data available test to pass |
|
|
342 ## below. |
|
|
343 data{data_idx} = 0; |
|
|
344 |
|
|
345 else |
|
|
346 ximg_data{++ximg_data_idx} = img_data; |
|
|
347 endif |
|
|
348 |
|
|
349 case "line" |
|
|
350 data_idx++; |
|
|
351 filespec{data_idx} = '-'; |
|
|
352 if (isempty (obj.keylabel)) |
|
|
353 titlespec{data_idx} = "title \"\""; |
|
|
354 else |
|
|
355 titlespec{data_idx} = strcat ("title \"", obj.keylabel, "\""); |
|
|
356 endif |
|
|
357 style = do_linestyle_command (obj, data_idx, plot_stream); |
|
|
358 usingclause{data_idx} = ""; |
|
|
359 withclause{data_idx} = sprintf ("with %s linestyle %d", |
|
|
360 style, data_idx); |
|
|
361 parametric(i) = true; |
|
|
362 if (! isempty (obj.zdata)) |
|
|
363 nd = 3; |
|
|
364 xdat = obj.xdata(:); |
|
|
365 ydat = obj.ydata(:); |
|
|
366 zdat = obj.zdata(:); |
|
|
367 if (xautoscale) |
|
|
368 xmin = min (xmin, min (xdat)); |
|
|
369 xmax = max (xmax, max (xdat)); |
|
|
370 xminp = min (xminp, min (xdat(xdat>0))); |
|
|
371 endif |
|
|
372 if (yautoscale) |
|
|
373 ymin = min (ymin, min (ydat)); |
|
|
374 ymax = max (ymax, max (ydat)); |
|
|
375 yminp = min (yminp, min (ydat(ydat>0))); |
|
|
376 endif |
|
|
377 if (zautoscale) |
|
|
378 zmin = min (zmin, min (zdat)); |
|
|
379 zmax = max (zmax, max (zdat)); |
|
|
380 zminp = min (zminp, min (zdat(zdat>0))); |
|
|
381 endif |
|
|
382 data{data_idx} = [xdat, ydat, zdat]'; |
|
|
383 usingclause{data_idx} = "using ($1):($2):($3)"; |
|
|
384 fputs (plot_stream, "set parametric;\n"); |
|
|
385 fputs (plot_stream, "unset hidden3d;\n"); |
|
|
386 fputs (plot_stream, "set style data lines;\n"); |
|
|
387 fputs (plot_stream, "set surface;\n"); |
|
|
388 fputs (plot_stream, "unset contour;\n"); |
|
|
389 else |
|
|
390 nd = 2; |
|
|
391 xdat = obj.xdata(:); |
|
|
392 ydat = obj.ydata(:); |
|
|
393 ldat = obj.ldata; |
|
|
394 yerr = xerr = false; |
|
|
395 if (! isempty (ldat)) |
|
|
396 yerr = true; |
|
|
397 ldat = ldat(:); |
|
|
398 endif |
|
|
399 udat = obj.udata; |
|
|
400 if (! isempty (udat)) |
|
|
401 udat = udat(:); |
|
|
402 endif |
|
|
403 xldat = obj.xldata; |
|
|
404 if (! isempty (xldat)) |
|
|
405 xerr = true; |
|
|
406 xldat = xldat(:); |
|
|
407 endif |
|
|
408 xudat = obj.xudata; |
|
|
409 if (! isempty (xudat)) |
|
|
410 xudat = xudat(:); |
|
|
411 endif |
|
|
412 if (yerr) |
|
|
413 ylo = ydat-ldat; |
|
|
414 yhi = ydat+udat; |
|
|
415 if (yautoscale) |
|
|
416 ty = [ydat; ylo; yhi]; |
|
|
417 ymin = min (ymin, min (ty)); |
|
|
418 ymax = max (ymax, max (ty)); |
|
|
419 yminp = min (yminp, min (ty(ty>0))); |
|
|
420 endif |
|
|
421 if (xerr) |
|
|
422 xlo = xdat-xldat; |
|
|
423 xhi = xdat+xudat; |
|
|
424 if (xautoscale) |
|
|
425 tx = [xdat; xlo; xhi]; |
|
|
426 xmin = min (xmin, min (tx)); |
|
|
427 xmax = max (xmax, max (tx)); |
|
|
428 xminp = min (xminp, min (tx(tx>0))); |
|
|
429 endif |
|
|
430 data{data_idx} = [xdat, ydat, xlo, xhi, ylo, yhi]'; |
|
|
431 usingclause{data_idx} = "using ($1):($2):($3):($4):($5):($6)"; |
|
|
432 withclause{data_idx} = sprintf ("with xyerrorbars linestyle %d", |
|
|
433 data_idx); |
|
|
434 else |
|
|
435 if (xautoscale) |
|
|
436 xmin = min (xmin, min (xdat)); |
|
|
437 xmax = max (xmax, max (xdat)); |
|
|
438 xminp = min (xminp, min (tx(tx>0))); |
|
|
439 endif |
|
|
440 data{data_idx} = [xdat, ydat, ylo, yhi]'; |
|
|
441 usingclause{data_idx} = "using ($1):($2):($3):($4)"; |
|
|
442 withclause{data_idx} = sprintf ("with yerrorbars linestyle %d", |
|
|
443 data_idx); |
|
|
444 endif |
|
|
445 elseif (xerr) |
|
|
446 xlo = xdat-xldat; |
|
|
447 xhi = xdat+xudat; |
|
|
448 if (xautoscale) |
|
|
449 tx = [xdat; xlo; xhi]; |
|
|
450 xmin = min (xmin, min (tx)); |
|
|
451 xmax = max (xmax, max (tx)); |
|
|
452 xminp = min (xminp, min (tx(tx>0))); |
|
|
453 endif |
|
|
454 if (yautoscale) |
|
|
455 ymin = min (ymin, min (ydat)); |
|
|
456 ymax = max (ymax, max (ydat)); |
|
|
457 yminp = min (yminp, min (ty(ty>0))); |
|
|
458 endif |
|
|
459 data{data_idx} = [xdat, ydat, xlo, xhi]'; |
|
|
460 usingclause{data_idx} = "using ($1):($2):($3):($4)"; |
|
|
461 withclause{data_idx} = sprintf ("with xerrorbars linestyle %d", |
|
|
462 data_idx); |
|
|
463 else |
|
|
464 if (xautoscale) |
|
|
465 xmin = min (xmin, min (xdat)); |
|
|
466 xmax = max (xmax, max (xdat)); |
|
|
467 xminp = min (xminp, min (xdat(xdat>0))); |
|
|
468 endif |
|
|
469 if (yautoscale) |
|
|
470 ymin = min (ymin, min (ydat)); |
|
|
471 ymax = max (ymax, max (ydat)); |
|
|
472 yminp = min (yminp, min (ydat(ydat>0))); |
|
|
473 endif |
|
|
474 data{data_idx} = [xdat, ydat]'; |
|
|
475 usingclause{data_idx} = "using ($1):($2)"; |
|
|
476 endif |
|
|
477 endif |
|
|
478 |
|
|
479 case "surface" |
|
|
480 data_idx++; |
|
|
481 style = do_linestyle_command (obj, data_idx, plot_stream); |
|
|
482 filespec{data_idx} = '-'; |
|
|
483 if (isempty (obj.keylabel)) |
|
|
484 titlespec{data_idx} = "title \"\""; |
|
|
485 else |
|
|
486 titlespec{data_idx} = strcat ("title \"", obj.keylabel, "\""); |
|
|
487 endif |
|
|
488 usingclause{data_idx} = ""; |
|
|
489 withclause{data_idx} = sprintf ("with %s linestyle %d", |
|
|
490 style, data_idx); |
|
|
491 parametric(i) = false; |
|
|
492 nd = 3; |
|
|
493 xdat = obj.xdata; |
|
|
494 ydat = obj.ydata; |
|
|
495 zdat = obj.zdata; |
|
|
496 if (xautoscale) |
|
|
497 tx = xdat(:); |
|
|
498 xmin = min (xmin, min (tx)); |
|
|
499 xmax = max (xmax, max (tx)); |
|
|
500 xminp = min (xminp, min (tx(tx>0))); |
|
|
501 endif |
|
|
502 if (yautoscale) |
|
|
503 ty = ydat(:); |
|
|
504 ymin = min (ymin, min (ty)); |
|
|
505 ymax = max (ymax, max (ty)); |
|
|
506 yminp = min (yminp, min (ty(ty>0))); |
|
|
507 endif |
|
|
508 if (zautoscale) |
|
|
509 tz = zdat(:); |
|
|
510 zmin = min (zmin, min (tz)); |
|
|
511 zmax = max (zmax, max (tz)); |
|
|
512 zminp = min (zminp, min (tz(tz>0))); |
|
|
513 endif |
|
|
514 err = false; |
|
|
515 if (isvector (xdat) && isvector (ydat) && ismatrix (zdat)) |
|
|
516 if (rows (zdat) == length (ydat) && columns (zdat) == length (xdat)) |
|
|
517 [xdat, ydat] = meshgrid (xdat, ydat); |
|
|
518 else |
|
|
519 err = true; |
|
|
520 endif |
|
|
521 elseif (ismatrix (xdat) && ismatrix (ydat) && ismatrix (zdat)) |
|
|
522 if (! (size_equal (xdat, ydat) && size_equal (xdat, zdat))) |
|
|
523 err = true; |
|
|
524 endif |
|
|
525 else |
|
|
526 err = true; |
|
|
527 endif |
|
|
528 if (err) |
|
|
529 error ("__go_draw_axes__: invalid grid data"); |
|
|
530 endif |
|
|
531 xlen = columns (zdat); |
|
|
532 ylen = rows (zdat); |
|
|
533 if (xlen == columns (xdat) && xlen == columns (ydat) |
|
|
534 && ylen == rows (xdat) && ylen == rows (ydat)) |
|
|
535 len = 3 * xlen; |
|
|
536 zz = zeros (ylen, len); |
|
|
537 k = 1; |
|
|
538 for kk = 1:3:len |
|
|
539 zz(:,kk) = xdat(:,k); |
|
|
540 zz(:,kk+1) = ydat(:,k); |
|
|
541 zz(:,kk+2) = zdat(:,k); |
|
|
542 k++; |
|
|
543 endfor |
|
|
544 data{data_idx} = zz; |
|
|
545 endif |
|
|
546 usingclause{data_idx} = "using ($1):($2):($3)"; |
|
|
547 withclause{data_idx} = "with line palette"; |
|
|
548 |
|
|
549 fputs (plot_stream, "unset parametric;\n"); |
|
|
550 fputs (plot_stream, "set hidden3d;\n"); |
|
|
551 fputs (plot_stream, "set style data lines;\n"); |
|
|
552 fputs (plot_stream, "set surface;\n"); |
|
|
553 fputs (plot_stream, "unset contour;\n"); |
|
|
554 fputs (plot_stream, "set palette defined (0 \"dark-blue\", 1 \"blue\", 2 \"cyan\", 3 \"yellow\", 4 \"red\" , 5 \"dark-red\");\n"); |
|
|
555 fputs (plot_stream, "unset colorbox;\n"); |
|
|
556 |
|
|
557 case "text" |
|
|
558 lpos = obj.position; |
|
|
559 label = obj.string; |
|
|
560 halign = obj.horizontalalignment; |
|
|
561 if (nd == 3) |
|
|
562 fprintf (plot_stream, "set label \"%s\" at %d,%d,%d %s;\n", |
|
|
563 undo_string_escapes (label), |
|
|
564 lpos(1), lpos(2), lpos(3), halign); |
|
|
565 else |
|
|
566 fprintf (plot_stream, "set label \"%s\" at %d,%d %s;\n", |
|
|
567 undo_string_escapes (label), |
|
|
568 lpos(1), lpos(2), halign); |
|
|
569 endif |
|
|
570 |
|
|
571 otherwise |
|
|
572 error ("__go_draw_axes__: unknown object class, %s", |
|
|
573 obj.type); |
|
|
574 endswitch |
|
|
575 |
|
|
576 endfor |
|
|
577 |
|
|
578 have_data = ! isempty (data); |
|
|
579 |
|
|
580 if (xautoscale && have_data) |
|
|
581 xlim = get_axis_limits (xmin, xmax, xminp, xlogscale); |
|
|
582 set (h, "xlim", xlim, "xlimmode", "auto"); |
|
|
583 else |
|
|
584 xlim = axis_obj.xlim; |
|
|
585 endif |
|
|
586 if (strcmp (axis_obj.xdir, "reverse")) |
|
|
587 xdir = "reverse"; |
|
|
588 else |
|
|
589 xdir = "noreverse"; |
|
|
590 endif |
|
|
591 fprintf (plot_stream, "set xrange [%g:%g] %s;\n", xlim, xdir); |
|
|
592 |
|
|
593 if (yautoscale && have_data) |
|
|
594 ylim = get_axis_limits (ymin, ymax, yminp, ylogscale); |
|
|
595 set (h, "ylim", ylim, "ylimmode", "auto"); |
|
|
596 else |
|
|
597 ylim = axis_obj.ylim; |
|
|
598 endif |
|
|
599 if (strcmp (axis_obj.ydir, "reverse")) |
|
|
600 ydir = "reverse"; |
|
|
601 else |
|
|
602 ydir = "noreverse"; |
|
|
603 endif |
|
|
604 fprintf (plot_stream, "set yrange [%g:%g] %s;\n", ylim, ydir); |
|
|
605 |
|
|
606 if (nd == 3) |
|
|
607 if (zautoscale && have_data) |
|
|
608 zlim = get_axis_limits (zmin, zmax, zminp, zlogscale); |
|
|
609 set (h, "zlim", zlim, "zlimmode", "auto"); |
|
|
610 else |
|
|
611 zlim = axis_obj.zlim; |
|
|
612 endif |
|
|
613 if (strcmp (axis_obj.zdir, "reverse")) |
|
|
614 zdir = "reverse"; |
|
|
615 else |
|
|
616 zdir = "noreverse"; |
|
|
617 endif |
|
|
618 fprintf (plot_stream, "set zrange [%g:%g] %s;\n", zlim, zdir); |
|
|
619 endif |
|
|
620 |
|
|
621 if (strcmp (axis_obj.box, "on")) |
|
|
622 if (nd == 3) |
|
|
623 fputs (plot_stream, "set border 4095;\n"); |
|
|
624 else |
|
|
625 fputs (plot_stream, "set border 431;\n"); |
|
|
626 endif |
|
|
627 else |
|
|
628 if (nd == 3) |
|
|
629 fputs (plot_stream, "set border 895;\n"); |
|
|
630 else |
|
|
631 fputs (plot_stream, "set border 3;\n"); |
|
|
632 fputs (plot_stream, "set xtics nomirror; set ytics nomirror;\n"); |
|
|
633 endif |
|
|
634 endif |
|
|
635 |
|
|
636 if (strcmp (axis_obj.key, "on")) |
|
|
637 if (strcmp (axis_obj.keybox, "on")) |
|
|
638 box = "box"; |
|
|
639 else |
|
|
640 box = "nobox"; |
|
|
641 endif |
|
|
642 inout = "inside"; |
|
|
643 switch (axis_obj.keypos) |
|
|
644 case -1 |
|
|
645 pos = "right bottom"; |
|
|
646 inout = "outside"; |
|
|
647 case 1 |
|
|
648 pos = "right top"; |
|
|
649 case 2 |
|
|
650 pos = "left top"; |
|
|
651 case 3 |
|
|
652 pos = "left bottom"; |
|
|
653 case 4 |
|
|
654 pos = "right bottom"; |
|
|
655 otherwise |
|
|
656 pos = ""; |
|
|
657 endswitch |
|
|
658 if (! have_newer_gnuplot) |
|
|
659 inout = ""; |
|
|
660 endif |
|
|
661 fprintf (plot_stream, "set key %s %s %s;\n", inout, pos, box); |
|
|
662 else |
|
|
663 fputs (plot_stream, "unset key;\n"); |
|
|
664 endif |
|
|
665 |
|
|
666 fputs (plot_stream, "set style data lines;\n"); |
|
|
667 |
|
|
668 if (! use_gnuplot_for_images) |
|
|
669 for i = 1:ximg_data_idx |
|
|
670 view_fcn (xlim, ylim, ximg_data{i}, view_zoom, view_cmd); |
|
|
671 endfor |
|
|
672 endif |
|
|
673 |
|
|
674 if (have_data) |
|
|
675 |
|
|
676 if (nd == 2) |
|
|
677 plot_cmd = "plot"; |
|
|
678 else |
|
|
679 plot_cmd = "splot"; |
|
|
680 rot_x = 90 - axis_obj.view(2); |
|
|
681 rot_z = axis_obj.view(1); |
|
|
682 while (rot_z < 0) |
|
|
683 rot_z += 360; |
|
|
684 endwhile |
|
|
685 fprintf (plot_stream, "set view %g, %g;\n", rot_x, rot_z); |
|
|
686 endif |
|
|
687 fprintf (plot_stream, "%s \"%s\" %s %s %s", plot_cmd, |
|
|
688 filespec{1}, usingclause{1}, titlespec{1}, withclause{1}); |
|
|
689 for i = 2:data_idx |
|
|
690 fprintf (plot_stream, ", \"%s\" %s %s %s", |
|
|
691 filespec{i}, usingclause{i}, titlespec{i}, withclause{i}); |
|
|
692 endfor |
|
|
693 fputs (plot_stream, ";\n"); |
|
|
694 for i = 1:data_idx |
|
|
695 if (strcmp (filespec{i}, "-")) |
|
|
696 if (nd == 2) |
|
|
697 fprintf (plot_stream, |
|
|
698 strcat (repmat ("%g ", 1, rows (data{i})), "\n"), |
|
|
699 data{i}); |
|
|
700 else |
|
|
701 if (parametric(i)) |
|
|
702 fprintf (plot_stream, "%g %g %g\n", data{i}); |
|
|
703 else |
|
|
704 tmp = data{i}; |
|
|
705 nc = columns (tmp); |
|
|
706 for j = 1:3:nc |
|
|
707 fprintf (plot_stream, "%g %g %g\n", tmp(:,j:j+2)'); |
|
|
708 fputs (plot_stream, "\n"); |
|
|
709 endfor |
|
|
710 endif |
|
|
711 endif |
|
|
712 fputs (plot_stream, "e\n"); |
|
|
713 endif |
|
|
714 endfor |
|
|
715 endif |
|
|
716 |
|
|
717 fflush (plot_stream); |
|
|
718 |
|
|
719 else |
|
|
720 print_usage (); |
|
|
721 endif |
|
|
722 |
|
|
723 endfunction |
|
|
724 |
|
|
725 ## Attempt to make "nice" limits from the actual max and min of the |
|
|
726 ## data. For log plots, we will also use the smallest strictly positive |
|
|
727 ## value. |
|
|
728 |
|
|
729 function lim = get_axis_limits (min_val, max_val, min_pos, logscale) |
|
|
730 |
|
|
731 if (logscale) |
|
|
732 if (isinf (min_pos)) |
|
|
733 warning ("axis: logscale with no positive values to plot"); |
|
|
734 endif |
|
|
735 if (min_val < 0) |
|
|
736 min_val = min_pos; |
|
|
737 if (max_val < 0) |
|
|
738 max_val = min_pos; |
|
|
739 endif |
|
|
740 warning ("axis: omitting negative data in log plot"); |
|
|
741 endif |
|
|
742 if (min_val == max_val) |
|
|
743 min_val = 0.9 * min_val; |
|
|
744 max_val = 1.1 * max_val; |
|
|
745 endif |
|
|
746 min_val = 10 ^ floor (log10 (min_val)); |
|
|
747 max_val = 10 ^ ceil (log10 (max_val)); |
|
|
748 else |
|
|
749 if (min_val == 0 && max_val == 0) |
|
|
750 min_val = -1; |
|
|
751 max_val = 1; |
|
|
752 elseif (min_val == max_val) |
|
|
753 min_val = 0.9 * min_val; |
|
|
754 max_val = 1.1 * max_val; |
|
|
755 endif |
|
|
756 ## FIXME -- to do a better job, we should consider the tic spacing. |
|
|
757 scale = 10 ^ floor (log10 (max_val - min_val) - 1); |
|
|
758 min_val = scale * floor (min_val / scale); |
|
|
759 max_val = scale * ceil (max_val / scale); |
|
|
760 endif |
|
|
761 |
|
|
762 lim = [min_val, max_val]; |
|
|
763 |
|
|
764 endfunction |
|
|
765 |
|
|
766 function style = do_linestyle_command (obj, idx, plot_stream) |
|
|
767 |
|
|
768 have_newer_gnuplot = compare_versions (__gnuplot_version__ (), "4.0", ">"); |
|
|
769 |
|
|
770 if (have_newer_gnuplot) |
|
|
771 fprintf (plot_stream, "set style line %d default;\n", idx); |
|
|
772 endif |
|
|
773 fprintf (plot_stream, "set style line %d", idx); |
|
|
774 |
|
|
775 found_style = false; |
|
|
776 |
|
|
777 if (isfield (obj, "color") && have_newer_gnuplot) |
|
|
778 color = obj.color; |
|
|
779 if (isnumeric (color)) |
|
|
780 fprintf (plot_stream, " linecolor rgb \"#%02x%02x%02x\"", |
|
|
781 round (255*color)); |
|
|
782 endif |
|
|
783 found_style = true; |
|
|
784 endif |
|
|
785 |
|
|
786 if (isfield (obj, "linestyle")) |
|
|
787 switch (obj.linestyle) |
|
|
788 case "-" |
|
|
789 lt = "lines"; |
|
|
790 case "--" |
|
|
791 lt = ""; |
|
|
792 case ":" |
|
|
793 lt = ""; |
|
|
794 case "-." |
|
|
795 lt = ""; |
|
|
796 case "none" |
|
|
797 lt = ""; |
|
|
798 otherwise |
|
|
799 lt = ""; |
|
|
800 endswitch |
|
|
801 else |
|
|
802 lt = ""; |
|
|
803 endif |
|
|
804 |
|
|
805 if (isfield (obj, "linewidth")) |
|
|
806 fprintf (plot_stream, " linewidth %f", obj.linewidth); |
|
|
807 found_style = true; |
|
|
808 endif |
|
|
809 |
|
|
810 if (isfield (obj, "marker")) |
|
|
811 switch (obj.marker) |
|
|
812 case "+" |
|
|
813 pt = "1"; |
|
|
814 case "o" |
|
|
815 pt = "7"; |
|
|
816 case "*" |
|
|
817 pt = "3"; |
|
|
818 case "." |
|
|
819 pt = "0"; |
|
|
820 case "x" |
|
|
821 pt = "2"; |
|
|
822 case {"square", "s"} |
|
|
823 pt = "5"; |
|
|
824 case {"diamond", "d"} |
|
|
825 pt = "13"; |
|
|
826 case "^" |
|
|
827 pt = "9"; |
|
|
828 case "v" |
|
|
829 pt = "11"; |
|
|
830 case ">" |
|
|
831 pt = "8"; |
|
|
832 case "<" |
|
|
833 pt = "10"; |
|
|
834 case {"pentagram", "p"} |
|
|
835 pt = "4"; |
|
|
836 case {"hexagram", "h"} |
|
|
837 pt = "6"; |
|
|
838 case "none" |
|
|
839 pt = ""; |
|
|
840 otherwise |
|
|
841 pt = ""; |
|
|
842 endswitch |
|
|
843 if (! isempty (pt)) |
|
|
844 fprintf (plot_stream, " pointtype %s", pt); |
|
|
845 found_style = true; |
|
|
846 endif |
|
|
847 else |
|
|
848 pt = ""; |
|
|
849 endif |
|
|
850 |
|
|
851 if (isfield (obj, "markersize")) |
|
|
852 fprintf (plot_stream, " pointsize %f", obj.markersize); |
|
|
853 found_style = true; |
|
|
854 endif |
|
|
855 |
|
|
856 style = "lines"; |
|
|
857 if (isempty (lt)) |
|
|
858 if (! isempty (pt)) |
|
|
859 style = "points"; |
|
|
860 endif |
|
|
861 elseif (! isempty (pt)) |
|
|
862 style = "linespoints"; |
|
|
863 endif |
|
|
864 |
|
|
865 if (have_newer_gnuplot && ! found_style) |
|
|
866 fputs (plot_stream, " default"); |
|
|
867 endif |
|
|
868 |
|
|
869 fputs (plot_stream, ";\n"); |
|
|
870 |
|
|
871 endfunction |
