使用OpenCV打开RAW文件

首先来介绍一下关于相机的基础知识：
1.CCD/CMOS相机的感光元件对波长（即颜色）不敏感，如果拿一个裸体的CCD/CMOS传感器去拍摄图像，只能得到灰度图；
2.因为上述这条，人们必须找到能够将波长区分开的方法，其中一种是使用三个滤光片（通常是RGB三色），在这三个滤光片之后放置三个CCD，这就是3CCD相机；

3.显然上一种方法的成本太高了，聪明的人类又想到了另一个方法：那就是在CCD/CMOS传感器矩阵之前放置一个滤色片矩阵，每个像素对应一个滤色片，将RGB三种颜色的滤色片均匀分布在这个矩阵中，拍摄到图像后将对应颜色的像素的值取出来并进行插值，就得到了三个通道的数据。

RAW文件存储的就是第三种方法拍摄的原始数据，可以使用photoshop打开，放大后能够很明显的观察到类似棋盘格一样的像素值。虽然没有白平衡设置，但最大的好处是真实的数据也没有被改变，基于这些数据操作者可以设置自己的通道权值，能够任意的调整色温和白平衡。

 

将图像放大：

局部放大

 

继续放大：

在一个项目中我用到了12bit的工业相机，通过调用该相机SDK中的函数可以将拍摄的图像数据保存为16bit单通道的raw文件，其中前12bit有效，末尾4bit是0。

下面来介绍如何用OpenCV打开raw文件。思路如下：
1.以二进制方式打开文件；
2.将每个像素对应的16bit的数据分成两个8bit，分别放入两个矩阵中；
3.对这两个矩阵进行色彩空间变换，转换为两个8bit三通道的图像；
4.将这两个图像合成为一个16bit三通道的图像。

废话少说，代码如下：

void CMy20120510readrawfileDlg::OnBnClickedButton2()
{
const int WIDTH = 1360;
const int HEIGHT = 1024;

CFile file;
file.Open(_T("aaa.raw"), CFile::modeRead | CFile::typeBinary);
file.SeekToBegin();

BYTE * pfilebuf = new BYTE[HEIGHT*WIDTH*2];

if (HEIGHT*WIDTH*2 != file.Read(pfilebuf, HEIGHT*WIDTH*2))
{
//提示文件读取错误
file.Close();
return;
}

file.Close();

//////////////////////////////////////////////////////////////////////////

CvMat* mat_a = cvCreateMat(1, HEIGHT*WIDTH, CV_8U); //单行矩阵便于赋值操作
CvMat* mat_b = cvCreateMat(1, HEIGHT*WIDTH, CV_8U); //同上

int i=0;

do {
CV_MAT_ELEM(*mat_a, unsigned char, 0, i) = pfilebuf[i*2]; //低8位信息
CV_MAT_ELEM(*mat_b, unsigned char, 0, i) = pfilebuf[i*2+1]; //高8位信息

i++;
}while(i<HEIGHT*WIDTH);

delete[] pfilebuf;

cvReshape(mat_a, mat_a, 0, HEIGHT); //把单行矩阵整形为二维矩阵
cvReshape(mat_b, mat_b, 0, HEIGHT);

IplImage* img_a = cvCreateImage(cvSize(WIDTH,HEIGHT), IPL_DEPTH_8U, 3);
IplImage* img_b = cvCreateImage(cvSize(WIDTH,HEIGHT), IPL_DEPTH_8U, 3);

cvCvtColor(mat_a, img_a, CV_BayerBG2RGB); //色彩空间转换，即Bayer模式转为RGB
cvCvtColor(mat_b, img_b, CV_BayerBG2RGB);

cvReleaseMat(&mat_a);
cvReleaseMat(&mat_b);

cvNamedWindow("img_a");
cvNamedWindow("img_b");
cvShowImage("img_a", img_a);
cvShowImage("img_b", img_b);

//////////////////////////////////////////////////////////////////////////

//因为cvAddWeighted需要参数矩阵都具有相同类型、相同大小
IplImage* img_a_16 = cvCreateImage(cvSize(WIDTH,HEIGHT), IPL_DEPTH_16U, 3);
IplImage* img_b_16 = cvCreateImage(cvSize(WIDTH,HEIGHT), IPL_DEPTH_16U, 3);

cvConvert(img_a, img_a_16);
cvConvert(img_b, img_b_16);

IplImage* img = cvCreateImage(cvSize(WIDTH,HEIGHT), IPL_DEPTH_16U, 3);

//高8位左移8位加上低8位合成一个16位图像
cvAddWeighted(img_a_16, 1, img_b_16, 256, 0, img);

cvNamedWindow("img");
cvShowImage("img", img);

cvWaitKey();
cvDestroyAllWindows();

cvReleaseImage(&img_a);
cvReleaseImage(&img_b);
cvReleaseImage(&img);
}

得到的图像如下：

可以看出img_a的结果是原16bit图像的低8位，直接显示的话是没什么意义的。
最后合成的16bit肉眼几乎观察不出有什么区别，但是对OpenCV来说已经是可以直接处理的16bit图像数据了！

接下来，要怎么处理就可以自由发挥啦。

顺便给出OpenCV帮助文档里cvCvtColor对Bayer模式转换的说明，请认真阅读哦：

 

 

http://www.opencv.org.cn/forum/viewtopic.php?f=1&t=18681&p=60197&hilit=%E5%8D%95%E5%BA%94%E7%9F%A9%E9%98%B5#p60197