拉文克劳10086 2019-01-28 23:33 采纳率: 100%
浏览 2403
已采纳

SimpleITK 做 三维矩阵插值

为什么如下代码里date1里面没有数据,注date是由一个三维矩阵的,
多谢大佬回答。感谢

import SimpleITK as sitk
import pandas as pd

file_path = "E:\data"  #记录: 此处不可以放叠加的文件夹,会扫描不到

#显示所有列
pd.set_option('display.max_columns', None)
#显示所有行
pd.set_option('display.max_rows', None)
reader = sitk.ImageSeriesReader()
series_IDs = sitk.ImageSeriesReader.GetGDCMSeriesIDs(file_path)
# 查看该文件夹下的序列数量
nb_series = len(series_IDs)
print(nb_series)
# 读取文件数据
series_file_names = sitk.ImageSeriesReader.GetGDCMSeriesFileNames(file_path, series_IDs[0])
reader.SetFileNames(series_file_names)
image3D = reader.Execute()
date = sitk.GetArrayFromImage(image3D)
#print(date)
#显示图像尺寸
print(image3D.GetSize())

# 做一个0.5 * 0.5 *0.5 的插值
resample = sitk.ResampleImageFilter()
resample.SetOutputDirection(image3D.GetDirection())
resample.SetOutputOrigin(image3D.GetOrigin())
newspacing =  [1,1,1]
resample.SetOutputSpacing(newspacing)
newimage = resample.Execute(image3D)
date1 = sitk.GetArrayFromImage(newimage)
print(date1)


  • 写回答

2条回答 默认 最新

  • 拉文克劳10086 2019-02-01 11:32
    关注
    本回答被题主选为最佳回答 , 对您是否有帮助呢?
    评论
查看更多回答(1条)

报告相同问题?

悬赏问题

  • ¥20 有关区间dp的问题求解
  • ¥15 多电路系统共用电源的串扰问题
  • ¥15 slam rangenet++配置
  • ¥15 有没有研究水声通信方面的帮我改俩matlab代码
  • ¥15 对于相关问题的求解与代码
  • ¥15 ubuntu子系统密码忘记
  • ¥15 信号傅里叶变换在matlab上遇到的小问题请求帮助
  • ¥15 保护模式-系统加载-段寄存器
  • ¥15 电脑桌面设定一个区域禁止鼠标操作
  • ¥15 求NPF226060磁芯的详细资料