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

tapas_physio_pca -> no voxel limitation #192

Open
wants to merge 9 commits into
base: development
Choose a base branch
from

Conversation

benoitberanger
Copy link
Contributor

This PR responds to #162

CHANGE : tapas_physio_pca

This function now uses the covariance matrix instead of the matrix itself.
This strategy is inpired by spm_regions function.

CHANGE : tapas_physio_create_noise_rois_regressors

The function have minor modifications due to the PCA rework.
However :

  • The "force coregister & reslice" warning is now a normal log message
  • Logs added to check the number of voxels remaining in the ROI after threshold & erosion.

NEW : function tapas_physio_interactive_noiseROI_mask

This function will load a mask and display it using spm_orthviews, apply a threshold + erosion then display the result using an overlay.
The noiseROI volume (the overlay) remains in memory, no volume is written. Maybe a feature to add ?
A panel is added on spm_orthviews so the threshold and erosion level can be changed and the result visualized dynamicaly.
exemple_interactive_ROI

DELETE : tapas_physio_pca_test

The function have been removed, since the approach have changed (covariance matrix)
This is a discussion the to have with the devs (@mrikasper). I'm opened to keep the test function, even though it needs to be complitly reworked to keep the old appreach and the new.

and added a description in the nifti heater
- this solves the error when 'nVoxels <= nVolumes'
- tapas_physio_pca now complete rewrite
- tapas_physio_create_noise_rois_regressors adapted for the new tapas_physio_pca
- more logs in tapas_physio_create_noise_rois_regressors
- 'different orientation" warning is now a normal log message
variables names correspond to spm_regions function
tapas_physio_create_noise_rois_regressors adapted
warnings -> simple logs
more help / comments
@mrikasper mrikasper self-assigned this Jul 18, 2022
@mrikasper mrikasper added the physio Issues related to PhysIO Toolbox label Jul 18, 2022
@mrikasper mrikasper added this to the PhysIO v9.0 milestone Jul 18, 2022
@benoitberanger
Copy link
Contributor Author

TODO

  • interactive : return the values for latter usage
  • unittest : provide dataset and scripts to check current and proposed PCA implementation
  • when "not enough voxels" still indicate how many before stopping execution

@mrikasper
Copy link
Member

Dear Benoit,

My apologies, but I didn't have enough time to include your PR into the upcoming TAPAS release. I am still working with Johanna Bayer on the summer project to improve PhysIO, but it will take a couple of weeks to overhaul everything. I will let you know once we integrate your PR.

All the best,
Lars

@benoitberanger
Copy link
Contributor Author

benoitberanger commented May 31, 2023

  • interactive : return the values for latter usage fc8a3a2
  • unittest : provide dataset and scripts to check current and proposed PCA implementation
  • when "not enough voxels" still indicate how many before stopping execution eff3868

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
physio Issues related to PhysIO Toolbox
Projects
Development

Successfully merging this pull request may close these issues.

None yet

2 participants