diff options
Diffstat (limited to 'xicc/xfit.h')
| -rw-r--r-- | xicc/xfit.h | 239 |
1 files changed, 239 insertions, 0 deletions
diff --git a/xicc/xfit.h b/xicc/xfit.h new file mode 100644 index 0000000..f9f3233 --- /dev/null +++ b/xicc/xfit.h @@ -0,0 +1,239 @@ + +#ifndef XFIT_H +#define XFIT_H + +/* + * Clut per channel curve fitting + * + * Author: Graeme W. Gill + * Date: 27/5/2007 + * Version: 1.00 + * + * Copyright 2000 - 2007 Graeme W. Gill + * All rights reserved. + * This material is licenced under the GNU AFFERO GENERAL PUBLIC LICENSE Version 3 :- + * see the License.txt file for licencing details. + * + * Based on the rspl.c and xlut.c code. + */ + +#define MXLUORD 60 /* Maximum possible shaper harmonic orders to use */ + +#define PCSD 3 /* PCS dimensions */ + +#define MXPARMS (2 * MXDI * MXLUORD + (1 << MXDI) * MXDO + MXDO * MXLUORD) + +#define MAX_NTP 11 /* Maximum test points within spans */ + +/* Optimisation mask legal combinations */ +typedef enum { + oc_i = 1, /* Input Shaper curve */ + oc_p = 2, /* Input Positioning curve */ + oc_m = 4, /* Temporary matrix */ + oc_o = 8, /* Output */ + + oc_io = oc_i + oc_o, + oc_im = oc_i + oc_m, + oc_mo = oc_m + oc_o, + oc_imo = oc_i + oc_m + oc_o, + + oc_pm = oc_p + oc_m, + oc_pmo = oc_p + oc_m + oc_o, + + oc_ip = oc_i + oc_p, + oc_ipo = oc_i + oc_p + oc_o +} optcomb; + +/* Flag values */ +#define XFIT_VERB 0x0001 /* Verbose output during fitting */ + +#define XFIT_FM_INPUT 0x0002 /* Use input space for fit metric (default is output) */ + +#define XFIT_IN_ZERO 0x0004 /* Adjust input curves 1 & 2 for zero (~~not impd. yet) */ +#define XFIT_OPTGRID_RANGE 0x0008 /* Optimize inner grid around used range (~~not impd. yet) */ + +#define XFIT_OUT_WP_REL 0x0010 /* Extract the white point and make output relative */ +#define XFIT_OUT_WP_REL_US 0x0030 /* Same as above but scale to avoid clipping above WP */ +#define XFIT_OUT_WP_REL_C 0x0050 /* Same as above but clip any cLUT values over D50 */ +#define XFIT_CLIP_WP 0x0080 /* Clip white point to have Y <= 1.0 (conflict with above) */ +#define XFIT_OUT_LAB 0x0100 /* Output space is LAB else XYZ for reading WP */ + +#define XFIT_OUT_ZERO 0x0200 /* Adjust output curves 1 & 2 for zero */ + +#define XFIT_MAKE_CLUT 0x1000 /* Create rspl clut, even if not needed */ + + +/* xfit reverse transform information for each test point */ +typedef struct { + double ide[MXDO][MXDI]; /* pseudo inverse: rout -> in */ +} xfit_piv; + +/* Context for optimising input and output luts */ +struct _xfit { + int verb; /* Verbose */ + int flags; /* Behaviour flags */ + int di, fdi; /* Dimensionaluty of input and output */ + optcomb tcomb; /* Target 1D curve elements to fit */ + + icxMatrixModel *skm; /* Optional skeleton model (used for input profiles) */ + + double *wp; /* Ref. to current white point if XFIT_OUT_WP_REL */ + double *dw; /* Ref. to device value that should map to D50 if XFIT_OUT_WP_REL */ + double fromAbs[3][3]; /* From abs to relative (used by caller) */ + double toAbs[3][3]; /* To abs from relative (used by caller) */ + int gres[MXDI]; /* clut resolutions being optimised for */ + rspl *clut; /* final rspl clut */ + + void *cntx2; /* Context of callback */ + double (*to_de2)(void *cntx, double *in1, double *in2); + /* callback to convert in or out value to fit metric squared */ + double (*to_dde2)(void *cntx, double dout[2][MXDIDO], double *in1, double *in2); + /* Same, but with partial derivatives */ + + int iluord[MXDI]; /* Input Shaper order actualy used (must be <= MXLUORD) */ + int sm_iluord; /* Smallest Input Shaper order used */ + int oluord[MXDO]; /* Output Shaper order actualy used (must be <= MXLUORD) */ + int sluord[MXDI]; /* Sub-grid shaper order */ + double in_min[MXDI]; /* Input value scaling minimum */ + double in_max[MXDI]; /* Input value scaling maximum */ + double out_min[MXDO]; /* Output value scaling minimum */ + double out_max[MXDO]; /* Output value scaling maximum */ + + int shp_off; /* Input parameters offset */ + int shp_offs[MXDI]; /* Input parameter offsets for each in channel from v[0] */ + int shp_cnt; /* Input parameters count */ + int mat_off; /* Matrix parameters offset from v[0] */ + int mat_offs[MXDO]; /* Matrix parameter offsets for each out channel from v[0] */ + int mat_cnt; /* Matrix parameters count */ + int out_off; /* Output parameters offset from v[0] */ + int out_offs[MXDO]; /* Output parameter offsets for each out channel from v[0] */ + int out_cnt; /* Output parameters count */ + int pos_off; /* Position parameters offset */ + int pos_offs[MXDI]; /* Position parameter offsets for each in channel from v[0] */ + int pos_cnt; /* Position parameters count */ + int tot_cnt; /* Total parameter count */ + + double *v; /* Holder for parameters */ + /* Optimisation parameters are layed out: */ + /* */ + /* Input pos or shape curves:, di groups of iluord[e] parameters */ + /* */ + /* (temp) Matrix: fdi groups of 2 ^ di parameters */ + /* */ + /* Output curves:, fdi groups of oluord[f] parameters */ + /* */ + /* Shaper curves:, di groups of iluord[e] parameters */ + /* */ + + int nodp; /* Number of data points */ + cow *ipoints; /* Reference to test points as in->out */ + cow *rpoints; /* Modified version of ipoints */ + xfit_piv *piv; /* Point inverse information for XFIT_FM_INPUT */ + double *uerrv; /* Array holding span width in DE for current opt chan */ + + double mat[3][3]; /* XYZ White point aprox relative to accurate relative matrix */ + double cmat[3][3]; /* Final rspl correction matrix */ + + double shp_smooth[MXDI];/* Smoothing factors for each input shape curve, nom = 1.0 */ + double out_smooth[MXDO]; + + /* Optimisation state */ + optcomb opt_msk; /* Optimisation mask: 3 = i+m, 2 = m, 6 = m+o, 7 = i+m+o */ + int opt_ssch; /* Single shaper channel mode flag */ + int opt_off; /* Optimisation parameters offset from v[0] */ + int opt_cnt; /* Optimisation parameters count */ + double *wv; /* Parameters being optimised */ + double *sa; /* Search area */ + int opt_ch; /* Channel being optimized */ + + /* Methods */ + void (*del)(struct _xfit *p); + + /* Do the fitting. Return nz on error */ + int (*fit)( + struct _xfit *p, + int flags, /* Xfit flag values */ + int di, /* Input dimensions */ + int fdi, /* Output dimensions */ + int rsplflags, /* clut rspl creation flags */ + double *wp, /* if flags & XFIT_OUT_WP_REL, */ + /* Initial white point, returns final wp */ + double *dw, /* Device white value to adjust to be D50 */ + double wpscale, /* If >= 0.0 scale final wp */ + double *dgw, /* Device space gamut boundary white for XFIT_OUT_WP_REL_US */ + /* (ie. RGB 1,1,1 CMYK 0,0,0,0, etc) */ + cow *ipoints, /* Array of data points to fit - referece taken */ + int nodp, /* Number of data points */ + icxMatrixModel *skm, /* Optional skeleton model (used for input profiles) */ + double in_min[MXDI], /* Input value scaling/domain minimum */ + double in_max[MXDI], /* Input value scaling/domain maximum */ + int gres[MXDI], /* clut resolutions being optimised for/returned */ + double out_min[MXDO], /* Output value scaling/range minimum */ + double out_max[MXDO], /* Output value scaling/range maximum */ + double smooth, /* clut rspl smoothing factor */ + double oavgdev[MXDO], /* Average output value deviation */ + int iord[], /* Order of input positioning/shaper curve for each dimension */ + int sord[], /* Order of input sub-grid shaper curve (not used) */ + int oord[], /* Order of output shaper curve for each dimension */ + double shp_smooth[MXDI],/* Smoothing factors for each curve, nom = 1.0 */ + double out_smooth[MXDO], + optcomb tcomb, /* Flag - target elements to fit. */ + void *cntx2, /* Context of callbacks */ + /* Callback to convert two fit values delta E squared */ + double (*to_de2)(void *cntx, double *in1, double *in2), + /* Same as above, with partial derivatives */ + double (*to_dde2)(void *cntx, double dout[2][MXDIDO], double *in1, double *in2) + ); + + /* Lookup a value though a combined input positioning and shaper curves */ + double (*incurve)(struct _xfit *p, double in, int chan); + + /* Inverse Lookup a value though a combined input positioning and shaper curves */ + double (*invincurve)(struct _xfit *p, double in, int chan); + + /* Lookup a value though an output curve */ + double (*outcurve)(struct _xfit *p, double in, int chan); + + /* Inverse Lookup a value though an output curve */ + double (*invoutcurve)(struct _xfit *p, double in, int chan); + +}; typedef struct _xfit xfit; + +xfit *new_xfit(); + +#endif /* XFIT_H */ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + |