|
Size: 11665
Comment:
|
Size: 11652
Comment:
|
| Deletions are marked like this. | Additions are marked like this. |
| Line 5: | Line 5: |
| ==== Data acquisition ==== * Improve the acquisition of the head points and the registration MEG / MRI * Real-time processing and display (based on !FieldTrip realtime module) * Support for intra-cranial recordings ==== Pre-processing ==== * Correction for head movements (using the continuous head localization coils) ==== Source modeling ==== * Computation of equivalent current dipoles * Display results of CTF SAM beamformer |
|
| Line 21: | Line 9: |
| ==== Data management ==== * Filter the display of the tree by type / name / tag |
==== EEG / epilepsy / intra-cranial recordings ==== * Editing the position of intracranial electrodes in the MRI viewer ==== Source modeling ==== * Computation of equivalent current dipoles * Beamformers |
| Line 33: | Line 25: |
| * Allow multiple RAW windows * Screen setups * Time scale: define in fixed s/mm (like the CTF tools) * Secondary windows: display length of time selection * If "Use SSP " option is selected, automatically select "Remove baseline" and "CTF compensations" |
|
| Line 39: | Line 26: |
| * Set the amplidute cale for the time series * CTRL+S : Save modifications |
* 2DLayout: Doesn't work when changing page => need refresh of !GlobalData.Preferences.!TopoLayoutOptions.!TimeWindow * EEG reference/storage: * Intracranial electrodes: Define in the MRI viewer * Bad channels that can be specified at the program level (for sites that have permanently bad channels) => AS Dubarry |
| Line 42: | Line 31: |
| * Make it work for all the file formats(at least bandpass filter) | * Make it work for all the file formats (at least bandpass filter + sin removal) |
| Line 44: | Line 33: |
| * Imported recordings: Offer the same interface as the RAW viewer: * Scroll bar * Events viewer / editor * Re-epoch (import from files in database) * SSP * Bad channels that can be specified at the program level (for sites that have permanently bad channels) => AS Dubarry |
|
| Line 51: | Line 34: |
| * Define manually minimum => 3 options: [0,max], [-max,max], [min,max] | |
| Line 53: | Line 35: |
| * Grey out the portion of the colorbars not displayed because of the threshold * bst_selections: * Add user defined combinations of sensors (eg. "double banana" for EEG) * Use this to produce "inversed polarity" displayes too (useful in EEG) |
|
| Line 62: | Line 40: |
| * Intracranial electrodes: * Define and display in the MRI viewer * Images of amplitude: [sensor x time], [trial x time], scout: [trial x time] (similaire to erpimage in eeglab) |
* Images of amplitude: [sensor x time], [trial x time], scout: [trial x time] * Can be done with Matrix > View as image: extract cluster, concatenate for all trials * 2D Layout for multiple conditions == Display == * Nicer 2D plots, standardized * Copy figures to clipboard (with the screencapture function) |
| Line 67: | Line 49: |
| * Figures: interaction with sensor selection / scouts selections * Calculation and use of p-values for the metrics * Trial analysis: Concatenations rather than averages for most metrics * Optimize bst_coherence * Display NxN as 1xN * Adapt colormaps for correlation (min and max properties) * PLV: Add a time integration |
* Tutorial coherence [1xN] * t-tests on connectivity measures * Graph view: * Fixed scales for intensity sliders * Fix zoom in one region * Text bigger * Too much data in appdata * Other metrics: * Coherence by bands: bst_coherence_band_welch.m * Granger by bands: bst_granger_band.m * Inter-trial coherence |
| Line 75: | Line 61: |
| * Circle plot: * Display Neuromag sensors * Use little squares * Little squares with scout color instead of point * Max distance slider: useless * Event-related coherence? |
|
| Line 84: | Line 64: |
| * Stat computed on time-frequency data * Validate all the operations on the time-freq files (check for mixed Measures) |
|
| Line 87: | Line 65: |
| * Scouts values for timefreq on surfaces | |
| Line 92: | Line 69: |
| * Parallel processing: Use parfor * Distributed processing: * Version of Brainstorm that can run without JAVA * Use the BIC/MEG cluster from the pipeline interface (PSOM, P.Bellec) |
* Smooth display of TF/PAC maps (option) * TF scouts: should display average of TF maps * Bandpass: Show warning when using inappropriate high-pass filter (precision too high) * Artifact detection: * Detection of bad segments in the RAW files * Artifact rejection like SPM: if bad in 20%, bad everywhere * Test difference between adjacent samples * Distributed processing: Brainstorm that can run without Java |
| Line 98: | Line 79: |
| * When processing multiple files: waitbar is not behaving well | * Show where the attenuation is projected:<<BR>>(sum(IK,2)-sum(SSP(k,:)*IK,2)./sum(IK,2) |
| Line 101: | Line 82: |
| * By subject AND condition | |
| Line 107: | Line 87: |
| * Apply to continuous recordings for correcting for head movements (using head position coils) | * Apply to continuous recordings for correcting head movements |
| Line 110: | Line 90: |
| * Detrending | |
| Line 116: | Line 97: |
| * Optical flow | |
| Line 118: | Line 100: |
| * Filter display of the database explorer (filename, file type, comment...) | |
| Line 120: | Line 101: |
| * GUI: Save configuration of windows (per protocol) | * Group matrix files => allow to process matrix files by trial types |
| Line 124: | Line 105: |
| * Screen captures: save straight to the database * Rename multiple files |
|
| Line 125: | Line 110: |
| * Stenroos 2014 paper: Include the following methods * Inner and outer skull surfaces generator from !FieldTrip (needs SPM, probably not so different from BST) * Nolte corrected-sphere model (good model re:Alex) * Fast BEM models |
|
| Line 126: | Line 115: |
| * Scouts: * Represent border as the middle of the adjacent triangles (to have a full segmentation) * Visualize Beamformer results (contact Zainab Fatima): |
* Visualize Beamformer results: |
| Line 130: | Line 117: |
| * Create new file type in the database | |
| Line 134: | Line 120: |
| * Define as default * Check all the processes |
|
| Line 137: | Line 121: |
| * Stat and connectivity: what to do? (re-send email John+Sylvain) | |
| Line 143: | Line 128: |
| * Optimize: grid_interp_mri * Magnetic extrapolation: * Do the same thing with EEG * Project sources: * Adapt smooth factor to the number of vertices * Number of neighbors to consider = average number of neighbors in the target mesh. * Compute by small time blocks |
* Magnetic extrapolation: Do the same thing with EEG |
| Line 159: | Line 138: |
| * EEG Source modeling: Manage references and bipolar montages properly (maybe not necessary) * MEG source modeling: Do reconstruction only for a subset of sensors for estimating dipoles? * Processes compute head model and sources: Additional option to set the file comment |
|
| Line 161: | Line 143: |
| * Use the FreeSurfer subject co-registration (downsample: keep the correspondance) | * Scouts:<<BR>> * Display edges in the middle of the faces instead of the vertices * Project scouts betweens subjects and between hemispheres * Display scouts in a tree: hemisphere, region, subregion * Downsample to atlas: allow on timefreq/connect files * Sort scouts by region in process options * Generate mixed density surfaces |
| Line 163: | Line 151: |
| * Improve ICP registration headpoints / scalp (chanfrein, multi-resolution, see with C Grova...) * Auto-reorientation of MRI after selected NAS / LPA / RPA |
|
| Line 166: | Line 152: |
| * ICBM brain * MINC MRI reader: EMMA, NIAK (Pierre Bellec), HDF5 directly read in Matlab * ICBM average surfaces + atlas * Using CIVET pipeline for extracting surfaces * Clustering cortex: Dimitrios, David, Yu-Teng |
* Use mid-gray instead of pial surface? |
| Line 174: | Line 155: |
| * Stat on scouts / clusters / "matrix" | |
| Line 180: | Line 160: |
| * http://www.adscience.fr/uploads/ckfiles/files/html_files/StatEL/statel_wilcoxon.htm * http://www.mathworks.fr/fr/help/stats/signrank.html * Less powerful than t-tests |
|
| Line 186: | Line 169: |
| * t-test on volume sources * Paired t-test on unconstrained sources: (convert to flat + Z-score) => !AnneSo * Question of Gaussianity of the samples: take a subset of samples + Kolmogorov-Smirnov / Shapiro-Wilk test * http://fr.wikipedia.org/wiki/Test_de_Shapiro-Wilk * http://stats.stackexchange.com/questions/362/what-is-the-difference-between-the-shapiro-wilk-test-of-normality-and-the-kolmog * http://www.mathworks.fr/fr/help/symbolic/mupad_ug/perform-shapiro-wilk-test.html * http://www.mathworks.fr/fr/help/symbolic/mupad_ref/stats-swgoft.html * http://stackoverflow.com/questions/14383115/shapiro-wilk-test-in-matlab * Create icons for Stat/PAC, Stat/Sprectrum, etc. * One sample t-test across subjects |
|
| Line 188: | Line 181: |
| * MRI: MINC format | |
| Line 191: | Line 183: |
| * NEUROFILE = COHERENCE EEG/video !LongTerm Monitoring | |
| Line 195: | Line 186: |
| * !FreeSurfer: orig.mgz | * Export TF maps to SPM / volumes * EEGLAB import: Selection of conditions in script mode |
| Line 198: | Line 190: |
| * Shortcuts: * Add Help buttons and menus (in popups, dialog windows...) => Links to the website. * List of all the keyboard and mouse shortcuts * Equivalents for MacOS |
* Add Help buttons and menus (in popups, dialog windows...) => Links to the website. |
| Line 203: | Line 192: |
| * Estimate time to complete each tutorial | * Processes: Describe all the processes |
| Line 205: | Line 194: |
| * Anatomy: Segmentation with !FreeSurfer | |
| Line 208: | Line 196: |
| * Coordinate sytems: How to convert between the different coordinates systems in scripts | |
| Line 211: | Line 198: |
| * Processes: Describe all the processes * Processes: How to write your own processes (user folder for processes) |
|
| Line 215: | Line 200: |
| * How to export sources for analysis in SPM * Advanced tutorials: * EEG (How to import an EEG cap) * Epilepsy / spike analysis * How to make and compress a movie (Brainstorm + !VirtualDub + XVid) * Ask users to send their channel files, align on Colin, distribute |
== Current bugs == * Import anatomy folder menu freezes Matlab on MacOS 10.9.3 / Matlab 2014a * Import anatomy folder menu crashes on MacOSX 10.8.5 / Matlab 2013b * Record tab: Text of epoch number is too big on MacOS * Menu "Use default EEG cap" doesn't work for a multiple selection (setting the same EEG cap for several subjects) * in_bst_data_multi: If trials have different sizes, output is random (the one of the first file) * tree_dependencies: sources files, reprojected on default anatomy; If based on data files that are bad trials, they should be ignored by tree_dependencies, and they are not * Image viewer has some bugs on some systems * Screen capture for reports never works: Find another solution * Screen capture when there is a fading effect in the window manager: captures the window * Close figure with coherence results should hide the "frequency" slider * Problems growing scouts on merged surfaces (Emily) * Edit scout in MRI: small modifications cause huge increase of the scout size |
| Line 223: | Line 216: |
| * Cleaning surfaces: should work with atlases * Send OpenGL bug report to the Mathworks |
* Interpolations: Use scatteredInterpolant, griddedInterpolant, triangulation.nearestNeighbor (Matlab 2014b) * Use Matlab Coder to optimize some processes: Bandpass filter, sinusoid removal |
| Line 227: | Line 220: |
| * Zip files created cannot be open with !WinZip * OpenGL options: {none, software, hardware} * Bug: Scout without overlay, adapt scale for each graph when "Uniform amplitude" option is unchecked (mixing sources + zscores) * Waitbars: * Replace old waitbars with java ones * Add a "Cancel" button on waitbars when the bounds are defined (ie. when bst controls the process) * Bug: Menu "Use default EEG cap" doesn't work for a multiple selection (setting the same EEG cap for several subjects) * Bug tree_dependencies: sources files, reprojected on default anatomy; If based on data files that are bad trials, they should be ignored by tree_dependencies, and they are not |
|
| Line 236: | Line 221: |
| * Bug in_bst_data_multi: If trials have different sizes, output is random (the one of the first file)... | |
| Line 238: | Line 222: |
| * Write shepards.m with new algorithm for nearest neighbors * Use Matlab Coder to compile / optimize some processes * Optimize calls to bst_get, now study and subject have necessarily the same folder name |
* 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 * Line smoothing / anti-aliasing (time series figures) * Fix all the 'todo' blocks in the code * Replace handle "0" with bst_get('groot') |
What's next
A roadmap to the future developments of Brainstorm.
Current topics
Functionnal connectivity
Integration of different metrics to study the brain connectivity:
Correlation, coherence, Granger causality, phase locking value- Development of new ways to represent the connectivity between sensors or brain regions
EEG / epilepsy / intra-cranial recordings
- Editing the position of intracranial electrodes in the MRI viewer
Source modeling
- Computation of equivalent current dipoles
- Beamformers
Large scale analysis
- Parallel processing: Reduce the computation times using the parallel processing toolbox
- Distributed processing: Integrate tools for sending Brainstorm processes on computation clusters
Recordings
- RAW file viewer:
- Pre-load next page of recordings
- Documentation: Add definition of bad segments
2DLayout: Doesn't work when changing page => need refresh of GlobalData.Preferences.TopoLayoutOptions.TimeWindow
- EEG reference/storage:
- Intracranial electrodes: Define in the MRI viewer
Bad channels that can be specified at the program level (for sites that have permanently bad channels) => AS Dubarry
- RAW processing:
- Make it work for all the file formats (at least bandpass filter + sin removal)
Events: advanced process for recombining. Example: http://www.erpinfo.org/erplab/erplab-documentation/manual/Binlister.html
- Colormaps:
- Create a colormap similar to MNE, where extrema are bright
- Import data:
- Save properties "baseline" and "resample" at the level of the protocol (to re-use for all the files)
- NIRS:
- Add new data type
- Display of sensors by pairs oxy/deoxy (red/blue), overlaid
- Images of amplitude: [sensor x time], [trial x time], scout: [trial x time]
Can be done with Matrix > View as image: extract cluster, concatenate for all trials
- 2D Layout for multiple conditions
Display
- Nicer 2D plots, standardized
- Copy figures to clipboard (with the screencapture function)
Connectivity
- Tutorial coherence [1xN]
- t-tests on connectivity measures
- Graph view:
- Fixed scales for intensity sliders
- Fix zoom in one region
- Text bigger
- Too much data in appdata
- Other metrics:
- Coherence by bands: bst_coherence_band_welch.m
- Granger by bands: bst_granger_band.m
- Inter-trial coherence
- Work on progress bars
Processes
- Time-frequency:
- Frequency bands: extended syntax (ex: [2 3 4], 10:5:90, ...)
- How to combine 3 orientations for unconstrained sources
- Display logs as negative
- 2D Layout in spectrum
- Make much faster and more memory efficient (C functions coded by Matti ?)
- Smooth display of TF/PAC maps (option)
- TF scouts: should display average of TF maps
- Bandpass: Show warning when using inappropriate high-pass filter (precision too high)
- Artifact detection:
- Detection of bad segments in the RAW files
- Artifact rejection like SPM: if bad in 20%, bad everywhere
- Test difference between adjacent samples
- Distributed processing: Brainstorm that can run without Java
- 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)
- Average:
- Remember how many trials were used per channel
- Save standard deviation
- Display standard deviation as a halo around the time series
- 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
Current Source Density (CSD) => Ghislaine
http://psychophysiology.cpmc.columbia.edu/software/CSDtoolbox/index.html- Other processes:
- Detrending
- Moving average
- Max
- Median
- Significance test (Dimitrios, Leo)
- Spatial smoothing: check / document parameters
Contact sheets & movies: use average of time windows instead of single instants, for each picture.
- Optical flow
Database
- MEG protocols: More flexible organization of the database; sub-conditions to allow different runs X different conditions.
Group matrix files => allow to process matrix files by trial types
- Add notes in the folders (text files, visible as nodes in the tree)
- Screen captures: save straight to the database
- Rename multiple files
Source modeling
- Stenroos 2014 paper: Include the following methods
Inner and outer skull surfaces generator from FieldTrip (needs SPM, probably not so different from BST)
- Nolte corrected-sphere model (good model re:Alex)
- Fast BEM models
- Dipole fitting
- Visualize Beamformer results:
- Read CTF SAM .svl
- Display as layers in the MRI viewer
- Unconstrained sources:
- Compute unconstrained and then project on the normal ?
- Difference and stat should be: norm(A) - norm(B)
- Stat and connectivity: what to do? (re-send email John+Sylvain)
- Overlapping spheres: improve the estimation of the spheres for the frontal lobes
- Volume grid:
- Scouts 3D
- Test volume sources with all the subsequent processes (timefreq, stat...)
- Optimize: 3D display (better that 9x9 cubes)
- Optimize: vol_dilate (with 26 neighbors)
- Magnetic extrapolation: Do the same thing with EEG
- Noise covariance matrix:
- 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
When deploying to other conditions: Apply destination SSP (NoiseCov = SSP . NoiseCov . SSP' )
Sources on surface: Display peak regions over time (time = color) => A.Gramfort
- Simulation: synthesize pseudo data-files from a cortex patch (duration, amplitude, noise)
Calculate ImagingKernel * Gain for a scout
- EEG Source modeling: Manage references and bipolar montages properly (maybe not necessary)
- MEG source modeling: Do reconstruction only for a subset of sensors for estimating dipoles?
- Processes compute head model and sources: Additional option to set the file comment
Anatomy
Scouts:<<BR>>
- Display edges in the middle of the faces instead of the vertices
- Project scouts betweens subjects and between hemispheres
- Display scouts in a tree: hemisphere, region, subregion
- Downsample to atlas: allow on timefreq/connect files
- Sort scouts by region in process options
- Generate mixed density surfaces
- Import / registration:
- Major bug when importing surfaces for an MRI that was re-oriented manually
- Use mid-gray instead of pial surface?
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
- Permutation tests:
- t-test only (wilcoxon? sign-test?): paired, equal var, unequal var
- nb permutations ~ 1000
- maximum statistic over "time" or "time and space"
- Permutations / clustering: cf fieldtrip
Threshold in time: keep only the regions that are significative for contiguous blocks of time, or over a certain number of time points
=> Process that creates a static representation of a temporal window- t-test on volume sources
Paired t-test on unconstrained sources: (convert to flat + Z-score) => AnneSo
- Question of Gaussianity of the samples: take a subset of samples + Kolmogorov-Smirnov / Shapiro-Wilk test
http://www.mathworks.fr/fr/help/symbolic/mupad_ug/perform-shapiro-wilk-test.html
http://www.mathworks.fr/fr/help/symbolic/mupad_ref/stats-swgoft.html
http://stackoverflow.com/questions/14383115/shapiro-wilk-test-in-matlab
- Create icons for Stat/PAC, Stat/Sprectrum, etc.
- One sample t-test across subjects
Input / output
- EEG File formats:
EEG CeeGraph
- EGI: Finish support for epoched files (formats 3,5,7)
FieldTrip structures: In / Out
- BCI2000 Input (via EEGLAB plugin)
- Export TF maps to SPM / volumes
- EEGLAB import: Selection of conditions in script mode
Distribution & documentation
Add Help buttons and menus (in popups, dialog windows...) => Links to the website.
- Introduction tutorials:
- Processes: Describe all the processes
- Clusters
- First steps: Brainstorm preferences
- Headmodel: explain the fields + how to get the constrained leadfield
- Sources: Modelized data
- Sources: theshold min. size (not documented yet)
- Import raw recordings: Add "detect bad trials/channels" in the pipeline
- Temporary folder
Current bugs
- Import anatomy folder menu freezes Matlab on MacOS 10.9.3 / Matlab 2014a
- Import anatomy folder menu crashes on MacOSX 10.8.5 / Matlab 2013b
- Record tab: Text of epoch number is too big on MacOS
- Menu "Use default EEG cap" doesn't work for a multiple selection (setting the same EEG cap for several subjects)
- in_bst_data_multi: If trials have different sizes, output is random (the one of the first file)
- tree_dependencies: sources files, reprojected on default anatomy; If based on data files that are bad trials, they should be ignored by tree_dependencies, and they are not
- Image viewer has some bugs on some systems
- Screen capture for reports never works: Find another solution
- Screen capture when there is a fading effect in the window manager: captures the window
- Close figure with coherence results should hide the "frequency" slider
- Problems growing scouts on merged surfaces (Emily)
- Edit scout in MRI: small modifications cause huge increase of the scout size
Geeky programming details
- Interpolations: Use scatteredInterpolant, griddedInterpolant, triangulation.nearestNeighbor (Matlab 2014b)
- Use Matlab Coder to optimize some processes: Bandpass filter, sinusoid removal
- Hide Java panels instead of deleting them
- mri2scs: convert arguments to meters
- bst_warp and channel_project: Use tess_parametrize_new instead of tess_parametrize
- Shared kernels: do the "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
- Line smoothing / anti-aliasing (time series figures)
- Fix all the 'todo' blocks in the code
- Replace handle "0" with bst_get('groot')