tif_color.c File Reference

#include <precomp.h>
#include <math.h>

Include dependency graph for tif_color.c:

Go to the source code of this file.

Macros
#define	RINT(R)   ((uint32)((R)>0?((R)+0.5):((R)-0.5)))
#define	SHIFT   16
#define	FIX(x)   ((int32)((x) * (1L<<SHIFT) + 0.5))
#define	ONE_HALF   ((int32)(1<<(SHIFT-1)))
#define	Code2V(c, RB, RW, CR)   ((((c)-(int32)(RB))*(float)(CR))/(float)(((RW)-(RB)!=0) ? ((RW)-(RB)) : 1))
#define	CLAMP(f, min, max)   ((f)<(min)?(min):(f)>(max)?(max):(f))
#define	HICLAMP(f, max)   ((f)>(max)?(max):(f))
#define	LumaRed   luma[0]
#define	LumaGreen   luma[1]
#define	LumaBlue   luma[2]

Functions
void	TIFFCIELabToXYZ (TIFFCIELabToRGB *cielab, uint32 l, int32 a, int32 b, float *X, float *Y, float *Z)
void	TIFFXYZToRGB (TIFFCIELabToRGB *cielab, float X, float Y, float Z, uint32 *r, uint32 *g, uint32 *b)
int	TIFFCIELabToRGBInit (TIFFCIELabToRGB *cielab, const TIFFDisplay *display, float *refWhite)
void	TIFFYCbCrtoRGB (TIFFYCbCrToRGB *ycbcr, uint32 Y, int32 Cb, int32 Cr, uint32 *r, uint32 *g, uint32 *b)
static float	CLAMPw (float v, float vmin, float vmax)
int	TIFFYCbCrToRGBInit (TIFFYCbCrToRGB *ycbcr, float *luma, float *refBlackWhite)

Macro Definition Documentation

CLAMP

#define CLAMP	(		f,
			min,
			max
	)		((f)<(min)?(min):(f)>(max)?(max):(f))

Definition at line 177 of file tif_color.c.

Code2V

#define Code2V	(		c,
			RB,
			RW,
			CR
	)		((((c)-(int32)(RB))*(float)(CR))/(float)(((RW)-(RB)!=0) ? ((RW)-(RB)) : 1))

Definition at line 176 of file tif_color.c.

FIX

#define FIX

(

x

)

((int32)((x) * (1L<<SHIFT) + 0.5))

Definition at line 174 of file tif_color.c.

HICLAMP

#define HICLAMP	(		f,
			max
	)		((f)>(max)?(max):(f))

Definition at line 178 of file tif_color.c.

LumaBlue

#define LumaBlue   luma[2]

LumaGreen

#define LumaGreen   luma[1]

LumaRed

#define LumaRed   luma[0]

ONE_HALF

#define ONE_HALF   ((int32)(1<<(SHIFT-1)))

Definition at line 175 of file tif_color.c.

RINT

#define RINT

(

R

)

((uint32)((R)>0?((R)+0.5):((R)-0.5)))

Definition at line 71 of file tif_color.c.

SHIFT

#define SHIFT   16

Definition at line 173 of file tif_color.c.

Function Documentation

CLAMPw()

static float CLAMPw ( float  v, float  vmin, float  vmax  )

static

Definition at line 202 of file tif_color.c.

{
    if( v < vmin )
    {
        /* printf("%f clamped to %f\n", v, vmin); */
        return vmin;
    }
    if( v > vmax )
    {
        /* printf("%f clamped to %f\n", v, vmax); */
        return vmax;
    }
    return v;
}

Referenced by TIFFYCbCrToRGBInit().

TIFFCIELabToRGBInit()

int TIFFCIELabToRGBInit	(	TIFFCIELabToRGB *	cielab,
		const TIFFDisplay *	display,
		float *	refWhite
	)

Definition at line 123 of file tif_color.c.

