TPCCLIB
Loading...
Searching...
No Matches
Functions
plotfit.c File Reference

Plot measured and fitted TACs in SVG format. More...

#include "libtpcmodext.h"

Go to the source code of this file.

Functions

int plot_fitrange_svg (DFT *dft1, DFT *dft2, char *main_title, double x1, double x2, double y1, double y2, char *fname, int verbose)
 
int plot_fit_svg (DFT *dft1, DFT *dft2, char *main_title, char *fname, int verbose)
 

Detailed Description

Plot measured and fitted TACs in SVG format.

Author
Vesa Oikonen

Definition in file plotfit.c.

Function Documentation

◆ plot_fit_svg()

int plot_fit_svg ( DFT * dft1,
DFT * dft2,
char * main_title,
char * fname,
int verbose )

Writes plots of original and fitted TACs in SVG 1.1 format. Data must not contain NaNs.

Returns
Returns 0 if successful, otherwise nonzero.
See also
plot_svg, plot_fitrange_svg
Parameters
dft1Measured data points.
dft2Fitted data points. Times can be different but unit must be the same.
main_titleString for plot main title, or "".
fnameSVG filename; existing file is backed up.
verboseVerbose level; if zero, then nothing is printed to stderr or stdout.

Definition at line 216 of file plotfit.c.

