RAWRGB to RGB32 ColoSpace

[wangjianhong0574]

i use your RAW to RGB32 the speed is very fast.Thank you for your great and selfless work. But in some cases the effect is not good, I have the following C language image quality is good. I hope it will be useful for you.

this is better
/*

• BayerConversion

• ---

• function: fairly basic interpolation scheme creating RGB values from surrounding R,G,B in bayer

•       image.bordercases are handled by coloring them black
  

• in: pointer to 8bpp Bayer buffer
• out: pointer to 32bpp RGB buffer
• special: using the following 3x3 kernel interpolation:

•       R G R     R      B G B     B
  

•        \|/.     |       \|/      |
  

•       G-B-G   B-G1-B   G-R-G   R-G2-R
  

•        /|\      |       /|\      |
  

•       R G R     R      B G B     B
  

/
/
G R
B G
*/

#define B_OFFSET 0
#define G_OFFSET 1
#define R_OFFSET 2

#define ODD(i) ((i)&1)
#define EVEN(i) (!((i)&1))
typedef enum {
Pattern_Recognition_Filter,
Laroche_Prescott_Filter,
Fast_Linear_Filter,
Anisotropic_Filter
}Bayer_Filter;

void BayerConversionGRBG(unsigned char *pIn, unsigned char pOut, int nRealWidth, int nRealHeight, Bayer_Filter BayerFilter)
{
int yWidth = 0, ym1Width = 0, yp1Width;
int x, y;
int pix;
/
* border cases (fill it in with black)
*/
/G R G R G R
B G B G B G
G R G R G R
B G B G B G/

for (y = 0; y<nRealHeight; y++)
{
	x = 0;
	pOut[((x + y*nRealWidth) << 2) + B_OFFSET] = 0;
	pOut[((x + y*nRealWidth) << 2) + G_OFFSET] = 0;
	pOut[((x + y*nRealWidth) << 2) + R_OFFSET] = 0;
	x = nRealWidth - 1;
	pOut[((x + y*nRealWidth) << 2) + B_OFFSET] = 0;
	pOut[((x + y*nRealWidth) << 2) + G_OFFSET] = 0;
	pOut[((x + y*nRealWidth) << 2) + R_OFFSET] = 0;
}

for (x = 0; x<nRealWidth; x++)
{
	y = 0;
	pOut[((x + y*nRealWidth) << 2) + B_OFFSET] = 0;
	pOut[((x + y*nRealWidth) << 2) + G_OFFSET] = 0;
	pOut[((x + y*nRealWidth) << 2) + R_OFFSET] = 0;
	y = nRealHeight - 1;
	pOut[((x + y*nRealWidth) << 2) + B_OFFSET] = 0;
	pOut[((x + y*nRealWidth) << 2) + G_OFFSET] = 0;
	pOut[((x + y*nRealWidth) << 2) + R_OFFSET] = 0;
}




if (BayerFilter == Anisotropic_Filter)
{
	for (y = 1; y<nRealHeight - 1; y++)
	{
		yWidth += nRealWidth;
		ym1Width = yWidth - nRealWidth;
		yp1Width = yWidth + nRealWidth;

		for (x = 1; x<nRealWidth - 1; x++)
		{
			pix = ((x + yWidth) << 2);
			if (ODD(y))
			if (ODD(x))
			{
				pOut[pix + B_OFFSET] = ((pIn[x - 1 + yWidth] +
					pIn[x + 1 + yWidth]) >> 1);
				pOut[pix + G_OFFSET] = ((pIn[x + yWidth] << 2) + (pIn[x - 1 + ym1Width] + pIn[x + 1 + ym1Width] + pIn[x - 1 + yp1Width] + pIn[x + 1 + yp1Width])) >> 3;
				pOut[pix + R_OFFSET] = ((pIn[x + ym1Width] +
					pIn[x + yp1Width]) >> 1);
			}
			else
			{   // ODD(y) EVEN(x)
				pOut[pix + B_OFFSET] = pIn[x + yWidth];
				if (abs(pIn[x - 1 + yWidth] - pIn[x + 1 + yWidth])<abs(pIn[x + ym1Width] - pIn[x + yp1Width]))
				{
					pOut[pix + G_OFFSET] = (pIn[x - 1 + yWidth] + pIn[x + 1 + yWidth]) >> 1;
				}
				else if (abs(pIn[x - 1 + yWidth] - pIn[x + 1 + yWidth])>abs(pIn[x + ym1Width] - pIn[x + yp1Width]))
				{
					pOut[pix + G_OFFSET] = (pIn[x + ym1Width] + pIn[x + yp1Width]) >> 1;
				}
				else
					pOut[pix + G_OFFSET] = ((pIn[x - 1 + yWidth] +
					pIn[x + 1 + yWidth] +
					pIn[x + ym1Width] +
					pIn[x + yp1Width]) >> 2);

				if (abs(pIn[x - 1 + ym1Width] - pIn[x + 1 + yp1Width])<abs(pIn[x + 1 + ym1Width] - pIn[x - 1 + yp1Width]))
				{
					pOut[pix + R_OFFSET] = (pIn[x - 1 + ym1Width] + pIn[x + 1 + yp1Width]) >> 1;
				}
				else if (abs(pIn[x - 1 + ym1Width] - pIn[x + 1 + yp1Width])>abs(pIn[x + 1 + ym1Width] - pIn[x - 1 + yp1Width]))
				{
					pOut[pix + R_OFFSET] = (pIn[x + 1 + ym1Width] + pIn[x - 1 + yp1Width]) >> 1;
				}
				else
					pOut[pix + R_OFFSET] = ((pIn[x - 1 + ym1Width] +
					pIn[x + 1 + ym1Width] +
					pIn[x - 1 + yp1Width] +
					pIn[x + 1 + yp1Width]) >> 2);
			}
			else
			if (ODD(x))
			{   // EVEN(y) ODD(x)
				if (abs(pIn[x - 1 + ym1Width] - pIn[x + 1 + yp1Width])<abs(pIn[x + 1 + ym1Width] - pIn[x - 1 + yp1Width]))
				{
					pOut[pix + B_OFFSET] = (pIn[x - 1 + ym1Width] + pIn[x + 1 + yp1Width]) >> 1;
				}
				else if (abs(pIn[x - 1 + ym1Width] - pIn[x + 1 + yp1Width])>abs(pIn[x + 1 + ym1Width] - pIn[x - 1 + yp1Width]))
				{
					pOut[pix + B_OFFSET] = (pIn[x + 1 + ym1Width] + pIn[x - 1 + yp1Width]) >> 1;
				}
				else
					pOut[pix + B_OFFSET] = ((pIn[x - 1 + ym1Width] +
					pIn[x + 1 + ym1Width] +
					pIn[x - 1 + yp1Width] +
					pIn[x + 1 + yp1Width]) >> 2);


				if (abs(pIn[x - 1 + yWidth] - pIn[x + 1 + yWidth])<abs(pIn[x + ym1Width] - pIn[x + yp1Width]))
				{
					pOut[pix + G_OFFSET] = (pIn[x - 1 + yWidth] + pIn[x + 1 + yWidth]) >> 1;
				}
				else if (abs(pIn[x - 1 + yWidth] - pIn[x + 1 + yWidth])>abs(pIn[x + ym1Width] - pIn[x + yp1Width]))
				{
					pOut[pix + G_OFFSET] = (pIn[x + ym1Width] + pIn[x + yp1Width]) >> 1;
				}
				else
					pOut[pix + G_OFFSET] = ((pIn[x - 1 + yWidth] +
					pIn[x + 1 + yWidth] +
					pIn[x + ym1Width] +
					pIn[x + yp1Width]) >> 2);

				pOut[pix + R_OFFSET] = pIn[x + yWidth];
			}
			else
			{   //EVEN(y) EVEN(x)
				pOut[pix + B_OFFSET] = ((pIn[x + ym1Width] +
					pIn[x + yp1Width]) >> 1);
				pOut[pix + G_OFFSET] = ((pIn[x + yWidth] << 2) + (pIn[x - 1 + ym1Width] + pIn[x + 1 + ym1Width] + pIn[x - 1 + yp1Width] + pIn[x + 1 + yp1Width])) >> 3;
				pOut[pix + R_OFFSET] = ((pIn[x - 1 + yWidth] +
					pIn[x + 1 + yWidth]) >> 1);
			}
		}
	}
}

}

[ermig1979]

I will add this issue to future plans.