dicom文件与bmp和jpg文件的相互转化

        前面工作需要，将dicom医学文件转化为普通图像，如bmp，jpg等，中间应用到了CxImage_x64和dcmtk包。实现过程中，遇到了不少麻烦，现将相关过程分享如下，希望能给需要的人提供参考，以及上传所有相关文档。

文档分为以下几部分：

（1）      工程文件，需要将支持库拷贝到对应位置，并设置工程库目录，包含目录之后才能运行；

（2）      可执行文件，拷贝到64位机器即可运行（其中包含医学影像文件image-000001.dcm）；

（3）      支持库：用于打开解析dicom文件，需要按照说明添加；

（4）      Readme.doc：描述文档结构和配置，代码，程序界面等信息。

1、支持库拷贝

先按照支持库里面说明，将两个包拷贝到常用位置，方便后面设置工程，具体可以参照我的设置测试方式。

 

2、工程设置

包含目录设置：

D:\Program Files (x86)\Microsoft VisualStudio 10.0\CxImage_x64\include（如果是32位自己下32位版本）

D:\Program Files\DCMTK\include

 


库目录设置：

D:\Program Files (x86)\Microsoft VisualStudio 10.0\CxImage_x64\lib

D:\Program Files\DCMTK\lib

 

连接器附加项设置：

cximage.lib

netapi32.lib

wsock32.lib

ofstd.lib

oflog.lib

dcmimgle.lib

ijg8.lib

ijg12.lib

ijg16.lib

dcmdata.lib

dcmimage.lib

dcmjpeg.lib

dcmnet.lib

3、工程字符集设为多字节

4、核心代码简单说明

4.1、选择文件

    TCHAR
szFilter[] = _T("文本À文件t(*.dcm)|*.dcm");  //选择dcm文件打开
    CFileDialog fileDlg(TRUE,_T("txt"),NULL, 0, szFilter,this);// 创建对话框 
    CStringstrFilePath;  
    if (IDOK ==fileDlg.DoModal())   // 点击确定保存文件路径
    {  
        strFilePath =fileDlg.GetPathName();   // 文件路径
        SetDlgItemText(IDC_EDIT_open,strFilePath);   // 设置文件路径显示与界面
         
        savepath = fileDlg.GetFolderPath(); // 生存图像保存路径
        filename= fileDlg.GetFileTitle();   // 生存图像保存文件名
 
    }  
 

4.2、打开并转化文件

if ((m_check_bmp.GetCheck()>0)||(m_check_jpg.GetCheck()>0)) // 指定需要转化的图片类型，目前只提供bmP和jpg
    {
        CString str,savestr;  //
        GetDlgItem(IDC_EDIT_open)->GetWindowText(str);  // 获取文件名
        DJEncoderRegistration::registerCodecs();//注册dicom文件封装类
        DcmFileFormat fileformat;
        if(fileformat.loadFile(str).good()) // 文件加载正常
        {      
            THU_STD_NAMESPACE::TDcmFileFormattdf(fileformat); // 转化为封装类
            if (m_check_bmp.GetCheck()>0)
            { // 如果选择了bmp转化
                tdf.saveToBmp(savepath+"\\"+filename+".bmp");   // 转化为bmp文件，具体代码可以点击函数定义进去查看，主要步骤包括将dcom文件的头文件信息以及数据索引指针赋值给bmp对象，在进行保存，jpg类似
            }
            if(m_check_jpg.GetCheck()>0)
            {
                tdf.saveToJpg(savepath+"\\"+filename+".jpg");//  转化为jpg文件
            }
 
                inth = tdf.getHeight();  // 获取影像高度
            int
w = tdf.getWidth(); //获取影像宽度
            CStringstr ;
            str.Format("%d",w);
            SetDlgItemText(IDC_ew,str);  //宽度显示于界面
            str.Format("%d",h);
            SetDlgItemText(IDC_eh,str);  //高度显示
    　　　　MessageBox("转化完毕！","提示");　//
信息提示
        }
        DJEncoderRegistration::cleanup();//注销注册类
    }
    else
        MessageBox("请选择转化类型","提示");

4.3、显示图像

