利用mask将2D序列转为4D数据

import os
import numpy as np
import nibabel as nib
from nilearn.input_data import NiftiMasker

os.environ["OMP_NUM_THREADS"]="2"
os.environ["CUDA_VISIBLE_DEVICES"]="0"

source='/path/'
X=np.load(source+'data.npy')
print("X.shape:",X.shape)
Y=X[:,:]#切片转化指定的中间序列
print("Y.shape:",Y.shape)


#%%%
mask_img =nib.load(source + 'mask_path.nii')
print(mask_img.shape)

masker = NiftiMasker(mask_img=mask_img, 
                     standardize=True,
                     detrend=1,
                     smoothing_fwhm=6,
                     memory_level=2,t_r=0.72)
masker = masker.fit()

components_img=masker.inverse_transform(Y)
print(components_img.shape)
nib.save(components_img,'1.nifti.nii.gz')#保存转换后的数据