227 {
228 int ret, n, ri, si, ei;
229 char x_title[64], y_title[64], tac_id[32], tac_title[64];
230 double minx, maxx, miny, maxy, tx1, tx2, ty1, ty2, f;
231 struct svg_viewports viewports; svg_init_viewports(&viewports);
232 int max_color_nr, color_nr;
233 int max_symbol_nr, symbol_nr;
234 SVG_LEGENDS legends; svg_init_legends(&legends);
235 FILE *fp_svg=NULL;
236
237
238 if(verbose>0) {
239 printf("%s(dft1, dft2, mt, fn, %d)\n", __func__, verbose);
240 }
241
242 /* Check data */
243 if(dft1==NULL || dft1->voiNr<1) return(1);
244 if(dft2==NULL) return(1);
245 if(dft2->voiNr!=dft1->voiNr) return(1);
246
247 int is_label=0; if(dft1->voiNr>1) is_label=1;
248
249 /* Check if file exists; backup, if necessary */
250 ret=backupExistingFile(fname, NULL, NULL); if(ret) return(2);
251
252 /* Determine the plot min and max values */
253 minx=maxx=miny=maxy=0;
254 ret=dftMinMax(dft1, &tx1, &tx2, &ty1, &ty2); if(ret) return(3);
255 minx=tx1; maxx=tx2; miny=ty1; maxy=ty2;
256 ret=dftMinMax(dft2, &tx1, &tx2, &ty1, &ty2); if(ret) return(3);
257 if(minx>tx1) minx=tx1;
258 if(maxx<tx2) maxx=tx2;
259 if(miny>ty1) miny=ty1;
260 if(maxy<ty2) maxy=ty2;
261 if(verbose>1) printf("minx:=%g\nmaxx:=%g\nminy:=%g\nmaxy:=%g\n", minx, maxx, miny, maxy);
262 if(miny>0.0) {f=maxy-miny; miny-=0.01*f;}
263
264 /* Calculate the axis ticks */
265 viewports.label_area_viewport.is=is_label; // needed for x axis ticks
266 viewports.x.fixed_min=0; viewports.y.fixed_min=0;
267 viewports.x.min=minx; viewports.x.max=maxx;
268 viewports.y.min=miny; viewports.y.max=maxy;
269 ret=svg_calculate_axes(&viewports, verbose-3); if(ret) return(4);
270
271 /* Set x and y axis titles based on activity and time units */
272 if(verbose>2) printf("set x title\n");
273 strcpy(x_title, "");
274 if(dft1->timeunit==DFTTIME_SEC || dft1->timeunit==DFTTIME_MIN)
275 sprintf(x_title, "Time (%s)", dftTimeunit(dft1->timeunit));
276 else if(dft1->timeunit!=DFTTIME_UNKNOWN)
277 strcpy(x_title, dftTimeunit(dft1->timeunit));
278 if(verbose>2) printf("set y title\n");
279 strcpy(y_title, "");
280 if(dftUnitId(dft1->unit)!=DFTUNIT_UNKNOWN)
281 strcpy(y_title, dftUnit(dftUnitId(dft1->unit)));
282
283 /* Set the plot window and window area sizes */
284 if(verbose>2) printf("set window sizes\n");
285 ret=svg_define_viewports(0, 0, strlen(main_title), strlen(y_title),
286 strlen(x_title), is_label, &viewports, verbose-3);
287 if(ret) return(5);
288
289 /* Initiate graphics file */
290 fp_svg=svg_initiate(fname, 0, 0, &viewports, NULL, verbose-3);
291 if(fp_svg==NULL) return(6);
292
293 /* Put the graph titles into their own viewports */
294 ret=svg_create_main_title(fp_svg, main_title, "", &viewports, NULL,verbose-3);
295 if(ret) return(7);
296 ret=svg_create_yaxis_title(fp_svg, y_title, &viewports, NULL, verbose-3);
297 if(ret) return(8);
298 ret=svg_create_xaxis_title(fp_svg, x_title, &viewports, NULL, verbose-3);
299 if(ret) return(9);
300
301 /* Put the plot into its own viewport */
302 ret=svg_start_plot_viewport(fp_svg, &viewports, NULL, verbose-3);
303 if(ret) return(10);
304
305 /* Start coordinate area viewport */
306 ret=svg_start_coordinate_viewport(fp_svg, &viewports, NULL, verbose-3);
307 if(ret) return(11);
308
309 /* Write plot axes */
310 ret=svg_write_axes(fp_svg, &viewports, NULL, verbose-3);
311 if(ret) {
312 if(verbose>0) printf("svg_write_axes() := %d\n", ret);
313 return(12);
314 }
315
316 /*
317 * Draw the plots
318 */
319 max_color_nr=0; while(svgColorName(max_color_nr)!=NULL) max_color_nr++;
320 if(max_color_nr<1) max_color_nr=1; // remainder works only if 2nd operator>0
321 if(verbose>3) printf("max_color_nr := %d\n", max_color_nr);
322 max_symbol_nr=0; while(svgSymbolName(max_symbol_nr)!=NULL) max_symbol_nr++;
323 if(max_symbol_nr<1) max_symbol_nr=1; // remainder works only if 2nd operator>0
324 if(verbose>3) printf("max_symbol_nr := %d\n", max_symbol_nr);
325 if(dft1->voiNr==1) color_nr=0; else color_nr=1;
326 symbol_nr=0;
327 for(ri=0, n=0; ri<dft1->voiNr; ri++) {
328 sprintf(tac_id, "plot_%d", n);
329/*
330 if(strlen(dft1->studynr)>0 && strcmp(dft1->studynr, ".")!=0)
331 sprintf(tac_title, "%s: %s", dft1->studynr, dft1->voi[ri].name);
332 else strcpy(tac_title, dft1->voi[ri].name);
333*/
334 rnameRmDots(dft1->voi[ri].name, tac_title);
335 /* Draw the fitted line */
336 for(si=0; si<dft2->frameNr; si++) if(!isnan(dft2->voi[ri].y[si])) break;
337 for(ei=dft2->frameNr-1; ei>=si; ei--) if(!isnan(dft2->voi[ri].y[ei])) break;
338 if((ei-si)>0) {
339 ret=svg_write_tac(fp_svg, &viewports, 1, tac_id, tac_title,
340 dft2->x+si, dft2->voi[ri].y+si, 1+ei-si,
341 svgColorName(color_nr%max_color_nr), symbol_nr%max_symbol_nr, SYMBOLFILLED,
342 NULL, verbose-3);
343 if(ret) {svg_legend_empty(&legends); return(21);}
344 }
345 /* Draw the measured points */
346 ret=svg_write_tac(fp_svg, &viewports, 2, tac_id, tac_title,
347 dft1->x, dft1->voi[ri].y, dft1->frameNr,
348 svgColorName(color_nr%max_color_nr), symbol_nr%max_symbol_nr, SYMBOLFILLED,
349 NULL, verbose-3);
350 if(ret) {svg_legend_empty(&legends); return(22);}
351 /* Set legend too, if requested */
352 if(is_label!=0) {
353 svg_legend_add(&legends, 0, symbol_nr%max_symbol_nr, SYMBOLFILLED,
354 color_nr%max_color_nr, tac_title);
355 }
356 /* Prepare for the next plot */
357 color_nr++; n++;
358 if(color_nr==max_color_nr) {symbol_nr++; color_nr=0;}
359 if(symbol_nr==max_symbol_nr) symbol_nr=0;
360 }
361
362 /* Close the coordinate viewport */
363 ret=svg_end_coordinate_viewport(fp_svg, NULL, verbose-3);
364 if(ret) {svg_legend_empty(&legends); return(91);}
365
366 /* Write the axis ticks */
367 if(svg_write_xticks(fp_svg, &viewports, NULL, verbose-3)!=0) {
368 svg_legend_empty(&legends); return(92);}
369 if(svg_write_yticks(fp_svg, &viewports, NULL, verbose-3)!=0) {
370 svg_legend_empty(&legends); return(93);}
371
372 /* Close the plot viewport */
373 ret=svg_end_plot_viewport(fp_svg, NULL, verbose-3);
374 if(ret) {svg_legend_empty(&legends); return(94);}
375
376 /* Make the plot legends into their own viewport */
377 if(viewports.label_area_viewport.is!=0) {
378 if(verbose>2) printf("creating plot legends\n");
379 ret=svg_create_legends(fp_svg, &viewports, &legends, NULL, verbose-3);
380 if(ret) {svg_legend_empty(&legends); return(95);}
381 }
382 svg_legend_empty(&legends);
383
384 /* Close the SVG file */
385 ret=svg_close(fp_svg, NULL, verbose-3); if(ret) return(101);
386
387 return(0);
388}
backupExistingFile
int backupExistingFile(char *filename, char *backup_ext, char *status)
Definition backup.c:14
dftMinMax
int dftMinMax(DFT *dft, double *minx, double *maxx, double *miny, double *maxy)
Definition dft.c:974
rnameRmDots
int rnameRmDots(char *rname1, char *rname2)
Definition rname.c:99
svg_legend_empty
void svg_legend_empty(SVG_LEGENDS *legends)
Definition svg_legend.c:26
svg_write_xticks
int svg_write_xticks(FILE *fp, struct svg_viewports *vp, char *errmsg, int verbose)
Definition svg_plot.c:555
svg_init_legends
void svg_init_legends(SVG_LEGENDS *legends)
Definition svg_legend.c:14
svg_create_xaxis_title
int svg_create_xaxis_title(FILE *fp, const char *title_text, struct svg_viewports *vp, char *errmsg, int verbose)
Definition svg_title.c:126
svgColorName
char * svgColorName(const svgColor index)
Definition svg_defs.c:38
svg_create_main_title
int svg_create_main_title(FILE *fp, const char *main_title_text, const char *sub_title_text, struct svg_viewports *vp, char *errmsg, int verbose)
Definition svg_title.c:17
svg_legend_add
int svg_legend_add(SVG_LEGENDS *legends, const int plot_type, const int symbol_type, const svgSymbolFill symbol_fill, const int color, const char *text)
Definition svg_legend.c:43
svgSymbolName
char * svgSymbolName(const svgSymbolType index)
Definition svg_defs.c:67
svg_define_viewports
int svg_define_viewports(const int main_viewport_width, const int main_viewport_height, const int is_main_title, const int is_yaxis_title, const int is_xaxis_title, const int is_label_area, struct svg_viewports *vp, int verbose)
Definition svg_vport.c:65
svg_end_plot_viewport
int svg_end_plot_viewport(FILE *fp, char *errmsg, int verbose)
Definition svg_plot.c:238
svg_start_coordinate_viewport
int svg_start_coordinate_viewport(FILE *fp, struct svg_viewports *vp, char *errmsg, int verbose)
Definition svg_plot.c:276
svg_create_legends
int svg_create_legends(FILE *fp, struct svg_viewports *vp, SVG_LEGENDS *legends, char *errmsg, int verbose)
Definition svg_legend.c:76
svg_init_viewports
void svg_init_viewports(struct svg_viewports *p)
Definition svg_vport.c:43
svg_initiate
FILE * svg_initiate(const char *filename, const double height, const double width, struct svg_viewports *vp, char *errmsg, int verbose)
Definition svg_file.c:22
svg_calculate_axes
int svg_calculate_axes(struct svg_viewports *vp, int verbose)
Definition svg_plot.c:364
svg_write_axes
int svg_write_axes(FILE *fp, struct svg_viewports *vp, char *errmsg, int verbose)
Definition svg_plot.c:437
svg_end_coordinate_viewport
int svg_end_coordinate_viewport(FILE *fp, char *errmsg, int verbose)
Definition svg_plot.c:327
svg_start_plot_viewport
int svg_start_plot_viewport(FILE *fp, struct svg_viewports *vp, char *errmsg, int verbose)
Definition svg_plot.c:181
svg_create_yaxis_title
int svg_create_yaxis_title(FILE *fp, const char *title_text, struct svg_viewports *vp, char *errmsg, int verbose)
Definition svg_title.c:197
svg_write_yticks
int svg_write_yticks(FILE *fp, struct svg_viewports *vp, char *errmsg, int verbose)
Definition svg_plot.c:655
svg_write_tac
int svg_write_tac(FILE *fp, struct svg_viewports *vp, const int plot_type, const char *tac_id, const char *tac_title, double *x, double *y, const int data_nr, const char *color, const svgSymbolType symbol_type, const svgSymbolFill symbol_fill, char *errmsg, int verbose)
Definition svg_plot.c:746
svg_close
int svg_close(FILE *fp, char *errmsg, int verbose)
Definition svg_file.c:107
DFT::voi
Voi * voi
Definition libtpccurveio.h:108
DFT::timeunit
int timeunit
Definition libtpccurveio.h:85
DFT::voiNr
int voiNr
Definition libtpccurveio.h:78
DFT::frameNr
int frameNr
Definition libtpccurveio.h:76
DFT::x
double * x
Definition libtpccurveio.h:102
DFT::unit
char unit[MAX_UNITS_LEN+1]
Definition libtpccurveio.h:83
SVG_LEGENDS
Definition libtpcsvg.h:138
Voi::y
double * y
Definition libtpccurveio.h:59
Voi::name
char name[MAX_REGIONNAME_LEN+1]
Definition libtpccurveio.h:48
svg_viewports
Definition libtpcsvg.h:102

