From 22f703cab05b7cd368f4de9e03991b7664dc5022 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?J=C3=B6rg=20Frings-F=C3=BCrst?= Date: Mon, 1 Sep 2014 13:56:46 +0200 Subject: Initial import of argyll version 1.5.1-8 --- xicc/icheck.c | 532 ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 532 insertions(+) create mode 100644 xicc/icheck.c (limited to 'xicc/icheck.c') diff --git a/xicc/icheck.c b/xicc/icheck.c new file mode 100644 index 0000000..ed3c6e5 --- /dev/null +++ b/xicc/icheck.c @@ -0,0 +1,532 @@ + +/* + * Argyll. + * + * Check for B2A table PCS->Device interpolation faults + * + * Author: Graeme W. Gill + * Date: 2000/12/11 + * Version: 1.00 + * + * Copyright 2000 Graeme W. Gill + * Please refer to License.txt file for details. + */ + +/* TTBD: + * + */ + +#include +#include +#include +#include +#include +#include +#include "numlib.h" +#include "copyright.h" +#include "aconfig.h" +#include "icc.h" + +void usage(void) { + fprintf(stderr,"Check PCS->Device Interpolation faults of ICC file, Version %s\n",ARGYLL_VERSION_STR); + fprintf(stderr,"Author: Graeme W. Gill, licensed under the AGPL Version 3\n"); + fprintf(stderr,"usage: icheck [-v] [-w] infile\n"); + fprintf(stderr," -v verbose\n"); + fprintf(stderr," -w create VRML visualisation\n"); + fprintf(stderr," -x Use VRML axies\n"); + exit(1); +} + +FILE *start_vrml(char *name, int doaxes); +void start_line_set(FILE *wrl); +void add_vertex(FILE *wrl, double pp[3]); +void make_lines(FILE *wrl, int ppset); +void end_vrml(FILE *wrl); + +int +main( + int argc, + char *argv[] +) { + int fa,nfa; /* argument we're looking at */ + int verb = 0; + int dovrml = 0; + int doaxes = 0; + char in_name[100]; + char out_name[100], *xl; + icmFile *rd_fp; + icc *rd_icco; + int rv = 0; + + /* Check variables */ + icmLuBase *luof, *luob; /* A2B and B2A table lookups */ + icmLuLut *lluof, *lluob; /* Lookup Lut type object */ + int gres; /* Grid resolution of B2A */ + icColorSpaceSignature ins, outs; /* Type of input and output spaces */ + int inn; /* Number of input chanels */ + icmLuAlgType alg; + FILE *wrl = NULL; + + error_program = argv[0]; + + if (argc < 2) + usage(); + + /* Process the arguments */ + for(fa = 1;fa < argc;fa++) { + nfa = fa; /* skip to nfa if next argument is used */ + if (argv[fa][0] == '-') { /* Look for any flags */ + char *na = NULL; /* next argument after flag, null if none */ + + if (argv[fa][2] != '\000') + na = &argv[fa][2]; /* next is directly after flag */ + else { + if ((fa+1) < argc) { + if (argv[fa+1][0] != '-') { + nfa = fa + 1; + na = argv[nfa]; /* next is seperate non-flag argument */ + } + } + } + + /* Verbosity */ + if (argv[fa][1] == 'v' || argv[fa][1] == 'V') { + verb = 1; + } + /* VRML */ + else if (argv[fa][1] == 'w' || argv[fa][1] == 'W') { + dovrml = 1; + } + /* Axes */ + else if (argv[fa][1] == 'x' || argv[fa][1] == 'X') { + doaxes = 1; + } + else if (argv[fa][1] == '?') + usage(); + else + usage(); + } + else + break; + } + + if (fa >= argc || argv[fa][0] == '-') usage(); + strcpy(in_name,argv[fa]); + + strcpy(out_name, in_name); + if ((xl = strrchr(out_name, '.')) == NULL) /* Figure where extention is */ + xl = out_name + strlen(out_name); + strcpy(xl,".wrl"); + + /* Open up the file for reading */ + if ((rd_fp = new_icmFileStd_name(in_name,"r")) == NULL) + error ("Read: Can't open file '%s'",in_name); + + if ((rd_icco = new_icc()) == NULL) + error ("Read: Creation of ICC object failed"); + + /* Read the header and tag list */ + if ((rv = rd_icco->read(rd_icco,rd_fp,0)) != 0) + error ("Read: %d, %s",rv,rd_icco->err); + + /* Get a Device to PCS conversion object */ + if ((luof = rd_icco->get_luobj(rd_icco, icmFwd, icAbsoluteColorimetric, icSigLabData, icmLuOrdNorm)) == NULL) { + if ((luof = rd_icco->get_luobj(rd_icco, icmFwd, icmDefaultIntent, icSigLabData, icmLuOrdNorm)) == NULL) + error ("%d, %s",rd_icco->errc, rd_icco->err); + } + + /* Get a PCS to Device conversion object */ + if ((luob = rd_icco->get_luobj(rd_icco, icmBwd, icAbsoluteColorimetric, icSigLabData, icmLuOrdNorm)) == NULL) { + if ((luob = rd_icco->get_luobj(rd_icco, icmBwd, icmDefaultIntent, icSigLabData, icmLuOrdNorm)) == NULL) + error ("%d, %s",rd_icco->errc, rd_icco->err); + } + + /* Get details of conversion (for B2A direction) */ + luob->spaces(luob, &outs, NULL, &ins, &inn, &alg, NULL, NULL, NULL, NULL); + + if (alg != icmLutType) { + error("Expecting Lut based profile"); + } + + if (outs != icSigLabData) { + error("Expecting Lab PCS"); + } + + lluof = (icmLuLut *)luof; /* Lookup Lut type object */ + lluob = (icmLuLut *)luob; /* Lookup Lut type object */ + + gres = lluob->lut->clutPoints; + + if (dovrml) { + wrl = start_vrml(out_name, doaxes); + start_line_set(wrl); + } + + { + double aerr = 0.0; + double ccount = 0.0; + double merr = 0.0; + double tcount = 0.0; + int co[3]; /* PCS grid counter */ + + /* Itterate throught the PCS clut grid cells */ + for (co[2] = 0; co[2] < (gres-1); co[2]++) { + for (co[1] = 0; co[1] < (gres-1); co[1]++) { + for (co[0] = 0; co[0] < (gres-1); co[0]++) { + int j, k, m; + int cc[3]; /* Cube corner offsets */ + double pcs[8][3], wpcsd; + double apcs[3]; + double adev[MAX_CHAN]; + double check[3]; /* Check PCS */ + double ier; /* Interpolation error */ + + apcs[0] = apcs[1] = apcs[2] = 0.0; + for (k = 0; k < inn; k++) + adev[k] = 0.0; + + /* For each corner of the PCS grid based at the current point, */ + /* average the PCS and Device values */ + m = 0; + for (cc[2] = 0; cc[2] < 2; cc[2]++, m++) { + for (cc[1] = 0; cc[1] < 2; cc[1]++) { + for (cc[0] = 0; cc[0] < 2; cc[0]++) { + double dev[MAX_CHAN]; + + pcs[m][0] = (co[0] + cc[0])/(gres - 1.0); + pcs[m][1] = (co[1] + cc[1])/(gres - 1.0); + pcs[m][2] = (co[2] + cc[2])/(gres - 1.0); + + /* Match icclib settable() range */ + pcs[m][0] = pcs[m][0] * 100.0; + pcs[m][1] = (pcs[m][1] * 254.0) - 127.0; + pcs[m][2] = (pcs[m][2] * 254.0) - 127.0; + +//printf("Input PCS %f %f %f\n", pcs[m][0], pcs[m][1], pcs[m][2]); + + /* PCS to (cliped) Device */ + if ((rv = lluob->clut(lluob, dev, pcs[m])) > 1) + error ("%d, %s",rd_icco->errc,rd_icco->err); + + /* (clipped) Device to (clipped) PCS */ + if ((rv = lluof->clut(lluof, pcs[m], dev)) > 1) + error ("%d, %s",rd_icco->errc,rd_icco->err); + + apcs[0] += pcs[m][0]; + apcs[1] += pcs[m][1]; + apcs[2] += pcs[m][2]; + +//printf("Corner PCS %f %f %f -> %f %f %f %f\n", +//pcs[m][0], pcs[m][1], pcs[m][2], dev[0], dev[1], dev[2], dev[3]); + + for (k = 0; k < inn; k++) + adev[k] += dev[k]; + } + } + } + + for (j = 0; j < 3; j++) + apcs[j] /= 8.0; + + for (k = 0; k < inn; k++) + adev[k] /= 8.0; + + /* Compute worst case distance of PCS corners to average PCS */ + wpcsd = 0.0; + for (m = 0; m < 8; m++) { + double ss; + for (ss = 0.0, j = 0; j < 3; j++) { + double tt = pcs[m][j] - apcs[j]; + ss += tt * tt; + } + ss = sqrt(ss); + if (ss > wpcsd) + wpcsd = ss; + } + wpcsd *= 0.75; /* Set threshold at 75% of most distant corner */ + /* Set a worst case */ + if (wpcsd < 1.0) + wpcsd = 1.0; + +// else if (wpcsd > 3.0) +// wpcsd = 3.0; + + +//printf("Average PCS %f %f %f, Average Device %f %f %f %f\n", +//apcs[0], apcs[1], apcs[2], adev[0], adev[1], adev[2], adev[3]); + + /* Average Device to PCS */ + if ((rv = lluof->clut(lluof, check, adev)) > 1) + error ("%d, %s",rd_icco->errc,rd_icco->err); + +//printf("Check PCS %f %f %f\n", +//check[0], check[1], check[2]); + + /* Compute error in PCS vs. Device interpolation */ + for (ier = 0.0, j = 0; j < 3; j++) { + double tt = apcs[j] - check[j]; + ier += tt * tt; + } + ier = sqrt(ier); + +//printf("Average PCS %f %f %f, Check PCS %f %f %f, error %f\n", +//apcs[0], apcs[1], apcs[2], check[0], check[1], check[2], ier); + + aerr += ier; + ccount++; + if (ier > merr) + merr = ier; + + if (ier > wpcsd) { + tcount++; + + printf("ier = %f, Dev = %f %f %f %f\n", + ier, adev[0], adev[1], adev[2], adev[3]); + if (dovrml) { + add_vertex(wrl, apcs); + add_vertex(wrl, check); + } + } + +//printf("~1 ier = %f\n",ier); +//printf("\n"); + + + if (verb) + printf("."), fflush(stdout); + } + } + } + + if (dovrml) { + make_lines(wrl, 2); + end_vrml(wrl); + } + + aerr /= ccount; + + printf("Average interpolation error %f, maximum %f\n",aerr, merr); + printf("Number outside corner radius = %f%%\n",tcount * 100.0/ccount); + } + + /* Done with lookup objects */ + luof->del(luof); + luob->del(luob); + + rd_icco->del(rd_icco); + rd_fp->del(rd_fp); + + return 0; +} + +/* ------------------------------------------------ */ +/* Some simple functions to do basix VRML work */ + +#define GAMUT_LCENT 50.0 +static int npoints = 0; +static int paloc = 0; +static struct { double pp[3]; } *pary; + +static void Lab2RGB(double *out, double *in); + +FILE *start_vrml(char *name, int doaxes) { + FILE *wrl; + struct { + double x, y, z; + double wx, wy, wz; + double r, g, b; + } axes[5] = { + { 0, 0, 50-GAMUT_LCENT, 2, 2, 100, .7, .7, .7 }, /* L axis */ + { 50, 0, 0-GAMUT_LCENT, 100, 2, 2, 1, 0, 0 }, /* +a (red) axis */ + { 0, -50, 0-GAMUT_LCENT, 2, 100, 2, 0, 0, 1 }, /* -b (blue) axis */ + { -50, 0, 0-GAMUT_LCENT, 100, 2, 2, 0, 1, 0 }, /* -a (green) axis */ + { 0, 50, 0-GAMUT_LCENT, 2, 100, 2, 1, 1, 0 }, /* +b (yellow) axis */ + }; + int i; + + if ((wrl = fopen(name,"w")) == NULL) + error("Error opening VRML file '%s'\n",name); + + npoints = 0; + + fprintf(wrl,"#VRML V2.0 utf8\n"); + fprintf(wrl,"\n"); + fprintf(wrl,"# Created by the Argyll CMS\n"); + fprintf(wrl,"Transform {\n"); + fprintf(wrl,"children [\n"); + fprintf(wrl," NavigationInfo {\n"); + fprintf(wrl," type \"EXAMINE\" # It's an object we examine\n"); + fprintf(wrl," } # We'll add our own light\n"); + fprintf(wrl,"\n"); + fprintf(wrl," DirectionalLight {\n"); + fprintf(wrl," direction 0 0 -1 # Light illuminating the scene\n"); + fprintf(wrl," direction 0 -1 0 # Light illuminating the scene\n"); + fprintf(wrl," }\n"); + fprintf(wrl,"\n"); + fprintf(wrl," Viewpoint {\n"); + fprintf(wrl," position 0 0 340 # Position we view from\n"); + fprintf(wrl," }\n"); + fprintf(wrl,"\n"); + if (doaxes != 0) { + fprintf(wrl,"# Lab axes as boxes:\n"); + for (i = 0; i < 5; i++) { + fprintf(wrl,"Transform { translation %f %f %f\n", axes[i].x, axes[i].y, axes[i].z); + fprintf(wrl,"\tchildren [\n"); + fprintf(wrl,"\t\tShape{\n"); + fprintf(wrl,"\t\t\tgeometry Box { size %f %f %f }\n", + axes[i].wx, axes[i].wy, axes[i].wz); + fprintf(wrl,"\t\t\tappearance Appearance { material Material "); + fprintf(wrl,"{ diffuseColor %f %f %f} }\n", axes[i].r, axes[i].g, axes[i].b); + fprintf(wrl,"\t\t}\n"); + fprintf(wrl,"\t]\n"); + fprintf(wrl,"}\n"); + } + fprintf(wrl,"\n"); + } + + return wrl; +} + +void +start_line_set(FILE *wrl) { + + fprintf(wrl,"\n"); + fprintf(wrl,"Shape {\n"); + fprintf(wrl," geometry IndexedLineSet { \n"); + fprintf(wrl," coord Coordinate { \n"); + fprintf(wrl," point [\n"); +} + +void add_vertex(FILE *wrl, double pp[3]) { + + fprintf(wrl,"%f %f %f,\n",pp[1], pp[2], pp[0]-GAMUT_LCENT); + + if (paloc < (npoints+1)) { + paloc = (paloc + 10) * 2; + if (pary == NULL) + pary = malloc(paloc * 3 * sizeof(double)); + else + pary = realloc(pary, paloc * 3 * sizeof(double)); + + if (pary == NULL) + error ("Malloc failed"); + } + pary[npoints].pp[0] = pp[0]; + pary[npoints].pp[1] = pp[1]; + pary[npoints].pp[2] = pp[2]; + npoints++; +} + + +void make_lines(FILE *wrl, int ppset) { + int i, j; + + fprintf(wrl," ]\n"); + fprintf(wrl," }\n"); + fprintf(wrl," coordIndex [\n"); + + for (i = 0; i < npoints;) { + for (j = 0; j < ppset; j++, i++) { + fprintf(wrl,"%d, ", i); + } + fprintf(wrl,"-1,\n"); + } + fprintf(wrl," ]\n"); + + /* Color */ + fprintf(wrl," colorPerVertex TRUE\n"); + fprintf(wrl," color Color {\n"); + fprintf(wrl," color [ # RGB colors of each vertex\n"); + + for (i = 0; i < npoints; i++) { + double rgb[3], Lab[3]; + Lab[0] = pary[i].pp[0]; + Lab[1] = pary[i].pp[1]; + Lab[2] = pary[i].pp[2]; + Lab2RGB(rgb, Lab); + fprintf(wrl," %f %f %f,\n", rgb[0], rgb[1], rgb[2]); + } + fprintf(wrl," ] \n"); + fprintf(wrl," }\n"); + /* End color */ + + fprintf(wrl," }\n"); + fprintf(wrl,"} # end shape\n"); + +} + +void end_vrml(FILE *wrl) { + + fprintf(wrl,"\n"); + fprintf(wrl," ] # end of children for world\n"); + fprintf(wrl,"}\n"); + + if (fclose(wrl) != 0) + error("Error closing VRML file\n"); +} + + +/* Convert a gamut Lab value to an RGB value for display purposes */ +static void +Lab2RGB(double *out, double *in) { + double L = in[0], a = in[1], b = in[2]; + double x,y,z,fx,fy,fz; + double R, G, B; + + /* Scale so that black is visible */ + L = L * (100 - 40.0)/100.0 + 40.0; + + /* First convert to XYZ using D50 white point */ + if (L > 8.0) { + fy = (L + 16.0)/116.0; + y = pow(fy,3.0); + } else { + y = L/903.2963058; + fy = 7.787036979 * y + 16.0/116.0; + } + + fx = a/500.0 + fy; + if (fx > 24.0/116.0) + x = pow(fx,3.0); + else + x = (fx - 16.0/116.0)/7.787036979; + + fz = fy - b/200.0; + if (fz > 24.0/116.0) + z = pow(fz,3.0); + else + z = (fz - 16.0/116.0)/7.787036979; + + x *= 0.9642; /* Multiply by white point, D50 */ + y *= 1.0; + z *= 0.8249; + + /* Now convert to sRGB values */ + R = x * 3.2410 + y * -1.5374 + z * -0.4986; + G = x * -0.9692 + y * 1.8760 + z * 0.0416; + B = x * 0.0556 + y * -0.2040 + z * 1.0570; + + if (R < 0.0) + R = 0.0; + else if (R > 1.0) + R = 1.0; + + if (G < 0.0) + G = 0.0; + else if (G > 1.0) + G = 1.0; + + if (B < 0.0) + B = 0.0; + else if (B > 1.0) + B = 1.0; + + R = pow(R, 1.0/2.2); + G = pow(G, 1.0/2.2); + B = pow(B, 1.0/2.2); + + out[0] = R; + out[1] = G; + out[2] = B; +} + -- cgit v1.2.3