我就废话不多说了,大家还是直接看代码吧~
## using simpleITK to load and save data. import SimpleITK as sitk itk_img = sitk.ReadImage('./nifti.nii.gz') img = sitk.GetArrayFromImage(itk_img) print("img shape:",img.shape) ## save out = sitk.GetImageFromArray(img) # # out.SetSpacing(itk_img.GetSpacing()) # # out.SetOrigin(itk_img.GetOrigin()) sitk.WriteImage(out,'simpleitk_save.nii.gz')
读取DICOM序列
import sys import SimpleITK as sitk import numpy as np dcm_directory = "./flair_dicom" series_ids = sitk.ImageSeriesReader.GetGDCMSeriesIDs(dcm_directory) print("series ids:",series_ids) if not series_ids: print("ERROR: given directory dose not a DICOM series.") sys.exit(1) series_file_names = sitk.ImageSeriesReader.GetGDCMSeriesFileNames(dcm_directory,series_ids[0]) series_reader = sitk.ImageSeriesReader() series_reader.SetFileNames(series_file_names) image3D = series_reader.Execute() size = image3D.GetSize() print( "Image size:", size[0], size[1], size[2] ) image_array = sitk.GetArrayFromImage(image3D) # SimpleITK读取的图像数据的坐标顺序为zyx,即从多少张切片到单张切片的宽和高; # 而据SimpleITK Image获取的origin和spacing的坐标顺序则是xyz origin = image3D.GetOrigin() spacing = image3D.GetSpacing() print("origin:",origin," spacing:",spacing)
读取DICOM单个文件
import SimpleITK as sitk import numpy as np image = sitk.ReadImage(slice_path) image_array = sitk.GetArrayFromImage(image) # z, y, x
补充知识:itk 读取图像,两种格式之间转换
我就废话不多说了,大家还是直接看代码吧~
/******************************************************* 控制台运行程序 输入: 程序名称 读取图像文件名 保存图像文件名 输出: 保存另一种类型的图像 功能: 图像在jpeg和bmp两种格式之间互换 ********************************************************/ #include "itkImageFileReader.h" #include "itkImageFileWriter.h" #include "itkImageIOFactory.h" #include "itkBMPImageIOFactory.h" #include "itkJPEGImageIOFactory.h" #include <iostream> int main(int argc , char *argv[]) { if (argc < 3) { std::cerr << "Usage: " << std::endl; std::cerr << argv[0] << " inputImageFile outputImageFile " << std::endl; return EXIT_FAILURE; } itk::BMPImageIOFactory::RegisterOneFactory(); itk::JPEGImageIOFactory::RegisterOneFactory(); typedef itk::Image< unsigned char, 2 > ImageType; typedef itk::ImageFileReader< ImageType > ReaderType; typedef itk::ImageFileWriter< ImageType > WriterType; ReaderType::Pointer reader = ReaderType::New(); WriterType::Pointer writer = WriterType::New(); const char * inputFilename = argv[1]; const char * outputFilename = argv[2]; reader->SetFileName( inputFilename ); writer->SetFileName( outputFilename ); writer->SetInput( reader->GetOutput() ); try { writer->Update(); } catch( itk::ExceptionObject & err ) { std::cerr << "ExceptionObject caught !" << std::endl; std::cerr << err << std::endl; return EXIT_FAILURE; } system("pause"); return 0; }
以上这篇使用SimpleITK读取和保存NIfTI/DICOM文件实例就是小编分享给大家的全部内容了,希望能给大家一个参考,也希望大家多多支持。
免责声明:本站文章均来自网站采集或用户投稿,网站不提供任何软件下载或自行开发的软件!
如有用户或公司发现本站内容信息存在侵权行为,请邮件告知! 858582#qq.com
白云城资源网 Copyright www.dyhadc.com
暂无“使用SimpleITK读取和保存NIfTI/DICOM文件实例”评论...
《魔兽世界》大逃杀!60人新游玩模式《强袭风暴》3月21日上线
暴雪近日发布了《魔兽世界》10.2.6 更新内容,新游玩模式《强袭风暴》即将于3月21 日在亚服上线,届时玩家将前往阿拉希高地展开一场 60 人大逃杀对战。
艾泽拉斯的冒险者已经征服了艾泽拉斯的大地及遥远的彼岸。他们在对抗世界上最致命的敌人时展现出过人的手腕,并且成功阻止终结宇宙等级的威胁。当他们在为即将于《魔兽世界》资料片《地心之战》中来袭的萨拉塔斯势力做战斗准备时,他们还需要在熟悉的阿拉希高地面对一个全新的敌人──那就是彼此。在《巨龙崛起》10.2.6 更新的《强袭风暴》中,玩家将会进入一个全新的海盗主题大逃杀式限时活动,其中包含极高的风险和史诗级的奖励。
《强袭风暴》不是普通的战场,作为一个独立于主游戏之外的活动,玩家可以用大逃杀的风格来体验《魔兽世界》,不分职业、不分装备(除了你在赛局中捡到的),光是技巧和战略的强弱之分就能决定出谁才是能坚持到最后的赢家。本次活动将会开放单人和双人模式,玩家在加入海盗主题的预赛大厅区域前,可以从强袭风暴角色画面新增好友。游玩游戏将可以累计名望轨迹,《巨龙崛起》和《魔兽世界:巫妖王之怒 经典版》的玩家都可以获得奖励。
更新日志
2024年10月06日
2024年10月06日
- 群星《前途海量 电影原声专辑》[FLAC/分轨][227.78MB]
- 张信哲.1992-知道新曲与精丫巨石】【WAV+CUE】
- 王翠玲.1995-ANGEL【新艺宝】【WAV+CUE】
- 景冈山.1996-我的眼里只有你【大地唱片】【WAV+CUE】
- 群星《八戒 电影原声带》[320K/MP3][188.97MB]
- 群星《我的阿勒泰 影视原声带》[320K/MP3][139.47MB]
- 纪钧瀚《胎教古典音乐 钢琴与大提琴的沉浸时光》[320K/MP3][148.91MB]
- 刘雅丽.2001-丽花皇后·EMI精选王【EMI百代】【FLAC分轨】
- 齐秦.1994-黄金十年1981-1990CHINA.TOUR.LIVE精丫上华】【WAV+CUE】
- 群星.2008-本色·百代音乐人创作专辑【EMI百代】【WAV+CUE】
- 群星.2001-同步过冬AVCD【环球】【WAV+CUE】
- 群星.2020-同步过冬2020冀待晴空【环球】【WAV+CUE】
- 沈雁.1986-四季(2012梦田复刻版)【白云唱片】【WAV+CUE】
- 纪钧瀚《胎教古典音乐 钢琴与大提琴的沉浸时光》[FLAC/分轨][257.88MB]
- 《国语老歌 怀旧篇 3CD》[WAV/分轨][1.6GB]