function voxellist = generateVoxelList(center, radius)
%     cx = ROICenter(1);
%     cy = ROICenter(2);
%     cz = ROICenter(3);
    
    voxellist = [];
    
    cx = center(1);
    cy = center(2);
    cz = center(3);

    for z=cz-radius+1:cz+radius-1
        for y=(cy-radius+1):(cy+radius-1)
            for x=(cx-radius+1):(cx+radius-1)
                voxellist = [voxellist [x;y;z]];
            end
        end
    end

end