程序博客网 > 海岛奇兵保险库及数据

基于VTK的MFC应用程序开发(2)

来源：互联网 发布：海岛奇兵保险库及数据 编辑：程序博客网 时间：2024/06/04 22:00

基于VTK的MFC应用程序开发(2)

分类： VTK应用示例 2013-03-29 13:03 6647人阅读 评论(18)收藏 举报
MFCVTK图像重采样

目录(?)[+]

现在基于VTK的MFC程序框架已经搭建起来。这一节我们来对上节的程序进行扩展，实现图像的重采样。重采样是图像处理中的一个常见功能，尤其是在医学图像处理中，一些三维图像数据量非常大，对于内存的要求比较高，在处理时为了时间和空间效率的平衡，对图像进行重采样（降采样）；另外还有的应用如配准中，为了保持图像的维数一致，还可能对图像进行升采样。但是在目前常使用到的医学图像处理软件中，并没有提供重采样的功能。这个也是本程序开发的初衷，在讲述重采样的原理的同时，提供一个方便的图像重采样软件。

关于图像重采样的内容在VTK在图像处理中的应用(5)一文中也已经介绍过。图像重采样实际上是一个插值的过程。原图像的像素可以看做一些规则分布的离散点，而重采样则意味在这些离散点的基础上通过插值的方式来计算新的离散点的灰度值。最简单的差值方式为最近邻法，即取原图像中距离当前采样点最近的像素值作为当前采样点的输出，这种方式计算简单，但是误差比较大，图像失真比较严重；第二种也是一个最常用的方式是线性差值方式，二维图像对应双线性插值（利用当前新采样点周围的四个原图像像素点进行线性插值），三维图像则是利用周围的8个像素点进行三次线性插值。采用线性插值方式计算相对也比较简单，但是同样会对图像具有平滑效果。另外还有三次多项式插值，B样条插值等方式，这种方式计算量比较大，但是精度比较高。这些方式可以根据不同的需要进行选择。

本篇内容将以线性插值方式为主，实现一个图像重采样应用程序，利用该程序可以方便的进行图像的重采样；同时，在本篇中还会介绍怎样利用CMake在MFC框架添加一个新类。关于如果搭建基于VTK的单文档程序，请参阅基于VTK的MFC应用程序开发(1)。

 

1.  添加一个对话框资源

 

 

添加相应的控件，设置对话框及控件ID。界面设计如上。用户首先需要选择改变图像的维数或者图像的像素间隔。这两个参数是相互关联的，即在每个方向上修改其中一个时，另一个也会随之改变，而两者的乘积保持不变。因此这里用户可以选择修改图像的维数还是像素间隔来进行重采样。采样参数共有6个，分别为X方向维数，Y方向维数，Z方向维数，X方向像素间隔，Y方向像素间隔和Z方向像素间隔。像素间隔在医学图像中是一个很重要的参数，表示像素的物理尺寸，普通图像较少使用，默认值为1。维数或者图像像素间隔，每次只能一项。当图像的某个方向的维数为1时，则不能进行重采样。

 

2.  将对话框添加为类

 

设置类的名字CResampleDialog，并选择基类为CDialog。这里需要注意的是，由于MFC工程通常代码文件与工程文件是在一个路径下，这里添加对话框类文件后，这些文件的路径将会是在工程路径下即bin目录中。而我们采用CMake来管理工程，代码文件与工程文件是分离的，因此这里不妨将.h文件和.cpp文件的路径设置为代码路径下，如下图所示红色区域显示：

 

 

点击完成后，可以发现在代码路径中多出了两个文件，即新添加的类文件。与此同时，在工程的解决方案中，同样会多出了两个文件，如下图中所示。

 

 

这样为了使用CMake进行管理，将两个文件从工程的解决方案中删除，并添加至CMakeLists.txt文件中，此时该工程的CMakeLists.txt文件内容如下。使用CMake对MFC工程进行管理时，添加类的过程相对比较麻烦。不知道是否还有更方便的方法，如果有，也欢迎各位赐教。

