export_volume.m

% Export brain volume images
%
% This script requires BrainDecoderToolbox2 (https://github.com/KamitaniLab/BrainDecoderToolbox2).
%
% Change Log:
%
% 2019-01-15  Shuntaro Aoki  <aoki@atr.jp>
%
%     *  Initial version.
%

clear all;

for i = 1:5
    bdata_file    = sprintf('./data/Subject%0d.mat', i);
    template_file = sprintf('./data/Subject%0d_SpaceTemplate.nii', i);
    output_file   = sprintf('./data/Subject%0d_Func.nii', i);
    
    [dataset, metadata] = load_data(bdata_file);

    % Voxel data
    voxel_data = get_dataset(dataset, metadata, 'VoxelData');

    % Voxel xyz
    voxel_x = get_metadata(metadata, 'voxel_x', 'RemoveNan', true);
    voxel_y = get_metadata(metadata, 'voxel_y', 'RemoveNan', true);
    voxel_z = get_metadata(metadata, 'voxel_z', 'RemoveNan', true);

    % Exporting to .nii image (1st volume)
    export_volumeimage(output_file, voxel_data(1, :), [voxel_x; voxel_y; voxel_z], template_file);
end