ITK和VTK读取DICOM图像文件

ITK和VTK读取DICOM图像文件

ITK读取DICOM图像

相比于VTK类库中vtkDICOMImageReader类读取DICOM序列图像,借助ITK类库实现对DICOM序列图像的读取要复杂许多。但是,使用ITK类库读取图像不像使用VTK类库那么功能局限,VTK类库对每种图像格式都有着专门的读写类,像vtkDICOMImageReader类只能读取普通的DICOM图像,对于压缩后的DICOM图像则只能通过ITK类库来实现读取。
在这里插入图片描述
当系统使用ITK来进行读取DICOM格式文件时,基于ITK类库中itkGDCMSeriesFileNames.h头文件,创建一个itk::GDCMSeriesFileNames类型的智能指针,该指针可用于读取并保存DICOM序列图像的文件夹路径,通过该类中的静态函数SetDirectory()实现。基于itkImageSeriesReader.h头文件,声明和实例化读写 DICOM 序列图像的itk::ImageSeriesReader对象,该指针主要负责在指定的文件夹路径下读取DICOM序列图像,通过该类中的静态函数SetFileNames()实现。基于itkImageToVTKImageFilter.h头文件,声明和实例化ITK到VTK转换通道的itk::ImageToVTKImageFilter对象,借助该类中的静态函数SetInput()和Update()实现。基于vtkImageFlip类,负责将从ITK转换到VTK的数据保存供下一步操作,借助该类中的静态函数SetInputData()和Update()实现。

代码实现

头文件

#include "itkImageToVTKImageFilter.h"
#include "itkImage.h"
#include "itkGDCMImageIO.h"
#include "itkGDCMSeriesFileNames.h"
#include "itkImageSeriesReader.h"

宏定义

using namespace std;
typedef short PixelType;
const unsigned int   Dimension = 2;
typedef itk::Image< PixelType, Dimension > Input2dImageType;
typedef itk::Image< PixelType, 3 > Input3dImageType;typedef itk::Image< PixelType, Dimension > Output2dImageType;
typedef itk::Image< PixelType, 3 > Output3dImageType;typedef itk::GDCMImageIO   ImageIOType;//GDCMImageIO读DICOM
ImageIOType::Pointer gdcmImageIO = ImageIOType::New();using ReaderType3d = itk::ImageSeriesReader< Input3dImageType >;
ReaderType3d::Pointer reader3d = ReaderType3d::New();typedef itk::ImageFileReader< Input2dImageType > ReaderType2d;
ReaderType2d::Pointer reader2d = ReaderType2d::New();using NamesGeneratorType = itk::GDCMSeriesFileNames;
NamesGeneratorType::Pointer nameGenerator = NamesGeneratorType::New();

主要实现

   reader3d->SetImageIO(gdcmImageIO);nameGenerator->SetUseSeriesDetails(true);nameGenerator->SetDirectory(file_path_str);using SeriesIdContainer = std::vector< std::string >;const SeriesIdContainer& seriesUID = nameGenerator->GetSeriesUIDs();auto seriesItr = seriesUID.begin();auto seriesEnd = seriesUID.end();using FileNamesContainer = std::vector< std::string >;FileNamesContainer fileNames;std::string seriesIdentifier;while (seriesItr != seriesEnd){seriesIdentifier = seriesItr->c_str();fileNames = nameGenerator->GetFileNames(seriesIdentifier);++seriesItr;}reader3d->SetFileNames(fileNames);reader3d->Update();typedef itk::ImageToVTKImageFilter< Input3dImageType> itkTovtkFilterType;itkTovtkFilterType::Pointer itkTovtkImageFilter = itkTovtkFilterType::New();itkTovtkImageFilter->SetInput(reader3d->GetOutput());//设置图像数据从ITK转向VTKitkTovtkImageFilter->Update();vtkSmartPointer< vtkImageFlip > reader1 = vtkSmartPointer< vtkImageFlip >::New();reader1->SetInputData(itkTovtkImageFilter->GetOutput());reader1->SetFilteredAxes(1);reader1->Update();

VTK读取DICOM

VTK作为适用于图像处理和三维重建的开源可视化软件包,在其类库中封装了很多种格式图像的读取类。但由于上节描述可知,使用VTK类库读取DICOM序列图像存在着一些功能局限。VTK类库对每种图像格式都有着专门的读写类,像vtkDICOMImageReader类只能读取普通的DICOM图像,针对压缩后的DICOM格式图像则无法有效读取,使用VTK类库实现图像读取实现简单。
基于VTK类库中图像的读取类实现对DICOM图像的读取。首先,通过vtkDICOMImageReader类中的静态函数SetDirectoryName()读取在磁盘中存储DICOM格式文件的相对路径。使用vtkImageViewer2类中的静态函数获取系统读取的DICOM图像给显示窗口,定义显示窗口的窗位和窗宽大小,进而实现对DICOM图像的可视化任务。对DICOM单张图像的读取和显示效果图如图所示,对DICOM序列图像的读取会在下章节三维重建算法实现中详细表述。
在这里插入图片描述

