VTK系列71_VTK对MHA格式文件三维感兴趣区域裁剪及MHA格式保存

今天药忘吃喽~ 2022-09-03 04:15 25阅读 0赞

# **实例12：基于VTK对MHA格式文件三维感兴趣区域裁剪及MHA格式保存** #

#include "vtkAutoInit.h" 
    VTK_MODULE_INIT(vtkRenderingOpenGL2);
    VTK_MODULE_INIT(vtkInteractionStyle);
    
    #include <vtkSmartPointer.h>
    #include <vtkImageData.h>
    #include <vtkRenderWindow.h>
    #include <vtkRenderWindowInteractor.h>
    #include <vtkInteractorStyleImage.h>
    #include <vtkRenderer.h>
    #include <vtkImageActor.h>
    #include <vtkBMPReader.h>
    #include <vtkExtractVOI.h>
    #include <vtkMetaImageReader.h>//mha、mad图像读取类
    #include <vtkMetaImageWriter.h>mha、mad图像写入类
    
    //测试图像：../data/lena.bmp
    int main(int argc, char* argv[])
    {
    	/*vtkSmartPointer<vtkBMPReader> reader =
    		vtkSmartPointer<vtkBMPReader>::New();
    	reader->SetFileName("lena.bmp");
    	reader->Update();*/
    
    	vtkSmartPointer<vtkMetaImageReader> reader = vtkSmartPointer<vtkMetaImageReader>::New();
    	//mhd与mha文件其实格式是一样的，记录mhd对应的raw文件应在同一目录
    	//mhd格式图像信息头与实际图像的存储分为两个文件（*.mhd文件记录图像信息头；*.raw或//*.zraw（zraw指有压缩）记录实际图像）
    	//mha格式将图像信息头与实际的像素值等数据写入到同一个文件中
    	//reader->SetFileName("test1.mhd");
    	  //.mha和.raw文件需要在同一个文件夹
    	reader->SetFileName("CT_6_spacing.mha");
    	reader->Update();
    
    	int dims[3];
    	reader->GetOutput()->GetDimensions(dims);
    
    	vtkSmartPointer<vtkExtractVOI> extractVOI =
    		vtkSmartPointer<vtkExtractVOI>::New();
    	extractVOI->SetInputConnection(reader->GetOutputPort());
    	//设置感兴趣区域：X_min、X_max、Y_min、Y_max、Z_min、Z_max
    	extractVOI->SetVOI(dims[0] / 4., 3. * dims[0] / 4., dims[1] / 4., 3. * dims[1] / 4., 0, 5);//这里设置0~5，表示裁剪六张切片
    	extractVOI->Update();
    
    	vtkSmartPointer<vtkImageActor> originalActor =
    		vtkSmartPointer<vtkImageActor>::New();
    	originalActor->SetInputData(reader->GetOutput());
    
    	vtkSmartPointer<vtkImageActor> voiActor =
    		vtkSmartPointer<vtkImageActor>::New();
    	voiActor->SetInputData(extractVOI->GetOutput());
    
    	double originalViewport[4] = { 0.0, 0.0, 0.5, 1.0 };
    	double voiviewport[4] = { 0.5, 0.0, 1.0, 1.0 };
    
    	vtkSmartPointer<vtkRenderer> originalRenderer =
    		vtkSmartPointer<vtkRenderer>::New();
    	originalRenderer->SetViewport(originalViewport);
    	originalRenderer->AddActor(originalActor);
    	originalRenderer->ResetCamera();
    	originalRenderer->SetBackground(1.0, 1.0, 1.0);
    
    	vtkSmartPointer<vtkRenderer> shiftscaleRenderer =
    		vtkSmartPointer<vtkRenderer>::New();
    	shiftscaleRenderer->SetViewport(voiviewport);
    	shiftscaleRenderer->AddActor(voiActor);
    	shiftscaleRenderer->ResetCamera();
    	shiftscaleRenderer->SetBackground(1.0, 1.0, 1.0);
    
    	vtkSmartPointer<vtkRenderWindow> renderWindow =
    		vtkSmartPointer<vtkRenderWindow>::New();
    	renderWindow->AddRenderer(originalRenderer);
    	renderWindow->AddRenderer(shiftscaleRenderer);
    	renderWindow->SetSize(900, 300);
    	renderWindow->Render();
    	renderWindow->SetWindowName("ExtractVOIExample");
    
    	vtkSmartPointer<vtkRenderWindowInteractor> renderWindowInteractor =
    		vtkSmartPointer<vtkRenderWindowInteractor>::New();
    	vtkSmartPointer<vtkInteractorStyleImage> style =
    		vtkSmartPointer<vtkInteractorStyleImage>::New();
    
    	//保存为mhd文件
    	//vtkMetaImageWriter* vtkWriter = vtkMetaImageWriter::New();
    	//vtkWriter->SetInputConnection(extractVOI->GetOutputPort());
    	以raw和mhd格式保存，去掉此句则以zraw保存
    	//vtkWriter->SetCompression(0);
    	//vtkWriter->SetFileName("BrainProtonDensity3Slices_1.mhd");
    	//vtkWriter->Write();
    
    	保存为mha文件
    	vtkMetaImageWriter* vtkWriter = vtkMetaImageWriter::New();
    	vtkWriter->SetInputConnection(extractVOI->GetOutputPort());
    	//以raw和mhd格式保存，去掉此句则以zraw保存
    	vtkWriter->SetCompression(0);
    	vtkWriter->SetFileName("CT_6_spacing_cut.mha");
    	vtkWriter->Write();
    
    	renderWindowInteractor->SetInteractorStyle(style);
    	renderWindowInteractor->SetRenderWindow(renderWindow);
    	renderWindowInteractor->Initialize();
    	renderWindowInteractor->Start();
    
    
    
    	return EXIT_SUCCESS;
    }

