gronod/argyllcms
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases 11 Wiki Activity
Files
development
argyllcms/link/pathplot.c
T
gronod 45bdf51201 Initial commit of v3.5.0
2026-08-20 20:28:40 +01:00

286 lines
7.5 KiB
C
Raw Permalink Blame History

/*
* Argyll Color Management System
*
* Plot the L in vs. L out for a path through the
* input device space, for a given link.
*
* We are assuming RGB device input space !!!!
*
* Author: Graeme W. Gill
* Date: 2002/1/28
* Version: 1.00
*
* Copyright 2002 Graeme W. Gill
*
* This material is licenced under the GNU AFFERO GENERAL PUBLIC LICENSE Version 3 :-
* see the License.txt file for licencing details.
*/
/* TTBD:
*
*/
#include <stdio.h>
#include <stdlib.h>
#include <stdarg.h>
#include <fcntl.h>
#include <string.h>
#include <math.h>
#include "copyright.h"
#include "aconfig.h"
#include "numlib.h"
#include "icc.h"
#include "conv.h"
#include "xicc.h"
#include "plot.h"
#include "ui.h"
#define PRES 100
static double start[13][4] = {
{ 0.0, 0.0, 0.0, 0.0},
{ 0.0, 0.0, 0.0, 0.0},
{ 0.0, 0.0, 0.0, 0.0},
{ 0.0, 0.0, 0.0, 0.0},
{ 0.0, 0.0, 0.0, 0.0},
{ 0.0, 0.0, 0.0, 0.0},
{ 0.0, 0.0, 0.0, 0.0},
{ 1.0, 0.0, 0.0, 0.0},
{ 0.0, 1.0, 0.0, 0.0},
{ 0.0, 0.0, 1.0, 0.0},
{ 0.0, 1.0, 1.0, 0.0},
{ 1.0, 0.0, 1.0, 0.0},
{ 1.0, 1.0, 0.0, 0.0}
};
static double end[13][4] = {
{ 1.0, 0.0, 0.0, 0.0},
{ 0.0, 1.0, 0.0, 0.0},
{ 0.0, 0.0, 1.0, 0.0},
{ 0.0, 1.0, 1.0, 0.0},
{ 1.0, 0.0, 1.0, 0.0},
{ 1.0, 1.0, 0.0, 0.0},
{ 1.0, 1.0, 1.0, 0.0},
{ 1.0, 1.0, 1.0, 0.0},
{ 1.0, 1.0, 1.0, 0.0},
{ 1.0, 1.0, 1.0, 0.0},
{ 1.0, 1.0, 1.0, 0.0},
{ 1.0, 1.0, 1.0, 0.0},
{ 1.0, 1.0, 1.0, 0.0}
};
static char *name[13] = {
"Black - Red", "Black - Green", "Black - Blue",
"Black - Cyan", "Black - Magenta", "Black - Yellow",
"Grey",
"Red - White", "Green - White", "Blue - White",
"Cyan - White", "Magenta - White", "Yellow - White",
};
#define USE_JAB /* Assume CRT in (2), print out (0) */
/* ---------------------------------------- */
void usage(void) {
fprintf(stderr,"Plot device space path L in/out curve from an ICC link file, Version %s\n",ARGYLL_VERSION_STR);
fprintf(stderr,"Author: Graeme W. Gill\n");
fprintf(stderr,"usage: pathplot inprof linkprof outprof\n");
fprintf(stderr," -v verbose\n");
exit(1);
}
int
main(
int argc,
char *argv[]
) {
int fa, nfa; /* argument we're looking at */
int verb = 0;
static char in_name[MAXNAMEL+1];
static char link_name[MAXNAMEL+1];
static char out_name[MAXNAMEL+1];
icmFile *in_fp, *link_fp, *out_fp;
icmErr err = { 0, { '\000'} };
icc *in_icco, *link_icco, *out_icco;
xicc *in_xicco, *out_xicco;
icmLuSpace *link_lu;
icxLuBase *in_lu, *out_lu;
icColorSpaceSignature pcsor; /* PCS to use */
icxViewCond ivc[1], ovc[1];
int rv = 0;
error_program = argv[0];
if (argc < 4)
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;
}
else if (argv[fa][1] == '?')
usage();
else
usage();
}
else
break;
}
if (fa >= argc || argv[fa][0] == '-') usage();
strncpy(in_name,argv[fa++],MAXNAMEL); in_name[MAXNAMEL] = '\000';
if (fa >= argc || argv[fa][0] == '-') usage();
strncpy(link_name,argv[fa++],MAXNAMEL); link_name[MAXNAMEL] = '\000';
if (fa >= argc || argv[fa][0] == '-') usage();
strncpy(out_name,argv[fa++],MAXNAMEL); out_name[MAXNAMEL] = '\000';
/* Open up the files for reading */
if ((in_fp = new_icmFileStd_name(&err,in_name,"r")) == NULL)
error ("Read: Can't open file '%s' (0x%x, '%s')",in_name,err.c,err.m);
if ((in_icco = new_icc(&err)) == NULL)
error ("Read: Creation of ICC object failed (0x%x, '%s')",err.c,err.m);
if ((rv = in_icco->read(in_icco,in_fp,0)) != 0)
error ("Read: %d, %s",rv,in_icco->e.m);
if ((in_xicco = new_xicc(in_icco)) == NULL)
error ("Creation of input profile xicc failed");
if ((link_fp = new_icmFileStd_name(&err,link_name,"r")) == NULL)
error ("Read: Can't open file '%s' (0x%x, '%s')",link_name,err.c,err.m);
if ((link_icco = new_icc(&err)) == NULL)
error ("Read: Creation of ICC object failed (0x%x, '%s')",err.c,err.m);
if ((rv = link_icco->read(link_icco,link_fp,0)) != 0)
error ("Read: %d, %s",rv,link_icco->e.m);
if ((out_fp = new_icmFileStd_name(&err,out_name,"r")) == NULL)
error ("Read: Can't open file '%s' (0x%x, '%s')",out_name,err.c,err.m);
if ((out_icco = new_icc(&err)) == NULL)
error ("Read: Creation of ICC object failed (0x%x, '%s')",err.c,err.m);
if ((rv = out_icco->read(out_icco,out_fp,0)) != 0)
error ("Read: %d, %s",rv,out_icco->e.m);
if ((out_xicco = new_xicc(out_icco)) == NULL)
error ("Creation of output profile xicc failed");
#ifdef USE_JAB
pcsor = icxSigJabData; /* Use CIECAM as PCS */
if (xicc_enum_viewcond(in_xicco, ivc, -2, "mt", 0, NULL) == -999) /* Set input at monitor in typical */
error ("%d, %s",in_xicco->e.c, in_xicco->e.m);
if (xicc_enum_viewcond(out_xicco, ovc, -2, "pp", 0, NULL) == -999) /* Set output at practical reflection print */
error ("%d, %s",out_xicco->e.c, out_xicco->e.m);
#else
pcsor = icSigLabData; /* Default use Lab as PCS */
#endif /* !USE_JAB */
/* Device to PCS conversion object */
if ((in_lu = in_xicco->get_luobj(in_xicco, ICX_CLIP_NEAREST, icmFwd, icAbsoluteColorimetric, pcsor, icmLuOrdNorm, ivc, NULL)) == NULL) {
if ((in_lu = in_xicco->get_luobj(in_xicco, ICX_CLIP_NEAREST, icmBwd, icmDefaultIntent, pcsor, icmLuOrdNorm, ivc, NULL)) == NULL)
error ("%d, %s",in_xicco->e.c, in_xicco->e.m);
}
/* Get a Device to Device conversion object */
if ((link_lu = (icmLuSpace *)link_icco->get_luobj(link_icco, icmFwd, icmDefaultIntent, pcsor, icmLuOrdNorm)) == NULL)
error ("%d, %s",link_icco->e.c, link_icco->e.m);
/* Get a Device to PCS conversion object */
if ((out_lu = out_xicco->get_luobj(out_xicco, ICX_CLIP_NEAREST, icmFwd, icAbsoluteColorimetric, pcsor, icmLuOrdNorm, ovc, NULL)) == NULL) {
if ((out_lu = out_xicco->get_luobj(out_xicco, ICX_CLIP_NEAREST, icmFwd, icmDefaultIntent, pcsor, icmLuOrdNorm, ovc, NULL)) == NULL)
error ("%d, %s",out_xicco->e.c, out_xicco->e.m);
}
{
double xx[PRES], yy[PRES];
int k, i, j;
double tt[10], tt2[10];
for (k = 0; k < 13; k++) {
printf("Doing %s\n",name[k]);
for (i = 0; i < PRES; i++) {
double frac = i/(PRES-1.0);
for (j = 0; j < 4; j++)
tt[j] = (end[k][j] - start[k][j]) * frac + start[k][j];
if (verb)
printf(" %f %f %f -> ",tt[0],tt[1],tt[2]);
/* input device space to PCS */
if ((rv = in_lu->lookup(in_lu, tt2, tt)) > 1)
error ("%d, %s",in_icco->e.c,in_icco->e.m);
xx[i] = tt2[0]; /* L value */
/* input device space to output device space */
if ((rv = link_lu->lookup_fwd(link_lu, tt, tt)) & icmPe_lurv_err)
error ("%d, %s",link_icco->e.c,link_icco->e.m);
/* output device space to PCS */
if ((rv = out_lu->lookup(out_lu, tt, tt)) > 1)
error ("%d, %s",out_icco->e.c,out_icco->e.m);
yy[i] = tt[0]; /* L value */
if (verb)
printf("%f %f %f\n",tt[0],tt[1],tt[2]);
}
if (do_plot(xx,yy,NULL,NULL,PRES) < 0)
error("do_plot returned -1!\n");
}
}
/* Done with lookup objects */
in_lu->del(in_lu);
link_lu->del(link_lu);
out_lu->del(out_lu);
in_icco->del(in_icco);
in_fp->del(in_fp);
link_icco->del(link_icco);
link_fp->del(link_fp);
out_icco->del(out_icco);
out_fp->del(out_fp);
return 0;
}
Reference in New Issue View Git Blame Copy Permalink