◆ plot_fitrange_svg()

int plot_fitrange_svg ( DFT * dft1,
DFT * dft2,
char * main_title,
double x1,
double x2,
double y1,
double y2,
char * fname,
int verbose )

Writes specified range of plots of original and fitted TACs in SVG 1.1 format.

Returns
Returns 0 if successful, otherwise nonzero.
See also
plot_fit_svg, plot_svg
Parameters
dft1Measured data points
dft2Fitted data points. Times can be different but unit must be the same.
main_titleString for plot main title, or NULL
x1Start time; NaN if determined from data
x2End time; NaN if determined from data
y1Minimum y value; NaN if determined from data
y2Maximum y value; NaN if determined from data
fnameSVG filename; existing file is backed up
verboseVerbose level; set to zero to not to print any comments

Definition at line 16 of file plotfit.c.

35 {
36 int ret, n, ri, si, ei;
37 char x_title[64], y_title[64], tac_id[32], tac_title[64];
38 double minx, maxx, miny, maxy, tx1, tx2, ty1, ty2, f;
39 struct svg_viewports viewports; svg_init_viewports(&viewports);
40 int max_color_nr, color_nr;
41 int max_symbol_nr, symbol_nr;
42 SVG_LEGENDS legends; svg_init_legends(&legends);
43 FILE *fp_svg=NULL;
44
45 if(verbose>0) {
46 printf("%s(dft1, dft2, mt, x1, x2, y1, y2, fn, %d)\n", __func__, verbose);
47 }
48
49 /* Check data */
50 if(dft1==NULL || dft1->voiNr<1) return(1);
51 if(dft2==NULL) return(1);
52 if(dft2->voiNr!=dft1->voiNr) return(1);
53
54 /* Check if file exists; backup, if necessary */
55 ret=backupExistingFile(fname, NULL, NULL); if(ret) return(2);
56
57 int is_label=0; if(dft1->voiNr>1) is_label=1;
58
59 /* Determine the plot min and max x values */
60 ret=dftMinMax(dft1, &tx1, &tx2, NULL, NULL); if(ret) return(3);
61 minx=tx1; maxx=tx2;
62 ret=dftMinMax(dft2, &tx1, &tx2, NULL, NULL); if(ret) return(3);
63 if(minx>tx1) minx=tx1;
64 if(maxx<tx2) maxx=tx2;
65 if(minx>0.0) {
66 f=maxx-minx; minx-=0.05*f;
67 if(minx<0.0) minx=0.0;
68 }
69 if(!isnan(x1)) minx=x1;
70 if(!isnan(x2)) maxx=x2;
71
72 /* Determine the plot min and max y values */
73 ret=dftMaxY(dft1, minx, maxx, &ty1, &ty2); if(ret) return(3);
74 miny=ty1; maxy=ty2;
75 ret=dftMaxY(dft2, minx, maxx, &ty1, &ty2); if(ret) return(3);
76 if(miny>ty1) miny=ty1;
77 if(maxy<ty2) maxy=ty2;
78 if(miny>0.0) {
79 f=maxy-miny; miny-=0.05*f;
80 if(miny<0.0) miny=0.0;
81 }
82 if(!isnan(y1)) miny=y1;
83 if(!isnan(y2)) maxy=y2;
84
85 if(verbose>1) printf("minx:=%g\nmaxx:=%g\nminy:=%g\nmaxy:=%g\n", minx, maxx, miny, maxy);
86
87 /* Calculate the axis ticks */
88 viewports.label_area_viewport.is=is_label; // needed for x axis ticks
89 viewports.x.min=minx; viewports.x.max=maxx;
90 viewports.y.min=miny; viewports.y.max=maxy;
91 if(isnan(x1) || isnan(x2)) viewports.x.fixed_min=0;
92 if(isnan(y1) || isnan(y2)) viewports.y.fixed_min=0;
93 else viewports.y.fixed_min=1;
94 ret=svg_calculate_axes(&viewports, verbose-3); if(ret) return(4);
95
96 /* Set x and y axis titles based on activity and time units */
97 strcpy(x_title, "");
98 if(dft1->timeunit==DFTTIME_SEC || dft1->timeunit==DFTTIME_MIN)
99 sprintf(x_title, "Time (%s)", dftTimeunit(dft1->timeunit));
100 else if(dft1->timeunit!=DFTTIME_UNKNOWN)
101 strcpy(x_title, dftTimeunit(dft1->timeunit));
102 strcpy(y_title, "");
103 if(dftUnitId(dft1->unit)!=DFTUNIT_UNKNOWN)
104 strcpy(y_title, dftUnit(dftUnitId(dft1->unit)));
105
106 /* Set the plot window and window area sizes */
107 ret=svg_define_viewports(0, 0, strlen(main_title), strlen(y_title),
108 strlen(x_title), is_label, &viewports, verbose-3);
109 if(ret) return(5);
110
111 /* Initiate graphics file */
112 fp_svg=svg_initiate(fname, 0, 0, &viewports, NULL, verbose-3);
113 if(fp_svg==NULL) return(6);
114
115 /* Put the graph titles into their own viewports */
116 ret=svg_create_main_title(fp_svg, main_title, "", &viewports, NULL,verbose-3);
117 if(ret) return(7);
118 ret=svg_create_yaxis_title(fp_svg, y_title, &viewports, NULL, verbose-3);
119 if(ret) return(8);
120 ret=svg_create_xaxis_title(fp_svg, x_title, &viewports, NULL, verbose-3);
121 if(ret) return(9);
122
123 /* Put the plot into its own viewport */
124 ret=svg_start_plot_viewport(fp_svg, &viewports, NULL, verbose-3);
125 if(ret) return(10);
126
127 /* Start coordinate area viewport */
128 ret=svg_start_coordinate_viewport(fp_svg, &viewports, NULL, verbose-3);
129 if(ret) return(11);
130
131 /* Write plot axes */
132 ret=svg_write_axes(fp_svg, &viewports, NULL, verbose-3);
133 if(ret) return(12);
134
135 /*
136 * Draw the plots
137 */
138 max_color_nr=0; while(svgColorName(max_color_nr)!=NULL) max_color_nr++;
139 if(max_color_nr==0) max_color_nr=1; // remainder works only if 2nd operator>0
140 if(verbose>3) printf("max_color_nr := %d\n", max_color_nr);
141 max_symbol_nr=0; while(svgSymbolName(max_symbol_nr)!=NULL) max_symbol_nr++;
142 if(max_symbol_nr==0) max_symbol_nr=1; // remainder works only if 2nd operator>0
143 if(verbose>3) printf("max_symbol_nr := %d\n", max_symbol_nr);
144 if(dft1->voiNr==1) color_nr=0; else color_nr=1;
145 symbol_nr=0;
146 for(ri=0, n=0; ri<dft1->voiNr; ri++) {
147 sprintf(tac_id, "plot_%d", n);
148/*
149 if(strlen(dft1->studynr)>0 && strcmp(dft1->studynr, ".")!=0)
150 sprintf(tac_title, "%s: %s", dft1->studynr, dft1->voi[ri].name);
151 else strcpy(tac_title, dft1->voi[ri].name);
152*/
153 rnameRmDots(dft1->voi[ri].name, tac_title);
154 /* Draw the fitted line */
155 for(si=0; si<dft2->frameNr; si++) if(!isnan(dft2->voi[ri].y[si])) break;
156 for(ei=dft2->frameNr-1; ei>=si; ei--) if(!isnan(dft2->voi[ri].y[ei])) break;
157 if((ei-si)>0) {
158 ret=svg_write_tac(fp_svg, &viewports, 1, tac_id, tac_title,
159 dft2->x+si, dft2->voi[ri].y+si, 1+ei-si,
160 svgColorName(color_nr%max_color_nr), symbol_nr%max_symbol_nr, SYMBOLFILLED,
161 NULL, verbose-3);
162 if(ret) {svg_legend_empty(&legends); return(21);}
163 }
164 /* Draw the measured points */
165 ret=svg_write_tac(fp_svg, &viewports, 2, tac_id, tac_title,
166 dft1->x, dft1->voi[ri].y, dft1->frameNr,
167 svgColorName(color_nr%max_color_nr), symbol_nr%max_symbol_nr, SYMBOLFILLED,
168 NULL, verbose-3);
169 if(ret) {svg_legend_empty(&legends); return(22);}
170 /* Set legend too, if requested */
171 if(is_label!=0) {
172 svg_legend_add(&legends, 0, symbol_nr%max_symbol_nr, SYMBOLFILLED,
173 color_nr%max_color_nr, tac_title);
174 }
175 /* Prepare for the next plot */
176 color_nr++; n++;
177 if(color_nr==max_color_nr) {symbol_nr++; color_nr=0;}
178 if(symbol_nr==max_symbol_nr) symbol_nr=0;
179 }
180
181 /* Close the coordinate viewport */
182 ret=svg_end_coordinate_viewport(fp_svg, NULL, verbose-3);
183 if(ret) {svg_legend_empty(&legends); return(91);}
184
185 /* Write the axis ticks */
186 if(svg_write_xticks(fp_svg, &viewports, NULL, verbose-3)!=0) {
187 svg_legend_empty(&legends); return(92);}
188 if(svg_write_yticks(fp_svg, &viewports, NULL, verbose-3)!=0) {
189 svg_legend_empty(&legends); return(93);}
190
191 /* Close the plot viewport */
192 ret=svg_end_plot_viewport(fp_svg, NULL, verbose-3);
193 if(ret) {svg_legend_empty(&legends); return(94);}
194
195 /* Make the plot legends into their own viewport */
196 if(viewports.label_area_viewport.is!=0) {
197 if(verbose>2) printf("creating plot legends\n");
198 ret=svg_create_legends(fp_svg, &viewports, &legends, NULL, verbose-3);
199 if(ret) {svg_legend_empty(&legends); return(95);}
200 }
201 svg_legend_empty(&legends);
202
203 /* Close the SVG file */
204 ret=svg_close(fp_svg, NULL, verbose-3); if(ret) return(101);
205
206 return(0);
207}
dftMaxY
int dftMaxY(DFT *dft, double t1, double t2, double *miny, double *maxy)
Definition dft.c:1090
Generated on Fri Sep 19 2025 13:32:07 for TPCCLIB by doxygen 1.12.0