2 * Copyright (c) 2011, Mickael Savinaud, Communications & Systemes <mickael.savinaud@c-s.fr>
5 * Redistribution and use in source and binary forms, with or without
6 * modification, are permitted provided that the following conditions
8 * 1. Redistributions of source code must retain the above copyright
9 * notice, this list of conditions and the following disclaimer.
10 * 2. Redistributions in binary form must reproduce the above copyright
11 * notice, this list of conditions and the following disclaimer in the
12 * documentation and/or other materials provided with the distribution.
14 * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS `AS IS'
15 * AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
16 * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
17 * ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
18 * LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
19 * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
20 * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
21 * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
22 * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
23 * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
24 * POSSIBILITY OF SUCH DAMAGE.
30 * Created on: 8 juil. 2011
40 #include "opj_config.h"
41 #include "opj_getopt.h"
44 #include "format_defs.h"
47 double* parseToleranceValues( char* inArg, const int nbcomp);
48 void comparePGXimages_help_display(void);
49 opj_image_t* readImageFromFilePGX(char* filename, int nbFilenamePGX, char *separator);
51 int imageToPNG(const opj_image_t* image, const char* filename, int num_comp_select);
54 typedef struct test_cmp_parameters
60 /** Number of components */
65 double* tabPEAKvalues;
69 char separator_base[2];
71 char separator_test[2];
73 } test_cmp_parameters;
75 /*******************************************************************************
76 * Command line help function
77 *******************************************************************************/
78 void comparePGXimages_help_display(void) {
79 fprintf(stdout,"\nList of parameters for the comparePGX function \n");
81 fprintf(stdout," -b \t REQUIRED \t filename to the reference/baseline PGX image \n");
82 fprintf(stdout," -t \t REQUIRED \t filename to the test PGX image\n");
83 fprintf(stdout," -n \t REQUIRED \t number of component of the image (used to generate correct filename)\n");
84 fprintf(stdout," -m \t OPTIONAL \t list of MSE tolerances, separated by : (size must correspond to the number of component) of \n");
85 fprintf(stdout," -p \t OPTIONAL \t list of PEAK tolerances, separated by : (size must correspond to the number of component) \n");
86 fprintf(stdout," -s \t OPTIONAL \t 1 or 2 filename separator to take into account PGX image with different components, "
87 "please indicate b or t before separator to indicate respectively the separator "
88 "for ref/base file and for test file. \n");
89 fprintf(stdout," -r \t OPTIONAL \t indicate if you want to run this function as conformance test or as non regression test\n");
93 /*******************************************************************************
95 *******************************************************************************/
96 int parse_cmdline_cmp(int argc, char **argv, test_cmp_parameters* param)
98 char *MSElistvalues = NULL; char *PEAKlistvalues= NULL;
99 char *separatorList = NULL;
100 int sizemembasefile, sizememtestfile;
101 int index, flagM=0, flagP=0;
102 const char optlist[] = "b:t:n:m:p:s:d";
106 param->base_filename = NULL;
107 param->test_filename = NULL;
109 param->tabMSEvalues = NULL;
110 param->tabPEAKvalues = NULL;
115 while ((c = opj_getopt(argc, argv, optlist)) != -1)
119 sizemembasefile = (int)strlen(opj_optarg)+1;
120 param->base_filename = (char*) malloc(sizemembasefile);
121 param->base_filename[0] = '\0';
122 strncpy(param->base_filename, opj_optarg, strlen(opj_optarg));
123 param->base_filename[strlen(opj_optarg)] = '\0';
124 /*printf("param->base_filename = %s [%d / %d]\n", param->base_filename, strlen(param->base_filename), sizemembasefile );*/
127 sizememtestfile = (int) strlen(opj_optarg) + 1;
128 param->test_filename = (char*) malloc(sizememtestfile);
129 param->test_filename[0] = '\0';
130 strncpy(param->test_filename, opj_optarg, strlen(opj_optarg));
131 param->test_filename[strlen(opj_optarg)] = '\0';
132 /*printf("param->test_filename = %s [%d / %d]\n", param->test_filename, strlen(param->test_filename), sizememtestfile);*/
135 param->nbcomp = atoi(opj_optarg);
138 MSElistvalues = opj_optarg;
142 PEAKlistvalues = opj_optarg;
149 separatorList = opj_optarg;
152 if ((opj_optopt == 'b') || (opj_optopt == 't') || (opj_optopt == 'n') || (opj_optopt == 'p') || (opj_optopt == 'm') || (opj_optopt
154 fprintf(stderr, "Option -%c requires an argument.\n", opj_optopt);
156 if (isprint(opj_optopt)) fprintf(stderr, "Unknown option `-%c'.\n", opj_optopt);
157 else fprintf(stderr, "Unknown option character `\\x%x'.\n", opj_optopt);
160 fprintf(stderr, "WARNING -> this option is not valid \"-%c %s\"\n", c, opj_optarg);
164 if (opj_optind != argc)
166 for (index = opj_optind; index < argc; index++)
167 fprintf(stderr,"Non-option argument %s\n", argv[index]);
171 if (param->nbcomp == 0)
173 fprintf(stderr,"Need to indicate the number of components !\n");
178 if ( flagM && flagP )
180 param->tabMSEvalues = parseToleranceValues( MSElistvalues, param->nbcomp);
181 param->tabPEAKvalues = parseToleranceValues( PEAKlistvalues, param->nbcomp);
182 if ( (param->tabMSEvalues == NULL) || (param->tabPEAKvalues == NULL))
184 fprintf(stderr,"MSE and PEAK values are not correct (respectively need %d values)\n",param->nbcomp);
194 /* Get separators after corresponding letter (b or t)*/
195 if (separatorList != NULL)
197 if( (strlen(separatorList) ==2) || (strlen(separatorList) ==4) )
199 /* keep original string*/
200 int sizeseplist = (int)strlen(separatorList)+1;
201 char* separatorList2 = (char*)malloc( sizeseplist );
202 separatorList2[0] = '\0';
203 strncpy(separatorList2, separatorList, strlen(separatorList));
204 separatorList2[strlen(separatorList)] = '\0';
205 /*printf("separatorList2 = %s [%d / %d]\n", separatorList2, strlen(separatorList2), sizeseplist);*/
207 if (strlen(separatorList) == 2) /* one separator behind b or t*/
209 char *resultT = NULL;
210 resultT = strtok(separatorList2, "t");
211 if (strlen(resultT) == strlen(separatorList)) /* didn't find t character, try to find b*/
213 char *resultB = NULL;
214 resultB = strtok(resultT, "b");
215 if (strlen(resultB) == 1)
217 param->separator_base[0] = separatorList[1];param->separator_base[1] = '\0';
218 param->separator_test[0] ='\0';
220 else /* not found b*/
222 free(separatorList2);
228 param->separator_base[0] ='\0';
229 param->separator_test[0] = separatorList[1];param->separator_test[1] = '\0';
231 /*printf("sep b = %s [%d] and sep t = %s [%d]\n",param->separator_base, strlen(param->separator_base), param->separator_test, strlen(param->separator_test) );*/
233 else /* == 4 characters we must found t and b*/
235 char *resultT = NULL;
236 resultT = strtok(separatorList2, "t");
237 if (strlen(resultT) == 3) /* found t in first place*/
239 char *resultB = NULL;
240 resultB = strtok(resultT, "b");
241 if (strlen(resultB) == 1) /* found b after t*/
243 param->separator_test[0] = separatorList[1];param->separator_test[1] = '\0';
244 param->separator_base[0] = separatorList[3];param->separator_base[1] = '\0';
246 else /* didn't find b after t*/
248 free(separatorList2);
252 else /* == 2, didn't find t in first place*/
254 char *resultB = NULL;
255 resultB = strtok(resultT, "b");
256 if (strlen(resultB) == 1) /* found b in first place*/
258 param->separator_base[0] = separatorList[1]; param->separator_base[1] = '\0';
259 param->separator_test[0] = separatorList[3]; param->separator_test[1] = '\0';
261 else /* didn't found b in first place => problem*/
263 free(separatorList2);
268 free(separatorList2);
270 else /* wrong number of argument after -s*/
277 if (param->nbcomp == 1)
279 param->separator_base[0] = '\0';
280 param->separator_test[0] = '\0';
284 fprintf(stderr,"If number of component is > 1, we need separator\n");
290 if ( (param->nr_flag) && (flagP || flagM) )
292 fprintf(stderr,"Wrong input parameters list: it is non-regression test or tolerance comparison\n");
295 if ( (!param->nr_flag) && (!flagP || !flagM) )
297 fprintf(stderr,"Wrong input parameters list: it is non-regression test or tolerance comparison\n");
304 /*******************************************************************************
305 * Parse MSE and PEAK input values (
307 *******************************************************************************/
308 double* parseToleranceValues( char* inArg, const int nbcomp)
310 double* outArgs= malloc(nbcomp * sizeof(double));
314 result = strtok( inArg, delims );
316 while( (result != NULL) && (it_comp < nbcomp ))
318 outArgs[it_comp] = atof(result);
319 result = strtok( NULL, delims );
323 if (it_comp != nbcomp)
331 /*******************************************************************************
332 * Create filenames from a filename by used separator and nb components
334 *******************************************************************************/
335 char* createMultiComponentsFilename(const char* inFilename, const int indexF, const char* separator)
338 char *outFilename, *ptr;
342 /*printf("inFilename = %s\n", inFilename);*/
343 if ((ptr = strrchr(inFilename, token)) != NULL)
345 posToken = (int) (strlen(inFilename) - strlen(ptr));
346 /*printf("Position of %c character inside inFilename = %d\n", token, posToken);*/
350 /*printf("Token %c not found\n", token);*/
351 outFilename = (char*)malloc(1);
352 outFilename[0] = '\0';
356 outFilename = (char*)malloc((posToken + 7) * sizeof(char)); /*6*/
358 strncpy(outFilename, inFilename, posToken);
360 outFilename[posToken] = '\0';
362 strcat(outFilename, separator);
364 sprintf(s, "%i", indexF);
365 strcat(outFilename, s);
367 strcat(outFilename, ".pgx");
369 /*printf("outfilename: %s\n", outFilename);*/
372 /*******************************************************************************
374 *******************************************************************************/
375 opj_image_t* readImageFromFilePGX(char* filename, int nbFilenamePGX, char *separator)
378 opj_image_t* image_read = NULL;
379 opj_image_t* image = NULL;
380 opj_cparameters_t parameters;
381 opj_image_cmptparm_t* param_image_read;
384 /* If separator is empty => nb file to read is equal to one*/
385 if ( strlen(separator) == 0 )
388 /* set encoding parameters to default values */
389 opj_set_default_encoder_parameters(¶meters);
390 parameters.decod_format = PGX_DFMT;
391 strncpy(parameters.infile, filename, sizeof(parameters.infile)-1);
394 param_image_read = malloc(nbFilenamePGX * sizeof(opj_image_cmptparm_t));
395 data = malloc(nbFilenamePGX * sizeof(*data));
398 for (it_file = 0; it_file < nbFilenamePGX; it_file++)
400 /* Create the right filename*/
401 char *filenameComponentPGX;
402 if (strlen(separator) == 0)
404 filenameComponentPGX = malloc((strlen(filename) + 1) * sizeof(*filenameComponentPGX));
405 strcpy(filenameComponentPGX, filename);
408 filenameComponentPGX = createMultiComponentsFilename(filename, it_file, separator);
410 /* Read the pgx file corresponding to the component */
411 image_read = pgxtoimage(filenameComponentPGX, ¶meters);
415 fprintf(stderr, "Unable to load pgx file\n");
417 free(param_image_read);
419 for (it_free_data = 0; it_free_data < it_file; it_free_data++) {
420 free(data[it_free_data]);
424 free(filenameComponentPGX);
429 /* Set the image_read parameters*/
430 param_image_read[it_file].x0 = 0;
431 param_image_read[it_file].y0 = 0;
432 param_image_read[it_file].dx = 0;
433 param_image_read[it_file].dy = 0;
434 param_image_read[it_file].h = image_read->comps->h;
435 param_image_read[it_file].w = image_read->comps->w;
436 param_image_read[it_file].bpp = image_read->comps->bpp;
437 param_image_read[it_file].prec = image_read->comps->prec;
438 param_image_read[it_file].sgnd = image_read->comps->sgnd;
441 data[it_file] = malloc(param_image_read[it_file].h * param_image_read[it_file].w * sizeof(int));
442 memcpy(data[it_file], image_read->comps->data, image_read->comps->h * image_read->comps->w * sizeof(int));
445 opj_image_destroy(image_read);
446 free(filenameComponentPGX);
449 image = opj_image_create(nbFilenamePGX, param_image_read, CLRSPC_UNSPECIFIED);
450 for (it_file = 0; it_file < nbFilenamePGX; it_file++)
452 /* Copy data into output image and free memory*/
453 memcpy(image->comps[it_file].data, data[it_file], image->comps[it_file].h * image->comps[it_file].w * sizeof(int));
458 free(param_image_read);
464 /*******************************************************************************
466 *******************************************************************************/
468 int imageToPNG(const opj_image_t* image, const char* filename, int num_comp_select)
470 opj_image_cmptparm_t param_image_write;
471 opj_image_t* image_write = NULL;
473 param_image_write.x0 = 0;
474 param_image_write.y0 = 0;
475 param_image_write.dx = 0;
476 param_image_write.dy = 0;
477 param_image_write.h = image->comps[num_comp_select].h;
478 param_image_write.w = image->comps[num_comp_select].w;
479 param_image_write.bpp = image->comps[num_comp_select].bpp;
480 param_image_write.prec = image->comps[num_comp_select].prec;
481 param_image_write.sgnd = image->comps[num_comp_select].sgnd;
483 image_write = opj_image_create(1, ¶m_image_write, CLRSPC_GRAY);
484 memcpy(image_write->comps->data, image->comps[num_comp_select].data, param_image_write.h * param_image_write.w * sizeof(int));
486 imagetopng(image_write, filename);
488 opj_image_destroy(image_write);
494 /*******************************************************************************
496 *******************************************************************************/
497 int main(int argc, char **argv)
499 test_cmp_parameters inParam;
502 int nbFilenamePGXbase, nbFilenamePGXtest;
503 char *filenamePNGtest= NULL, *filenamePNGbase = NULL, *filenamePNGdiff = NULL;
504 int memsizebasefilename, memsizetestfilename, memsizedifffilename;
505 int valueDiff = 0, nbPixelDiff = 0;
506 double sumDiff = 0.0;
507 /* Structures to store image parameters and data*/
508 opj_image_t *imageBase = NULL, *imageTest = NULL, *imageDiff = NULL;
509 opj_image_cmptparm_t* param_image_diff;
511 /* Get parameters from command line*/
512 if( parse_cmdline_cmp(argc, argv, &inParam) == EXIT_FAILURE )
514 comparePGXimages_help_display();
515 if (inParam.tabMSEvalues) free(inParam.tabMSEvalues);
516 if (inParam.tabPEAKvalues) free(inParam.tabPEAKvalues);
517 if (inParam.base_filename) free(inParam.base_filename);
518 if (inParam.test_filename) free(inParam.test_filename);
522 /* Display Parameters*/
523 printf("******Parameters********* \n");
524 printf(" base_filename = %s\n"
525 " test_filename = %s\n"
526 " nb of Components = %d\n"
527 " Non regression test = %d\n"
528 " separator Base = %s\n"
529 " separator Test = %s\n",
530 inParam.base_filename, inParam.test_filename, inParam.nbcomp,
531 inParam.nr_flag, inParam.separator_base, inParam.separator_test);
533 if ( (inParam.tabMSEvalues != NULL) && (inParam.tabPEAKvalues != NULL))
535 printf(" MSE values = [");
536 for (it_comp = 0; it_comp < inParam.nbcomp; it_comp++)
537 printf(" %f ", inParam.tabMSEvalues[it_comp]);
539 printf(" PEAK values = [");
540 for (it_comp = 0; it_comp < inParam.nbcomp; it_comp++)
541 printf(" %f ", inParam.tabPEAKvalues[it_comp]);
543 printf(" Non-regression test = %d\n", inParam.nr_flag);
546 if (strlen(inParam.separator_base) == 0)
547 nbFilenamePGXbase = 0;
549 nbFilenamePGXbase = inParam.nbcomp;
551 if (strlen(inParam.separator_test) == 0)
552 nbFilenamePGXtest = 0;
554 nbFilenamePGXtest = inParam.nbcomp;
556 printf(" NbFilename to generate from base filename = %d\n", nbFilenamePGXbase);
557 printf(" NbFilename to generate from test filename = %d\n", nbFilenamePGXtest);
558 printf("************************* \n");
560 /*----------BASELINE IMAGE--------*/
561 memsizebasefilename = (int)strlen(inParam.test_filename) + 1 + 5 + 2 + 4;
562 memsizetestfilename = (int)strlen(inParam.test_filename) + 1 + 5 + 2 + 4;
564 imageBase = readImageFromFilePGX( inParam.base_filename, nbFilenamePGXbase, inParam.separator_base);
565 if ( imageBase != NULL)
567 filenamePNGbase = (char*) malloc(memsizebasefilename);
568 filenamePNGbase[0] = '\0';
569 strncpy(filenamePNGbase, inParam.test_filename, strlen(inParam.test_filename));
570 filenamePNGbase[strlen(inParam.test_filename)] = '\0';
571 strcat(filenamePNGbase, ".base");
572 /*printf("filenamePNGbase = %s [%d / %d octets]\n",filenamePNGbase, strlen(filenamePNGbase),memsizebasefilename );*/
576 if (inParam.tabMSEvalues) free(inParam.tabMSEvalues);
577 if (inParam.tabPEAKvalues) free(inParam.tabPEAKvalues);
578 if (inParam.base_filename) free(inParam.base_filename);
579 if (inParam.test_filename) free(inParam.test_filename);
583 /*----------TEST IMAGE--------*/
585 imageTest = readImageFromFilePGX(inParam.test_filename, nbFilenamePGXtest, inParam.separator_test);
586 if ( imageTest != NULL)
588 filenamePNGtest = (char*) malloc(memsizetestfilename);
589 filenamePNGtest[0] = '\0';
590 strncpy(filenamePNGtest, inParam.test_filename, strlen(inParam.test_filename));
591 filenamePNGtest[strlen(inParam.test_filename)] = '\0';
592 strcat(filenamePNGtest, ".test");
593 /*printf("filenamePNGtest = %s [%d / %d octets]\n",filenamePNGtest, strlen(filenamePNGtest),memsizetestfilename );*/
597 if (imageBase) opj_image_destroy(imageBase);
598 if (inParam.tabMSEvalues) free(inParam.tabMSEvalues);
599 if (inParam.tabPEAKvalues) free(inParam.tabPEAKvalues);
600 if (inParam.base_filename) free(inParam.base_filename);
601 if (inParam.test_filename) free(inParam.test_filename);
602 free(filenamePNGbase);
606 /*----------DIFF IMAGE--------*/
609 param_image_diff = malloc( imageBase->numcomps * sizeof(opj_image_cmptparm_t));
611 /* Comparison of header parameters*/
612 printf("Step 1 -> Header comparison\n");
614 for (it_comp = 0; it_comp < imageBase->numcomps; it_comp++)
616 param_image_diff[it_comp].x0 = 0;
617 param_image_diff[it_comp].y0 = 0;
618 param_image_diff[it_comp].dx = 0;
619 param_image_diff[it_comp].dy = 0;
621 if (imageBase->comps[it_comp].sgnd != imageTest->comps[it_comp].sgnd)
623 printf("ERROR: sign mismatch [comp %d] (%d><%d)\n", it_comp, ((imageBase->comps)[it_comp]).sgnd, ((imageTest->comps)[it_comp]).sgnd);
627 param_image_diff[it_comp].sgnd = 0 ;
629 if (((imageBase->comps)[it_comp]).prec != ((imageTest->comps)[it_comp]).prec)
631 printf("ERROR: prec mismatch [comp %d] (%d><%d)\n", it_comp, ((imageBase->comps)[it_comp]).prec, ((imageTest->comps)[it_comp]).prec);
635 param_image_diff[it_comp].prec = 8 ;
637 if (((imageBase->comps)[it_comp]).bpp != ((imageTest->comps)[it_comp]).bpp)
639 printf("ERROR: byte per pixel mismatch [comp %d] (%d><%d)\n", it_comp, ((imageBase->comps)[it_comp]).bpp, ((imageTest->comps)[it_comp]).bpp);
643 param_image_diff[it_comp].bpp = 1 ;
645 if (((imageBase->comps)[it_comp]).h != ((imageTest->comps)[it_comp]).h)
647 printf("ERROR: height mismatch [comp %d] (%d><%d)\n", it_comp, ((imageBase->comps)[it_comp]).h, ((imageTest->comps)[it_comp]).h);
651 param_image_diff[it_comp].h = imageBase->comps[it_comp].h ;
653 if (((imageBase->comps)[it_comp]).w != ((imageTest->comps)[it_comp]).w)
655 printf("ERROR: width mismatch [comp %d] (%d><%d)\n", it_comp, ((imageBase->comps)[it_comp]).w, ((imageTest->comps)[it_comp]).w);
659 param_image_diff[it_comp].w = imageBase->comps[it_comp].w ;
662 /* If only one parameter is different, we stop the test*/
665 free(inParam.tabMSEvalues);
666 free(inParam.tabPEAKvalues);
667 free(inParam.base_filename);
668 free(inParam.test_filename);
670 free(filenamePNGbase);
671 free(filenamePNGtest);
673 opj_image_destroy(imageBase);
674 opj_image_destroy(imageTest);
676 free(param_image_diff);
681 imageDiff = opj_image_create(imageBase->numcomps, param_image_diff, CLRSPC_UNSPECIFIED);
683 free(param_image_diff);
685 /* Measurement computation*/
686 printf("Step 2 -> measurement comparison\n");
688 memsizedifffilename = strlen(inParam.test_filename) + 1 + 5 + 2 + 4;
689 filenamePNGdiff = (char*) malloc(memsizedifffilename);
690 filenamePNGdiff[0] = '\0';
691 strncpy(filenamePNGdiff, inParam.test_filename, strlen(inParam.test_filename));
692 filenamePNGdiff[strlen(inParam.test_filename)] = '\0';
693 strcat(filenamePNGdiff, ".diff");
694 /*printf("filenamePNGdiff = %s [%d / %d octets]\n",filenamePNGdiff, strlen(filenamePNGdiff),memsizedifffilename );*/
696 /* Compute pixel diff*/
697 for (it_comp = 0; it_comp < imageDiff->numcomps; it_comp++)
701 char *filenamePNGbase_it_comp, *filenamePNGtest_it_comp, *filenamePNGdiff_it_comp;
703 filenamePNGbase_it_comp = (char*) malloc(memsizebasefilename);
704 filenamePNGbase_it_comp[0] = '\0';
705 strncpy(filenamePNGbase_it_comp,filenamePNGbase,strlen(filenamePNGbase));
706 filenamePNGbase_it_comp[strlen(filenamePNGbase)] = '\0';
708 filenamePNGtest_it_comp = (char*) malloc(memsizetestfilename);
709 filenamePNGtest_it_comp[0] = '\0';
710 strncpy(filenamePNGtest_it_comp,filenamePNGtest,strlen(filenamePNGtest));
711 filenamePNGtest_it_comp[strlen(filenamePNGtest)] = '\0';
713 filenamePNGdiff_it_comp = (char*) malloc(memsizedifffilename);
714 filenamePNGdiff_it_comp[0] = '\0';
715 strncpy(filenamePNGdiff_it_comp,filenamePNGdiff,strlen(filenamePNGdiff));
716 filenamePNGdiff_it_comp[strlen(filenamePNGdiff)] = '\0';
718 for (itpxl = 0; itpxl < ((imageDiff->comps)[it_comp]).w * ((imageDiff->comps)[it_comp]).h; itpxl++)
720 if (abs( ((imageBase->comps)[it_comp]).data[itpxl] - ((imageTest->comps)[it_comp]).data[itpxl] ) > 0)
722 valueDiff = ((imageBase->comps)[it_comp]).data[itpxl] - ((imageTest->comps)[it_comp]).data[itpxl];
723 ((imageDiff->comps)[it_comp]).data[itpxl] = abs(valueDiff);
724 sumDiff += (double)valueDiff;
727 SE += (double)(valueDiff * valueDiff);
728 PEAK = (PEAK > abs(valueDiff)) ? PEAK : abs(valueDiff);
731 ((imageDiff->comps)[it_comp]).data[itpxl] = 0;
734 MSE = SE / ( ((imageDiff->comps)[it_comp]).w * ((imageDiff->comps)[it_comp]).h );
736 if (!inParam.nr_flag && (inParam.tabMSEvalues != NULL) && (inParam.tabPEAKvalues != NULL))
737 { /* Conformance test*/
738 printf("<DartMeasurement name=\"PEAK_%d\" type=\"numeric/double\"> %f </DartMeasurement> \n", it_comp, PEAK);
739 printf("<DartMeasurement name=\"MSE_%d\" type=\"numeric/double\"> %f </DartMeasurement> \n", it_comp, MSE);
741 if ( (MSE > inParam.tabMSEvalues[it_comp]) || (PEAK > inParam.tabPEAKvalues[it_comp]) )
743 printf("ERROR: MSE (%f) or PEAK (%f) values produced by the decoded file are greater "
744 "than the allowable error (respectively %f and %f) \n",
745 MSE, PEAK, inParam.tabMSEvalues[it_comp], inParam.tabPEAKvalues[it_comp]);
749 else /* Non regression-test */
751 if ( nbPixelDiff > 0)
756 printf("<DartMeasurement name=\"NumberOfPixelsWithDifferences_%d\" type=\"numeric/int\"> %d </DartMeasurement> \n", it_comp, nbPixelDiff);
757 printf("<DartMeasurement name=\"ComponentError_%d\" type=\"numeric/double\"> %f </DartMeasurement> \n", it_comp, sumDiff);
760 sprintf(it_compc, "_%i", it_comp);
761 strcat(it_compc,".png");
762 strcat(filenamePNGbase_it_comp, it_compc);
763 /*printf("filenamePNGbase_it = %s [%d / %d octets]\n",filenamePNGbase_it_comp, strlen(filenamePNGbase_it_comp),memsizebasefilename );*/
764 strcat(filenamePNGtest_it_comp, it_compc);
765 /*printf("filenamePNGtest_it = %s [%d / %d octets]\n",filenamePNGtest_it_comp, strlen(filenamePNGtest_it_comp),memsizetestfilename );*/
766 strcat(filenamePNGdiff_it_comp, it_compc);
767 /*printf("filenamePNGdiff_it = %s [%d / %d octets]\n",filenamePNGdiff_it_comp, strlen(filenamePNGdiff_it_comp),memsizedifffilename );*/
769 if ( imageToPNG(imageBase, filenamePNGbase_it_comp, it_comp) == EXIT_SUCCESS )
771 printf("<DartMeasurementFile name=\"BaselineImage_%d\" type=\"image/png\"> %s </DartMeasurementFile> \n", it_comp, filenamePNGbase_it_comp);
774 if ( imageToPNG(imageTest, filenamePNGtest_it_comp, it_comp) == EXIT_SUCCESS )
776 printf("<DartMeasurementFile name=\"TestImage_%d\" type=\"image/png\"> %s </DartMeasurementFile> \n", it_comp, filenamePNGtest_it_comp);
779 if ( imageToPNG(imageDiff, filenamePNGdiff_it_comp, it_comp) == EXIT_SUCCESS )
781 printf("<DartMeasurementFile name=\"DiffferenceImage_%d\" type=\"image/png\"> %s </DartMeasurementFile> \n", it_comp, filenamePNGdiff_it_comp);
787 free(filenamePNGbase_it_comp);
788 free(filenamePNGtest_it_comp);
789 free(filenamePNGdiff_it_comp);
792 /*-----------------------------*/
794 opj_image_destroy(imageBase);
795 opj_image_destroy(imageTest);
796 opj_image_destroy(imageDiff);
798 free(filenamePNGbase);
799 free(filenamePNGtest);
800 free(filenamePNGdiff);
802 free(inParam.tabMSEvalues);
803 free(inParam.tabPEAKvalues);
804 free(inParam.base_filename);
805 free(inParam.test_filename);
811 printf("---- TEST SUCCEED ----\n");