dicom文件的显示(基于dcmtk实现)

转载 2012年03月26日 09:11:51

采用VC++MFC实现,我把主要的代码给贴上来了。虽然网上也可以搜到部分代码,但是没有具体的实现过程。由于我们采用dcmtk能够将dib读取出来,但是没有调色板,所以我们需要自己构造我们的调试板,具体的代码如下:

DcmFileFormat * pDicomFile=new DcmFileFormat();
 OFCondition Fcond=pDicomFile->loadFile("txt.DCM");//打开文件,此处可以根据具体路径进行更改
 if(Fcond.good())
  pDC->TextOut(50,50,(CString )"load file success");
 else
  pDC->TextOut(50,50,(CString)"load file error");


 DcmDataset * pDataset=pDicomFile->getDataset();

//这个pDataset里面包含了dicom文件都属性的信息,包括病人的id,姓名等以及其他,可以通过
//OFString  ofstr;

//pDataset->findAndGetOFString(DCM_PatientID,ofstr);

//来得到具体的属性,其中DCM_PatientID为dctmk内部定义的宏定义,具体请看dcmtk帮助文档

 E_TransferSyntax xfer=pDataset->getOriginalXfer();

/////////////////////////////////
 DicomImage * pDicomImg=new DicomImage(pDataset,xfer,0,0,0);

 pDicomImg=new DicomImage(FileName,0,0,0);

//上面两个函数都可以实现产生DicomImage 的对象,网上搜到的程序用的基本上都是第一个函数,但是我通过比较,发现第二个函数耗的内存少一些,毕竟少了两个参数,直接把文件名做为参数更直接一些。这在读取大文件的时候直接可以在任务管理器中查看这两个函数的内存消耗显著差异。

 pDicomImg->rotateImage(180);//实现图像翻转的,可以去掉
 void * pDicomDibits;//dib设备无关位图的指针
 unsigned long iDataSize= pDicomImg->createWindowsDIB(pDicomDibits,0,0,8,1,1);

//这个函数比较重要,产生dib的指针,后面的参数为表示产生第一帧的8位的dib位图。

//有了这个dib指针,我们就可以使用VC++的函数直接来显示图像了,但是我们没有调色板,所以,首先还得构造调色板。下面是构造调试板的内容。构造8位的调色板,主要关注的内容是BitmapInfoHeader.biHeight,BitmapInfoHeader.biWidth这两个参数,其它的都是固定的,下面的代码可以直接写到view类都OnDraw()函数中使用就可以显示了。

BITMAPINFOHEADER   BitmapInfoHeader;  
 BitmapInfoHeader.biBitCount=8;
 BitmapInfoHeader.biClrImportant=0;
 BitmapInfoHeader.biClrUsed=0;
 BitmapInfoHeader.biCompression=BI_RGB;
 BitmapInfoHeader.biPlanes=1;
 BitmapInfoHeader.biHeight=pDicomImg->getHeight();//得到图像的高度
 BitmapInfoHeader.biWidth=pDicomImg->getWidth();//得到图像的宽度
 BitmapInfoHeader.biXPelsPerMeter=0;
 BitmapInfoHeader.biYPelsPerMeter=0;
 BitmapInfoHeader.biSize=sizeof(BITMAPINFOHEADER);
 BitmapInfoHeader.biSizeImage=0;
 //////////////////////////////////////
 
 ////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
 RGBQUAD pColorTable[256];
 for(int i=0;i<256;i++)
 {

  pColorTable[i].rgbBlue=i;
  pColorTable[i].rgbGreen=i;
  pColorTable[i].rgbRed=i;
  pColorTable[i].rgbReserved=0;
 }
 //定义一个bmpinfo需要内存大小的内存
 char p[sizeof(BITMAPINFOHEADER)+256*4];
 memcpy(p,&BitmapInfoHeader,sizeof(BITMAPINFOHEADER));
 memcpy(p+sizeof(BITMAPINFOHEADER),pColorTable,1024);
 PBITMAPINFO BitmapInfo = (PBITMAPINFO)p;

 //颜色表
 int colorTableLng;
 colorTableLng=256;

 // 创建调色板
 HPALETTE hPalette=0,hOldPal;
 if (colorTableLng!=0)
 {
  //定义颜色表指针pColorTable,指向DIB的颜色表
  //申请缓冲区,生成LOGPALETTE结构
  LPLOGPALETTE pLogPal = (LPLOGPALETTE)new char[2*sizeof(WORD)
   +colorTableLng * sizeof(PALETTEENTRY)];
  pLogPal->palVersion = 0x300;
  pLogPal->palNumEntries =colorTableLng;
  for(int i = 0; i < colorTableLng; i++)
  {
   pLogPal->palPalEntry[i].peRed= pColorTable[i].rgbRed;
   pLogPal->palPalEntry[i].peGreen =pColorTable[i].rgbGreen;
   pLogPal->palPalEntry[i].peBlue = pColorTable[i].rgbBlue;
   pLogPal->palPalEntry[i].peFlags = 0;
  }
  //
  //创建逻辑调色板
  hPalette =::CreatePalette(pLogPal);
  // 将调色板选入系统
  hOldPal=::SelectPalette(pDC->GetSafeHdc(), hPalette, TRUE);
  //实现调色板
  pDC->RealizePalette();
  //释放缓冲区
  delete []pLogPal;
 }

 //
 //DIB显示所需要的模式
 pDC->SetStretchBltMode(COLORONCOLOR);

//显示DIB到显示器,
 ::StretchDIBits(pDC->GetSafeHdc(), 0, 0,BitmapInfoHeader.biWidth ,BitmapInfoHeader.biHeight,
  0,0, BitmapInfoHeader.biWidth,BitmapInfoHeader.biHeight,pDicomDibits,BitmapInfo, DIB_RGB_COLORS, SRCCOPY);

相关文章推荐

关于DICOM与Qt的关系

最近Qt群有人问到如何使用C++编写代码实现DICOM文件与其他格式的转换工作,于是乎就查阅了下Qt下面的事情,google到了QtDcm这个基于C++的Qt的开源库,于是下载了其源代码,发现源代码还...
  • Esonpo
  • Esonpo
  • 2013-10-08 17:03
  • 2186

用VC++和DCMTK显示DICOM医学图像

用VC++和DCMTK显示DICOM医学图像 摘 要 DICOM是医学影像存储和传输的国际标准,它的制定为不同的医学影像设备和用户提供了统一的接口标准和交互协议。解读 DICOM 的医学图像文件格式以...

DCMTK开源库的学习笔记3:dcmtk文件中数据元的修改

dcm文件是医学领域DICOM3.0标准所对应的主要的文件格式。前两篇学习笔记中,学习了读取dcm文件的相关信息,如信息头MetaInformation元素、像素数据元素,只停留在了读取和显示dcm文...

GDI+ 显示DICOM 图像测试

GDI+ 显示 DICOM 图像测试 主要使用:VS2008  MFC 、 GDI+ 、  DCMTK。      一、MFC下使用GDI+需要稍微配置一下,过程很简单,网上也很多,不详细介绍。 ...

DCMTK:DCM文件转QImage

代码如下,暂做记录QImage *image = 0; ///register RLE decompression codecs DcmRLEDecoderRegistration::register...

【学习DCMTK】学习资料

搭建开发环境
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:深度学习:神经网络中的前向传播和反向传播算法推导
举报原因:
原因补充:

(最多只允许输入30个字)