代码实现

VTK关于DICOM图像文件的读取实现相比于ITK要简单很多。

vtkSmartPointer<vtkDICOMImageReader> DICOMreader = vtkSmartPointer<vtkDICOMImageReader>::New();DICOMreader->SetDirectoryName(file_path_str);//DICOM所在文件夹DICOMreader->Update();

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.rhkb.cn/news/52931.html

如若内容造成侵权/违法违规/事实不符,请联系长河编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

VTK使用

入门参考的这篇&#xff08;修改了很多&#xff09;&#xff1a;QT加载显示DICOM文件浏览 但还是必须在main()中先读文件再show&#xff0c;先show了再打开文件是黑的显示不出来&#xff08;后续再查&#xff09;。 //.h中要加入 #include <vtkAutoInit.h> VTK_MODULE_I…

VTK教程3--------打开vtk文件

下面&#xff0c;本教程将介绍如何在MFC下打开vtk文件&#xff0c;最终的效果如图所示&#xff1a; 如果你看过一些三维重建相关的论文&#xff0c;那么一定对这个图像很熟悉。很多论文都会以这个图像作为例子讲解&#xff0c;好了&#xff0c;闲话少说&#xff0c;直接讲实现…

Android中读取电话本Contacts联系人的所有电话号信息

1.首先&#xff0c;要知道android 的contacts里的电话信息有多类&#xff1a;moblie,家庭&#xff0c;工作&#xff0c;传真等。如图&#xff1a; 2.android的Contacts是通过ContentProvider来提供的&#xff0c;其实android把contacts和SMS给组织成数据库文件了&#xff0c;你…

虚拟化 : VT-x VT-d VT-c的分别

虚拟化 &#xff1a; VT-x VT-d VT-c的分别: VT-x [运行ESXI上的64bit Guest OS基本指令] Intel运用Virtualization虚拟化技术中的一个指令集。VT-x有助于提高基于软件的虚拟化解决方案的灵活性与稳定性。通过按照纯软件虚拟化的要求消除虚拟机监视器(VMM&#xff09;代表…

VTK-vtkPolyDataNormals

前言&#xff1a;本博文主要记录vtkPolyDataNormals的用途&#xff0c;使用方法以及实现原理&#xff0c;帮助更多的小伙伴更好的理解和应用vtk。 目录 1. vtkPolyDataNormals 2. 使用示例 3. 实现原理 1. vtkPolyDataNormals vtkPolyDataNormals可以计算多边形网格数据的点…

手机联系人批量导入(字符编码、xlsx与vcf互转)

一、前言 单位同事每人都办理了一个192新手机号&#xff0c;一打过来不知道是谁&#xff0c;又懒的一个个保存姓名。一想是不是可以批量导入呢?电子表格号码名单我有啊&#xff0c;试试吧。 二、实践 1.先下载手机联系人文件看看吧。在手机联系人设置中有个导出选项&#x…

信创干部人事档案管理系统单机版 - 人力资源档案管理系统软件

信创干部人事档案管理系统单机版v2.0&#xff08;以下简称系统&#xff09;&#xff0c;是一套具有先进性、安全性、前瞻性的人力资源档案管理系统&#xff0c;是在总结近二十年为万余家单位档案信息化建设实践经验的基础上&#xff0c;遵循ISO15489、ISO23081、ISO14721等国际…

VTK-vtkFieldData

欢迎大家加入社区&#xff0c;雪易VTK社区-CSDN社区云 前言&#xff1a;为区分vtkPoints和vtkPointData的区别&#xff0c;了解vtkFieldData在VTK中的存在意义&#xff0c;从而系统的掌握vtk中关于数据的表达方式。 vtk中通过vtkDataArray进行数据的存储&#xff0c;通过vtkD…

虚拟内存统计----Vmstat命令

目录 一、Vmstat命令 概述 1.1 物理内存和虚拟内存区别 1.1.1 物理内存 1.1.2 虚拟内存 1.2 vmstat 命令语法 1.2.1 格式 1.2.2 命令参数 二、 vmstat 示例 2.1 显示虚拟内存使用情况 2.1.1 查看 2.1.2 表示在1秒时间内进行2次采样 2.1.3 指定的MB 单位输…

票据机器人-OCR自动识别助力智能财务

