代码参考:将一个MRI(nii格式)数据进行切片
想要批量抽取MRI数据的几个切片,借鉴了上述代码的Python版本,原始代码抽取了全部切片。本次实验主要是抽取sMRI的横断面(axial)方向的其中30张切片,数量根据自己需求可自行设置。
import scipy, numpy, shutil, os, nibabel
import sys, getopt
import imageio
def niito2D(filepath):
inputfiles = os.listdir(filepath) #遍历文件夹数据
outputfile = 'result\\' #输出文件夹
print('Input file is ', inputfiles)
print('Output folder is ', outputfile)
for inputfile in inputfiles:
image_array = nibabel.load(filepath+inputfile).get_data() #数据读取
print(len(image_array.shape))
# set destination folder
if not os.path.exists(outputfile):
os.makedirs(outputfile) #不存在输出文件夹则新建
print("Created ouput directory: " + outputfile)
print('Reading NIfTI file...')
total_slices = 30 #总切片数
slice_counter = 70 #从第几个切片开始
# iterate through slices
for current_slice in ra
346



&spm=1001.2101.3001.12003&articleId=114547207&d=1&t=3&u=33ad49c426724f1184b35ac80b1cffbe)

被折叠的 条评论
为什么被折叠?