[plain] view plaincopy
  1. #----------------------------------------------------------------------------------  
  2.   
  3. cmake_minimum_required( VERSION 2.8 )  
  4.   
  5. project( vtkSDI )  
  6.   
  7. #----------------------------------------------------------------------------------  
  8.   
  9. # 查找并包含VTK工具包  
  10.   
  11. find_package( VTK )  
  12.   
  13. if (VTK_FOUND)  
  14.   
  15.   include (${VTK_USE_FILE})  
  16.   
  17. else (VTK_FOUND)  
  18.   
  19.   message (FATAL_ERROR "Cannot build without VTK. Please set VTK_DIR")  
  20.   
  21. endif( VTK_FOUND )  
  22.   
  23. #----------------------------------------------------------------------------------  
  24.   
  25. # 这里添加本工程的文件  
  26.   
  27. # 主要分为两部分:  
  28.   
  29. # 一是新建的单文档程序中的非工程文件  
  30.   
  31. # 二是用户后续添加的类文件  
  32.   
  33. SET( PROJECT_SRCS  
  34.   
  35.     MainFrm.h  
  36.   
  37.     MainFrm.cpp  
  38.   
  39.    
  40.     stdafx.h  
  41.   
  42.     stdafx.cpp  
  43.   
  44.    
  45.     vtkSDI.h  
  46.   
  47.     vtkSDI.cpp  
  48.   
  49.    
  50.     vtkSDIDoc.h  
  51.   
  52.     vtkSDIDoc.cpp  
  53.   
  54.    
  55.     vtkSDIView.h  
  56.   
  57.     vtkSDIView.cpp  
  58.   
  59.    
  60.     ResampleDialog.h  
  61.   
  62.     ResampleDialog.cpp  
  63.   
  64.    
  65.     targetver.h  
  66.   
  67.     Resource.h  
  68.   
  69.     vtkSDI.rc  
  70.   
  71.       
  72.     res/vtkSDI.rc2  
  73.   
  74.     res/vtkSDI.ico  
  75.   
  76.     res/vtkSDIDoc.ico  
  77.   
  78.     res/Toolbar.bmp  
  79.   
  80.     res/Toolbar256.bmp  
  81.   
  82.     )  
  83. #----------------------------------------------------------------------------------  
  84.   
  85. # 设置工程包含的vtk模块，这里根据需要加载对应的模块  
  86.   
  87. include("${VTK_DIR}/GUISupport/MFC/VTKMFCSettings.cmake")  
  88.   
  89. set( VTK_LIBS ${vtk_libraries}   
  90.   
  91.               vtkMFC  
  92.   
  93.               vtkIO  
  94.   
  95.               vtkRendering   
  96.   
  97.               vtkGraphics   
  98.   
  99.               vtkHybrid   
  100.   
  101.               vtkFiltering   
  102.   
  103.               vtkCommon   
  104.   
  105.               vtkImaging   
  106.   
  107.    )  
  108.   
  109. VTK_MFC_ADD_DELAYLOAD_FLAGS(CMAKE_EXE_LINKER_FLAGS  
  110.   
  111.     vtkMFC.dll  
  112.   
  113.     vtkIO.dll  
  114.   
  115.     vtkRendering.dll  
  116.   
  117.     vtkGraphics.dll  
  118.   
  119.     vtkHybrid.dll  
  120.   
  121.     vtkFiltering.dll  
  122.   
  123.     vtkCommon.dll  
  124.   
  125.     vtkImaging.dll  
  126.   
  127.     )  
  128. #----------------------------------------------------------------------------------  
  129.   
  130. ADD_EXECUTABLE( vtkSDI WIN32  ${PROJECT_SRCS} )  
  131.   
  132. TARGET_LINK_LIBRARIES ( vtkSDI ${VTK_LIBS} )  
  133. #----------------------------------------------------------------------------------  


 

3.  实现图像重采样对话框类

 

