想做无忧无虑的咸鱼
2019-01-28 23:33
采纳率: 100%
浏览 2.3k

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条回答 默认 最新

相关推荐 更多相似问题