nii.gz格式四维转三维(给爷写笑了)

把这里当一下垃圾篓(别管我啦 发疯.jpg

import nibabel as nib

def convert_4d_to_3d(input_file, output_file, drop_axis=2):
    # 加载四维的nii.gz文件
    img = nib.load(input_file)

    # 获取四维数据
    data = img.get_fdata()

    # 丢弃指定的维度,这里假设drop_axis为3表示丢弃第四个维度(即时间维度)
    # 如果drop_axis为0,则丢弃第一个维度,以此类推
    new_data = data.squeeze(axis=drop_axis)

    # 创建一个新的affine矩阵,因为丢弃维度后空间坐标可能会改变
    # 这里我们简单地复制原始affine,但根据具体情况可能需要调整
    new_affine = img.affine.copy()

    # 创建新的header,并更新维度信息
    header = img.header.copy()
    header.set_data_shape(new_data.shape)

    # 使用新的数据和header创建一个新的nibabel图像对象
    new_img = nib.Nifti1Image(new_data, new_affine, header)

    # 保存为新的三维nii.gz文件
    nib.save(new_img, output_file)

input_file = 'E:\\zxxzl\\srtp02\\four_jiqixuexi\\pngtransform\\pngto3d_test\\output\\1_4d.nii.gz'  # 替换为您的四维nii.gz文件的路径
output_file = 'E:\\zxxzl\\srtp02\\four_jiqixuexi\\pngtransform\\pngto3d_test\\output\\1_4d_3d.nii.gz'  # 您想要保存的三维nii.gz文件的路径
convert_4d_to_3d(input_file, output_file)

output_path = 'E:\\zxxzl\\srtp02\\four_jiqixuexi\\pngtransform\\pngto3d_test\\output\\1_4d_3d.nii.gz'
# 加载保存的NIfTI文件以验证维度
loaded_nii = nib.load(output_path)
print("Loaded image dimensions:", loaded_nii.shape)
print("Loaded image affine:", loaded_nii.affine)

批量四维转三维


import os
import numpy as np
from PIL import Image
import nibabel as nib

# 读取PNG文件夹路径
png_folder = "E:/zxxzl/srtp02/four_jiqixuexi/pngtransform/nii4dto3d/4ddeck"
# 创建保存NIfTI文件的目标文件夹
output_folder = "E:/zxxzl/srtp02/four_jiqixuexi/pngtransform/nii4dto3d/output"
os.makedirs(output_folder, exist_ok=True)

def convert_4d_to_3d(input_file, output_file, drop_axis=2):
    # 加载四维的nii.gz文件
    img = nib.load(input_file)

    # 获取四维数据
    data = img.get_fdata()

    # 丢弃指定的维度,这里假设drop_axis为3表示丢弃第四个维度(即时间维度)
    # 如果drop_axis为0,则丢弃第一个维度,以此类推
    new_data = data.squeeze(axis=drop_axis)

    # 创建一个新的affine矩阵,因为丢弃维度后空间坐标可能会改变
    # 这里我们简单地复制原始affine,但根据具体情况可能需要调整
    new_affine = img.affine.copy()

    # 创建新的header,并更新维度信息
    header = img.header.copy()
    header.set_data_shape(new_data.shape)

    # 使用新的数据和header创建一个新的nibabel图像对象
    new_img = nib.Nifti1Image(new_data, new_affine, header)

    # 保存为新的三维nii.gz文件
    nib.save(new_img, output_file)

for filename in os.listdir(png_folder):
    if filename.endswith(".nii.gz"):
        input_file = os.path.join(png_folder, filename)
        # 构建目标文件的完整路径和名称
        nii_filename = os.path.splitext(filename)[0] + '.gz'
        output_file = os.path.join(output_folder, nii_filename)
        convert_4d_to_3d(input_file, output_file)


output_path1 = "E:/zxxzl/srtp02/four_jiqixuexi/pngtransform/nii4dto3d/output/1_mask.nii.gz"
loaded_nii = nib.load(output_path1)
print("Loaded image dimensions:", loaded_nii.shape)
print("Loaded image affine:", loaded_nii.affine)


 

  • 5
    点赞
  • 9
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
对于将nii.gz格式三维图像换为dicom格式,可以使用以下Python代码: ```python import SimpleITK as sitk # 读取nii.gz文件 image = sitk.ReadImage('example.nii.gz') # 将nii.gz文件换为dicom格式 sitk.WriteImage(image, 'example.dcm') ``` 上述代码使用了SimpleITK库,可以从nii.gz文件中读取三维图像数据,并将其换为dicom格式。需要注意的是,dicom格式需要指定文件名后缀为".dcm"。 对于将dicom三维图像换为三维重建VTK的输入格式,可以使用以下Python代码: ```python import SimpleITK as sitk import vtk # 读取dicom文件 reader = vtk.vtkDICOMImageReader() reader.SetDirectoryName('example.dcm') reader.Update() # 将dicom三维图像换为VTK格式 sitk_img = sitk.ReadImage('example.dcm') vtk_img = sitk.Cast(sitk_img, sitk.sitkUInt8) vtk_img = sitk.GetArrayFromImage(vtk_img) vtk_img = vtk.vtkImageData() vtk_img.SetDimensions(sitk_img.GetSize()[::-1]) vtk_img.SetSpacing(sitk_img.GetSpacing()[::-1]) vtk_img.SetOrigin(sitk_img.GetOrigin()[::-1]) vtk_img.GetPointData().SetScalars(vtk.vtkDataArray.CreateDataArray(vtk.VTK_UNSIGNED_CHAR, array=vtk_img.ravel(), numComponents=1)) # 保存VTK格式三维图像数据 writer = vtk.vtkXMLImageDataWriter() writer.SetFileName('example.vti') writer.SetInputData(vtk_img) writer.Write() ``` 上述代码使用了VTK库和SimpleITK库,可以从dicom文件中读取三维图像数据,并将其换为VTK格式。需要注意的是,VTK格式需要指定文件名后缀为".vti"。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值