基于VTK和ITK的MRI图像（头皮）三维重建代码（包括各项异性滤波+ostu+形态学+移动立方体算法）

最新推荐文章于 2024-05-16 12:31:22 发布

HIT_KyleChen

最新推荐文章于 2024-05-16 12:31:22 发布

阅读量2.1k

点赞数 3

本文链接：https://blog.csdn.net/vjhghjghj/article/details/84335463

版权

本文介绍了使用VTK和ITK库进行MRI图像的三维重建过程，具体涉及Marching Cubes算法，包括了各项异性滤波、Otsu阈值处理、形态学操作以及移动立方体算法的应用，详细解析了源代码`MarchingCubes.cxx`的实现步骤，并提供了`CMakeLists.txt`配置文件以供项目构建。

摘要由CSDN通过智能技术生成

MarchingCubes

build

MarchingCubes.cxx

/*
头皮面绘制：插值+各项异性滤波+ostu+形态学+移动立方体
处理后的文件存放在processed_mri_data02
*/
#include "itkGDCMImageIO.h"
#include "itkGDCMSeriesFileNames.h"
#include "itkImageSeriesReader.h"
#include "itkImageSeriesWriter.h"
#include "itkImageFileReader.h"
#include "itkImageFileWriter.h"
#include "itkImage.h"
#include "itkImageToVTKImageFilter.h"
#include "itkGDCMImageIO.h"
#include "itkGDCMSeriesFileNames.h"
#include "itkImageSeriesReader.h"
#include <vtkSmartPointer.h>
#include "itkOtsuThresholdImageFilter.h"
#include <vtkDICOMImageReader.h>
#include <vtkVersion.h>
#include <vtkSmartPointer.h>
#include <vtkMarchingCubes.h>
#include <vtkVoxelModeller.h>
#include <vtkSphereSource.h>
#include <vtkImageData.h>
#include <vtkBMPReader.h>
#include <vtkDICOMImageReader.h>
#include <vtkActor.h>
#include <vtkPolyDataMapper.h>
#include <vtkRenderWindowInteractor.h>
#include <vtkRenderWindow.h>
#include <vtkRenderer.h>
#include <vtkStringArray.h>
#include<vtkImageAppend.h>
#include <vtkSmartPointer.h>
#include <vtkObjectFactory.h>
#include <vtkRenderWindow.h>
#include <vtkRenderWindowInteractor.h>
#include <vtkRenderer.h>
#include <vtkActor.h>
#include <vtkImageViewer2.h>
#include <vtkDICOMImageReader.h>
#include <vtkInteractorStyleImage.h>
#include <vtkActor2D.h>
#include <vtkTextProperty.h>
#include <vtkTextMapper.h>
#include <sstream>
#include <iostream>
#include <vtkStringArray.h>
#include<vtkImageAppend.h>
#include "itkGradientAnisotropicDiffusionImageFilter.h"
#include "itkRescaleIntensityImageFilter.h"
#include "itkBinaryMorphologicalClosingImageFilter.h"
#include "itkBinaryBallStructuringElement.h"
using namespace std;
/*在显示图像时会用到窗口信息类*/
class  StatusMessage
{
public:
    static std::string Format(int slice, int maxSlice)
    {
        std::stringstream tmp;
        tmp << "Slice Number " << slice + 1 << "/" << maxSlice + 1;
        return tmp.str();
    }

};

class myvtkInteractorStyleImage: public vtkInteractorStyleImage
{
public:
    static myvtkInteractorStyleImage* New();
    vtkTypeMacro(myvtkInteractorStyleImage, vtkInteractorStyleImage);

protected:
    vtkImageViewer2* _ImageViewer;
    vtkTextMapper* _StatusMapper;
    int _Slice;
    int _MinSlice;
    int _MaxSlice;

public:
    void SetImageViewer(vtkImageViewer2* imageViewer) {
        _ImageViewer = imageViewer;
        _MinSlice = imageViewer->GetSliceMin();
        _MaxSlice = imageViewer->GetSliceMax();
        _Slice = _MinSlice;
        cout << "Slice: Min = " << _MinSlice << ", Max = " << _MaxSlice << endl;
    }

    void SetStatusMapper(vtkTextMapper* statusMapper)
    {
        _StatusMapper = statu