主要实现构造函数，初始化函数，Radio消息响应，文本框消息响应，确定和取消按钮消息响应。由于维数与间隔两个参数是互斥的，在修改一个同时另外一个要随之改变。例如一个图像X方向的维数为200，像素间隔为1.0，当用户改变其维数为100时，像素间隔就要相应的变为2.0。

 

4.  实现图像动态加载

 

为了尽可能多的支持图像类型，而不是单一的支持一种图像，对图像加载函数进行修改。代码中出现了一个新的类vtkImageReader2Factory，该类能够根据设置的图像路径名字来返回一个可用的图像读取类对象，找不到则返回NULL。这样的话便能够根据不同的图像类型自适应读取图像。

 

[cpp] view plaincopy
  1. void CvtkSDIDoc::OnFileOpen()  
  2. {  
  3.     // TODO: Add your command handler code here  
  4.     CString FilePathName;  
  5.     m_pImageData = NULL;  
  6.   
  7.     CFileDialog dlg(TRUE, NULL, NULL,  
  8.                 OFN_HIDEREADONLY | OFN_OVERWRITEPROMPT,  
  9.                 (LPCTSTR)_TEXT("All Files (*.*)|*.*||"),  
  10.                     NULL);  
  11.   
  12.     if(dlg.DoModal()==IDOK)  
  13.     {  
  14.         FilePathName=dlg.GetPathName();  
  15.           
  16.         // Read file  
  17.         vtkSmartPointer<vtkImageReader2Factory> readerFactory =  
  18.             vtkSmartPointer<vtkImageReader2Factory>::New();  
  19.         vtkImageReader2 * reader =   
  20.             readerFactory->CreateImageReader2(FilePathName.GetBuffer(0));  
  21.         if (reader)  
  22.         {  
  23.             reader->SetFileName(FilePathName.GetBuffer(0));  
  24.             reader->Update();  
  25.   
  26.             m_pImageData = reader->GetOutput();  
  27.             reader->Delete();  
  28.   
  29.             UpdateAllViews(0);  
  30.         }  
  31.         else  
  32.         {  
  33.             AfxMessageBox("Image Format is not supported now!");  
  34.         }  
  35.     }  
  36.     else  
  37.     {  
  38.         return;  
  39.     }  
  40. }  

 

5.  添加图像重采样菜单及响应

 

 

将菜单的响应函数响应到CvtkSDIDoc类中。当用户调用该菜单时，弹出重采样设置对话框，由用户设置参数。根据用户设置的参数对图像进行重采样，代码如下。

 

[cpp] view plaincopy
  1. void CvtkSDIDoc::OnSizeResample()  
  2.   
  3. {  
  4.   
  5.     // TODO: Add your command handler code here  
  6.   
  7.    
  8.   
  9.     //get dimension and space of original image  
  10.   
  11.     int dims[3];  
  12.   
  13.     m_pImageData->GetDimensions(dims);  
  14.   
  15.    
  16.   
  17.     double spaces[3];  
  18.   
  19.     m_pImageData->GetSpacing(spaces);  
  20.   
  21.    
  22.   
  23.     CResampleDialog dlg(dims[0], dims[1], dims[2], spaces[0], spaces[1], spaces[2]);  
  24.   
  25.     if(dlg.DoModal() == IDOK)  
  26.   
  27.     {  
  28.   
  29.         //get the new dimension and space information  
  30.   
  31.         int dx = dlg.m_iDx;  
  32.   
  33.         int dy = dlg.m_iDy;  
  34.   
  35.         int dz = dlg.m_iDz;  
  36.   
  37.    
  38.   
  39.         float sx = dlg.m_fSx;  
  40.   
  41.         float sy = dlg.m_fSy;  
  42.   
  43.         float sz = dlg.m_fSz;  
  44.   
  45.    
  46.   
  47.         vtkSmartPointer<vtkImageResample> resample =   
  48.   
  49.             vtkSmartPointer<vtkImageResample>::New();  
  50.   
  51.    
  52.   
  53.         resample->SetInput(m_pImageData);  
  54.   
  55.         resample->SetAxisOutputSpacing(0, sx);  
  56.   
  57.         resample->SetAxisOutputSpacing(1, sy);  
  58.   
  59.         resample->SetAxisOutputSpacing(2, sz);  
  60.   
  61.         resample->SetInterpolationModeToLinear();  
  62.   
  63.         resample->Update();  
  64.   
  65.    
  66.   
  67.         m_pImageData = NULL;  
  68.   
  69.         m_pImageData = resample->GetOutput();  
  70.   
  71.         UpdateAllViews(0);  
  72.   
  73.     }  
  74.   
  75. }  


