source: trunk/src/Plotter2.cpp@ 2895

Last change on this file since 2895 was 2895, checked in by WataruKawasaki, 11 years ago

New Development: Yes

JIRA Issue: Yes CAS-6168

Ready for Test: Yes

Interface Changes:

What Interface Changed:

Test Programs:

Put in Release Notes: Yes

Module(s): sd.plotter2

Description: added methods to sd.plotter2 to enable changing the size and shape of the entire region to be plotted (which is referred to as 'viewsurface' in PGPLOT).
File size: 35.7 KB
Line 
1#include "Plotter2.h"
2
3namespace asap {
4
5Plotter2RectInfo::Plotter2RectInfo() {
6 xmin = 0.0;
7 xmax = 1.0;
8 ymin = 0.0;
9 ymax = 1.0;
10 color = 15; // gray
11 fill = 4; // hatch
12 width = 1;
13}
14
15Plotter2RectInfo::~Plotter2RectInfo() {
16}
17
18Plotter2DataInfo::Plotter2DataInfo() {
19 xData.clear();
20 yData.clear();
21
22 drawLine = true;
23 lineStyle = 1; // solid line
24 lineWidth = 1;
25 lineColor = -1; // undefined (default color should be assigned)
26
27 drawMarker = false;
28 markerType = 20; // small circle
29 markerSize = 1.0;
30 markerColor = 1; // default foreground color (black)
31
32 hasData = false; // has no data
33}
34
35Plotter2DataInfo::~Plotter2DataInfo() {
36}
37
38Plotter2ViewportInfo::Plotter2ViewportInfo() {
39 showViewport = true;
40
41 vpPosXMin = 0.1;
42 vpPosXMax = 0.9;
43 vpPosYMin = 0.1;
44 vpPosYMax = 0.9;
45
46 vpRangeXMin = 0.0;
47 vpRangeXMax = 1.0;
48 vpRangeYMin = 0.0;
49 vpRangeYMax = 1.0;
50 isAutoRangeX = true;
51 isAutoRangeY = true;
52 autoRangeMarginX = 0.0;
53 autoRangeMarginY = 0.1;
54
55 hasDataRange = false;
56
57 isAutoTickIntervalX = true;
58 isAutoTickIntervalY = true;
59 majorTickIntervalX = 0.1;
60 majorTickIntervalY = 0.1;
61 nMajorTickWithinTickNumsX = 2;
62 nMajorTickWithinTickNumsY = 2;
63 nMinorTickWithinMajorTicksX = 5;
64 nMinorTickWithinMajorTicksY = 5;
65
66 numLocationX = "b";
67 numLocationY = "l";
68
69 fontSizeDef = 1.0;
70
71 vData.clear();
72 vRect.clear();
73
74 labelXString = "";
75 labelXPosX = 0.5;
76 labelXPosY = 0.05;
77 labelXAngle = 0.0;
78 labelXFJust = 0.5;
79 labelXSize = fontSizeDef * 1.1;
80 labelXColor = 1;
81 labelXBColor = 0;
82
83 labelYString = "";
84 labelYPosX = 0.05;
85 labelYPosY = 0.5;
86 labelYAngle = 90.0;
87 labelYFJust = 0.5;
88 labelYSize = fontSizeDef * 1.1;
89 labelYColor = 1;
90 labelYBColor = 0;
91
92 titleString = "";
93 titlePosX = 0.5;
94 titlePosY = 0.95;
95 titleAngle = 0.0;
96 titleFJust = 0.5;
97 titleSize = fontSizeDef * 1.5;
98 titleColor = 1;
99 titleBColor = 0;
100
101 vpBColor = -1; // transparent (<0)
102}
103
104Plotter2ViewportInfo::~Plotter2ViewportInfo() {
105 vData.clear();
106}
107
108void Plotter2ViewportInfo::adjustRange() {
109 if (hasDataRange) {
110 if (isAutoRangeX) {
111 adjustRangeX(&vpRangeXMin, &vpRangeXMax);
112 }
113 if (isAutoRangeY) {
114 adjustRangeY(&vpRangeYMin, &vpRangeYMax);
115 }
116 }
117}
118
119void Plotter2ViewportInfo::adjustRangeX(float* xmin, float* xmax) {
120 float xmargin = (maxXData - minXData) * autoRangeMarginX;
121 *xmin = minXData - xmargin;
122 *xmax = maxXData + xmargin;
123}
124
125void Plotter2ViewportInfo::adjustRangeY(float* ymin, float* ymax) {
126 float ymargin = (maxYData - minYData) * autoRangeMarginY;
127 *ymin = minYData - ymargin;
128 *ymax = maxYData + ymargin;
129}
130
131std::vector<float> Plotter2ViewportInfo::getRangeX() {
132 float minX, maxX;
133
134 if (isAutoRangeX) {
135 adjustRangeX(&minX, &maxX);
136 } else {
137 minX = vpRangeXMin;
138 maxX = vpRangeXMax;
139 }
140
141 std::vector<float> res;
142 res.clear();
143 res.push_back(minX);
144 res.push_back(maxX);
145
146 return res;
147}
148
149std::vector<float> Plotter2ViewportInfo::getRangeY() {
150 float minY, maxY;
151
152 if (isAutoRangeY) {
153 adjustRangeY(&minY, &maxY);
154 } else {
155 minY = vpRangeYMin;
156 maxY = vpRangeYMax;
157 }
158
159 std::vector<float> res;
160 res.clear();
161 res.push_back(minY);
162 res.push_back(maxY);
163
164 return res;
165}
166
167void Plotter2ViewportInfo::adjustTickInterval() {
168 if (hasDataRange) {
169 if (isAutoTickIntervalX) {
170 adjustTickIntervalX(vpRangeXMin, vpRangeXMax);
171 }
172 if (isAutoTickIntervalY) {
173 adjustTickIntervalY(vpRangeYMin, vpRangeYMax);
174 }
175 }
176}
177
178void Plotter2ViewportInfo::adjustTickIntervalX(const float xmin, const float xmax) {
179 majorTickIntervalX = (float)pow(10.0, ceil(log10((xmax - xmin)/10.0)));
180 if ((xmax - xmin) / majorTickIntervalX < 4.0) {
181 majorTickIntervalX /= 2.0;
182 }
183 if ((xmax - xmin) / majorTickIntervalX < 4.0) {
184 majorTickIntervalX /= 2.0;
185 }
186}
187
188void Plotter2ViewportInfo::adjustTickIntervalY(const float ymin, const float ymax) {
189 majorTickIntervalY = (float)pow(10.0, ceil(log10((ymax - ymin)/10.0)));
190 if ((ymax - ymin) / majorTickIntervalY < 4.0) {
191 majorTickIntervalY /= 2.0;
192 }
193 if ((ymax - ymin) / majorTickIntervalY < 4.0) {
194 majorTickIntervalY /= 2.0;
195 }
196}
197
198void Plotter2ViewportInfo::setData(const std::vector<float>& inXData, const std::vector<float>& inYData, const int id) {
199 if (!hasDataRange) {
200 minXData = inXData[0];
201 maxXData = inXData[0];
202 minYData = inYData[0];
203 maxYData = inYData[0];
204
205 hasDataRange = true;
206 }
207
208 Plotter2DataInfo* info = &vData[id];
209
210 info->xData.clear();
211 info->xData.reserve(inXData.size());
212 for (unsigned int i = 0; i < inXData.size(); ++i) {
213 info->xData.push_back(inXData[i]);
214
215 if (!info->hasData) {
216 updateXDataRange(inXData[i]);
217 }
218 }
219
220 info->yData.clear();
221 info->yData.reserve(inYData.size());
222 for (unsigned int i = 0; i < inYData.size(); ++i) {
223 info->yData.push_back(inYData[i]);
224
225 if (!info->hasData) {
226 updateYDataRange(inYData[i]);
227 }
228 }
229
230 if (info->hasData) {
231 updateAllDataRanges();
232 } else {
233 info->hasData = true;
234 }
235
236 adjustRange();
237 adjustTickInterval();
238}
239
240void Plotter2ViewportInfo::updateXDataRange(const float data) {
241 if (data < minXData) {
242 minXData = data;
243 }
244 if (maxXData < data) {
245 maxXData = data;
246 }
247}
248
249void Plotter2ViewportInfo::updateYDataRange(const float data) {
250 if (data < minYData) {
251 minYData = data;
252 }
253 if (maxYData < data) {
254 maxYData = data;
255 }
256}
257
258void Plotter2ViewportInfo::updateAllDataRanges() {
259 minXData = vData[0].xData[0];
260 maxXData = minXData;
261 minYData = vData[0].yData[0];
262 maxYData = minYData;
263
264 for (unsigned int i = 0; i < vData.size(); ++i) {
265 for (unsigned int j = 0; j < vData[i].xData.size(); ++j) {
266 updateXDataRange(vData[i].xData[j]);
267 updateYDataRange(vData[i].yData[j]);
268 }
269 }
270}
271
272void Plotter2ViewportInfo::getWorldCoordByWindowCoord(const float winX, const float winY, float* worldX, float* worldY) {
273 float xratio = (winX - vpPosXMin) / (vpPosXMax - vpPosXMin);
274 if (winX < 0.0) {
275 xratio = 0.5;
276 }
277 *worldX = vpRangeXMin + xratio * (vpRangeXMax - vpRangeXMin);
278 float yratio = (winY - vpPosYMin) / (vpPosYMax - vpPosYMin);
279 if (winY < 0.0) {
280 yratio = 0.5;
281 }
282 *worldY = vpRangeYMin + yratio * (vpRangeYMax - vpRangeYMin);
283}
284
285Plotter2::Plotter2() {
286 filename = "";
287 device = "xwindow";
288 hasDevice = false;
289
290 vInfo.clear();
291 Plotter2ViewportInfo vi;
292 vInfo.push_back(vi);
293
294 hasDefaultViewport = true;
295 currentViewportId = 0;
296
297 width = 8.82796; // default viewsurface width seems to be this value.
298 aspect = 0.75; // default viewsurface aspect
299}
300
301Plotter2::~Plotter2() {
302 close();
303 vInfo.clear();
304}
305
306std::string Plotter2::getFileName() {
307 return filename;
308}
309
310void Plotter2::setFileName(const std::string& inFilename) {
311 filename = inFilename;
312}
313
314std::string Plotter2::getDevice() {
315 return device;
316}
317
318void Plotter2::setDevice(const std::string& inDevice) {
319 device = inDevice;
320}
321
322void Plotter2::open() {
323 cpgopen((filename + "/" + device).c_str());
324 hasDevice = true;
325}
326
327float Plotter2::getViewSurfaceWidth() {
328 return width;
329}
330
331float Plotter2::getViewSurfaceAspect() {
332 return aspect;
333}
334
335void Plotter2::setViewSurface(const float inWidth, const float inAspect) {
336 width = inWidth;
337 aspect = inAspect;
338}
339
340int Plotter2::addViewport(const float xmin, const float xmax, const float ymin, const float ymax) {
341 Plotter2ViewportInfo vi;
342
343 vi.vpPosXMin = xmin;
344 vi.vpPosXMax = xmax;
345 vi.vpPosYMin = ymin;
346 vi.vpPosYMax = ymax;
347
348 vInfo.push_back(vi);
349 currentViewportId = vInfo.size() - 1;
350
351 return currentViewportId;
352}
353
354void Plotter2::setViewport(const float xmin, const float xmax, const float ymin, const float ymax, const int id) {
355 Plotter2ViewportInfo* vi = &vInfo[id];
356
357 vi->vpPosXMin = xmin;
358 vi->vpPosXMax = xmax;
359 vi->vpPosYMin = ymin;
360 vi->vpPosYMax = ymax;
361
362 hasDefaultViewport = false;
363}
364
365void Plotter2::showViewport(const int inVpid) {
366 int vpid = inVpid;
367 if (vpid < 0) {
368 vpid = vInfo.size() - 1;
369 }
370 if (vpid < 0) {
371 exit(0);
372 }
373
374 Plotter2ViewportInfo* vi = &vInfo[vpid];
375 vi->showViewport = true;
376}
377
378void Plotter2::hideViewport(const int inVpid) {
379 int vpid = inVpid;
380 if (vpid < 0) {
381 vpid = vInfo.size() - 1;
382 }
383 if (vpid < 0) {
384 exit(0);
385 }
386
387 Plotter2ViewportInfo* vi = &vInfo[vpid];
388 vi->showViewport = false;
389}
390
391bool Plotter2::getHasDefaultViewport() {
392 return hasDefaultViewport;
393}
394
395int Plotter2::getCurrentViewportId() {
396 return currentViewportId;
397}
398
399void Plotter2::getViewInfo() {
400 std::string colorNames[16] = {"white", "black", "red", "green",
401 "blue", "cyan", "magenta", "yellow",
402 "orange", "yellowgreen", "emerald", "skyblue",
403 "purple", "pink", "gray", "lightgray"};
404 std::string lstyleNames[6] = {"", "solid", "dashed", "dash-dotted", "dotted", "dash-dot-dot-dotted"};
405 std::string fstyleNames[5] = {"", "solid", "outline", "hatched", "crosshatched"};
406
407 for (unsigned int i = 0; i < vInfo.size(); ++i) {
408 std::cout << "Viewport [ID = " << i << "] (" << (i+1) << "/" << vInfo.size() << ") ";
409 std::cout << "=============================================================" << std::endl;
410
411 std::cout << " Visible: " << (vInfo[i].showViewport ? "Yes" : "No") << std::endl;
412
413 std::cout << " Position in Window Coordinate: ";
414 std::cout << "X(left=" << vInfo[i].vpPosXMin << ",right=" << vInfo[i].vpPosXMax << "), ";
415 std::cout << "Y(bottom=" << vInfo[i].vpPosYMin << ",top=" << vInfo[i].vpPosYMax << ")" << std::endl;
416
417 std::cout << " Display Range: ";
418 std::cout << "X(" << (vInfo[i].isAutoRangeX ? "automatic" : "manual") << ", ";
419 std::cout << "min=" << vInfo[i].vpRangeXMin << ", max=" << vInfo[i].vpRangeXMax << "), ";
420 std::cout << "Y(" << (vInfo[i].isAutoRangeY ? "automatic" : "manual") << ", ";
421 std::cout << "min=" << vInfo[i].vpRangeYMin << ", max=" << vInfo[i].vpRangeYMax << ")" << std::endl;
422
423 std::cout << " Numbering/Ticks:" << std::endl;
424 std::cout << " X(" << ((vInfo[i].numLocationX == "b") ? "bottom" : ((vInfo[i].numLocationX == "t") ? "top" : (vInfo[i].numLocationX + " [*INVAILD*]"))) << ", interval[" << (vInfo[i].isAutoTickIntervalX ? "automatic" : "manual") << ", ";
425 std::cout << "numbering=" << (vInfo[i].majorTickIntervalX * vInfo[i].nMajorTickWithinTickNumsX) << ", ";
426 std::cout << "majortick=" << vInfo[i].majorTickIntervalX << ", ";
427 std::cout << "minortick=" << (vInfo[i].majorTickIntervalX / vInfo[i].nMinorTickWithinMajorTicksX) << "])" << std::endl;
428 std::cout << " Y(" << ((vInfo[i].numLocationY == "l") ? "left" : ((vInfo[i].numLocationY == "r") ? "right" : (vInfo[i].numLocationY + " [*INVAILD*]"))) << ", interval[" << (vInfo[i].isAutoTickIntervalY ? "automatic" : "manual") << ", ";
429 std::cout << "numbering=" << (vInfo[i].majorTickIntervalY * vInfo[i].nMajorTickWithinTickNumsY) << ", ";
430 std::cout << "majortick=" << vInfo[i].majorTickIntervalY << ", ";
431 std::cout << "minortick=" << (vInfo[i].majorTickIntervalY / vInfo[i].nMinorTickWithinMajorTicksY) << "])" << std::endl;
432
433 std::cout << " X Label: ";
434 if (vInfo[i].labelXString != "") {
435 std::cout << "\"" << vInfo[i].labelXString << "\"" << std::endl;
436 std::cout << " position=(" << vInfo[i].labelXPosX << "," << vInfo[i].labelXPosY << "), ";
437 std::cout << "justification=" << vInfo[i].labelXFJust << ", ";
438 std::cout << "angle=" << vInfo[i].labelXAngle << "deg, ";
439 std::cout << "fontsize=" << vInfo[i].labelXSize << std::endl;
440 } else {
441 std::cout << "No" << std::endl;
442 }
443
444 std::cout << " Y Label: ";
445 if (vInfo[i].labelYString != "") {
446 std::cout << "\"" << vInfo[i].labelYString << "\"" << std::endl;
447 std::cout << " position=(" << vInfo[i].labelYPosX << "," << vInfo[i].labelYPosY << "), ";
448 std::cout << "justification=" << vInfo[i].labelYFJust << ", ";
449 std::cout << "angle=" << vInfo[i].labelYAngle << "deg, ";
450 std::cout << "fontsize=" << vInfo[i].labelYSize << std::endl;
451 } else {
452 std::cout << "No" << std::endl;
453 }
454
455 std::cout << " Title: ";
456 if (vInfo[i].titleString != "") {
457 std::cout << "\"" << vInfo[i].titleString << "\"" << std::endl;
458 std::cout << " position=(" << vInfo[i].titlePosX << "," << vInfo[i].titlePosY << "), ";
459 std::cout << "justification=" << vInfo[i].titleFJust << ", ";
460 std::cout << "angle=" << vInfo[i].titleAngle << "deg, ";
461 std::cout << "fontsize=" << vInfo[i].titleSize << std::endl;
462 } else {
463 std::cout << "No" << std::endl;
464 }
465
466 std::cout << " Background Color = ";
467 if (vInfo[i].vpBColor < 0) {
468 std::cout << "transparent" << std::endl;
469 } else if (vInfo[i].vpBColor < 16) {
470 std::cout << colorNames[vInfo[i].vpBColor] << std::endl;
471 } else {
472 std::cout << "index:" << vInfo[i].vpBColor << " [INVAILD VALUE!]" << std::endl;
473 }
474
475 for (unsigned int j = 0; j < vInfo[i].vData.size(); ++j) {
476 std::cout << " Dataset [ID = " << j << "] (" << (j+1) << "/" << vInfo[i].vData.size() << ") ";
477 std::cout << "----------------------------------------------" << std::endl;
478
479 bool showDataset = (vInfo[i].vData[j].drawLine || vInfo[i].vData[j].drawMarker);
480 bool showViewport = vInfo[i].showViewport;
481 std::cout << " Visible: " << (showViewport ? "" : "(");
482 std::cout << (showDataset ? "Yes" : "No") << (showViewport ? "" : ")") << std::endl;
483 std::cout << " Number of Data Points: " << vInfo[i].vData[j].xData.size() << std::endl;
484 if (vInfo[i].vData[j].drawLine) {
485 std::cout << " Line: color=";
486 if (vInfo[i].vData[j].lineColor < 0) {
487 int defaultColorIdx = (j + 1) % 15 + 1;
488 std::cout << "default(" << colorNames[defaultColorIdx] << ")";
489 } else if (vInfo[i].vData[j].lineColor < 16) {
490 std::cout << colorNames[vInfo[i].vData[j].lineColor];
491 } else {
492 std::cout << "index:" << vInfo[i].vData[j].lineColor << " [*INVAILD*]";
493 }
494 std::cout << ", width=" << vInfo[i].vData[j].lineWidth;
495 std::cout << ", style=" << lstyleNames[vInfo[i].vData[j].lineStyle] << std::endl;
496 }
497 if (vInfo[i].vData[j].drawMarker) {
498 std::cout << " Marker: color=";
499 if (vInfo[i].vData[j].markerColor < 0) {
500 std::cout << "default";
501 } else if (vInfo[i].vData[j].markerColor < 16) {
502 std::cout << colorNames[vInfo[i].vData[j].markerColor];
503 } else {
504 std::cout << "index:" << vInfo[i].vData[j].markerColor << " [*INVAILD*]";
505 }
506 std::cout << ", shape=" << vInfo[i].vData[j].markerType;
507 std::cout << ", size=" << vInfo[i].vData[j].markerSize << std::endl;
508 }
509
510 }
511
512 for (unsigned int j = 0; j < vInfo[i].vRect.size(); ++j) {
513 std::cout << " XMask [ID = " << j << "] (" << (j+1) << "/" << vInfo[i].vRect.size() << ") ";
514 std::cout << "----------------------------------------------" << std::endl;
515 std::cout << " Range: (min=" << vInfo[i].vRect[j].xmin << ", max=" << vInfo[i].vRect[j].xmax << ")" << std::endl;
516 std::cout << " Attributes: (color=" << colorNames[vInfo[i].vRect[j].color];
517 std::string fstyle = fstyleNames[vInfo[i].vRect[j].fill];
518 std::cout << ", fillstyle=" << fstyle;
519 if (fstyle == "outline") {
520 std::cout << ", outlinewidth=" << vInfo[i].vRect[j].width;
521 }
522 if ((fstyle == "hatched")||(fstyle == "crosshatched")) {
523 std::cout << ", hatchspacing=" << vInfo[i].vRect[j].hsep;
524 }
525 std::cout << ")" << std::endl;
526 }
527 }
528 std::cout << "=====================================================================================" << std::endl << std::flush;
529}
530
531void Plotter2::setRange(const float xmin, const float xmax, const float ymin, const float ymax, const int inVpid) {
532 setRangeX(xmin, xmax, inVpid);
533 setRangeY(ymin, ymax, inVpid);
534}
535
536void Plotter2::setRangeX(const float xmin, const float xmax, const int inVpid) {
537 int vpid = inVpid;
538 if (vpid < 0) {
539 vpid = vInfo.size() - 1;
540 }
541 if (vpid < 0) {
542 Plotter2ViewportInfo vi;
543 vInfo.push_back(vi);
544 vpid = 0;
545 }
546
547 Plotter2ViewportInfo* vi = &vInfo[vpid];
548 vi->vpRangeXMin = xmin;
549 vi->vpRangeXMax = xmax;
550 vi->isAutoRangeX = false;
551}
552
553void Plotter2::setRangeY(const float ymin, const float ymax, const int inVpid) {
554 int vpid = inVpid;
555 if (vpid < 0) {
556 vpid = vInfo.size() - 1;
557 }
558 if (vpid < 0) {
559 Plotter2ViewportInfo vi;
560 vInfo.push_back(vi);
561 vpid = 0;
562 }
563
564 Plotter2ViewportInfo* vi = &vInfo[vpid];
565 vi->vpRangeYMin = ymin;
566 vi->vpRangeYMax = ymax;
567 vi->isAutoRangeY = false;
568}
569
570std::vector<float> Plotter2::getRangeX(const int inVpid) {
571 int vpid = inVpid;
572 if (vpid < 0) {
573 vpid = vInfo.size() - 1;
574 }
575 if (vpid < 0) {
576 exit(0);
577 }
578
579 return vInfo[vpid].getRangeX();
580}
581
582std::vector<float> Plotter2::getRangeY(const int inVpid) {
583 int vpid = inVpid;
584 if (vpid < 0) {
585 vpid = vInfo.size() - 1;
586 }
587 if (vpid < 0) {
588 exit(0);
589 }
590
591 return vInfo[vpid].getRangeY();
592}
593
594void Plotter2::setAutoRange(const int inVpid) {
595 setAutoRangeX(inVpid);
596 setAutoRangeY(inVpid);
597}
598
599void Plotter2::setAutoRangeX(const int inVpid) {
600 int vpid = inVpid;
601 if (vpid < 0) {
602 vpid = vInfo.size() - 1;
603 }
604 if (vpid < 0) {
605 Plotter2ViewportInfo vi;
606 vInfo.push_back(vi);
607 vpid = 0;
608 }
609
610 Plotter2ViewportInfo* vi = &vInfo[vpid];
611 vi->isAutoRangeX = true;
612}
613
614void Plotter2::setAutoRangeY(const int inVpid) {
615 int vpid = inVpid;
616 if (vpid < 0) {
617 vpid = vInfo.size() - 1;
618 }
619 if (vpid < 0) {
620 Plotter2ViewportInfo vi;
621 vInfo.push_back(vi);
622 vpid = 0;
623 }
624
625 Plotter2ViewportInfo* vi = &vInfo[vpid];
626 vi->isAutoRangeY = true;
627}
628
629void Plotter2::setFontSizeDef(const float size, const int inVpid) {
630 int vpid = inVpid;
631 if (vpid < 0) {
632 vpid = vInfo.size() - 1;
633 }
634 if (vpid < 0) {
635 Plotter2ViewportInfo vi;
636 vInfo.push_back(vi);
637 vpid = 0;
638 }
639
640 Plotter2ViewportInfo* vi = &vInfo[vpid];
641 vi->fontSizeDef = size;
642}
643
644void Plotter2::setTicksX(const float interval, const int num, const int inVpid) {
645 int vpid = inVpid;
646 if (vpid < 0) {
647 vpid = vInfo.size() - 1;
648 }
649 if (vpid < 0) {
650 Plotter2ViewportInfo vi;
651 vInfo.push_back(vi);
652 vpid = 0;
653 }
654
655 Plotter2ViewportInfo* vi = &vInfo[vpid];
656 vi->majorTickIntervalX = interval;
657 vi->nMinorTickWithinMajorTicksX = num;
658 vi->isAutoTickIntervalX = false;
659}
660
661void Plotter2::setTicksY(const float interval, const int num, const int inVpid) {
662 int vpid = inVpid;
663 if (vpid < 0) {
664 vpid = vInfo.size() - 1;
665 }
666 if (vpid < 0) {
667 Plotter2ViewportInfo vi;
668 vInfo.push_back(vi);
669 vpid = 0;
670 }
671
672 Plotter2ViewportInfo* vi = &vInfo[vpid];
673 vi->majorTickIntervalY = interval;
674 vi->nMinorTickWithinMajorTicksY = num;
675 vi->isAutoTickIntervalY = false;
676}
677
678void Plotter2::setAutoTicks(const int inVpid) {
679 setAutoTicksX(inVpid);
680 setAutoTicksY(inVpid);
681}
682
683void Plotter2::setAutoTicksX(const int inVpid) {
684 int vpid = inVpid;
685 if (vpid < 0) {
686 vpid = vInfo.size() - 1;
687 }
688 if (vpid < 0) {
689 Plotter2ViewportInfo vi;
690 vInfo.push_back(vi);
691 vpid = 0;
692 }
693
694 Plotter2ViewportInfo* vi = &vInfo[vpid];
695 vi->isAutoTickIntervalX = true;
696}
697
698void Plotter2::setAutoTicksY(const int inVpid) {
699 int vpid = inVpid;
700 if (vpid < 0) {
701 vpid = vInfo.size() - 1;
702 }
703 if (vpid < 0) {
704 Plotter2ViewportInfo vi;
705 vInfo.push_back(vi);
706 vpid = 0;
707 }
708
709 Plotter2ViewportInfo* vi = &vInfo[vpid];
710 vi->isAutoTickIntervalY = true;
711}
712
713void Plotter2::setNumIntervalX(const float interval, const int inVpid) {
714 int vpid = inVpid;
715 if (vpid < 0) {
716 vpid = vInfo.size() - 1;
717 }
718 if (vpid < 0) {
719 Plotter2ViewportInfo vi;
720 vInfo.push_back(vi);
721 vpid = 0;
722 }
723
724 Plotter2ViewportInfo* vi = &vInfo[vpid];
725 vi->nMajorTickWithinTickNumsX = (int)(interval / vi->majorTickIntervalX);
726}
727
728void Plotter2::setNumIntervalY(const float interval, const int inVpid) {
729 int vpid = inVpid;
730 if (vpid < 0) {
731 vpid = vInfo.size() - 1;
732 }
733 if (vpid < 0) {
734 Plotter2ViewportInfo vi;
735 vInfo.push_back(vi);
736 vpid = 0;
737 }
738
739 Plotter2ViewportInfo* vi = &vInfo[vpid];
740 vi->nMajorTickWithinTickNumsY = (int)(interval / vi->majorTickIntervalY);
741}
742
743void Plotter2::setNumLocationX(const std::string& side, const int inVpid) {
744 int vpid = inVpid;
745 if (vpid < 0) {
746 vpid = vInfo.size() - 1;
747 }
748 if (vpid < 0) {
749 Plotter2ViewportInfo vi;
750 vInfo.push_back(vi);
751 vpid = 0;
752 }
753
754 Plotter2ViewportInfo* vi = &vInfo[vpid];
755 vi->numLocationX = side;
756}
757
758void Plotter2::setNumLocationY(const std::string& side, const int inVpid) {
759 int vpid = inVpid;
760 if (vpid < 0) {
761 vpid = vInfo.size() - 1;
762 }
763 if (vpid < 0) {
764 Plotter2ViewportInfo vi;
765 vInfo.push_back(vi);
766 vpid = 0;
767 }
768
769 Plotter2ViewportInfo* vi = &vInfo[vpid];
770 vi->numLocationY = side;
771}
772
773void Plotter2::setData(const std::vector<float>& xdata, const std::vector<float>& ydata, const int inVpid, const int inDataid) {
774 int vpid = inVpid;
775 if (vpid < 0) {
776 vpid = vInfo.size() - 1;
777 }
778 if (vpid < 0) {
779 Plotter2ViewportInfo vi;
780 vInfo.push_back(vi);
781 vpid = 0;
782 }
783
784 Plotter2ViewportInfo* vi = &vInfo[vpid];
785
786 int dataid = inDataid;
787 if (dataid < 0) {
788 Plotter2DataInfo di;
789 vi->vData.push_back(di);
790 dataid = vi->vData.size() - 1;
791 }
792
793 vi->setData(xdata, ydata, dataid);
794}
795
796void Plotter2::setLine(const int color, const int width, const int style, const int inVpid, const int inDataid) {
797 int vpid = inVpid;
798 if (vpid < 0) {
799 vpid = vInfo.size() - 1;
800 }
801 if (vpid < 0) {
802 Plotter2ViewportInfo vi;
803 vInfo.push_back(vi);
804 vpid = 0;
805 }
806
807 int dataid = inDataid;
808 if (dataid < 0) {
809 dataid = vInfo[vpid].vData.size() - 1;
810 }
811
812 Plotter2ViewportInfo* vi = &vInfo[vpid];
813 vi->vData[dataid].drawLine = true;
814 vi->vData[dataid].lineColor = color;
815 vi->vData[dataid].lineWidth = width;
816 vi->vData[dataid].lineStyle = style;
817}
818
819void Plotter2::showLine(const int inVpid, const int inDataid) {
820 int vpid = inVpid;
821 if (vpid < 0) {
822 vpid = vInfo.size() - 1;
823 }
824 if (vpid < 0) {
825 exit(0);
826 }
827
828 int dataid = inDataid;
829 if (dataid < 0) {
830 dataid = vInfo[vpid].vData.size() - 1;
831 }
832 if (dataid < 0) {
833 exit(0);
834 }
835
836 Plotter2ViewportInfo* vi = &vInfo[vpid];
837 vi->vData[dataid].drawLine = true;
838}
839
840void Plotter2::hideLine(const int inVpid, const int inDataid) {
841 int vpid = inVpid;
842 if (vpid < 0) {
843 vpid = vInfo.size() - 1;
844 }
845 if (vpid < 0) {
846 exit(0);
847 }
848
849 int dataid = inDataid;
850 if (dataid < 0) {
851 dataid = vInfo[vpid].vData.size() - 1;
852 }
853 if (dataid < 0) {
854 exit(0);
855 }
856
857 Plotter2ViewportInfo* vi = &vInfo[vpid];
858 vi->vData[dataid].drawLine = false;
859}
860
861void Plotter2::setPoint(const int type, const float size, const int color, const int inVpid, const int inDataid) {
862 int vpid = inVpid;
863 if (vpid < 0) {
864 vpid = vInfo.size() - 1;
865 }
866 if (vpid < 0) {
867 Plotter2ViewportInfo vi;
868 vInfo.push_back(vi);
869 vpid = 0;
870 }
871
872 int dataid = inDataid;
873 if (dataid < 0) {
874 dataid = vInfo[vpid].vData.size() - 1;
875 }
876
877 Plotter2ViewportInfo* vi = &vInfo[vpid];
878 vi->vData[dataid].drawMarker = true;
879 vi->vData[dataid].markerType = type;
880 vi->vData[dataid].markerSize = size;
881 vi->vData[dataid].markerColor = color;
882}
883
884void Plotter2::showPoint(const int inVpid, const int inDataid) {
885 int vpid = inVpid;
886 if (vpid < 0) {
887 vpid = vInfo.size() - 1;
888 }
889 if (vpid < 0) {
890 exit(0);
891 }
892
893 int dataid = inDataid;
894 if (dataid < 0) {
895 dataid = vInfo[vpid].vData.size() - 1;
896 }
897 if (dataid < 0) {
898 exit(0);
899 }
900
901 Plotter2ViewportInfo* vi = &vInfo[vpid];
902 vi->vData[dataid].drawMarker = true;
903}
904
905void Plotter2::hidePoint(const int inVpid, const int inDataid) {
906 int vpid = inVpid;
907 if (vpid < 0) {
908 vpid = vInfo.size() - 1;
909 }
910 if (vpid < 0) {
911 exit(0);
912 }
913
914 int dataid = inDataid;
915 if (dataid < 0) {
916 dataid = vInfo[vpid].vData.size() - 1;
917 }
918 if (dataid < 0) {
919 exit(0);
920 }
921
922 Plotter2ViewportInfo* vi = &vInfo[vpid];
923 vi->vData[dataid].drawMarker = false;
924}
925
926void Plotter2::setMaskX(const float xmin, const float xmax, const int color, const int fill, const int width, const float hsep, const int inVpid) {
927 int vpid = inVpid;
928 if (vpid < 0) {
929 vpid = vInfo.size() - 1;
930 }
931 if (vpid < 0) {
932 Plotter2ViewportInfo vi;
933 vInfo.push_back(vi);
934 vpid = 0;
935 }
936
937 Plotter2ViewportInfo* vi = &vInfo[vpid];
938
939 Plotter2RectInfo ri;
940 ri.xmin = xmin;
941 ri.xmax = xmax;
942 //y-range of xmask should be calculated in plot().
943 //std::vector<float> yrange = vi->getRangeY();
944 //float yexcess = 0.1*(yrange[1] - yrange[0]);
945 //ri.ymin = yrange[0] - yexcess;
946 //ri.ymax = yrange[1] + yexcess;
947 ri.color = color;
948 ri.fill = fill;
949 ri.width = width;
950 ri.hsep = hsep;
951
952 vi->vRect.push_back(ri);
953}
954
955void Plotter2::setLabelX(const std::string& label, const float inPosx, const float inPosy, const float size, const std::string& style, const int color, const int bgcolor, const int inVpid) {
956 int vpid = inVpid;
957 if (vpid < 0) {
958 vpid = vInfo.size() - 1;
959 }
960 if (vpid < 0) {
961 Plotter2ViewportInfo vi;
962 vInfo.push_back(vi);
963 vpid = 0;
964 }
965
966 Plotter2ViewportInfo* vi = &vInfo[vpid];
967
968 std::string styleString;
969 if (style == "") {
970 styleString = "";
971 } else if (style == "roman") {
972 styleString = "\\fr";
973 } else if (style == "italic") {
974 styleString = "\\fi";
975 } else if (style == "script") {
976 styleString = "\\fs";
977 }
978 vi->labelXString = styleString + label;
979
980 float posx = inPosx;
981 if (posx < 0.0) {
982 posx = 0.5*(vi->vpPosXMin + vi->vpPosXMax);
983 }
984 vi->labelXPosX = posx;
985
986 float posy = inPosy;
987 if (posy < 0.0) {
988 posy = 0.35*vi->vpPosYMin;
989 }
990 vi->labelXPosY = posy;
991
992 vi->labelXAngle = 0.0;
993 vi->labelXFJust = 0.5;
994 vi->labelXSize = size;
995 vi->labelXColor = color;
996 vi->labelXBColor = bgcolor;
997}
998
999void Plotter2::setLabelY(const std::string& label, const float inPosx, const float inPosy, const float size, const std::string& style, const int color, const int bgcolor, const int inVpid) {
1000 int vpid = inVpid;
1001 if (vpid < 0) {
1002 vpid = vInfo.size() - 1;
1003 }
1004 if (vpid < 0) {
1005 Plotter2ViewportInfo vi;
1006 vInfo.push_back(vi);
1007 vpid = 0;
1008 }
1009
1010 Plotter2ViewportInfo* vi = &vInfo[vpid];
1011
1012 std::string styleString;
1013 if (style == "") {
1014 styleString = "";
1015 } else if (style == "roman") {
1016 styleString = "\\fr";
1017 } else if (style == "italic") {
1018 styleString = "\\fi";
1019 } else if (style == "script") {
1020 styleString = "\\fs";
1021 }
1022 vi->labelYString = styleString + label;
1023
1024 float posx = inPosx;
1025 if (posx < 0.0) {
1026 posx = 0.35*vi->vpPosXMin;
1027 }
1028 vi->labelYPosX = posx;
1029
1030 float posy = inPosy;
1031 if (posy < 0.0) {
1032 posy = 0.5*(vi->vpPosYMin + vi->vpPosYMax);
1033 }
1034 vi->labelYPosY = posy;
1035
1036 vi->labelYAngle = 90.0;
1037 vi->labelYFJust = 0.5;
1038 vi->labelYSize = size;
1039 vi->labelYColor = color;
1040 vi->labelYBColor = bgcolor;
1041}
1042
1043void Plotter2::setTitle(const std::string& label, const float inPosx, const float inPosy, const float size, const std::string& style, const int color, const int bgcolor, const int inVpid) {
1044 int vpid = inVpid;
1045 if (vpid < 0) {
1046 vpid = vInfo.size() - 1;
1047 }
1048 if (vpid < 0) {
1049 Plotter2ViewportInfo vi;
1050 vInfo.push_back(vi);
1051 vpid = 0;
1052 }
1053
1054 Plotter2ViewportInfo* vi = &vInfo[vpid];
1055
1056 std::string styleString;
1057 if (style == "") {
1058 styleString = "";
1059 } else if (style == "roman") {
1060 styleString = "\\fr";
1061 } else if (style == "italic") {
1062 styleString = "\\fi";
1063 } else if (style == "script") {
1064 styleString = "\\fs";
1065 }
1066 vi->titleString = styleString + label;
1067
1068 float posx = inPosx;
1069 if (posx < 0.0) {
1070 posx = 0.5*(vi->vpPosXMin + vi->vpPosXMax);
1071 }
1072 vi->titlePosX = posx;
1073
1074 float posy = inPosy;
1075 if (posy < 0.0) {
1076 posy = vi->vpPosYMax + 0.25*(1.0 - vi->vpPosYMax);
1077 }
1078 vi->titlePosY = posy;
1079
1080 vi->titleAngle = 0.0;
1081 vi->titleFJust = 0.5;
1082 vi->titleSize = size;
1083 vi->titleColor = color;
1084 vi->titleBColor = bgcolor;
1085}
1086
1087void Plotter2::setViewportBackgroundColor(const int bgcolor, const int inVpid) {
1088 int vpid = inVpid;
1089 if (vpid < 0) {
1090 vpid = vInfo.size() - 1;
1091 }
1092 if (vpid < 0) {
1093 exit(0);
1094 }
1095
1096 Plotter2ViewportInfo* vi = &vInfo[vpid];
1097 vi->vpBColor = bgcolor;
1098}
1099
1100/*
1101void Plotter2::setAnnotation(const std::string& label, const float posx, const float posy, const float angle, const float fjust, const float size, const std::string& style, const int color, const int bgcolor, const int inVpid) {
1102 int vpid = inVpid;
1103 if (vpid < 0) {
1104 vpid = vInfo.size() - 1;
1105 }
1106 if (vpid < 0) {
1107 Plotter2ViewportInfo vi;
1108 vInfo.push_back(vi);
1109 vpid = 0;
1110 }
1111
1112 std::string styleString;
1113 if (style == "") {
1114 styleString = "";
1115 } else if (style == "roman") {
1116 styleString = "\\fr";
1117 } else if (style == "italic") {
1118 styleString = "\\fi";
1119 } else if (style == "script") {
1120 styleString = "\\fs";
1121 }
1122
1123 Plotter2ViewportInfo* vi = &vInfo[vpid];
1124 //vi->titleString = styleString + label;
1125 //vi->titlePosX = posx;
1126 //vi->titlePosY = posy;
1127 //vi->titleAngle = angle;
1128 //vi->titleFJust = fjust;
1129 //vi->titleSize = size;
1130 //vi->titleColor = color;
1131 //vi->titleBColor = bgcolor;
1132}
1133*/
1134
1135void Plotter2::close() {
1136 if (hasDevice) {
1137 cpgclos();
1138 hasDevice = false;
1139 }
1140}
1141
1142void Plotter2::plot() {
1143 open();
1144
1145 if ((width > 0.0) && (aspect > 0.0)) {
1146 cpgpap(width, aspect);
1147 }
1148
1149 cpgscr(0, 1.0, 1.0, 1.0); // set background color white
1150 cpgscr(1, 0.0, 0.0, 0.0); // set foreground color black
1151
1152 for (unsigned int i = 0; i < vInfo.size(); ++i) {
1153 Plotter2ViewportInfo vi = vInfo[i];
1154
1155 if (vi.showViewport) {
1156 cpgstbg(0); // reset background colour to the initial one (white)
1157 cpgsci(1); // reset foreground colour to the initial one (black)
1158 cpgsls(1); // reset line style to solid
1159 cpgslw(1); // reset line width to 1
1160 cpgscf(1); // reset font style to normal
1161 cpgsch(vi.fontSizeDef);// reset font size
1162 cpgsfs(1); // reset fill style (solid)
1163
1164 // setup viewport
1165 cpgsvp(vi.vpPosXMin, vi.vpPosXMax, vi.vpPosYMin, vi.vpPosYMax);
1166 cpgswin(vi.vpRangeXMin, vi.vpRangeXMax, vi.vpRangeYMin, vi.vpRangeYMax);
1167
1168 // background color (default is transparent)
1169 if (vi.vpBColor >= 0) {
1170 cpgsci(vi.vpBColor);
1171 cpgrect(vi.vpRangeXMin, vi.vpRangeXMax, vi.vpRangeYMin, vi.vpRangeYMax);
1172 cpgsci(1); // reset foreground colour to the initial one (black)
1173 }
1174
1175 // data
1176 for (unsigned int j = 0; j < vi.vData.size(); ++j) {
1177 cpgstbg(0); // reset background colour to the initial one (white)
1178 cpgsci(1); // reset foreground colour to the initial one (black)
1179 cpgsls(1); // reset line style to solid
1180 cpgslw(1); // reset line width to 1
1181 cpgscf(1); // reset font style to normal
1182 cpgsch(vi.fontSizeDef);// reset font size
1183
1184 Plotter2DataInfo di = vi.vData[j];
1185 std::vector<float> vxdata = di.xData;
1186 int ndata = vxdata.size();
1187 float* pxdata = new float[ndata];
1188 float* pydata = new float[ndata];
1189 for (int k = 0; k < ndata; ++k) {
1190 pxdata[k] = di.xData[k];
1191 pydata[k] = di.yData[k];
1192 }
1193
1194 if (di.drawLine) {
1195 cpgsls(di.lineStyle);
1196 cpgslw(di.lineWidth);
1197 int colorIdx = di.lineColor;
1198 if (colorIdx < 0) {
1199 colorIdx = (j + 1) % 15 + 1;
1200 }
1201 cpgsci(colorIdx);
1202 cpgline(ndata, pxdata, pydata);
1203 }
1204
1205 if (di.drawMarker) {
1206 cpgsch(di.markerSize);
1207 cpgsci(di.markerColor);
1208 cpgpt(ndata, pxdata, pydata, di.markerType);
1209 }
1210
1211 delete [] pxdata;
1212 delete [] pydata;
1213 }
1214
1215 //calculate y-range of xmasks
1216 std::vector<float> yrange = vi.getRangeY();
1217 float yexcess = 0.1*(yrange[1] - yrange[0]);
1218 float xmaskymin = yrange[0] - yexcess;
1219 float xmaskymax = yrange[1] + yexcess;
1220
1221 // masks
1222 for (unsigned int j = 0; j < vi.vRect.size(); ++j) {
1223 cpgstbg(0); // reset background colour to the initial one (white)
1224 cpgsci(1); // reset foreground colour to the initial one (black)
1225 cpgsls(1); // reset line style to solid
1226 cpgslw(1); // reset line width to 1
1227 cpgscf(1); // reset font style to normal
1228 cpgsch(vi.fontSizeDef);// reset font size
1229 cpgsfs(1); // reset fill style (solid)
1230
1231 Plotter2RectInfo ri = vi.vRect[j];
1232 cpgsci(ri.color);
1233 cpgsfs(ri.fill);
1234 cpgslw(ri.width);
1235 cpgshs(45.0, ri.hsep, 0.0);
1236 float* mxdata = new float[4];
1237 float* mydata = new float[4];
1238 mxdata[0] = ri.xmin;
1239 mxdata[1] = ri.xmax;
1240 mxdata[2] = ri.xmax;
1241 mxdata[3] = ri.xmin;
1242 mydata[0] = xmaskymin;
1243 mydata[1] = xmaskymin;
1244 mydata[2] = xmaskymax;
1245 mydata[3] = xmaskymax;
1246 cpgpoly(4, mxdata, mydata);
1247 }
1248
1249 cpgstbg(0); // reset background colour to the initial one (white)
1250 cpgsci(1); // reset foreground colour to the initial one (black)
1251 cpgsls(1); // reset line style to solid
1252 cpgslw(1); // reset line width to 1
1253 cpgscf(1); // reset font style to normal
1254 cpgsch(vi.fontSizeDef);// reset font size
1255 cpgsfs(1); // reset fill style (solid)
1256
1257 cpgbox("BCTS", vi.majorTickIntervalX, vi.nMinorTickWithinMajorTicksX,
1258 "BCTSV", vi.majorTickIntervalY, vi.nMinorTickWithinMajorTicksY);
1259
1260 // viewport outline, ticks and number labels
1261 std::string numformatx, numformaty;
1262 if (vi.numLocationX == "b") {
1263 numformatx = "N";
1264 } else if (vi.numLocationX == "t") {
1265 numformatx = "M";
1266 } else if (vi.numLocationX == "") {
1267 numformatx = "";
1268 }
1269 if (vi.numLocationY == "l") {
1270 numformaty = "NV";
1271 } else if (vi.numLocationY == "r") {
1272 numformaty = "MV";
1273 } else if (vi.numLocationY == "") {
1274 numformaty = "";
1275 }
1276
1277 cpgbox(numformatx.c_str(), vi.majorTickIntervalX * vi.nMajorTickWithinTickNumsX, 0,
1278 numformaty.c_str(), vi.majorTickIntervalY * vi.nMajorTickWithinTickNumsY, 0);
1279
1280 float xpos, ypos;
1281
1282 // x-label
1283 vi.getWorldCoordByWindowCoord(vi.labelXPosX, vi.labelXPosY, &xpos, &ypos);
1284 cpgsch(vi.labelXSize);
1285 cpgsci(vi.labelXColor);
1286 cpgstbg(vi.labelXBColor); //outside viewports, works ONLY with /xwindow
1287 cpgptxt(xpos, ypos, vi.labelXAngle, vi.labelXFJust, vi.labelXString.c_str());
1288
1289 // y-label
1290 vi.getWorldCoordByWindowCoord(vi.labelYPosX, vi.labelYPosY, &xpos, &ypos);
1291 cpgsch(vi.labelYSize);
1292 cpgsci(vi.labelYColor);
1293 cpgstbg(vi.labelYBColor); //outside viewports, works ONLY with /xwindow
1294 cpgptxt(xpos, ypos, vi.labelYAngle, vi.labelYFJust, vi.labelYString.c_str());
1295
1296 // title
1297 vi.getWorldCoordByWindowCoord(vi.titlePosX, vi.titlePosY, &xpos, &ypos);
1298 cpgsch(vi.titleSize);
1299 cpgsci(vi.titleColor);
1300 cpgstbg(vi.titleBColor); //outside viewports, works ONLY with /xwindow
1301 cpgptxt(xpos, ypos, vi.titleAngle, vi.titleFJust, vi.titleString.c_str());
1302 }
1303
1304 }
1305
1306 close();
1307}
1308
1309} // namespace asap
Note: See TracBrowser for help on using the repository browser.