计算mask的体素数量
import numpy as np
import nibabel as nib # 用于处理神经影像数据的库
# 从文件中加载mask图像
mask_image = nib.load(r'E:\mask.nii.gz')
# 获取图像数据
mask_data = mask_image.get_fdata()
# 计算非零像素的数量，即白质骨架的体素总数
voxel_count = np.count_nonzero(mask_data)
print(f"白质骨架的体素总数为：{voxel_count}")
65611