{
    int i;
    double dfGamma;
 
    cielab->range = CIELABTORGB_TABLE_RANGE;
 
    _TIFFmemcpy(&cielab->display, display, sizeof(TIFFDisplay));
 
    /* Red */
    dfGamma = 1.0 / cielab->display.d_gammaR ;
    cielab->rstep =
        (cielab->display.d_YCR - cielab->display.d_Y0R) / cielab->range;
    for(i = 0; i <= cielab->range; i++) {
        cielab->Yr2r[i] = cielab->display.d_Vrwr
            * ((float)pow((double)i / cielab->range, dfGamma));
    }
 
    /* Green */
    dfGamma = 1.0 / cielab->display.d_gammaG ;
    cielab->gstep =
        (cielab->display.d_YCR - cielab->display.d_Y0R) / cielab->range;
    for(i = 0; i <= cielab->range; i++) {
        cielab->Yg2g[i] = cielab->display.d_Vrwg
            * ((float)pow((double)i / cielab->range, dfGamma));
    }
 
    /* Blue */
    dfGamma = 1.0 / cielab->display.d_gammaB ;
    cielab->bstep =
        (cielab->display.d_YCR - cielab->display.d_Y0R) / cielab->range;
    for(i = 0; i <= cielab->range; i++) {
        cielab->Yb2b[i] = cielab->display.d_Vrwb
            * ((float)pow((double)i / cielab->range, dfGamma));
    }
 
    /* Init reference white point */
    cielab->X0 = refWhite[0];
    cielab->Y0 = refWhite[1];
    cielab->Z0 = refWhite[2];
 
    return 0;
}

Referenced by initCIELabConversion().

TIFFCIELabToXYZ()

void TIFFCIELabToXYZ	(	TIFFCIELabToRGB *	cielab,
		uint32	l,
		int32	a,
		int32	b,
		float *	X,
		float *	Y,
		float *	Z
	)

Definition at line 44 of file tif_color.c.

{
    float L = (float)l * 100.0F / 255.0F;
    float cby, tmp;
 
    if( L < 8.856F ) {
        *Y = (L * cielab->Y0) / 903.292F;
        cby = 7.787F * (*Y / cielab->Y0) + 16.0F / 116.0F;
    } else {
        cby = (L + 16.0F) / 116.0F;
        *Y = cielab->Y0 * cby * cby * cby;
    }
 
    tmp = (float)a / 500.0F + cby;
    if( tmp < 0.2069F )
        *X = cielab->X0 * (tmp - 0.13793F) / 7.787F;
    else    
        *X = cielab->X0 * tmp * tmp * tmp;
 
    tmp = cby - (float)b / 200.0F;
    if( tmp < 0.2069F )
        *Z = cielab->Z0 * (tmp - 0.13793F) / 7.787F;
    else    
        *Z = cielab->Z0 * tmp * tmp * tmp;
}

Referenced by DECLAREContigPutFunc().

TIFFXYZToRGB()

void TIFFXYZToRGB	(	TIFFCIELabToRGB *	cielab,
		float	X,
		float	Y,
		float	Z,
		uint32 *	r,
		uint32 *	g,
		uint32 *	b
	)

Definition at line 76 of file tif_color.c.

{
    int i;
    float Yr, Yg, Yb;
    float *matrix = &cielab->display.d_mat[0][0];
 
    /* Multiply through the matrix to get luminosity values. */
    Yr =  matrix[0] * X + matrix[1] * Y + matrix[2] * Z;
    Yg =  matrix[3] * X + matrix[4] * Y + matrix[5] * Z;
    Yb =  matrix[6] * X + matrix[7] * Y + matrix[8] * Z;
 
    /* Clip input */
    Yr = TIFFmax(Yr, cielab->display.d_Y0R);
    Yg = TIFFmax(Yg, cielab->display.d_Y0G);
    Yb = TIFFmax(Yb, cielab->display.d_Y0B);
 
    /* Avoid overflow in case of wrong input values */
    Yr = TIFFmin(Yr, cielab->display.d_YCR);
    Yg = TIFFmin(Yg, cielab->display.d_YCG);
    Yb = TIFFmin(Yb, cielab->display.d_YCB);
 
    /* Turn luminosity to colour value. */
    i = (int)((Yr - cielab->display.d_Y0R) / cielab->rstep);
    i = TIFFmin(cielab->range, i);
    *r = RINT(cielab->Yr2r[i]);
 
    i = (int)((Yg - cielab->display.d_Y0G) / cielab->gstep);
    i = TIFFmin(cielab->range, i);
    *g = RINT(cielab->Yg2g[i]);
 
    i = (int)((Yb - cielab->display.d_Y0B) / cielab->bstep);
    i = TIFFmin(cielab->range, i);
    *b = RINT(cielab->Yb2b[i]);
 
    /* Clip output. */
    *r = TIFFmin(*r, cielab->display.d_Vrwr);
    *g = TIFFmin(*g, cielab->display.d_Vrwg);
    *b = TIFFmin(*b, cielab->display.d_Vrwb);
}

Referenced by DECLAREContigPutFunc().

TIFFYCbCrtoRGB()