首先根据当前图像的维数和像素间隔信息初始化重采样对话框，当用户设置完毕新的图像维数和间隔并确定后，开始执行重采样。重采样采用vtkImageResample，该类继承自vtkImageReslice。vtkImageReslice在图像处理章节中介绍过，可以提取图像的任意切面，是一个功能非常强大的类。vtkImageResample中通过SetAxisOutputSpacing()函数设置坐标轴的像素宽度，如0对应x轴，1对应y轴，2对应z轴。而函数SetInterpolationModeToLinear()则可以设置采样的插值函数类型，这里采用线性插值方式。最后将执行的结果进行保存，并更新视图显示新的采样图像。但是，当更新视图后，我们发现图像显示的大小并没有发生变化。这是因为，在改变图像维数的时候，像素间隔也随之改变。维数增大，像素间隔变小；反之，像素间隔增大。而图像的物理大小（维数*像素间隔）是不变的，所以在显示的时候，图像是没有改变的。这与一般的图像采样有所区别。不过，当我们把图像保存为一般的图像格式时，便可以看出变化了。

 

6.  保存图像

 

重新响应菜单File中的Save菜单，实现保存当前图像。这里对保存图像进行分析，获取文件的后缀，根据不同的格式调用相应的类进行保存，主要支持bmp,jpg,png,tif,mhd和mha。

 

[cpp] view plaincopy
  1. void CvtkSDIDoc::OnFileSave()  
  2. {  
  3.     // TODO: Add your command handler code here  
  4.   
  5.     if(m_pImageData == NULL) return;  
  6.   
  7.     CString FilePathName;  
  8.   
  9.     CFileDialog dlg(FALSE, NULL, NULL,  
  10.         OFN_HIDEREADONLY | OFN_OVERWRITEPROMPT,  
  11.         (LPCTSTR)_TEXT("All Files (*.*)|*.*||"),  
  12.         NULL);  
  13.   
  14.     if(dlg.DoModal()==IDOK)  
  15.     {  
  16.         FilePathName=dlg.GetPathName();  
  17.         m_ImageFormat =   
  18. FilePathName.Right(FilePathName.GetLength()-FilePathName.ReverseFind('.')-1);  
  19.         if(m_ImageFormat == "") m_ImageFormat = "mhd";  
  20.   
  21.         if(m_ImageFormat.Compare("bmp") == 0)  
  22.         {  
  23.             vtkSmartPointer<vtkBMPWriter> writer   
  24. = vtkSmartPointer<vtkBMPWriter>::New();  
  25.             writer->SetFileName(FilePathName.GetBuffer(0));  
  26.             writer->SetInput(m_pImageData);  
  27.             writer->Write();  
  28.         }  
  29.   
  30.         if(m_ImageFormat.Compare("jpg") == 0)  
  31.         {  
  32.             vtkSmartPointer<vtkJPEGWriter> writer   
  33. = vtkSmartPointer<vtkJPEGWriter>::New();  
  34.             writer->SetFileName(FilePathName.GetBuffer(0));  
  35.             writer->SetInput(m_pImageData);  
  36.             writer-> Write ();  
  37.         }  
  38.   
  39.         if(m_ImageFormat.Compare("tif") == 0)  
  40.         {  
  41.             vtkSmartPointer<vtkTIFFWriter> writer   
  42. = vtkSmartPointer<vtkTIFFWriter>::New();  
  43.             writer->SetFileName(FilePathName.GetBuffer(0));  
  44.             writer->SetInput(m_pImageData);  
  45.             writer-> Write ();  
  46.         }  
  47.   
  48.         if(m_ImageFormat.Compare("png") == 0)  
  49.         {  
  50.             vtkSmartPointer<vtkPNGWriter> writer   
  51. = vtkSmartPointer<vtkPNGWriter>::New();  
  52.             writer->SetFileName(FilePathName.GetBuffer(0));  
  53.             writer->SetInput(m_pImageData);  
  54.             writer->Update();  
  55.         }  
  56.   
  57.         if(m_ImageFormat.Compare("mhd") == 0)  
  58.         {  
  59.             vtkSmartPointer<vtkMetaImageWriter> writer   
  60. = vtkSmartPointer<vtkMetaImageWriter>::New();  
  61.             writer->SetFileName(FilePathName.GetBuffer(0));  
  62.             writer->SetInput(m_pImageData);  
  63.             writer-> Write ();  
  64.         }  
  65.   
  66.     }  
  67.     else  
  68.     {  
  69.         return;  
  70.     }  
  71. }  


