在這裏推薦一個生成磁共振mask的工具包,Bart工具箱https://mrirecon.github.io/bart/。可以用python調用,代碼如下
import sys
import os
os.environ['TOOLBOX_PATH'] = '/home/txx/program/software/bart-0.4.04'
sys.path.append('/home/txx/program/software/bart-0.4.04/python')
import numpy as np
import h5py
import time
import cfl
from bart import bart
import matplotlib.pyplot as plt
from pylab import rcParams
rcParams['figure.figsize'] = 7, 7
import scipy.io as sio
import skimage
import skimage.measure as meas
mask=bart(1,'poisson -Y 256 -y 2 -Z 256 -z 2 -C 32 -e')
plt.imshow(np.abs(mask[0,:,:]),cmap='gray')
plt.show()
結果圖:
BART:https://github.com/mrirecon/bart-workshop/blob/master/demos/intro/intro.ipynb