atlas = ft_read_atlas('/Users/anaskhan/Desktop/Apps/fieldtrip/template/atlas/yeo/Yeo2011_7Networks_MNI152_FreeSurferConformed1mm_LiberalMask_colin27.nii');
atlas.coordsys = 'mni';

colorMap = containers.Map;
colorMap('tissue 1') = [120, 18, 134] / 255;  % Visual
colorMap('tissue 2') = [70, 130, 180] / 255;  % Somatomotor
colorMap('tissue 3') = [0, 118, 14] / 255;    % Dorsal Attention
colorMap('tissue 4') = [196, 58, 250] / 255;  % Ventral Attention/Salience
colorMap('tissue 5') = [220, 248, 164] / 255; % Limbic
colorMap('tissue 6') = [230, 148, 34] / 255;  % Frontoparietal (FPCN)
colorMap('tissue 7') = [205, 62, 78] / 255;   % Default Mode

for lab = 1:numel(atlas.tissuelabel)
    
    cfg = [];
    cfg.atlas      = atlas;
    cfg.roi        = atlas.tissuelabel{lab};
  
    mask = ft_volumelookup(cfg, atlas);

    seg = keepfields(atlas, {'dim','unit','coordsys','transform'});
    seg.brain = mask;

    cfg = [];
    cfg.method      = 'iso2mesh';   
    cfg.radbound    = 2;            
    cfg.maxsurf     = 0;
    cfg.tissue      = 'brain';
    cfg.numvertices = 100000;
    cfg.smooth      = 3;
    cfg.spmversion  = 'spm12';
    roi_mesh = ft_prepare_mesh(cfg, seg);

    ft_plot_mesh(roi_mesh, ...
        'facecolor', colorMap(atlas.tissuelabel{lab}), ...
        'edgecolor', 'none', ...
        'facealpha', 1.0);
    
end