当前票据的数字化问题尚未完全解决&#xff0c;人工处理海量票据的工作量较大&#xff0c;把票据录入系统、校验真伪的成本较高。 票据机器人对财务报销环节涉及的增票、火车票、飞机行程单、出租车票、定额发票、汽车票等票据进行自动分类、识别录入、验真及归档。 票据机器人…

第一次申请去美国面签,需要注意哪些事项提高成功率?

Teresa &#xff0c;PM 296 人赞同 哈哈哈&#xff0c;作为一个刚毕业的无房无车无钱妹纸&#xff0c;刚顺利拿到美国旅游十年签证&#xff01;我要去西雅图看TI5啦&#xff01;整个过程耗时两周~~~热乎乎的攻略奉上~~ 所有问题都可在官网进行查询&#xff0c;请熟读此网站&am…

java行程单解析获取内容_java如何抓取网页上的动态信息,获取源代码后如何分析JS?...

首先明确我指的动态数据是什么。 名词定义&#xff1a;动态数据在这里指的是网页中由Javascript动态生成的页面内容&#xff0c;即网页源文件中没有&#xff0c;在页面加载到浏览器后动态生成的。 下面进入正题。 抓取静态页面很简单&#xff0c;通过Java获取到html源码&#x…

系统采集--行程单、火车票、汽车票方案

方案一&#xff1a;进项旅客运输采集&#xff08;行程单、火车票、汽车票&#xff09; 实现步骤&#xff1a; 建立视图&#xff08;view_trip、view_trip、view_car&#xff09;&#xff0c;其中视图中税号需要关联组织表获取&#xff1b;后台&#xff1a;根据当前组织id去组…

百望云入选电子凭证会计数据标准深化试点服务保障单位 打通数字化“最后一公里”

电子凭证试点正在加速推进中&#xff01; 近日&#xff0c;财政部会同税务总局、人民银行、国务院国资委、国家档案局、标准委、国电联办、民航局、国铁集团9部委&#xff0c;组织召开电子凭证会计数据标准深化试点启动会。 之后&#xff0c;财政部会计司相继发布《关于公布电子…

【云驻共创】华为云文字识别服务的体验之旅

文章目录 前言一、初识OCR技术1.文字识别技术的概念1.1 文字识别&#xff08;OCR&#xff09;介绍1.2 文字识别&#xff08;OCR&#xff09;应用场景 2.文字识别技术的发展历程2.1 OCR识别领域发展历程2.1.1 西文OCR2.1.1 Tesseract2.1.3 LeNet2.1.4 深度学习OCR 2.2 OCR识别流…

Python自动化办公(一) :滴滴行程单信息提取存入excel表格 2021-04-13

ython自动化办公&#xff08;一&#xff09;:滴滴行程单信息提取存入excel表格 提示&#xff1a;这里可以添加系列文章的所有文章的目录&#xff0c;目录需要自己手动添加 例如&#xff1a;第一章 Python 机器学习入门之pandas的使用 提示&#xff1a;写完文章后&#xff0c;目…

智能客服“听不懂人话”?消费者很“闹心”

智能客服与人工客服并非互相取代的关系。人工客服不能缺位&#xff0c;应通过人机协同更好地回应消费者诉求。在金融消费领域&#xff0c;建议与资金安全等密切相关的业务能够设置人工服务“一键转接”或者“一键冻结”功能&#xff0c;保障消费者财产安全&#xff1b;为“银发…

让ChatGPT调用10万+开源AI模型!HuggingFace新功能爆火:大模型可随取随用多模态AI工具...

萧箫 发自 凹非寺量子位 | 公众号 QbitAI 只需和ChatGPT聊聊天&#xff0c;它就能帮你调用10万个HuggingFace模型&#xff01; 这是抱抱脸最新上线的功能HuggingFace Transformers Agents&#xff0c;一经推出就获得极大关注&#xff1a; 这个功能&#xff0c;相当于给ChatGPT等…

chatgpt赋能python:用Python编写录音代码:打造高效音频处理工具

用Python编写录音代码&#xff1a;打造高效音频处理工具 Python已经成为数据科学、人工智能和机器学习的首选编程语言&#xff0c;因此它也是制作音频处理工具的完美选择。Python编程语言简单易学&#xff0c;并能够快速实现高效的音频处理。通过使用Python&#xff0c;您可以…

如何把图片去水印?来试试这三种图片去水印的方法

如何把图片去水印&#xff1f;来试试这三种图片去水印的方法 相信大家在日常上网的时候&#xff0c;看到一些漂亮的图片时&#xff0c;想将它们保存下来&#xff0c;但又碍于图片带有水印不得不放弃。但其实我们是可以通过一些方法将图片中的水印给去除掉&#xff0c;那我们该…