void Cdicomview20170310Dlg::OnBnClickedshowpic()
{
    int
cx, cy; 
    CImage 
image; 
    CRect  
rect; 
    if (PathFileExists(savepath+"\\"+filename+".bmp"))// 判断图像是否存在
    {
        image.Load(savepath+"\\"+filename+".bmp"); // 加载图像
    } else
if ( PathFileExists(savepath+"\\"+filename+".jpg" ))
    {
        image.Load(savepath+"\\"+filename+".jpg");   // 加载图像
    }
    else
        return;
 
    cx  =image.GetWidth();   // 获取图像宽高
    cy  =
image.GetHeight(); 
 
    GetDlgItem(IDC_pic)->GetWindowRect(&rect);  // 获取显示控件区域
    ScreenToClient(&rect); 
    GetDlgItem(IDC_pic)->MoveWindow(rect.left, rect.top,cx,cy,TRUE); 
// 移动到指定区域
    CWnd *pWnd=
NULL; 
    pWnd    =
GetDlgItem(IDC_pic);//获取显示控件句柄
    pWnd->GetClientRect(&rect);//获取句柄对应区域
 
    CDC *pDc = NULL; 
    pDc =
pWnd->GetDC();//显示图像
 
    image.Draw(pDc->m_hDC,rect);//绘制图像
    ReleaseDC(pDc);  // 释放资源
}
 

4.4、bmp转dcm

OFConditionstatus;  // 定义变量
    DcmFileFormat fileformat; 
    DcmDataset*mydatasete=fileformat.getDataset();
    DicomUtils du;
    du.AddDicomElements((DcmDataset*&)mydatasete); 
    Uint16 rows,cols,samplePerPixel,bitsAlloc,bitsStored,highBit,pixelRpr,planConf,pixAspectH,pixAspectV; 
    OFString photoMetrInt; 
    Uint32 length; 
    E_TransferSyntax ts; 
    //循-环¡¤添¬¨ª加¨®4张?图ª?片? 
    //for(inti=0;i<4;++i) 
    //{ 
        charnumtmp[500];  // 申明文件名字符串
        CString tpath = savepath+"\\"+filename+".bmp"  ;
        memset(numtmp,0,sizeof(tpath));
        sprintf(numtmp,"%s",tpath);//转化文件名字符串
        I2DBmpSource *bmpSource = new I2DBmpSource();//I2DBmpSource();
        bmpSource->setImageFile(numtmp);  ／／　读取bmp图像数据
 
        char*pixData=NULL; 
bmpSource->readPixelData(rows,cols,samplePerPixel,photoMetrInt,bitsAlloc,bitsStored,highBit,pixelRpr,planConf,pixAspectH,pixAspectV,pixData,length,ts);  // 解析图像数据
        //deletebmpSource; 
    //}; 
    mydatasete->putAndInsertUint16(DCM_SamplesPerPixel,samplePerPixel);  // 添加bmp图像数据到dcm文件
    mydatasete->putAndInsertString(DCM_NumberOfFrames,"1"); 
    mydatasete->putAndInsertUint16(DCM_Rows,rows); 
    mydatasete->putAndInsertUint16(DCM_Columns,cols); 
    mydatasete->putAndInsertUint16(DCM_BitsAllocated,bitsAlloc); 
    mydatasete->putAndInsertUint16(DCM_BitsStored,bitsStored); 
    mydatasete->putAndInsertUint16(DCM_HighBit,highBit); 
    mydatasete->putAndInsertUint8Array(DCM_PixelData,(Uint8*)pixData,1*length); 
    mydatasete->putAndInsertOFStringArray(DCM_PhotometricInterpretation,photoMetrInt); 
    mydatasete->putAndInsertString(DCM_PlanarConfiguration,"1");
    tpath = savepath+"\\"+filename+"new.dcm"  ;
    memset(numtmp,0,sizeof(tpath));
    sprintf(numtmp,"%s",tpath);
    status=fileformat.saveFile(numtmp,ts);  //保存ｄｃｍ文件
    if(status.bad())    // 保存状态
    { 
        MessageBox("转Áa换?失º¡ì败㨹！ê?","提¬¨¢示º?");
    }  else {
        MessageBox("转Áa换?完ª¨º毕À?！ê?","提¬¨¢示º?");
    }
 

5、程序效果

步骤分3步：

（1）点击打开文件，选择待转化文件；

（2）选择jpg或者bmp待转化格式；

（3）点击转换；

（4）显示转化图像的效果；

（5）退出；

 

转化生成的图片保存在原始影像文件所在目录，如下图所示：