1、保存.mat文件
from scipy.io import savemat
savemat(out_dir / “xxx.mat”, {‘dwi’:’dwi’, ‘t1’:’t1’})
2、保存.h5文件
import h5py
with h5py.File(out_dir / ‘gwh.h5’, ‘w’) as hf:
hf.create_dataset(‘t1’, data=t1)
hf.create_dataset(‘dwi’, data=dwi)
3、访问.mat文件变量
import scipy.io as sio
#读取.mat文件
data = sio.loadmat(‘xxx.mat’)
#访问键为‘variable’的变量
var1 = data[‘variable’]
3、画图函数,方便变量展示
import matplotlib.pyplot as plt
Plt.figure()
Plt.subplot(121)
Plt.imshow(var1, cmap=cmap, vmin=np.min(var1), vmax=np.max(var1))
plt.subplot(122)
Plt.imshow(var2, camp=camp, vmin=np.min(var2), vmax=np.max(var2))
Plt.show()
4、For函数、If函数
for循环
for num in range(slice_num):
print(‘slice_num: ’, slice_num)
if分支
if fname_1.name == fname_2.name:
print(fname_1.name)
1、load(xxx.mat) 函数
2、size()读出数组的维度
3、squeeze()函数,将维度值为1的维度去除
4、reshape()函数
5、画图调试
figure(); plot(x, y, ‘r-’, ‘Linewidth’,2); hold on;
6、disp()
7、h5read(‘xxx路径’, ‘/变量名’)
8、for循环
for i=1:xxx
end
9、imshow(xxx,[0, 5000])
显示多行多列:
figure(1);imshow([output(:, :, 1) output(:, :, 2); output(:, :, 3) output(:, :, 4)], [0, bar]); colormap jet;
10、save(‘filename.mat’, ‘var1’);
11、if语句
if b3_max > b2_max
b2=rand_small_map;
else
b2=rand_map1;
end
12、可以通过列表来拼凑字符串
save_name = [ outdir, ‘mask.mat’];
Save(sava_name, ‘mask’);
13、转换维度
permute(output, [2, 3, 1]);
14、将数组中non置为0
temp = nonzeros(temp);