将DCM文件保存成BMP文件或数据流

徘徊边缘 提交于 2020-03-17 06:19:18

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(&mDibImage[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浏览软件中的结果一致。
 

易学教程内所有资源均来自网络或用户发布的内容,如有违反法律规定的内容欢迎反馈
该文章没有解决你所遇到的问题?点击提问,说说你的问题,让更多的人一起探讨吧!