addpath(genpath('.../spm12/'))% unzip
input_file_name ='/home/data/IXI012-HH-1211-MRA_001.nii.gz';gunzip(input_file_name)% path
input_file_name ='/home/data/IXI012-HH-1211-MRA_001.nii';
output_file_name ='/home/data/IXI012-HH-1211-MRA_001_skeleton.nii';% read nii
img =spm_vol(input_file_name);% read nii (not support nii.gz)
testvol =spm_read_vols(img);% get mat data
testvol =logical(testvol);% logical type for analysis% process data (extract skeleton of the vessel segmentation for example)
skel =Skeleton3D(testvol);% extract skeleton (process mat data)% save nii with changed filename and data
skel =uint8(skel);% data format (uint8) for saving nii
img.fname = output_file_name;%更改输出文件名,否则会覆盖原文件spm_write_vol(img, skel);% write image