from neurosynth.base.dataset import Dataset
from neurosynth.analysis import network

dataset = Dataset.load('dataset.pkl')

roi = 'my_roi_mask.nii.gz'

network.coactivation(dataset, roi, threshold=0.1, outroot='coactivation_from_image')

network.coactivation(dataset, [[0, 20, 28]], threshold=0.1, outroot='coactivation_from_coords', r=10)