ITK对DICOM数据的读与写

最新推荐文章于 2024-06-16 08:00:00 发布

cuihaolong

最新推荐文章于 2024-06-16 08:00:00 发布

阅读量8.1k

点赞数 5

分类专栏： 3D Print

本文链接：https://blog.csdn.net/cuihaolong/article/details/66968866

版权

本文介绍了如何使用ITK库进行DICOM数据的读写，包括处理DICOM序列和单张数据，同时探讨了在操作过程中输入输出类型转换的注意事项，以及fastMarching算法的应用问题。

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

1.读写DICOM序列

将一系列DICOM数据读取形成三维体数据，经过所需的处理之后再将其写回成二维DICOM数据，其中需要GDCMImageIO作为DICOM格式的ImageIO。

读写过程如下：

//IO相关
#include "itkGDCMImageIO.h"
#include "itkGDCMSeriesFileNames.h"
#include "itkImageSeriesReader.h"
#include "itkImageSeriesWriter.h"

int main(int argc, char *argv[])
{
	
	typedef   float           InternalPixelType;
	const     unsigned int    Dimension = 3;//输入为三维
	typedef itk::Image< InternalPixelType, Dimension >  InternalImageType;
	
	typedef  float          OutputPixelType;
	const unsigned int      OutputDimension = 2;//输出为二维
	typedef itk::Image< OutputPixelType, OutputDimension > OutputImageType;

	typedef  itk::ImageSeriesReader< InternalImageType > ReaderType;
	typedef  itk::ImageSeriesWriter<  InternalImageType,OutputImageType  > WriterType;//自动在输入和输出之间转换

	ReaderType::Pointer reader = ReaderType::New();
	WriterType::Pointer writer = WriterType::New();

	//DICOM数据的ImageIO
	typedef itk::GDCMImageIO             ImageIOType;
	ImageIOType::Pointer gdcmIO = ImageIOType::New();

	//文件生成器
	typedef itk::GDCMSeriesFileNames     NamesGeneratorType;
	NamesGeneratorType::Pointer namesGenerator = NamesGeneratorType::New();
	namesGenerator->SetInputDirectory("D:/VTKCode/Data/Threshold");

	//保存文件名
	const ReaderType::FileNamesContainer & filenames =
		na