DCMTK开源库的学习笔记1：将DCM文件保存成BMP文件或数据流（即数组）

DCMTK开源库介绍：

DCMTK是目前最全面实现DICOM3.0标准的开源库，通过结合DCMTK开源库和CxImage图像开源库，能够很方便的开发属于自己的DCM文件编辑浏览软件。

DCMTK与CxImage的结合：

下面记录”DCM文件到BMP文件格式的转换“的学习和探索的过程。

1）BMP文件的格式这里不细讲（可自行参阅MSDN相关资料）。BMP属于设备无关图（DIB），内存中只要依次包含BMP文件头（BITMAPFILEHEADER）、BMP信息头（BITMAPINFOHEADER）、调色板（Palette），以及真实像素信息，就可以在所有显示设备中进行显示。因此将DCM文件转换到BMP格式的主要工作是：从DCM众多数据元中挑选出BMP文件头、信息头所必须的信息（如图像宽度、图像高度、每个像素所占的空间、像素数据等）。

2）此次采用了dcmtk中的dcmimgle开发包，利用DicomImage类提供的createWindowsDIB接口以及writeBMP来进行格式转化。具体代码如下：

DcmFileFormat *mImage=new DcmFileFormat();
DJDecoderRegistration::registerCodecs(); // register JPEG codecs
DcmDataset *dataset = mImage->getDataset();
E_TransferSyntax xfer=dataset->getOriginalXfer();
dataset->chooseRepresentation(EXS_LittleEndianExplicit,NULL);
DicomImage *pdcmImage=new DicomImage((DcmFileFormat*)mImage,EXS_LittleEndianExplicit,CIF_AcrNemaCompatibility,0,1);
if(pdcmImage->getStatus()!=EIS_Normal)
exit(0);
void* pDIB=NULL;
int size=0;
pdcmImage->createWindowsDIB(pDIB,size,0,8,1,1);

至此，pDIB内存中就获得了”c:\test.dcm“文件像素数据多对应的BMP格式的数据。如果想将pDIB保存成BMP文件，还需要补充BMP文件头（BITMAPFILEHEADER）、BMP信息头（BITMAPINFOHEADER）、调色板（Palette）。参照MSDN中BITMAPFILEHEADER和BITMAPINFOHEADER的结构，逐个成员赋值即可。假设我们生成的文件头变量为mBmpFileHeader，信息头变量为mBmpInfoHeader。由于createWindowsDIB函数中的BMP像素位数设定的是8，所以BMP图像属于索引图像，其文件中需要调色板，调色板的生成如下所示（最简单的生成方法）：

RGBQUAD* mPalette=new RGBQUAD[256];
memset(mPalette,0,sizeof(RGBQUAD)*256);
for(int i=0;i<256;++i)
{
mPalette[i].rgbBlue=i;
mPalette[i].rgbGreen=i;
mPalette[i].rgbRed=i;
}

随后将mBmpFileHeader、mBmpInfoHeader、mPalette添加到pDIB的前面，就构成了一幅BMP图像，直接将此数据保存成bmp文件即可。具体的链接步骤如下：

<span style="white-space:pre">	</span>int imgsize=width*height;
int headsize=sizeof(BITMAPFILEHEADER)+sizeof(BITMAPINFOHEADER)+sizeof(RGBQUAD)*256;
int size=headsize+imgsize;
BYTE* mDibImage=new BYTE[size];
memset(mDibImage,0,size);
memcpy(mDibImage,mBmpFileHeader,sizeof(BITMAPFILEHEADER));
memcpy(&mDibIma
3ff0
ge[sizeof(BITMAPFILEHEADER)],mBmpInfoHeader,sizeof(BITMAPINFOHEADER));
memcpy(&mDibImage[sizeof(BITMAPFILEHEADER)+sizeof(BITMAPINFOHEADER)],mPalette,sizeof(RGBQUAD)*256);
memcpy(&mDibImage[headsize],pDIB,imgsize);

至此mDibImage内存存储的数据就是一幅完整的BMP文件。将其保存后的结果如下图（左）所示：





相比dcm浏览软件中看到的结果（上图右），自己生成的图像模糊很多。思考原因，可能是窗宽窗位的关系。DicomImage中也提供了调整窗宽窗位的接口setWindow。在调用createWindowsDIB函数之前，调用setW函数，随后生成的BMP图像结果如下：



与常用的dcm浏览软件中的结果一致。

源代码下载：

鉴于有些网友向我询问完整源码的问题，以前都是直接邮箱发送的，最近嫌麻烦，就将完整的.h和.cpp文件上传到了CSDN中，资源链接是：http://download.csdn.net/detail/zssureqh/7860671，需要2个积分奥。哈哈，如果不是CSDN的朋友，或者没有积分的朋友也不要着急，也可以到我的Github主页进行免费下载，地址是：https://github.com/zssure-thu/DcmtkAndCxImage，有问题及时反馈，多多交流，共同进步。

作者：zssure@163.com

时间：2014-09-03（最后更新时间）