[pymvpa] Using searchlight for quick spherical ROIs? (Michael Hanke)

Yaroslav Halchenko yoh at psychology.rutgers.edu
Tue Mar 17 20:43:36 UTC 2009


yeah there should be more transparent way ;)

smth like ... (again not tested)

feat_ids = [tds.mapper.getOutId(coord[::-1]) for coord in ROI_coords]
sl = Searchlight(cv, radius=slr, center_ids=feat_ids) #Restrict to
sl_map = sl(tds) # RUN
out = N.zeros(tds.nfeatures)
out[N.array(feat_ids)] = sl[:]
tds.map2Nifti(out).save('bugi.nii.gz')

But I wonder -- what value of slr do you take? imho it might happen that
you will touch some other ROI? or they are distant enough?

-- 
Yaroslav Halchenko
Research Assistant, Psychology Department, Rutgers-Newark
Student  Ph.D. @ CS Dept. NJIT
Office: (973) 353-1412 | FWD: 82823 | Fax: (973) 353-1171
        101 Warren Str, Smith Hall, Rm 4-105, Newark NJ 07102
WWW:     http://www.linkedin.com/in/yarik        



More information about the Pkg-ExpPsy-PyMVPA mailing list