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

时间:2020-12-02 10:04:19

采用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);

相关文章