{学海无涯}【基于VTK的三维建模】——1.Dicom文件的读取与显示

时间:2022-08-06 10:05:57

学习中,记笔记是个好习惯。

似乎有点晚了,不过没关系,从现在开始!


这次要做的项目是老师给的,要从一批Dicom文件中读取头部模型的数据并进行三维建模。先是读取头部CT图片,再是读取头部MIR(核磁共振成像)图片。不过医学图像都是存储为Dicom格式。


今天完成的效果图如下:

{学海无涯}【基于VTK的三维建模】——1.Dicom文件的读取与显示


Dicom 文件不仅包含图像,还包含医院信息,患者信息等许多信息,格式等内容不再赘述,我用了一个写好的类CDicom,在文档中创建指针CDicom* m_pDicomCreate。利用该指针,完成相应的操作。


一.读取

代码如下

void CControlView::GetPictures(CString FolderDir){
CMainFrame* pMainFrame=(CMainFrame*)::AfxGetMainWnd();
if (pMainFrame ==NULL)
{
return;
}
CImageView* pView= static_cast<CImageView*>(pMainFrame->m_wndSplitter.GetPane(0,0));
CDicom002Doc* pDoc=(CDicom002Doc*)pMainFrame->GetActiveDocument();

if (pView == NULL)
return;
if (pDoc == NULL)
return;

CFileFind FFind;
CString szDir=FolderDir;
if(szDir.Right(1)!='\\')
szDir+='\\';
szDir+="001.dcm";
BOOL Ressult=FFind.FindFile(szDir);

if(Ressult)
{
Ressult=FFind.FindNextFile();

pDoc->m_pDicomCreate=new CDicom;
ASSERT(pDoc->m_pDicomCreate!= NULL);
pDoc->m_pDicomCreate->read_dicom(szDir);//以下都是通过指针得到Dicom文件中的各种信息

pDoc->m_dRowPixelSpacing = pDoc->m_pDicomCreate->pixel_spacing.x;//x space
pDoc->m_dColPixelSpacing = pDoc->m_pDicomCreate->pixel_spacing.y;//y space
pDoc->m_dThickness = pDoc->m_pDicomCreate->slice_thickness;

pDoc->m_nRow = pDoc->m_pDicomCreate->rows;//
pDoc->m_nColumn = pDoc->m_pDicomCreate->columns;//
pDoc->m_nLayer = (int)( pDoc->m_nFileNumberCreate* pDoc->m_dThickness/ pDoc->m_dRowPixelSpacing);
pDoc->m_nPixels = pDoc->m_nRow* pDoc->m_nColumn;//一层图像的像素点数

pDoc->m_nWW = pDoc->m_pDicomCreate->window_width;//从m_pDicomCreate中得到窗宽窗位
pDoc->m_nWC = pDoc->m_pDicomCreate->window_center;

pDoc->m_pDicomBuffer=pDoc->m_pDicomCreate->ct_value;//你妹的,找半天,这是ct图像缓冲存放的地方,非常关键

pDoc->m_strModality = pDoc->m_pDicomCreate->modality;

if(pDoc->m_nRow > pDoc->m_nColumn)
pDoc->m_nNumN = pDoc->m_nRow;//m_nNumNCreate 导航视图缓冲的长和宽
else
pDoc->m_nNumN = pDoc->m_nColumn;


UpdateData(TRUE);
m_id=pDoc->m_pDicomCreate->patient_id;
m_name=pDoc->m_pDicomCreate->patient_name;
m_date=pDoc->m_pDicomCreate->patient_birth_date;
m_sex=pDoc->m_pDicomCreate->patient_sex;
UpdateData(FALSE);
pDoc->m_flagPic=true;
}

}

学到了很多东西,你必须进入程序仔细研读。比如说CDicom这个类中,缓冲图像是存在ct_value中的。


二.数据转化

我们得到的数据是存放在缓冲区的,不能够直接显示。必须转化为位图格式,才能够被操作识别。代码如下:

void CDicom002Doc::MakeDicomBMP(){//if(m_pBMPN!=NULL)//{//delete [] m_pBMPN;//m_pBMPN=NULL;//}m_pBMPN=new unsigned char[m_nRow*m_nColumn];memset(m_pBMPN, 0, m_nRow*m_nColumn*sizeof(unsigned char));TransferToBMP();//if(m_pBMPNs!=NULL)//存储经简单处理过的bmp图  如旋转一定角度  纯数据//{//delete [] m_pBMPNs;//m_pBMPNs=NULL;//}double b=double(m_nColumn)/m_nRow;//int bytoa;//int ay;if (b == 1)//htow==1{//ay=m_nRow;m_pBMPNs=new unsigned char[m_nColumn* m_nRow];memset(m_pBMPNs, 0, m_nColumn* m_nRow* sizeof(unsigned char));//预置为黑色for(int x=0; x< m_nRow; x++)for(int y=0; y< m_nColumn; y++)m_pBMPNs[y* m_nRow+ x]=m_pBMPN[y* m_nRow+ x];//m_pBMPTs=m_pBMPT;}}void CDicom002Doc::TransferToBMP(){int i,j;double ratio,wMin,wMax;ratio=255.0/m_nWW;wMin=m_nWC-(m_nWW>>1);wMax=m_nWC+(m_nWW>>1);for ( i= 0; i< m_nRow; i++)   for ( j= 0; j< m_nColumn; j++)   {      if (m_pDicomBuffer[j* m_nRow+ i] >= (wMax))   m_pBMPN[(m_nColumn-1-j)* m_nRow+ i] = 255;      else if (m_pDicomBuffer[j* m_nRow + i] < (wMin))   m_pBMPN[(m_nColumn-1-j)* m_nRow+ i] = 0;      else   m_pBMPN[(m_nColumn-1-j)* m_nRow+ i] =unsigned char( ratio*(m_pDicomBuffer[j* m_nRow + i]-(wMin)) );   }//m_nRow方向数据颠倒}

其实m_pBMPN[(m_nColumn-1-j)* m_nRow+ i] =unsigned char( ratio*(m_pDicomBuffer[j* m_nRow + i]-(wMin)) );

这个原理,我还是不太懂,如果有哪位同学碰巧知道,请不吝赐教。


三.显示

显示同样利用了现成的类CSTDIB,用来对位图进行操作的类。


void CImageView::DrawDicom(CDC *pDC){CRect rect;CDicom002Doc* pDoc = static_cast<CDicom002Doc* >(GetDocument());CBrush m_Brush(RGB(0,0,0));pDoc->MakeDicomBMP();//if(pDoc->m_pdib!=NULL)//{//delete pDoc->m_pdib;//pDoc->m_pdib=NULL;//}pDoc->m_pdib=new CSTDIB;rect=CRect(10,10,360,360);pDoc->m_pdib->CreateDib256(pDoc->m_pBMPNs,pDoc->m_nNumN,pDoc->m_nNumN);////根据外部图像数据创建256色(灰度)DIB (宽度可以不是4的倍数)pDoc->m_pdib->DrawEx(pDC,rect.Width(),rect.Height(),rect.left,rect.top);}


恩,大体是这样,后面还有很多工作要做。

大概列个提纲:

1.读取一个文件夹中的所有图片数据,存放在数组里

2.利用VTK对数据进行建模。形成3D图像