void TIFFYCbCrtoRGB	(	TIFFYCbCrToRGB *	ycbcr,
		uint32	Y,
		int32	Cb,
		int32	Cr,
		uint32 *	r,
		uint32 *	g,
		uint32 *	b
	)

Definition at line 181 of file tif_color.c.

{
    int32 i;
 
    /* XXX: Only 8-bit YCbCr input supported for now */
    Y = HICLAMP(Y, 255);
    Cb = CLAMP(Cb, 0, 255);
    Cr = CLAMP(Cr, 0, 255);
 
    i = ycbcr->Y_tab[Y] + ycbcr->Cr_r_tab[Cr];
    *r = CLAMP(i, 0, 255);
    i = ycbcr->Y_tab[Y]
        + (int)((ycbcr->Cb_g_tab[Cb] + ycbcr->Cr_g_tab[Cr]) >> SHIFT);
    *g = CLAMP(i, 0, 255);
    i = ycbcr->Y_tab[Y] + ycbcr->Cb_b_tab[Cb];
    *b = CLAMP(i, 0, 255);
}

Referenced by DECLARESepPutFunc().

TIFFYCbCrToRGBInit()

int TIFFYCbCrToRGBInit	(	TIFFYCbCrToRGB *	ycbcr,
		float *	luma,
		float *	refBlackWhite
	)

Definition at line 234 of file tif_color.c.

{
    TIFFRGBValue* clamptab;
    int i;
    
#define LumaRed     luma[0]
#define LumaGreen   luma[1]
#define LumaBlue    luma[2]
 
    clamptab = (TIFFRGBValue*)(
    (uint8*) ycbcr+TIFFroundup_32(sizeof (TIFFYCbCrToRGB), sizeof (long)));  
    _TIFFmemset(clamptab, 0, 256);      /* v < 0 => 0 */
    ycbcr->clamptab = (clamptab += 256);
    for (i = 0; i < 256; i++)
    clamptab[i] = (TIFFRGBValue) i;
    _TIFFmemset(clamptab+256, 255, 2*256);  /* v > 255 => 255 */
    ycbcr->Cr_r_tab = (int*) (clamptab + 3*256);
    ycbcr->Cb_b_tab = ycbcr->Cr_r_tab + 256;
    ycbcr->Cr_g_tab = (int32*) (ycbcr->Cb_b_tab + 256);
    ycbcr->Cb_g_tab = ycbcr->Cr_g_tab + 256;
    ycbcr->Y_tab = ycbcr->Cb_g_tab + 256;
 
    { float f1 = 2-2*LumaRed;       int32 D1 = FIX(CLAMP(f1,0.0F,2.0F));
      float f2 = LumaRed*f1/LumaGreen;  int32 D2 = -FIX(CLAMP(f2,0.0F,2.0F));
      float f3 = 2-2*LumaBlue;      int32 D3 = FIX(CLAMP(f3,0.0F,2.0F));
      float f4 = LumaBlue*f3/LumaGreen; int32 D4 = -FIX(CLAMP(f4,0.0F,2.0F));
      int x;
 
#undef LumaBlue
#undef LumaGreen
#undef LumaRed
      
      /*
       * i is the actual input pixel value in the range 0..255
       * Cb and Cr values are in the range -128..127 (actually
       * they are in a range defined by the ReferenceBlackWhite
       * tag) so there is some range shifting to do here when
       * constructing tables indexed by the raw pixel data.
       */
      for (i = 0, x = -128; i < 256; i++, x++) {
        int32 Cr = (int32)CLAMPw(Code2V(x, refBlackWhite[4] - 128.0F,
                refBlackWhite[5] - 128.0F, 127),
                            -128.0F * 32, 128.0F * 32);
        int32 Cb = (int32)CLAMPw(Code2V(x, refBlackWhite[2] - 128.0F,
                refBlackWhite[3] - 128.0F, 127),
                            -128.0F * 32, 128.0F * 32);
 
        ycbcr->Cr_r_tab[i] = (int32)((D1*Cr + ONE_HALF)>>SHIFT);
        ycbcr->Cb_b_tab[i] = (int32)((D3*Cb + ONE_HALF)>>SHIFT);
        ycbcr->Cr_g_tab[i] = D2*Cr;
        ycbcr->Cb_g_tab[i] = D4*Cb + ONE_HALF;
        ycbcr->Y_tab[i] =
            (int32)CLAMPw(Code2V(x + 128, refBlackWhite[0], refBlackWhite[1], 255),
                                  -128.0F * 32, 128.0F * 32);
      }
    }
 
    return 0;
}

Referenced by initYCbCrConversion().