bayer类型转RGB

通过把bayer类型的图片转换成bgr类型的图片（在matlab里图像读取的数据方式是rgb）,bayer类型是单通道的数据图像，像素排列格式为奇数行为grgrgr...偶数行为bgbgbg......

首先建立一个三通道的图像cvCreateImage(cvSize(img->width,img->height),IPL_DEPTH_8U,3)把bayer类型数据grgrgr.../bgbgbg......放到三通道BGR相应的位置处。也就是说把bayer类型的第一个像素点G（此时一个像素是一个字节）放到三通道第一个像素的G位置（此时一个像素是三个字节）。然后通双线性插值公式进行插值。如果要看分离后的像素值的大小，分离后保存的图片最好为bmp格式，jpg格式的图片在matlab里所显示的像素值可能会和opencv里所保存的像素值不一样（我也不知道为什么）。还有 我的双线性插值所得的像素值没有取小数，取得都是整数，显示出来的效果很差 不知道是不是这个原因（欢迎大神指点），不知道opencv像素值怎么可以保存小数到图片里，因为创建的是PL_DEPTH_8U 无符号八位整形，其他格式看了 没有找到无符号八位浮点型（希望搞过的可以交流下）。 

#include<highgui.h>
#include<cv.h>
#include<iostream>
#include<stdio.h>
using namespace std;
using namespace cv;


int main()
{
int k=1;
IplImage *img= cvLoadImage("bayer_test.tif",-1);




cvNamedWindow("org");
cvShowImage("org",img);
//IplImage *imgdst=cvLoadImage("test.jpg");
IplImage *imgdst = cvCreateImage(cvSize(img->width,img->height),IPL_DEPTH_8U,3);
把单通道像素分配到三通道里///
for(int i =0;i <720; i=i+2)
{

for(int j = 0;j <1280; j=j+2)  
{

((uchar*)(imgdst->imageData+i*imgdst->widthStep))[j*imgdst->nChannels+1]=((uchar*)(img->imageData+i*img->widthStep))[j];
}

for(int j=1;j<1280;j=j+2)     
{
((uchar*)(imgdst->imageData+i*imgdst->widthStep))[j*imgdst->nChannels+2]=((uchar*)(img->imageData+i*img->widthStep))[j];
}


}


for(int i = 1;i < 720; i=i+2)     
{
for(int j = 0; j< 1280;j =j+2 )
{
((uchar*)(imgdst->imageData+i*imgdst->widthStep))[j*imgdst->nChannels+0]=((uchar*)(img->imageData+i*img->widthStep))[j];
}
for(int j =1;j < 1280;j =j+2)   
{
((uchar*)(imgdst->imageData+i*imgdst->widthStep))[j*imgdst->nChannels+1]=((uchar*)(img->imageData+i*img->widthStep))[j];
}
}


/*FILE *pFile;
pFile= fopen("aa.txt","a+");
for(int i=0;i<720;i++)
{
for(int j=0;j<1280;j++)
{
((uchar*)(imgdst->imageData+i*imgdst->widthStep))[j*imgdst->nChannels+0]=120;
((uchar*)(imgdst->imageData+i*imgdst->widthStep))[j*imgdst->nChannels+1]=121;
((uchar*)(imgdst->imageData+i*imgdst->widthStep))[j*imgdst->nChannels+2]=122;
fprintf(pFile,"%d%d%d\n",((uchar*)(imgdst->imageData+i*imgdst->widthStep))[j*imgdst->nChannels+0],((uchar*)(imgdst->imageData+i*imgdst->widthStep))[j*imgdst->nChannels+1],((uchar*)(imgdst->imageData+i*imgdst->widthStep))[j*imgdst->nChannels+2]);
}
}fclose(pFile);*/




cvNamedWindow("org1");
cvShowImage("org1",img);
cvSaveImage("org1.tif",img);


cvNamedWindow("dst0");
cvShowImage("dst0",imgdst);
cvSaveImage("dst0.bmp",imgdst);

/进行GR排列进行插值
for(int i =2;i<720; i=i+2)
{
for(int j=2;j<1280;j++)
{    ///在G处插值R
((uchar*)(imgdst->imageData+i*imgdst->widthStep))[j*imgdst->nChannels+2]=((img->imageData+i*img->widthStep)[j-1]+(img->imageData+i*img->widthStep)[j+1])/2;
   //在G处插值B
((uchar*)(imgdst->imageData+i*imgdst->widthStep))[j*imgdst->nChannels+0]=((img->imageData+(i-1)*img->widthStep)[j]+(img->imageData+(i+1)*img->widthStep)[j])/2;

}
}




for(int i=2;i <720;i=i+2)
{
for(int j=1; j<1280;j++)
{
//在R处插值G
((uchar*)(imgdst->imageData+i*imgdst->widthStep))[j*imgdst->nChannels+1]=((img->imageData+i*img->widthStep)[j-1]+(img->imageData+i*img->widthStep)[j+1]+(img->imageData+(i+1)*img->widthStep)[j]+(img->imageData+(i-1)*img->widthStep)[j+1])/4;
   //在R处插值B
            ((uchar*)(imgdst->imageData+i*imgdst->widthStep))[j*imgdst->nChannels+0]=((img->imageData+(i-1)*img->widthStep)[j-1]+(img->imageData+(i+1)*img->widthStep)[j+1]+(img->imageData+(i+1)*img->widthStep)[j-1]+(img->imageData+(i-1)*img->widthStep)[j+1])/4;
}
}

//对BG排列进行插值///
for(int i =1;i<720; i=i+2)
{
for(int j=2;j<1280;j++)
{
///在B处插值R
((uchar*)(imgdst->imageData+i*imgdst->widthStep))[j*imgdst->nChannels+0]=((img->imageData+(i-1)*img->widthStep)[j-1]+(img->imageData+(i+1)*img->widthStep)[j+1]+(img->imageData+(i+1)*img->widthStep)[j-1]+(img->imageData+(i-1)*img->widthStep)[j+1])/4;
    ///在B处插值G
((uchar*)(imgdst->imageData+i*imgdst->widthStep))[j*imgdst->nChannels+1]=((img->imageData+i*img->widthStep)[j-1]+(img->imageData+i*img->widthStep)[j+1]+(img->imageData+(i-1)*img->widthStep)[j]+(img->imageData+i+1*img->widthStep)[j])/4;
}


}



  for(int i =1;i<720; i=i+2)
{
for(int j=1;j<1280;j++)
{
///在G处插值R
((uchar*)(imgdst->imageData+i*imgdst->widthStep))[j*imgdst->nChannels+2]=((img->imageData+i*img->widthStep)[j+1]+(img->imageData+i*img->widthStep)[j-1])/2;
  ///在G处插值B
      ((uchar*)(imgdst->imageData+i*imgdst->widthStep))[j*imgdst->nChannels+0]=((img->imageData+(i+1)*img->widthStep)[j]+(img->imageData+(i-1)*img->widthStep)[j])/2;
}


}
 
cvNamedWindow("dst");
cvShowImage("dst",imgdst);
cvSaveImage("dst.bmp",imgdst);


cvWaitKey();

}