昨天看书的心得

昨天做图像算法，总觉得YIQ色彩和RGB之间通过相互转换就能得到颜色的灰度值。

如果得到一个RGB的灰度值，就能进行灰度变换。图像增强基本都是基于灰度做出来的。

 void RGB2YIQ(REAL *y, REAL *i, REAL*q, BYTE r, BYTE g, BYTE b);

 void YIQ2RGB(BYTE *r, BYTE *g, BYTE *b, REAL y, REAL i, REAL q);

REAL RGB2YIQMat[3][3] =

        {{ 0.30f,  0.59f,  0.11f},

         { 0.60f, -0.28f, -0.32f},

         { 0.21f, -0.52f,  0.31f}};

REAL YIQ2RGBMat[3][3] = 

        {{ 1.0000f,  0.9483f,  0.6240f},

         { 1.0000f, -0.2761f, -0.6398f},

         { 1.0000f, -1.1055f,  1.7299f}};

/*Converts r, g, b, into y, i, q*/

void CImpactPlus::RGB2YIQ(REAL *y, REAL *i, REAL*q, BYTE r, BYTE g, BYTE b)

{

    REAL rT, gT, bT;

    rT = r / 255.0f; gT = g / 255.0f; bT = b / 255.0f;

    *y = RGB2YIQMat[0][0] * rT + RGB2YIQMat[0][1] * gT + RGB2YIQMat[0][2] * bT; 

    *i = RGB2YIQMat[1][0] * rT + RGB2YIQMat[1][1] * gT + RGB2YIQMat[1][2] * bT; 

    *q = RGB2YIQMat[2][0] * rT + RGB2YIQMat[2][1] * gT + RGB2YIQMat[2][2] * bT;

    *i = (*i + 0.6f) / 1.2f;

    *q = (*q + 0.52f) / 1.04f;

    if (*y < 0.0) *y = 0.0f;

    if (*y > 1.0) *y = 1.0f;

    if (*i < 0.0) *i = 0.0f;

    if (*i > 1.0) *i = 1.0f;

    if (*q < 0.0) *q = 0.0f;

    if (*q > 1.0) *q = 1.0f;

}

/*Converts y, i, q, into r, g, b*/

void CImpactPlus::YIQ2RGB(BYTE *r, BYTE *g, BYTE *b, REAL y, REAL i, REAL q)

{

    REAL rT, gT, bT;

    i = i * 1.2f - 0.6f;

    q = q * 1.04f - 0.52f;

    rT = YIQ2RGBMat[0][0] * y + YIQ2RGBMat[0][1] * i + YIQ2RGBMat[0][2] * q; 

    gT = YIQ2RGBMat[1][0] * y + YIQ2RGBMat[1][1] * i + YIQ2RGBMat[1][2] * q; 

    bT = YIQ2RGBMat[2][0] * y + YIQ2RGBMat[2][1] * i + YIQ2RGBMat[2][2] * q;

    if (rT < 0.0f) rT = 0.0f;

    if (rT > 1.0f) rT = 1.0f;

    if (gT < 0.0f) gT = 0.0f;

    if (gT > 1.0f) gT = 1.0f;

    if (bT < 0.0f) bT = 0.0f;

    if (bT > 1.0f) bT = 1.0f;

    *r = max(0, min(255, (BYTE)(rT * 255.0f)));

    *g = max(0, min(255, (BYTE)(gT * 255.0f)));

    *b = max(0, min(255, (BYTE)(bT * 255.0f)));

}

在数字图像处理的书中，经常有这样的程序：

  PixelValue = (BYTE)(0.299*(float)lpDIBBits[p+2] + 0.587*(float)lpDIBBits[p+1]+0.114*(float)lpDIBBits[p]+0.1;

这行程序的意思是：得到当前像素的YIQ之Y值。

而lpDIBBits[p+2]就是红色成分， lpDIBBits[p+1]是绿色成分， lpDIBBits[p]是蓝色成分。

唯一不明白的是他为什么要加0.1？

下面是我依照他的原理来实现图像的线性灰度变换。

Status CImpactPlus::LinerTrans(Bitmap* pBitmap, float fa, float fb)

{

 Status s;

 TestGdiPointer(_T("HistogramEqual"), pBitmap);

 

 
 INT nWidth = pBitmap->GetWidth();

 INT nHeight = pBitmap->GetHeight();
 BitmapData bitmapData;

 Rect rt(0, 0, nWidth, nHeight);

 if(StatusBad(s = pBitmap->LockBits(&rt,

  ImageLockModeRead | ImageLockModeWrite,

  PixelFormat32bppARGB, &bitmapData)))

 

 return SGTRACE_ERR(_T("LockBits"), s);

 ARGB* pixels = (ARGB*)bitmapData.Scan0;

 

 float fy, fi, fq;

 fy = fi = fq = 0.0f;

 BYTE a, r, g, b;

 for(int row=rt.Y; row<rt.Height; row++)

 {

  for(int col=rt.X; col<rt.Width; col++)

  {

   ARGB c = pixels[row * bitmapData.Stride / 4 + col];

   Color color(c);

   a = color.GetA();

   r = color.GetR();

   g = color.GetG();

   b = color.GetB();

   

   RGB2YIQ(&fy, &fi, &fq, r, g, b);

   float fTemp = fa * fy + fb;

   YIQ2RGB(&r, &g, &b, fTemp, fi, fq);

   

   pixels[row * bitmapData.Stride / 4 + col] = color.MakeARGB(a,r,g,b);

  }

 }

 if(StatusBad(s = pBitmap->UnlockBits(&bitmapData)))

  return SGTRACE_ERR(_T("UnlockBits"), s);

 return Ok;

}

这样方法简单明快，即使是很复杂的算法，实现起来也很清晰。