Hi @Wang,
Thanks for reaching out.
If you want to use the osp_plotVoxelOverlap
function you have to run the following steps:
MRSCont = OspreyJob('PathToYourJob.m');
MRSCont = OspreyLoad(MRSCont);
MRSCont = OspreyProcess(MRSCont);
MRSCont = OspreyCoreg(MRSCont);
MRSCont = OspreySeg(MRSCont);
or the equivalent GUI buttons.
In the jobfile you should include the MRS data and T1 images in subject space. And you have to make sure that there are no SPM segmentation results from a prior SPM/Osprey analysis in the output folder.
Afterward, you should be able to run the osp_plotVoxelOverlap
function. See more details about options for this function in this post (Osprey Voxel Overlap Heatmap Problem).
Best,
Helge