下图中左侧图像的大小为1024*768，采用本篇的程序进行重采样为512*768，保存的结果如右图所示。可以看出图像在X方向上维数减少一半。

 

     

 

到此，基于VTK的图像重采样程序已经实现完毕。该程序支持基本的图像格式，并可以将当前图像根据需要进行存储，存储格式目前支持bmp,jpg,png,tif,mhd和mha。对图像进行重采样时，需要选择改变维数或者改变像素间隔，满足用户不同的需求。程序执行效果如下。

 

 

==========欢迎转载，转载时请保留该声明信息==========

版权归@东灵工作室所有，更多信息请访问东灵工作室


教程系列导航：http://blog.csdn.net/www_doling_net/article/details/8763686

================================================

0 0
海岛奇兵保险库及数据 海岛奇兵保险库及数据
原创粉丝点击
热门IT博客
手机淘宝店铺发货流程手机淘宝店铺发货流程
女生下体有味道 知乎女生下体有味道 知乎
ti6双败赛制 知乎
网络侦探杰斯兽
淘宝客 完全开源
mac怎么打开html文件
海淀区单片机竞赛成绩
淘宝怎么看店铺排名
网络英语培训
快手绝世的容颜淘宝店
阿里云幕布邮寄要多久
北京北大软件上班
哈利波特怪味豆 淘宝
淘宝下架再上架
微商管理系统源码
java用户注册界面代码
美工个人总结
算法导论 第三版
c语言现代编程方法
日本化妆品推荐知乎
热门问题 老师的惩罚 人脸识别 我在镇武司摸鱼那些年 重生之率土为王 我在大康的咸鱼生活 盘龙之生命进化 天生仙种 凡人之先天五行 春回大明朝 姑娘不必设防，我是瞎子 桂枝香王安石 桂枝香 王安石桂枝香 桂枝香 金陵怀古 青桐树介绍 青桐树 桂技茯苓丸的功效与作用 参桂百补丸 大沥桂江三鸟批发市场 涂桂洪 桂浩明 桂源铺奶茶 桂源铺奶茶怎么样 桂源铺奶茶店 爽洋洋 满陇桂雨公园 满觉陇赏桂 满陇桂雨公园门票 满觉陇桂花 下满觉陇 上海桂满陇 桂满陇团购 桂珍紫砂壶价格 桂由美婚纱摄影 桂由美 西安清桂电子科技有限公司 桂皮树 桂皮醛 桂盛卡 广西农村信用社桂盛卡 广西农村信用社手机银行下载 桂盛卡是什么卡 广西农村信用社网上银行 桂单竹 绥中新桂竹宾馆 桂竹笋 桂的组词 桂字组词 桂的组词是什么 桂组词有哪些 渠组词

程序博客网，程序员的互联网技术博客家园。csdn论坛精品 msdn技术资料都在这里