Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Error Reading MGH Files #3

Open
smile0227 opened this issue May 6, 2024 · 1 comment
Open

Error Reading MGH Files #3

smile0227 opened this issue May 6, 2024 · 1 comment

Comments

@smile0227
Copy link

smile0227 commented May 6, 2024

Dear Sebenius,
I apologize for any inconvenience caused by reaching out to you. I am currently encountering an issue with retrieving cortical thickness values for each vertex within every brain region after smoothing.
Below is the code I am using:

surf_dir = 'K:/NC11'
features = [('K:/NC11/surf/lh.thickness.fwhm10.fsaverage.mgh','K:/NC11/surf/rh.thickness.fwhm10.fsaverage.mgh')]
parcellation = 'aparc'
vertex_data, regions, features_used = get_vertex_df(surf_dir, features, parcellation)

However, I am encountering the following error:
hemi_data[i + 1] = hemi_data_dict[feature]
ValueError: could not broadcast input array from shape (163842,) into shape (138345,)

It seems that this error arises from the dimensions of .thickness.fwhm10.fsaverage.mgh being 1638421, while lh.aparc.annot has dimensions of 138345*1.

Would you happen to have any insights on how to resolve this issue?
Looking forward to your response.
Best regards,
Yanxi
2024.05.07

@smile0227
Copy link
Author

Dear Sebenius, I apologize for any inconvenience caused by reaching out to you. I am currently encountering an issue with retrieving cortical thickness values for each vertex within every brain region after smoothing. Below is the code I am using:

surf_dir = 'K:/NC11' features = [('K:/NC11/surf/lh.thickness.fwhm10.fsaverage.mgh','K:/NC11/surf/rh.thickness.fwhm10.fsaverage.mgh')] parcellation = 'aparc' vertex_data, regions, features_used = get_vertex_df(surf_dir, features, parcellation)

However, I am encountering the following error: hemi_data[i + 1] = hemi_data_dict[feature] ValueError: could not broadcast input array from shape (163842,) into shape (138345,) It seems that this error arises from the dimensions of _.thickness.fwhm10.fsaverage.mgh being 163842_1, while lh.aparc.annot has dimensions of 138345*1.

Would you happen to have any insights on how to resolve this issue? Looking forward to your response. Best regards, Yanxi 2024.05.07

I think I understand now. I should read lh.aparc.annot in the fsaverage space instead of individual space.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant