拉文克劳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条)

报告相同问题?

悬赏问题

  • ¥15 #MATLAB仿真#车辆换道路径规划
  • ¥15 java 操作 elasticsearch 8.1 实现 索引的重建
  • ¥15 数据可视化Python
  • ¥15 要给毕业设计添加扫码登录的功能!!有偿
  • ¥15 kafka 分区副本增加会导致消息丢失或者不可用吗?
  • ¥15 微信公众号自制会员卡没有收款渠道啊
  • ¥100 Jenkins自动化部署—悬赏100元
  • ¥15 关于#python#的问题:求帮写python代码
  • ¥20 MATLAB画图图形出现上下震荡的线条
  • ¥15 关于#windows#的问题:怎么用WIN 11系统的电脑 克隆WIN NT3.51-4.0系统的硬盘