ITK 实现多张图像转成单个nii.gz或mha文件案例
作者:zhimingf 时间:2023-02-18 20:39:51
主要实现的部分是利用NameGeneratorType读入系列图像,见头文件#include "itkNumericSeriesFileNames.h"。
需要包含的头文件有:
#include "itkImage.h"
#include "itkImageSeriesReader.h"
#include "itkImageFileWriter.h"
#include "itkNumericSeriesFileNames.h"
#include "itkPNGImageIO.h"//转成JPG格式,将PNG替换成JPEG就可以。
int main( int argc, char ** argv )
{
// 需要四个参数,分别是程序起点,第一张图像的编号和最后一张图像的变化,输出文件的名称(包含路径)
if( argc < 4 )
{
std::cerr << "Usage: " << std::endl;
std::cerr << argv[0] << " firstSliceValue lastSliceValue outputImageFile " << std::endl;
return EXIT_FAILURE;
}
//定义读入图像类型,创建对应的reader
typedef unsigned char PixelType;
const unsigned int Dimension = 3;
typedef itk::Image< PixelType, Dimension > ImageType;
typedef itk::ImageSeriesReader< ImageType > ReaderType;
typedef itk::ImageFileWriter< ImageType > WriterType;
ReaderType::Pointer reader = ReaderType::New();
WriterType::Pointer writer = WriterType::New();
//输入参数定义
const unsigned int first = atoi( argv[1] );
const unsigned int last = atoi( argv[2] );
const char * outputFilename = argv[3];//输出的文件名加上对应格式的后缀即可,如mha或nii.gz
//系列图像读入
typedef itk::NumericSeriesFileNames NameGeneratorType;
NameGeneratorType::Pointer nameGenerator = NameGeneratorType::New();
nameGenerator->SetSeriesFormat( "vwe%03d.png" );
nameGenerator->SetStartIndex( first );
nameGenerator->SetEndIndex( last );
nameGenerator->SetIncrementIndex( 1 );//张数的增长间距
//读入图像,写出图像,进行Update
reader->SetImageIO( itk::PNGImageIO::New() );
reader->SetFileNames( nameGenerator->GetFileNames() );
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;
}
return EXIT_SUCCESS;
}
补充知识:将一组png图片转为nii.gz
主要之前使用matlab 对numpy数组存放方式不是很了解.应该是[z,x,y]这样在itksnamp上看就对了
import SimpleITK as sitk
import glob
import numpy as np
from PIL import Image
import cv2
import matplotlib.pyplot as plt # plt 用于显示图片
def save_array_as_nii_volume(data, filename, reference_name = None):
"""
save a numpy array as nifty image
inputs:
data: a numpy array with shape [Depth, Height, Width]
filename: the ouput file name
reference_name: file name of the reference image of which affine and header are used
outputs: None
"""
img = sitk.GetImageFromArray(data)
if(reference_name is not None):
img_ref = sitk.ReadImage(reference_name)
img.CopyInformation(img_ref)
sitk.WriteImage(img, filename)
image_path = './oriCvLab/testCvlab/img/'
image_arr = glob.glob(str(image_path) + str("/*"))
image_arr.sort()
print(image_arr, len(image_arr))
allImg = []
allImg = np.zeros([165, 768,1024], dtype='uint8')
for i in range(len(image_arr)):
single_image_name = image_arr[i]
img_as_img = Image.open(single_image_name)
# img_as_img.show()
img_as_np = np.asarray(img_as_img)
allImg[i, :, :] = img_as_np
# np.transpose(allImg,[2,0,1])
save_array_as_nii_volume(allImg, './testImg.nii.gz')
print(np.shape(allImg))
img = allImg[:, :, 55]
# plt.imshow(img, cmap='gray')
# plt.show()
来源:https://blog.csdn.net/zhimingf/article/details/71404557
标签:ITK,图像,nii.gz,mha
0
投稿
猜你喜欢
javascript禁用Tab键脚本实例
2023-10-13 10:53:58
PassWord输入框代码分享
2024-04-18 09:34:36
mysql启动时出现ERROR 2003 (HY000)问题的解决方法
2024-01-22 02:53:13
MySQL-group-replication 配置步骤(推荐)
2024-01-17 20:51:36
实例详解JavaScript中setTimeout函数的执行顺序
2024-04-22 13:25:09
C#命名空间System.ComponentModel属性方法汇总
2024-06-05 09:24:02
js处理括弧配对替换的方法
2008-01-16 13:48:00
Golang实现简易的rpc调用
2024-04-25 15:07:55
MySQL过滤重复数据的两种方法示例
2024-01-16 12:54:27
TensorFlow平台下Python实现神经网络
2023-10-11 22:31:31
Matplotlib绘制雷达图和三维图的示例代码
2022-09-04 03:31:33
python实现邮件自动发送
2023-06-10 16:24:44
MySQL8.0中的窗口函数的示例代码
2024-01-14 12:30:26
对python创建及引用动态变量名的示例讲解
2023-02-03 07:10:49
匹配 IP 地址与域名的正则表达式
2023-06-17 05:55:48
基于Python预测一下世界杯最后赢家
2021-01-30 16:16:48
利用Python多线程实现图片下载器
2022-02-12 06:07:41
Python利用前序和中序遍历结果重建二叉树的方法
2023-03-01 14:56:35
python列表推导式操作解析
2022-08-26 04:07:39
Python defaultdict方法使用分析
2023-01-02 18:05:35