Does anyone know of a good solution for extracting ROIs from cifti data in python? Specifically:
I have 7T HCP functional data (fsLR 59k, surface+volume)
I have the MMP parcellation (fsLR 59k, surface) it is not released in this resolution but someone projected it for now.
Say I want to apply a mask of a few different MMP parcels to extract an ROI time series for all of the grayordinates there from my functional data (without averaging over this roi). If I get_fdata() and create an index from the MMP data to apply to the functional data it does not work because they are different lengths and do not line up.