23626
Comment:
|
28790
|
Deletions are marked like this. | Additions are marked like this. |
Line 6: | Line 6: |
== Current topics == ==== Source modeling ==== * Implementation of a new unified minimum norm/beamformer framework (almost finished) ==== Documentation ==== * Cross-validation with MNE and FieldTrip (work in progress): <<BR>>http://martinos.org/mne/dev/auto_tutorials/plot_brainstorm_auditory.html * New tutorial: MEG resting-state / OMEGA database ==== Functional connectivity ==== * Significance thresholding of the connectivity matrices (not started) ==== Computation ==== * Removing the dependence to the Java interface to run in headless mode (not started) * Call MNE-Python and FieldTrip functions from Brainstorm GUI (work in progress) <<BR>><<BR>> |
|
Line 25: | Line 7: |
* RAW file viewer: * '''Faster display''': Downsample time series before plotting |
* RAW file viewer:<<BR>> * Bottom bar: display extend events as segments instead of dots (to see their extension) * Downsample before filtering? (attention to the filter design) * Add parameter to make the visual downsampling more or less aggressive |
Line 32: | Line 16: |
* EEG: Laplacian montage (see doc sent by Jeremy) * Bad trials: When changing the status of bad to good: remove the bad segments as well, otherwise it is not processed by processes like the PSD. |
|
Line 35: | Line 21: |
* Overlay multiple conditions | |
Line 37: | Line 22: |
* Add a proper amplitude scale that gets updated when shift+scroll, to compare figures * Create heat maps: Maybe with matlab function heatmap? |
|
Line 39: | Line 26: |
* Start Brainstorm without Java (-nodesktop) * Option to increase the size of all the elements of the interface (for high res laptop screens) |
* Add a warning when computing a forward model with > 100000 sources (check selection) * Multiple screens: Add option to set on which display the main bst window should be * File filters: Add boolean logic: https://github.com/brainstorm-tools/brainstorm3/issues/68 * Snapshot: Save as image / all figures (similar to Movie/all figure) |
Line 46: | Line 35: |
* Set parula (or others) as the default, not jet: https://bids.github.io/colormap/ * https://jakevdp.github.io/blog/2014/10/16/how-bad-is-your-colormap/ * https://www.youtube.com/watch?v=xAoljeRJ3lU |
|
Line 48: | Line 40: |
* Removing all the CTRL and SHIFT in the keyboard shortcuts | |
Line 58: | Line 49: |
* Connectivity on unconstrained sources: how to group the three orientations? * Display as straight lines |
* Connectivity on unconstrained sources: "Default signal extraction for volume grids should be the time series of the first principal component of the triplet signals after each has been zero-meaned" (SB) * Display of connectivity graphs: * Display as straight lines * Recode 2D graphs * 3D display with anatomical constrains * Display using real position of EEG electrodes |
Line 66: | Line 61: |
* Matrix view of NxN graphs: Add legend of the elements along X and Y axis | |
Line 90: | Line 86: |
* PAC:<<BR>> | * PAC: |
Line 110: | Line 106: |
* Explore Matlab-Python bridge:<<BR>> https://mail.python.org/pipermail/neuroimaging/2016-June/001001.html | * Use CUDA for speeding up some operations (filtering, wavelets, etc) |
Line 112: | Line 108: |
* Add MNE-Python functions: * scikit-learn classifiers: S Marti / G Dehaene * Implement data exchange with MNE-Python: write FIF files from Brainstorm and/or pass python objects in memory instead of FIF files * SSS/tSSS cleaning * Reproduce other tutorials / examples * Change the graphic renderer from Matlab |
|
Line 113: | Line 115: |
* ft_sourceanalysis: * Check noise covariance * Check all the options of all the methods * Single trial reconstructions + noise covariance? * Filters?? http://www.fieldtriptoolbox.org/example/common_filters_in_beamforming * Beamformers: Save ftSource.avg.mom <<BR>>http://www.fieldtriptoolbox.org/workshop/meg-uk-2015/fieldtrip-beamformer-demo * http://www.natmeg.se/ft_beamformer/beamformer.html * http://www.fieldtriptoolbox.org/tutorial/beamformingextended * Baseline? Two inputs? * ft_prepare_sourcemodel: Compute MNI transformation (linear and non-linear) => Peter * Freqanalysis: ITC |
|
Line 117: | Line 130: |
* ft_prepare_sourcemodel: * http://www.fieldtriptoolbox.org/tutorial/minimumnormestimate * Compute MNI transformation (linear and non-linear) => Peter * DICS beamformer: * http://www.fieldtriptoolbox.org/workshop/meg-uk-2015/fieldtrip-beamformer-demo * http://www.natmeg.se/ft_beamformer/beamformer.html * http://www.fieldtriptoolbox.org/tutorial/beamformingextended |
|
Line 125: | Line 131: |
* Use CUDA for speeding up some operations (filtering, wavelets, etc) | * Decoding/Classifiers: Faster algorithms (MNE-Python?) |
Line 127: | Line 133: |
* Bug: After "convert to continuous", the time of the following processes should change | |
Line 128: | Line 135: |
* ITC: Inter-trial coherence (see MNE reports for group tutorial) * PLS: * https://www.rotman-baycrest.on.ca/index.php?section=84 * meg-pls dot weebly dot com * Krishnan 2011: http://www.ncbi.nlm.nih.gov/pubmed/20656037 * Cheung 2015: http://www.sciencedirect.com/science/article/pii/S1053811915007648 * McIntosh 2012: http://www.ncbi.nlm.nih.gov/pubmed/22804773 * McIntosh 2004: http://www.ncbi.nlm.nih.gov/pubmed/15501095 |
* Events: Allow selection from a drop-down list (similar to option "channelname" in panel_process_selection) * ITC: Inter-trial coherence (see MNE reports for group tutorial)<<BR>>http://www.sciencedirect.com/science/article/pii/S1053811916304232 |
Line 138: | Line 139: |
* Why doesn't the ICA process converge when using 25 components in the EEG tutorial? * Add an option to resample the signals before computing the ICA decomposition * Add a stand-alone tutorial |
|
Line 141: | Line 145: |
* Use faster methods (MNE-Python?) | |
Line 146: | Line 151: |
* Understand why EEG/Epilepsy tutorial data crashes if we don't limit the number of components | |
Line 150: | Line 154: |
* ICA+machine learning: https://www.ncbi.nlm.nih.gov/pubmed/28497769 * Automated artifact rejection: https://arxiv.org/abs/1612.08194 * Save IC time series in database * Use EYE-EEG: EEGLAB toolbox for eye-tracker guided ICA (Olaf Dimigen): http://www2.hu-berlin.de/eyetracking-eeg/ * Other EEGLAB functions: * Step function detection: https://github.com/lucklab/erplab/wiki/Artifact-Detection:-Tutorial |
|
Line 156: | Line 166: |
* SSS/tSSS: Get implementation from MNE * Bandpass: Use new filters in all the functions using a bandpass ('bst-hfilter' instead of 'bst-fft-fir') * process_evt_detect_threshold * bst_connectivity * bst_pac |
* Bandpass: * Use new filters in all the functions using a bandpass ('bst-hfilter' instead of 'bst-fft-fir'): process_evt_detect_threshold * Weird bug: Filter(import) != Import(Filter) in the HCP tutorial... to investigate * Bandpass |
Line 174: | Line 183: |
* Multi-tapers | * Optimization: bst_timefreq (around l.136), remove evoked in source space: Average should be computed in sensor space instead of source space (requested by Dimitrios) |
Line 188: | Line 197: |
* Allow baseline normalization of files computed with time bands | |
Line 199: | Line 209: |
* Other processes: * Moving average * Max * Median |
|
Line 211: | Line 217: |
* Faster DB searches (for Emily) | |
Line 213: | Line 220: |
* Navigator: Use F2/F3 to explore the entire database (right now F2 works in a weird way) | |
Line 224: | Line 230: |
* protocol.mat can be too big: do not store the results links in it (and recreate when loading) | * Auto-save: * protocol.mat can be too big: do not store the results links in it (and recreate when loading)- http://neuroimage.usc.edu/forums/t/abnormally-slow-behavior/2065/10 * Improve auto-save: add tracking file next to protocol.mat, do not save all the time, only when closing app, and reload protocol at stratup if tracking file is still there |
Line 237: | Line 246: |
* Mixed head models: * Create scout form maximum doesn't work (menu Sources > Max value) |
* Project sources: Very poor algorithm to project sub-cortical regions and cerebellum (algorithm to fit surfaces should be imrpoved) * Menu head model > Copy to other conditions/subjects (check if applicable first) * Menu Sources > Maximum value: Doesn't work with volume or mixed head models * Mixed head models: <<BR>> |
Line 281: | Line 292: |
* Add eyes models to attract eye activity | |
Line 283: | Line 295: |
* FreeSurfer: * Do not save "mid" surface unless specified explicitely with the process version * Add cerebellum to default model generated with "Import FS anatomy" * Import volume atlases in MNI space in subjects (Jeremy) * Parcelating volume grids: scikit-learn.cluster.Ward * Edit fiducials: Replace 6 text boxes with 1 for easy copy-paste (see fiducials.m) * BrainVISA: Add support for MarsAtlas (Guillaume A) |
* BrainVISA: Add support for MarsAtlas * MNI transformation: Use non-linear MNI transformation y_... * Registration:<<BR>> * Getting electrode positions from 3D scanners: https://sccn.ucsd.edu/wiki/Get_chanlocs * GARDEL: http://meg.univ-amu.fr/wiki/GARDEL:presentation * Use the same registration for multiple recording sessions that have already re-registered previously (eg. with MaxFilter) * When linking multiple EEG recordings including 3D positions, do the registration only once and copy it to all the runs * Check templates: MNI transformation and volume atlases for ICBM152 vs Colin27 (loading the AAL atlas as surface or voume scouts do not align well on the ICBM152) * Compute non-linear MNI registration instead of linear * Select and remove bad digitized head points before automatic coregistration * Load the MNE -transf.fif: http://neuroimage.usc.edu/forums/showthread.php?2830 * MRI Viewer: * Add display of "world coordinates" when "vox2ras" is available * Pan in zoomed view (shift + click + move?) * Zoom in/out with mouse (shift + scroll?) * Ruler tool to measure distances * Add keyboard shortcuts to scroll in the three orientations (same in MRI 3D) * Display scouts as additional volumes * Render surface envelope in the MRI as a thin line instead of the full interpolation matrix * Edit fiducials: Replace 6 text boxes with 1 for easy copy-paste (see fiducials.m) * Optimize computation interpolation MRI-surface (tess_tri_interp) => spm_mesh_to_grid |
Line 292: | Line 317: |
* Add new labels to all BrainSuite anatomy templates | |
Line 294: | Line 320: |
* Warping: Scale option has to be fixed, it is currently very unstable | |
Line 298: | Line 323: |
* Sort scouts by region in process options | |
Line 299: | Line 325: |
* Sort scouts by region in process options * Menu head model > Copy to other conditions/subjects (check if applicable first) * Generate mixed density surfaces * Optimize computation interpolation MRI-surface (tess_tri_interp) * Render surface envelope in the MRI as a thin line instead of the full interpolation matrix |
* Project from one hemisphere to the other using registered spheres/squares (http://neuroimage.usc.edu/forums/t/how-to-create-mirror-roi-in-the-other-hemisphere/5910/8) * Parcellating volume grids: scikit-learn.cluster.Ward |
Line 306: | Line 329: |
* Add eyes models to attract eye activity * tess_mrimask: Not robust enough to major holes in the brain |
* Surface>Volume interpolation: Use spm_mesh_to_grid |
Line 309: | Line 331: |
* MRI Viewer: Display scouts as additional volumes | |
Line 312: | Line 333: |
* Contact positions: Import / set / detect * Import/export electrodes positions in MNI/SCS/MRI coordinates * Co-register MRI and CT for electrodes marking in the MRI Viewer<<BR>>=> Mutual information: BrainVISA > Outils > Recalage rigide * Automatic segmentation of CT scans: See project of Rodrigo Paz (Marseille/Timone) * Import from FreeSurfer/FreeView: Problem of difference between RAS and TkRegRAS: http://neuroimage.usc.edu/forums/showthread.php?1958-SEEG-electrodes-and-subject-s-anatomy-are-not-alligned |
* Electrode models: Better way for managing/updating/adding electrode models * SEEG: Project contacts to SEEG / Compare with IntrAnat * Contact positions: Import / set / detect<<BR>> * New option: Align on none|inner|cortex to replace ECOG-mid |
Line 320: | Line 340: |
* ECOG: Default names of contacts? * Add menu to import implantation channel file in imported recordings * Automatic segmentation of CT: * GARDEL: http://meg.univ-amu.fr/wiki/GARDEL:presentation * Arnulfo: https://bmcbioinformatics.biomedcentral.com/articles/10.1186/s12859-015-0511-6 * MAP07 / SPM: https://www.epi.ch/_files/Artikel_Epileptologie/Huppertz_2_13.pdf * ECOG: * Project and display contacts on cortex surface should consider the rigidity of the grids: Contacts cannot rotate, and distance between contacts should remain constant across runs * Method for contacts projection: https://pdfs.semanticscholar.org/f10d/6b899d851f3c4b115404298d7b997cf1d5ab.pdf * ECOG: Brain shift: When creating contact positions on a post-implantation image, the brain shift should be taken into account for creating images of the ECOG contacts on the pre-op brain => iELVis (http://ielvis.pbworks.com/w/page/116347253/FrontPage) |
|
Line 321: | Line 351: |
* Display SEEG+ECOG contacts at the same time * Better way to represent ECOG strips that are not in contact with the skull (new type?) * SEEG length: count from the most superficial contact * SEEG default view: transparent cortex instead of MRI viewer (Marcel) * ECOG: Project on cortex instead of inner skull, just for display (Marcel) * Color each strip/grid differently |
|
Line 328: | Line 352: |
* ECOG: Topography display without the 3D coordinates (Marcel) * 2DLayout: Should not try to use the actual position of the sensors, just the grids/strips. |
|
Line 331: | Line 353: |
* Display curved SEEG electrodes * Display wires for ECOG? * Re-referencing: * Create new average reference montages with a specific list of channels, with the possibility to edit the order of the channels (for Jeremy) * Closest white reference (Arnulfo) * Detection CEEP stim artifacts: Use ImaGIN code ImaGIN_StimDetect |
|
Line 332: | Line 360: |
* Prepare tutorial with public dataset (Marcel) | |
Line 341: | Line 368: |
* PLS: Partial Least Squares | |
Line 348: | Line 374: |
* Output .nii have incorrect sform/qform when using the options to downsample the volume of cut the empty slices * DICOM converter: * Add dcm2nii (MRICron) * Add MRIConvert * FieldTrip: Import/Export time-frequency: * Export: http://neuroimage.usc.edu/forums/t/export-time-frequency-to-fieldtrip/1968 * Import: http://neuroimage.usc.edu/forums/t/import-time-frequency-data-from-fieldtrip/2644 |
|
Line 356: | Line 389: |
* Export to EDF | |
Line 358: | Line 390: |
* Nihon Koden * Nicolet EEG: Get from from MNE |
|
Line 362: | Line 392: |
* XLTEK: https://github.com/danielmhanover/OpenXLT * Persyst .lay: https://github.com/ieeg-portal/Persyst-Reader * Nervus .eeg: https://github.com/ieeg-portal/Nervus-Reader * Biopac .acq: https://github.com/ieeg-portal/Biopac-Reader * gTec EEG recordings: Read directly from the HDF5 files instead of the Matlab exports. |
|
Line 367: | Line 402: |
* BIDS-MEG: Automatic loading of BIDS-structured databases | |
Line 369: | Line 403: |
* gTec EEG recordings: Read directly from the HDF5 files instead of the Matlab exports. * Use new Matlab functions readtable/writetable (2014a): for Excel and text files |
* SPM .mat/.dat: Fix the import of the EEG/SEEG coordinates * Get acquisition date from files: Missing for 4D |
Line 373: | Line 407: |
* Account creation: Send authenticated emails: * https://support.google.com/mail/answer/81126 * https://support.google.com/mail/answer/180707 |
* Compile with FieldTrip and SPM * Add tags to the forum posts for easier listing by topic * Tutorial OMEGA/BIDS: * Add review of literature for the resting state MEG * Download example datasets directly from the OMEGA repository * New tutorials: <<BR>> * Other public datasets: [[https://github.com/INCF/BIDS-examples/tree/bep008_meg|https://github.com/INCF/BIDS-examples/tree/bep008_meg/]] * Rat PAC + high gamma (Soheila) * EEG/research * FieldTrip ECOG tutorial: http://www.fieldtriptoolbox.org/tutorial/human_ecog * FieldTrip cortico-muscular coherence tutorial: http://www.fieldtriptoolbox.org/tutorial/coherence * Reproduce tutorials from MNE-Python: https://martinos.org/mne/stable/tutorials.html * Cam-CAN database: https://camcan-archive.mrc-cbu.cam.ac.uk/dataaccess/<<BR>>(download new datasets, including maxfiltered files and manual fiducial placements) * MEG steady-state / high-gamma visual / frequency tagging * BIDS-EEG example datasets * Stand-alone ICA tutorial * Co-register MEG runs (Beth) * Move all the files to download on the cloud for faster download everywhere in the world |
Line 381: | Line 432: |
* Reference tutorials on Google scholar + ResearchGate | |
Line 384: | Line 434: |
* New tutorials: * MEG steady-state / high-gamma visual * MEG resting state + OMEGA (Guoimar, Jeremy, Beth) * Intra-cranial recordings (Marcel) * Connectivity: Cortico-muscular coherence * Co-register MEG runs (Beth) |
|
Line 394: | Line 438: |
* Clean up the wiki: * Remove all the wiki pages that are not used * Check all the links in all the pages * Check that all the TODO blocks have been properly handled * Remove useless images from all tutorials * Update page count on the main tutorials page |
|
Line 421: | Line 471: |
* Interface scaling: Rewrite class IconLoader to scale only once the icons at startup instead of at each request of an icon (might improve the speed of the rendering of the tree) |
What's next
A roadmap to the future developments of Brainstorm.
Contents
Recordings
RAW file viewer:
- Bottom bar: display extend events as segments instead of dots (to see their extension)
- Downsample before filtering? (attention to the filter design)
- Add parameter to make the visual downsampling more or less aggressive
- Pre-load next page of recordings
- Keep the filter specifications in memory instead of recomputing for every page
- Add field "comment" to markers: For clinicians to add notes (Marcel)
- Events: Change the category of a selected event easily, instead of deleting/marking new
Events: Advanced process for recombining.
Example: http://www.erpinfo.org/erplab/erplab-documentation/manual/Binlister.html
- EEG: Laplacian montage (see doc sent by Jeremy)
- Bad trials: When changing the status of bad to good: remove the bad segments as well, otherwise it is not processed by processes like the PSD.
- Review clinical recordings: Reduce the dimensionality of the data with a simple inverse problem, similar to what we do for the magnetic extrapolation ("Regional sources" in BESA, cf S Rampp)
- MEG/EEG registration: Apply the same transformation to multiple runs
- 2DLayout:
- Does not work when DC offset is not removed
- Add a proper amplitude scale that gets updated when shift+scroll, to compare figures
- Create heat maps: Maybe with matlab function heatmap?
Interface
Add a warning when computing a forward model with > 100000 sources (check selection)
- Multiple screens: Add option to set on which display the main bst window should be
File filters: Add boolean logic: https://github.com/brainstorm-tools/brainstorm3/issues/68
- Snapshot: Save as image / all figures (similar to Movie/all figure)
Generalize the use of the units (field .DisplayUnits): Rewrite processes to save the units correctly
- Colormaps:
- Allow brightness/contrast manipulations on the custom colormaps
- Create a colormap similar to MNE, where extrema are bright
- Global colormap max: Should get the maximum across all the open files
Set parula (or others) as the default, not jet: https://bids.github.io/colormap/
https://jakevdp.github.io/blog/2014/10/16/how-bad-is-your-colormap/
- Open new figures as tab (docked in the Figures window)
- Copy figures to clipboard (with the screencapture function)
- Display warning before opening files that are too big
- Smooth display from figure_image (ERPimage, raster plot...)
Contact sheets & movies: use average of time windows instead of single instants, for each picture.
- Contact sheets: Allow explicit list of times in input (+ display as in MNE-Python with TS)
- Display CTF coils: Show discs instead of squares
- Use boundary() instead of conhull() in all the display functions (ie. 2DDisc)
Connectivity
- Thresholding and stat tests the connectivity matrices
- Connectivity on unconstrained sources: "Default signal extraction for volume grids should be the time series of the first principal component of the triplet signals after each has been zero-meaned" (SB)
- Display of connectivity graphs:
- Display as straight lines
- Recode 2D graphs
- 3D display with anatomical constrains
- Display using real position of EEG electrodes
- Use new band-pass filters in bst_connectivity ('bst-hfilter' instead of 'bst-fft-fir')
Review by Jan-Mathijs: http://journal.frontiersin.org/article/10.3389/fnsys.2015.00175/full
- Connectivity based on band limited power (Sylvain):
- Compute Hilbert/Bandpass + correlation of the envelopes
- Bandpass envelopes before computing correlations?
- Compute Hilbert(sensors) and then project to source space?
- Matrix view of NxN graphs: Add legend of the elements along X and Y axis
- Graph view:
- Does not display negative values correctly (correlation or difference of coherence)
- Re-write using pure Matlab code and smoothed graphics
- Fixed scales for intensity sliders
- Text bigger
- Too much data in appdata
- Fixed scales for intensity sliders
- Add "=" shortcut for having graphs with similar configurations
- Disable zoom in one region (serious bugs)
- NxN on sensors: does not place the sensors correctly in space
- Coherence:
Average cross-spectra instead of concatenating epochs (to avoid discontinuities)
Explore inter-trial approaches (Esther refers to chronux toolbox)
- Granger:
- Crashes sometimes: improve stability
- Check for minimum time window (Esther: min around 500-1000 data points)
- Re-write and optimize code
- Add progress bar
- PLV:
- Add p-values
- Remove evoked
- Optimize code
- Add time integration
- Unconstrained sources
- Add warning when running of short windows (because of filters)
- PAC:
- Add input TF , to disconnect TF decomposition and PAC computation (Peter)
- Refine frequency vector of low frequencies
- How many central frequencies to use in bst_pac?
- Change filters: no chirplet functions
- bst_freqfilter: Use nfcomponents like in bst_pac
- Esther recommended a larger frequency binning of the PAC estimation
- PAC maps: Display all sensors at once (like TF and DynamicPAC)
- Hui-Ling's PAC:
- Time-resolved correlation/coherence: Display as time bands
- Other metrics:
- Coherence by bands: bst_coherence_band_welch.m
- Granger by bands: bst_granger_band.m
- Inter-trial coherence
Tutorial coherence [1xN] : Reproduce FieldTrip results?
Connect NxN: Display as time series > Display warning before trying to open too many signals
Processes
- Use CUDA for speeding up some operations (filtering, wavelets, etc)
- Allow processes in Python and Java
- Add MNE-Python functions:
- scikit-learn classifiers: S Marti / G Dehaene
- Implement data exchange with MNE-Python: write FIF files from Brainstorm and/or pass python objects in memory instead of FIF files
- SSS/tSSS cleaning
- Reproduce other tutorials / examples
- Change the graphic renderer from Matlab
Add FieldTrip functions:
- ft_sourceanalysis:
- Check noise covariance
- Check all the options of all the methods
- Single trial reconstructions + noise covariance?
Filters?? http://www.fieldtriptoolbox.org/example/common_filters_in_beamforming
Beamformers: Save ftSource.avg.mom
http://www.fieldtriptoolbox.org/workshop/meg-uk-2015/fieldtrip-beamformer-demohttp://www.fieldtriptoolbox.org/tutorial/beamformingextended
- Baseline? Two inputs?
ft_prepare_sourcemodel: Compute MNI transformation (linear and non-linear) => Peter
- Freqanalysis: ITC
- ft_read_atlas('TTatlas+tlrc.BRICK');
ft_volumereslice: http://www.fieldtriptoolbox.org/faq/how_change_mri_orientation_size_fov
- ft_freqanalysis
- ft_combineplanar
Allow FieldTrip functions in compiled version
- ft_sourceanalysis:
- Decoding/Classifiers: Faster algorithms (MNE-Python?)
- Pipeline editor:
- Bug: After "convert to continuous", the time of the following processes should change
- Add loops over subjects/conditions/trial groups
- Events: Allow selection from a drop-down list (similar to option "channelname" in panel_process_selection)
ITC: Inter-trial coherence (see MNE reports for group tutorial)
http://www.sciencedirect.com/science/article/pii/S1053811916304232- ICA:
- Why doesn't the ICA process converge when using 25 components in the EEG tutorial?
- Add an option to resample the signals before computing the ICA decomposition
- Add a stand-alone tutorial
- Exploration: Add window with spectral decomposition (useful for muscle artifacts)
- Export IC time series (and then compute their spectrum): solves the problem above
Comparison JADE/Infomax:
http://journals.plos.org/plosone/article?id=10.1371/journal.pone.0030135- Use faster methods (MNE-Python?)
- Add methods: SOBI, Fastica, AMICA/CUDICA (recommended by S Makeig)
Dimension reduction with PCA adds artifacts: Not done by default in EEGLAB
Contact: Stephen Shall Jones ( shall-jones@infoscience.otago.ac.nz )
Student Carl Leichter detailed this in his thesis- S Makeig: Use ICA to select the IC of interest instead of only removing artifacts
- Display of spectrum for components (PSD/FFT)
- Use FastICA (algo crashing)
- Add components preselection: Correlation with EOG/ECG
- Import ICA matrices available in EEGLAB .set files
EEGLAB recommends ICA + trial rejection + ICA again: Impossible right now with Brainstorm
(http://sccn.ucsd.edu/wiki/Chapter_09:_Decomposing_Data_Using_ICA)ICA+machine learning: https://www.ncbi.nlm.nih.gov/pubmed/28497769
Automated artifact rejection: https://arxiv.org/abs/1612.08194
- Save IC time series in database
Use EYE-EEG: EEGLAB toolbox for eye-tracker guided ICA (Olaf Dimigen): http://www2.hu-berlin.de/eyetracking-eeg/
- Other EEGLAB functions:
Step function detection: https://github.com/lucklab/erplab/wiki/Artifact-Detection:-Tutorial
- SSP:
Display warning if changing the ChannelFlag while there is a Projector applied
Show where the attenuation is projected:
(sum(IK,2)-sum(SSP(k,:)*IK,2)./sum(IK,2)
- Pipeline editor:
- When computing sources from the pipeline editor: doesn't reselect the options if you click twice on "edit" (works for minnorm, but not for lcmv)
- When computing time-frequency/hilbert/psd: Find a way not to force the user to click on Edit
- Bandpass:
- Use new filters in all the functions using a bandpass ('bst-hfilter' instead of 'bst-fft-fir'): process_evt_detect_threshold
- Weird bug: Filter(import) != Import(Filter) in the HCP tutorial... to investigate
- Bandpass
- Spectral flattening (John):
- ARIMA(5,0,1): Apply on the signal before any frequency/connectivity/PAC analysis
- PSD:
- Rewrite to have the same input as coherence (frequency resolution instead of window length)
- Use the progress bar
Allow display of Avg+StdErr
Remove line noise: http://www.nitrc.org/projects/cleanline
- Use Matlab Coder to optimize some processes: Wavelets, bandpass filter, sinusoid removal
- Reports: Click on link reopens exactly the figure
- Time-frequency:
- Optimization: bst_timefreq (around l.136), remove evoked in source space: Average should be computed in sensor space instead of source space (requested by Dimitrios)
Short-time Fourier transform: http://www.mikexcohen.com/lectures.html
Matching pursuit: http://m.jneurosci.org/content/36/12/3399.abstract?etoc
- Bug: Display logs as negative
- Bug: 3D figures: Colormaps with "log" option doesn't work
- Bug: Difference of power displayed in log: problems (Soheila)
- 2D Layout in spectrum
- Make much faster and more memory efficient (C functions coded by Matti ?)
- TF scouts: should display average of TF maps
- Impossible to keep complex values for unconstrained sources
- Pad short epochs with zero values for getting lower frequencies
- Hilbert with time bands very slow on very long files (eg. 3600s at 1000Hz) because the time vector is still full (10^7 values): save compressed time vector instead.
- Extend clusters tab to display of TF to overlay TF signals (Svet)
- When normalizing with baseline: Propagate with the edge effects marked in TFmask
- Allow baseline normalization of files computed with time bands
- Allow running TF on montages
- Review continuous files in time-frequency space (for epilepsy)
- Bug when computing TF on constrained and unconstrained scouts at the same time (in mixed head models for instance): uses only the constrained information and doesn't sum the 3 orientations for the unconstrained regions.
- Artifact detection:
- Artifact rejection like SPM: if bad in 20%, bad everywhere
- Test difference between adjacent samples
- Events detection: Add option "std" vs "amplitude"
- Co-registration of MEG runs:
- SSP: Group projectors coming from different files
- Finish validation of the method
- Apply to continuous recordings for correcting head movements
- Simulation:
Fix units in simulation processes => no *1e-9 in "simulate recordings"
- Use "add noise" process from Hui-Ling (in Work/Dev/Divers)
- Use field process field "Group" to separate Input/Processing/Output options
- Use new Matlab functions: movmean, movsum, movmedian, movmax, movmin, movvar, movstd
- Process ft_prepare_heamodel: Add support from BEM surfaces from the Brainstorm database
Database
- Faster DB searches (for Emily)
- Add buttons to sort files: by name, by comment, by date
- MEG protocols: More flexible organization of the database; sub-conditions to allow different runs X different conditions.
- Matrix files: Allow to be dependent from other files
- Add notes in the folders (text files, visible as nodes in the tree)
- Screen captures: save straight to the database
- Rename multiple files
- Default headmodel lost when reloaded: Keep selection on the hard drive (in brainstormstudy.mat)
- Auto-save:
protocol.mat can be too big: do not store the results links in it (and recreate when loading)- http://neuroimage.usc.edu/forums/t/abnormally-slow-behavior/2065/10
- Improve auto-save: add tracking file next to protocol.mat, do not save all the time, only when closing app, and reload protocol at stratup if tracking file is still there
Distributed computing
Options from FieldTrip:
Loose collection of computers: https://github.com/fieldtrip/fieldtrip/tree/master/peer
Alternative, with less limitations: http://research.cs.wisc.edu/htcondor/
Single multicore machine: https://github.com/fieldtrip/fieldtrip/tree/master/engine
Batch system: https://github.com/fieldtrip/fieldtrip/tree/master/qsub
Documentation: http://fieldtrip.fcdonders.nl/faq#distributed_computing_with_fieldtrip_and_matlab
Various initiatives: http://samirdas.github.io/Data_sharing.html#/
Source modeling
- Project sources: Very poor algorithm to project sub-cortical regions and cerebellum (algorithm to fit surfaces should be imrpoved)
Menu head model > Copy to other conditions/subjects (check if applicable first)
Menu Sources > Maximum value: Doesn't work with volume or mixed head models
Mixed head models:
- Set loose parameter from the interface
- Volume grid:
- Optimize: 3D display (better than 9x9 cubes)
- Optimize: vol_dilate (with 26 neighbors)
Menu Sources > Simulate recordings:
- Do not close the 3D figures after generating a new file
- Add a process equivalent to this menu
- Panel Get coordinates: Add button "find maximum"
- Dipoles:
- Project individual dipoles files on a template
- panel_dipoles: Doesn't work with multiple figures
- BEM single sphere: Get implementation from MNE
- Visualize Beamformer results:
- Read CTF SAM .svl
- Display as layers in the MRI viewer
- Unconstrained sources:
- Stat and connectivity: what to do? (re-send email John+Sylvain)
- Overlapping spheres: improve the estimation of the spheres for the frontal lobes
- Magnetic extrapolation: Do the same thing with EEG
- Noise covariance matrix:
- Display with figure_image()
- Storage of multiple noise covariance matrices (just like the head models)
- Always save as full, then at inversion time, we can decide between full, heteroskedastic (diagonal) or homoskedastic (i.i.d, scalar)
Problem of having inividual trials + averages in the condition => Display warning or not?
- Save nAvg in noisecov file, to make it easier to scale to other recordings
Sources on surface: Display peak regions over time (time = color) => A.Gramfort
Calculate ImagingKernel * Gain for a scout
- Time-frequency beamformers:
- Band-pass everything in different frequency bands + Source estimation + TF
- Ask data to Sarang where he sees effects that cannot be extracted with MN followed by TF
- Process "Extract scouts time series": Add PCA option (replace isnorm with choice PCA/Norm)
- BEM: Fix unstable results when one vertex is too close from the layers (5mm ?)
- Hui-Ling beamformers:
- More explanations about what is in NAI and Spatial filters
Explain that is this is better to study effects extended in time (Ntime > Nsensors)
- Group LCMV+MCB
- Condition LEFT median nerve: very bad results
- Keep options for inverse computation
- OpenMEEG: Post example datasets for the remaining issues:
- Example protocol ECOG: doesn't work
- Add eyes models to attract eye activity
Anatomy
BrainVISA: Add support for MarsAtlas
- MNI transformation: Use non-linear MNI transformation y_...
Registration:<<BR>>
Getting electrode positions from 3D scanners: https://sccn.ucsd.edu/wiki/Get_chanlocs
Use the same registration for multiple recording sessions that have already re-registered previously (eg. with MaxFilter)
- When linking multiple EEG recordings including 3D positions, do the registration only once and copy it to all the runs
- Check templates: MNI transformation and volume atlases for ICBM152 vs Colin27 (loading the AAL atlas as surface or voume scouts do not align well on the ICBM152)
- Compute non-linear MNI registration instead of linear
- Select and remove bad digitized head points before automatic coregistration
Load the MNE -transf.fif: http://neuroimage.usc.edu/forums/showthread.php?2830
- MRI Viewer:
- Add display of "world coordinates" when "vox2ras" is available
- Pan in zoomed view (shift + click + move?)
- Zoom in/out with mouse (shift + scroll?)
- Ruler tool to measure distances
- Add keyboard shortcuts to scroll in the three orientations (same in MRI 3D)
- Display scouts as additional volumes
- Render surface envelope in the MRI as a thin line instead of the full interpolation matrix
- Edit fiducials: Replace 6 text boxes with 1 for easy copy-paste (see fiducials.m)
Optimize computation interpolation MRI-surface (tess_tri_interp) => spm_mesh_to_grid
BrainSuite:
Add new labels to all BrainSuite anatomy templates
Use BrainSuite inner skull for surface generation
- Use same colors for left and right for anatomical atlases
- Scouts:
- Display edges in the middle of the faces instead of the vertices
- Display scouts in a tree: hemisphere, region, subregion
- Sort scouts by region in process options
- Downsample to atlas: allow on timefreq/connect files
Project from one hemisphere to the other using registered spheres/squares (http://neuroimage.usc.edu/forums/t/how-to-create-mirror-roi-in-the-other-hemisphere/5910/8)
- Parcellating volume grids: scikit-learn.cluster.Ward
- Major bug when importing surfaces for an MRI that was re-oriented manually
- Smooth surface: Fix little spikes due to irregularities in the mesh
Surface>Volume interpolation: Use spm_mesh_to_grid
- Bug: Hide scouts in the preview of the grid for volume head models
ECOG/SEEG
- Electrode models: Better way for managing/updating/adding electrode models
SEEG: Project contacts to SEEG / Compare with IntrAnat
Contact positions: Import / set / detect
- New option: Align on none|inner|cortex to replace ECOG-mid
- Add history: Save modifications and transformations applied to the channel files (Marcel)
- ECOG: How to handle cases where not all the grid contacts are in the channel file? (Marcel)
- Project contact positions across subjects or templates (Marcel)
- ECOG: Default names of contacts?
- Add menu to import implantation channel file in imported recordings
- Automatic segmentation of CT:
- ECOG:
- Project and display contacts on cortex surface should consider the rigidity of the grids: Contacts cannot rotate, and distance between contacts should remain constant across runs
Method for contacts projection: https://pdfs.semanticscholar.org/f10d/6b899d851f3c4b115404298d7b997cf1d5ab.pdf
ECOG: Brain shift: When creating contact positions on a post-implantation image, the brain shift should be taken into account for creating images of the ECOG contacts on the pre-op brain => iELVis (http://ielvis.pbworks.com/w/page/116347253/FrontPage)
- Display:
- Bad channels: Contacts greyed out instead of ignored (Marcel)
- Display time in H:M:S
- Display curved SEEG electrodes
- Display wires for ECOG?
- Re-referencing:
- Create new average reference montages with a specific list of channels, with the possibility to edit the order of the channels (for Jeremy)
- Closest white reference (Arnulfo)
Detection CEEP stim artifacts: Use ImaGIN code ImaGIN_StimDetect
Alternatives to OpenMEEG: SimBio/FieldTrip? Matti Stenroos? NFT/NIST?
Statistics
- ANOVA:
- Use LENA functions(?)
- Output = 1 file per effect, all grouped in a node "ANOVA"
- Display several ANOVA maps (from several files) on one single figure, using a "graphic accumulator", towards which one can send any type of graphic object
- New process to test for Gaussianity using swtest
- Simulate recordings with specific properties, for stat validation
Quality control before statistics, on condition averages across subjects:
mean(baseline)/std(baseline): shows bad subject quickly.Use SurfStat: Impements interesting things, like an analytical cluster-based p-value correction (Random-field theory which is used in SPM) - Peter
Input / output
XDF import: Use the EEGLAB plugin, contact Martin Bleichner (Oldenburg)
- Output .nii have incorrect sform/qform when using the options to downsample the volume of cut the empty slices
- DICOM converter:
- Add dcm2nii (MRICron)
- Add MRIConvert
FieldTrip: Import/Export time-frequency:
- 4D file format:
- Use reader from MNE-Python: mne.io.read_raw_kit (doesn't require Yokogawa slow library)
- Reference gradiometers: Keep the orientation of the first or second coil?
- Reference gradiometers: Add the sensor definition from coil_def.dat
- Validate with phantom recordings that noise compensation is properly taken into account
- References at too far from the head sensors in Marseille 4D system
- The noise compensation is considered to be always applied on the recordings, not sure this assumption is always correct
- 4D phantom tutorial (JM Badier?)
- EEG File formats:
EEG CeeGraph
- EGI: Finish support for epoched files (formats 3,5,7)
Persyst .lay: https://github.com/ieeg-portal/Persyst-Reader
Nervus .eeg: https://github.com/ieeg-portal/Nervus-Reader
Biopac .acq: https://github.com/ieeg-portal/Biopac-Reader
- gTec EEG recordings: Read directly from the HDF5 files instead of the Matlab exports.
- BCI2000 Input (via EEGLAB plugin)
- EEGLAB import:
- Support for binary AND epoched files (now it's one or the other)
- Allow epoched files with recordings saved in external files
- BST-BIN: Add compressionto .bst
- Review raw on all the file formats (ASCII EEG and Cartool missing)
- SPM .mat/.dat: Fix the import of the EEG/SEEG coordinates
- Get acquisition date from files: Missing for 4D
Distribution & documentation
Compile with FieldTrip and SPM
- Add tags to the forum posts for easier listing by topic
- Tutorial OMEGA/BIDS:
- Add review of literature for the resting state MEG
- Download example datasets directly from the OMEGA repository
New tutorials:
Other public datasets: https://github.com/INCF/BIDS-examples/tree/bep008_meg/
- Rat PAC + high gamma (Soheila)
- EEG/research
FieldTrip ECOG tutorial: http://www.fieldtriptoolbox.org/tutorial/human_ecog
FieldTrip cortico-muscular coherence tutorial: http://www.fieldtriptoolbox.org/tutorial/coherence
Reproduce tutorials from MNE-Python: https://martinos.org/mne/stable/tutorials.html
Cam-CAN database: https://camcan-archive.mrc-cbu.cam.ac.uk/dataaccess/<<BR>>(download new datasets, including maxfiltered files and manual fiducial placements)
- MEG steady-state / high-gamma visual / frequency tagging
- BIDS-EEG example datasets
- Stand-alone ICA tutorial
- Co-register MEG runs (Beth)
- Move all the files to download on the cloud for faster download everywhere in the world
- Provide secure way of sending password over HTTPS for:
- Account creation
- Forum exchanges
org.brainstorm.dialog.CloneControl
Workflows FieldTrip: http://www.fieldtriptoolbox.org/faq/what_types_of_datasets_and_their_respective_analyses_are_used_on_fieldtrip
Count GitHub clones in the the download stats
- Google Analytics: Create template and update the section of the Community page
- Missing in page "Cite Brainstorm": Add all the methods used in the software
- Deface the MRIs of all the tutorials
- Compiled R2016b: Color picker doesn't work (for changing surface color for instance)
- Clean up the wiki:
- Remove all the wiki pages that are not used
- Check all the links in all the pages
- Check that all the TODO blocks have been properly handled
- Remove useless images from all tutorials
- Update page count on the main tutorials page
Current bugs
- Screen capture:
- Bug on Win8/Win10: doesn't capture the correct part of the screen
- Window managers with fading effect: captures the top window
- Image viewer:
- Difficult to get to 100%
- Buggy on some systems
- 2DLayout:
- (TF) Units are weird with % values
- (TF) Difficult to navigate in frequencies: Scaling+changing frequency resets the scaling
- Progress bar:
- Doesn't close properly on some Linux systems
- Focus requests change workspace when processing constantly (Linux systems)
- MacOS bugs:
- Buttons {Yes,No,Cancel} listed backwards
- Record tab: Text of epoch number is too big
- Colormap menus: Do not work well on compiled MacOSX 10.9.5 and 10.10
- Matlab bugs:
- Interface looks small on screens with very high resolutions: Reduce the resolution
- Event markers are not visible anymore with the sequence: Open MEG, open EOG, close MEG
- in_bst_data_multi: If trials have different sizes, output is random (the one of the first file)
- Edit scout in MRI: small modifications cause huge increase of the scout size
- Canolty maps computation: Fix progress bar
Geeky programming details
- bst_bsxfun: After 2016b, we can use directly the scalar operators (./ .* ...) instead of bsxfun. Update bst_bsxfun to skip the use of bsxfun when possible.
Interface scaling: Rewrite class IconLoader to scale only once the icons at startup instead of at each request of an icon (might improve the speed of the rendering of the tree)
- Hide Java panels instead of deleting them
- Processes with "radio" and "radio_line" options: Replace with "radio_label" and "radio_linelabel"
- Interpolations: Use scatteredInterpolant, griddedInterpolant, triangulation.nearestNeighbor (2014b)
- bst_warp and channel_project: Use tess_parametrize_new instead of tess_parametrize
- Shared kernels: "get bad channels" operation in a different way (reading all the files is too slow)
- Optimize bst_get:
- Now study and subject have necessarily the same folder name
- Replace big switch with separate functions
- Progress bar:
- Add different levels (to handle sub-processes)
- Make work correctly with RAW on resting tutorial
- Uniformize calls in bst_process/Run
- Add a "Cancel" button
- Fix all the 'todo' blocks in the code
- Error message: Add a link to report directly the bug on the forum