python处理BrainWeb20数据,20个真实头模型T1加权(181、256、256)、标签图(362、434、362);模拟生成常规图
#/author: guohui time:2020/6/3
#/E/python file
'''读取rawb文件'''
import numpy as np
from PIL import Image
from matplotlib import pyplot as plt
import os
import scipy.ndimage
import cv2
def crop_center(img,cropx,cropy):
y,x = img.shape
startx = x//2-(cropx//2)
starty = y//2-(cropy//2)
return img[starty:starty+cropy,startx:startx+cropx]
path = r'E:\guohui\data_set\brainweb\T1_1m_0_0\脑脊液\phantom_1.0mm_normal_csf.rawb'
img_data = np.fromfile(path, dtype='uint8')
# data_new_shape = img_data.reshape(181, 256, 256) #20个真实头模型T1加权
# data_new_shape = img_data.reshape(362