# Overlay another hopeful beauty (functional overlay)
        brain.add_overlay(surf_data, name="overlap", sign=sw0.sign)

        # Update colorbar
        tmp = brain.overlays["overlap"]
        lut = tmp.pos_bar.lut.table.to_array()
        lut[:, 0:3] = [230, 171, 2]
        tmp.pos_bar.lut.table = lut

        # Refresh
        brain.show_view("lat")
        brain.hide_colorbar()

    # Save the beauts
    brain.save_imageset("%s_%s" % (sw0.outprefix, hemi),
                        sw0.views,
                        'jpg',
                        colorbar=None)

    # End a great journey, till another life
    brain.close()

sw0.cropify()
sw0.montage("box")
sw0.montage("stick")

####
## Setup
#
## 3dcalc -a /home2/data/Projects/CWAS/nki/cwas/short/compcor_kvoxs_smoothed_to_kvoxs_smoothed/iq_age+sex+meanFD.mdmr/cluster_correct_v05_c05/clust_logp_FSIQ.nii.gz -b /home2/data/Projects/CWAS/nki/cwas/medium/compcor_kvoxs_smoothed_to_kvoxs_smoothed/iq_age+sex+meanFD.mdmr/cluster_correct_v05_c05/clust_logp_FSIQ.nii.gz -expr 'step(a)+(2*step(b))' -prefix /home2/data/Projects/CWAS/figures/fig_03/C_overlap.nii.gz
#
#basedir = "/home2/data/Projects/CWAS"
    if (sum(abs(surf_data)) > 0):
        # Overlay another hopeful beauty (functional overlay)
        brain.add_overlay(surf_data, name="overlap", sign=sw0.sign)
    
        # Update colorbar
        tmp = brain.overlays["overlap"]
        lut = tmp.pos_bar.lut.table.to_array()
        lut[:,0:3] = [230,171,2]
        tmp.pos_bar.lut.table = lut
    
        # Refresh
        brain.show_view("lat")
        brain.hide_colorbar()
    
    # Save the beauts
    brain.save_imageset("%s_%s" % (sw0.outprefix, hemi), sw0.views, 
                        'jpg', colorbar=None)
    
    # End a great journey, till another life
    brain.close()

sw0.cropify()
sw0.montage("box")
sw0.montage("stick")



####
## Setup
#
## 3dcalc -a /home2/data/Projects/CWAS/nki/cwas/short/compcor_kvoxs_smoothed_to_kvoxs_smoothed/iq_age+sex+meanFD.mdmr/cluster_correct_v05_c05/clust_logp_FSIQ.nii.gz -b /home2/data/Projects/CWAS/nki/cwas/medium/compcor_kvoxs_smoothed_to_kvoxs_smoothed/iq_age+sex+meanFD.mdmr/cluster_correct_v05_c05/clust_logp_FSIQ.nii.gz -expr 'step(a)+(2*step(b))' -prefix /home2/data/Projects/CWAS/figures/fig_03/C_overlap.nii.gz
#