读取的MHA源文件:

![20210804105649130.png][]

选定区域剪切后保存的MHA文件：

![20210804105804604.png][]![watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzMyODA5MDkz_size_16_color_FFFFFF_t_70][]

其中前四张剪切前后效果对比：

![watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzMyODA5MDkz_size_16_color_FFFFFF_t_70 1][] ![watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzMyODA5MDkz_size_16_color_FFFFFF_t_70 2][] ![watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzMyODA5MDkz_size_16_color_FFFFFF_t_70 3][] ![watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzMyODA5MDkz_size_16_color_FFFFFF_t_70 4][]

![watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzMyODA5MDkz_size_16_color_FFFFFF_t_70 5][] ![watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzMyODA5MDkz_size_16_color_FFFFFF_t_70 6][] ![watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzMyODA5MDkz_size_16_color_FFFFFF_t_70 7][] ![watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzMyODA5MDkz_size_16_color_FFFFFF_t_70 8][]

注：此例程MHD和NII格式文件同样适用。

**[本例程配套素材见源码整理文章下载（点击进入）][Link 1]**

## **VTK系列目录：** ##

## **[1 VTK基本概念][1 VTK]** ##

## **[2 VTK图像处理][2 VTK]** ##

## **[3 VTK图形处理][3 VTK]** ##

## **[4 VTK体绘制][4 VTK]** ##

[20210804105649130.png]: /images/20220829/4676d09e6ca94a458c4198caad65865f.png
[20210804105804604.png]: https://img-blog.csdnimg.cn/20210804105804604.png
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzMyODA5MDkz_size_16_color_FFFFFF_t_70]: /images/20220829/5b07d187bed547319396a604fe13272e.png
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzMyODA5MDkz_size_16_color_FFFFFF_t_70 1]: https://img-blog.csdnimg.cn/20210804110911139.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzMyODA5MDkz,size_16,color_FFFFFF,t_70
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzMyODA5MDkz_size_16_color_FFFFFF_t_70 2]: https://img-blog.csdnimg.cn/20210804110956613.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzMyODA5MDkz,size_16,color_FFFFFF,t_70
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzMyODA5MDkz_size_16_color_FFFFFF_t_70 3]: https://img-blog.csdnimg.cn/20210804111014136.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzMyODA5MDkz,size_16,color_FFFFFF,t_70
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzMyODA5MDkz_size_16_color_FFFFFF_t_70 4]: /images/20220829/881f5b6eeb0a4fa888d00c356d6f0fee.png
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzMyODA5MDkz_size_16_color_FFFFFF_t_70 5]: https://img-blog.csdnimg.cn/20210804110002205.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzMyODA5MDkz,size_16,color_FFFFFF,t_70
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzMyODA5MDkz_size_16_color_FFFFFF_t_70 6]: https://img-blog.csdnimg.cn/20210804110021171.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzMyODA5MDkz,size_16,color_FFFFFF,t_70
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzMyODA5MDkz_size_16_color_FFFFFF_t_70 7]: https://img-blog.csdnimg.cn/20210804105928749.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzMyODA5MDkz,size_16,color_FFFFFF,t_70
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzMyODA5MDkz_size_16_color_FFFFFF_t_70 8]: /images/20220829/51537ff3cda94ee79a8b77ea4f4a6306.png
[Link 1]: https://blog.csdn.net/qq_32809093/article/details/115689942
[1 VTK]: https://blog.csdn.net/qq_32809093/article/details/105380666
[2 VTK]: https://blog.csdn.net/qq_32809093/article/details/107083519
[3 VTK]: https://blog.csdn.net/qq_32809093/article/details/107083554
[4 VTK]: https://blog.csdn.net/qq_32809093/article/details/107126102