sng-newy / eigenstrapping Goto Github PK

AttributeError                            Traceback (most recent call last)
Input In [82], in <module>
      1 from eigenstrapping import SurfaceEigenstrapping
----> 2 t = SurfaceEigenstrapping(surface='../tpl-avg_space-canonical_den-1mm_label-hipp_midthickness.surf.gii'
      3                               ,data='../thickness_midsurf_L_mean.shape.gii', num_modes=200)

File /export03/data/opt/venv/lib/python3.8/site-packages/eigenstrapping/base.py:329, in SurfaceEigenstrapping.__init__(self, data, surface, evals, emodes, num_modes, save_surface, seed, decomp_method, medial, randomize, resample, n_jobs, use_cholmod, permute, add_res, save_rotations, parcellation, normalize)
    326     self.surface = load_surface(self.surface_file)
    328 self.data = copy.deepcopy(data)
--> 329 self.n_vertices = self.data.shape[0]
    331 if emodes is not None and evals is not None:
    332     self.evals = evals[:num_modes]

AttributeError: 'str' object has no attribute 'shape'

src = mne.setup_source_space(subject, spacing='oct6', subjects_dir=subjects_dir, add_dist=False, n_jobs=10)
vox16_to_vox8 = {ivox:i for i, ivox in enumerate(src[0]['vertno'])}
used_tris_dsmp = [[vox16_to_vox8[tri[0]], vox16_to_vox8[tri[1]], vox16_to_vox8[tri[2]]] for tri in src[0]['use_tris']]
surface=(src[0]['rr'][src[0]['inuse'].astype(bool)].T, used_tris_dsmp)
tria = TriaMesh(surface[0], surface[1])
evals, emodes = calc_eig(tria, 1000, False)

eigen = SurfaceEigenstrapping(
    emodes=emodes,
    evals=evals,
    data=data[:4098],
    num_modes=1000,
    save_rotations=True,
    normalize=False
)