/* $Xorg: LabMxC.c,v 1.3 2000/08/17 19:44:40 cpqbld Exp $ */ /* * Code and supporting documentation (c) Copyright 1990 1991 Tektronix, Inc. * All Rights Reserved * * This file is a component of an X Window System-specific implementation * of XCMS based on the TekColor Color Management System. Permission is * hereby granted to use, copy, modify, sell, and otherwise distribute this * software and its documentation for any purpose and without fee, provided * that this copyright, permission, and disclaimer notice is reproduced in * all copies of this software and in supporting documentation. TekColor * is a trademark of Tektronix, Inc. * * Tektronix makes no representation about the suitability of this software * for any purpose. It is provided "as is" and with all faults. * * TEKTRONIX DISCLAIMS ALL WARRANTIES APPLICABLE TO THIS SOFTWARE, * INCLUDING THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A * PARTICULAR PURPOSE. IN NO EVENT SHALL TEKTRONIX BE LIABLE FOR ANY * SPECIAL, INDIRECT OR CONSEQUENTIAL DAMAGES OR ANY DAMAGES WHATSOEVER * RESULTING FROM LOSS OF USE, DATA, OR PROFITS, WHETHER IN AN ACTION OF * CONTRACT, NEGLIGENCE, OR OTHER TORTIOUS ACTION, ARISING OUT OF OR IN * CONNECTION WITH THE USE OR THE PERFORMANCE OF THIS SOFTWARE. * * * * NAME * CIELabMxC.c * * DESCRIPTION * Source for the XcmsCIELabQueryMaxC() gamut boundary * querying routine. * */ /* $XFree86: xc/lib/X11/LabMxC.c,v 1.3 2001/01/17 19:41:39 dawes Exp $ */ #ifdef HAVE_CONFIG_H #include <config.h> #endif #include "Xlibint.h" #include "Xcmsint.h" #include <math.h> #include "Cv.h" /* * DEFINES */ #define MAXBISECTCOUNT 100 #define EPS (XcmsFloat)0.001 #define START_CHROMA (XcmsFloat)3.6 #define TOPL (XcmsFloat)100.0 /************************************************************************ * * * PUBLIC ROUTINES * * * ************************************************************************/ /* * NAME * XcmsCIELabQueryMaxC - max chroma for a hue_angle and L_star * * SYNOPSIS */ Status XcmsCIELabQueryMaxC(ccc, hue_angle, L_star, pColor_return) XcmsCCC ccc; XcmsFloat hue_angle; /* hue angle in degrees */ XcmsFloat L_star; XcmsColor *pColor_return; /* * DESCRIPTION * Return the maximum chroma for a specific hue_angle and L_star. * The returned format is in XcmsCIELabFormat. * * * ASSUMPTIONS * This routine assumes that the white point associated with * the color specification is the Screen White Point. The * Screen White Point will also be associated with the * returned color specification. * * RETURNS * XcmsFailure - Failure * XcmsSuccess - Succeeded * */ { XcmsCCCRec myCCC; XcmsColor tmp; XcmsColor max_lc; XcmsFloat n_L_star, last_L_star, prev_L_star; XcmsFloat hue, lastaStar, lastbStar, /*lastChroma,*/ maxDist, nT, rFactor; XcmsRGBi rgb_saved; int nCount, nMaxCount; /* * Check Arguments */ if (ccc == NULL || pColor_return == NULL) { return(XcmsFailure); } /* Use my own CCC and inherit screen white Pt */ memcpy ((char *)&myCCC, (char *)ccc, sizeof(XcmsCCCRec)); myCCC.clientWhitePt.format = XcmsUndefinedFormat; myCCC.gamutCompProc = (XcmsCompressionProc)NULL;/* no gamut comp func */ while (hue_angle < 0.0) { hue_angle += 360.0; } while (hue_angle >= 360.0) { hue_angle -= 360.0; } hue = radians(hue_angle); tmp.spec.CIELab.L_star = L_star; tmp.spec.CIELab.a_star = XCMS_CIEASTAROFHUE(hue, START_CHROMA); tmp.spec.CIELab.b_star = XCMS_CIEBSTAROFHUE(hue, START_CHROMA); tmp.pixel = pColor_return->pixel; tmp.format = XcmsCIELabFormat; /* Step 1: compute the maximum L_star and chroma for this hue. */ memcpy((char *)&max_lc, (char *)&tmp, sizeof(XcmsColor)); if (_XcmsCIELabQueryMaxLCRGB(&myCCC, hue, &max_lc, &rgb_saved) == XcmsFailure) { return(XcmsFailure); } /* * Step 2: Do a bisection here to compute the maximum chroma * Note the differences between when the point to be found * is above the maximum LC point and when it is below. */ if (L_star <= max_lc.spec.CIELab.L_star) { maxDist = max_lc.spec.CIELab.L_star; } else { maxDist = TOPL - max_lc.spec.CIELab.L_star; } n_L_star = L_star; last_L_star = -1.0; nMaxCount = MAXBISECTCOUNT; rFactor = 1.0; for (nCount = 0; nCount < nMaxCount; nCount++) { prev_L_star = last_L_star; last_L_star = tmp.spec.CIELab.L_star; /* lastChroma = XCMS_CIELAB_PMETRIC_CHROMA(tmp.spec.CIELab.a_star, */ /* tmp.spec.CIELab.b_star); */ lastaStar = tmp.spec.CIELab.a_star; lastbStar = tmp.spec.CIELab.b_star; nT = (n_L_star - max_lc.spec.CIELab.L_star) / maxDist * rFactor; if (nT > 0) { tmp.spec.RGBi.red = rgb_saved.red * (1.0 - nT) + nT; tmp.spec.RGBi.green = rgb_saved.green * (1.0 - nT) + nT; tmp.spec.RGBi.blue = rgb_saved.blue * (1.0 - nT) + nT; } else { tmp.spec.RGBi.red = rgb_saved.red + (rgb_saved.red * nT); tmp.spec.RGBi.green = rgb_saved.green + (rgb_saved.green * nT); tmp.spec.RGBi.blue = rgb_saved.blue + (rgb_saved.blue * nT); } tmp.format = XcmsRGBiFormat; /* convert from RGB to CIELab */ if (_XcmsConvertColorsWithWhitePt(&myCCC, &tmp, ScreenWhitePointOfCCC(&myCCC), 1, XcmsCIELabFormat, (Bool *) NULL) == XcmsFailure) { return(XcmsFailure); } /* * Now check if we've reached the target L_star */ /* printf("result Lstar = %lf\n", tmp.spec.CIELab.L_star); */ if (tmp.spec.CIELab.L_star <= L_star + EPS && tmp.spec.CIELab.L_star >= L_star - EPS) { memcpy((char *)pColor_return, (char *)&tmp, sizeof(XcmsColor)); return(XcmsSuccess); } if (nT > 0) { n_L_star += ((TOPL - n_L_star) * (L_star - tmp.spec.CIELab.L_star)) / (TOPL - L_star); } else { n_L_star *= L_star / tmp.spec.CIELuv.L_star; } if (tmp.spec.CIELab.L_star <= prev_L_star + EPS && tmp.spec.CIELab.L_star >= prev_L_star - EPS) { rFactor *= 0.5; /* selective relaxation employed */ /* printf("rFactor = %lf\n", rFactor); */ } } if (XCMS_FABS(last_L_star - L_star) < XCMS_FABS(tmp.spec.CIELab.L_star - L_star)) { tmp.spec.CIELab.a_star = lastaStar; tmp.spec.CIELab.b_star = lastbStar; /* tmp.spec.CIELab.a_star = XCMS_CIEASTAROFHUE(hue, lastChroma); */ /* tmp.spec.CIELab.b_star = XCMS_CIEBSTAROFHUE(hue, lastChroma); */ } tmp.spec.CIELab.L_star = L_star; memcpy((char *)pColor_return, (char *)&tmp, sizeof(XcmsColor)); return(XcmsSuccess); }