I’m using get_clusters_table function from nilearn. I outputs cluster coordinates in xyz. Do you know any python packages/functions that can take this coordinates and output a corresponding atlas label (e.g. Talairach Brodmann areas)?
I saw that one can do it in MRIcro manually, but I would really like to find a more automatizw solution.
An alternative approach that may interest you is to create z-maps (stat images) rather than cluster tables, and then use atlasreader to obtain coordinates + atlas labels.
AFNI has whereami in the GUI and on the command line for this purpose. You can feed it tables of coordinates for a particular template space and for atlas or multiple atlases. See video link below.
you can simply transform coordinates to voxel indices by multiplying the coordinates by the inverse affine, and then index the atlas image with the result