10073
Comment:
|
11262
|
Deletions are marked like this. | Additions are marked like this. |
Line 2: | Line 2: |
A roadmap for the future developments of Brainstorm. == Next six months == ==== Pre-processing ==== * Extend processing of continuous CTF files to all file formats * Continuous viewer: Many small improvements to the usability of its interface * Improve detection and correction of artifacts with SSP * Co-registration of several MEG runs on one single head position * Make all the main operations available in the pipeline editor ==== File formats ==== * MINC MRI * Stellate EEG * Brain Products / BrainAmp EEG * Read and display NIRS recordings * SAM Beamformer results |
A roadmap to the future developments of Brainstorm. == Current topics == |
Line 20: | Line 6: |
* Implementation of methods developed at USC <<BR>><<BR>> == Details of those features == * SSP: * Make SSP projections dynamic, and keep the full list instead of always them combining them * Take the bad channels in account in the application of the SSP * Refresh display after changing the list of bad channels (imported files + raw) so the SSP are applied correctly * Write documentation |
* Integration of different metrics to study the brain connectivity: <<BR>>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 ==== * Full EEG/epilepsy online tutorial * 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 <<BR>><<BR>><<BR>> == Recordings == |
Line 31: | Line 25: |
* Buttons for changing scales (time + amplitude) * Add tooltips to all the buttons / shortcuts * Adding events using configurable shortcuts (CTRL+number) |
* Pre-load next page of recordings * Allow multiple RAW windows(columns display) * Screen setups: save configuration of windows * Time scale: define in fixed s/mm (like the CTF tools) * Secondary windows: display length of time selection |
Line 36: | Line 32: |
* Homogenize a selection of several subjects/conditions * Popup menu when more than one study selected * Creation of a common channel file (match channels by names, not by order) * Register MEG runs (recompute fields for a different set of sensors, MEGCoregister from old brainstorm) * Complete processing pipeline: * Import + pre-process * Sources / head model / noise covariance * Project sources * bst_selections: * Add user defined combinations of sensors (eg. "double banana" for EEG) * Use this to produce "inversed polarity" displayes too (useful in EEG) * Standard setups for al the EEG caps |
* Set the amplitude scale for the time series * 2DLayout: Doesn't work when changing page => need refresh of GlobalData.Preferences.TopoLayoutOptions.TimeWindow * EEG reference/storage:<<BR>> * Intracranial electrodes: Define and display in the MRI viewer + 3D figures * Re-reference recordings * 2D topo plots: Display bipolar montages with electrodes positions = average of the two electrodes * Bad locations: indicate with [NaN NaN NaN] instead of [0 0 0] * Location of bipolar electrodes: Save the 2 positions in the channel file, with weights=[1 -1] * Bad channels that can be specified at the program level (for sites that have permanently bad channels) => AS Dubarry |
Line 49: | Line 42: |
* Allow to overwrite RAW files * Update file definition + events if time changes (ex: resample) * Make it work for the file formats == Other == * Command-line Brainstorm: for working on clusters (make sure that there are no interface interruptions) |
* 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: <<BR>> * 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 == Connectivity == * Inter-trial coherence * Figures: interaction with sensor selection / scouts selections * Display NxN as 1xN * PLV: Add a time integration * Correlation with delays? * Work on progress bars * Circle plot: Display Neuromag sensors == Processes == * Time-frequency: * Frequency bands: extended syntax (ex: [2 3 4], 10:5:90, ...) * Scouts values for timefreq on surfaces * 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 * When processing multiple files: waitbar is not behaving well * Improve interface for SSP on imported recordings * Define SSP from one time point * Show where the attenuation is projected:<<BR>>(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<<BR>>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 * Spectrum figure: use bst_selection == Database == |
Line 57: | Line 108: |
* Sources on surface: Display peak regions over time (time = color) => A.Gramfort * Images of amplitude: [sensor x time], [trial x time], scout: [trial x time] (similaire to erpimage in eeglab) * Simulation: synthesize pseudo data-files from a cortex patch (duration, amplitude, noise) |
* 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 |
Line 63: | Line 115: |
* Inverse: * Stable LCMV Beamformer * MUSIC * |
* 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 (contact Zainab Fatima): * Read CTF SAM .svl * Create new file type in the database * Display as layers in the MRI viewer |
Line 68: | Line 125: |
* Define as default * Check all the processes |
* Compute unconstrained and then project on the normal ? |
Line 71: | Line 127: |
* Stat and connectivity: what to do? (re-send email John+Sylvain) | |
Line 77: | Line 134: |
* Optimize: grid_interp_mri | |
Line 80: | Line 136: |
* Optimize bst_extrapm.m, add waitbar * Use the noise covariance from the database instead of recomputing it * 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 |
|
Line 87: | Line 137: |
* 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? |
|
Line 89: | Line 142: |
* 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? |
|
Line 91: | Line 149: |
* MRI import: auto-reorientation of MRI after selected NAS / LPA / RPA. * Major bug when importing surfaces for an MRI that was re-oriented manually * Finalize Brodmann scouts * Remove NCS/Talairach coordinate system, or fix it => Sylvain, Karim ND * Clustering cortex: Dimitrios, David, Yu-Teng == Processes == * Other processes: * Moving average * Remove linear trend * Power line removal * Bug: gradnorm crashes with bad channels * Spatial smoothing: check / document parameters * Sinusoid removal: fix new function * Time-frequency: * Write script for timefreq tutorial * Make much faster and more memory efficient (C functions coded by Matti ?) * Power spectrum: computation (FFT/welch, or average of TF) + display(f/Power, 2DLayout) * Display stat computed on time-frequency data * Display TF maps separately for the two gradiometers (if not: overlap) * Source reconstruction by frequency bands * Scouts on surface / time-freq * Process selection interface: * Bug when redimensioning window (with more than one process) * isAvgRef: warning quand process necessite des donnees en AVG REF en entree * Save "freqband" option when edited from custom processes * John's noise cancellation filters * Chain calls of similar processes (ex. bandpass+sin_remove), to avoid reading the full raw files several times |
* Scouts: * Mix constrained/unconstrained/volume sources, using the "Source model" atlas * 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? |
Line 132: | Line 173: |
== Display == * OpenGL options: {none, software, hardware} * Colormap: Set colormap max with right-click + move on the colorbar * Bug: Mixing 2 views MRI/3d, white cuts appear after !SetCurrentFigure * Bug: Scout without overlay, adapt scale for each graph when "Uniformize" 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) * Contact sheets & movies: use average of time windows instead of single instants, for each picture. == I/O == |
* t-test on volume sources * Create icons for Stat/PAC, Stat/Sprectrum, etc. * One sample t-test across subjects == Input / output == |
Line 145: | Line 179: |
* Nihon Kohden * EDF+ * !BrainVision / !BrainAmp: Get functions from EEGLAB |
|
Line 149: | Line 180: |
* NEUROFILE = COHERENCE EEG/video !LongTerm Monitoring => Manfred Spueler | |
Line 151: | Line 181: |
* Other file formats * MEG160 (KIT) * CTF: * Read continuous RAW files split in several .meg4 files ( > 2 Go) * Read STIM channel and generate !MarkerFile * EEGLAB: Apply ICA matrices, get number of trials for AVG files * !FieldTrip structures: In / Out * Gzipped Nifti (new BrainVISA standard) * Analyze / Nifti:<<BR>> * Fix output function (output MRIs cannot be imported in SPM) * Use correctly field "orient" to align automatically MRI when importing * Define scouts from SPM / Analyze 3D masks * Files > 2Gb: display warnings == Programming details == * Double-click doesn't work well on some Linux workstations * Bug: Menu "Use default EEG cap" doesn't work for a multiple selection (setting the same EEG cap for several subjects) * Bug node selection: click on sources > TF: select node-source, not node-condition * 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 * bst_warp and channel_project: Use tess_parametrize_new instead of tess_parametrize * Testing functions: test_ctf "folder" / test_neuromag "folder" / test "folder" => Alex * Bug in_bst_data_multi: If trials have different sizes, output is random (the one of the first file)... * sLORETA: Values are now multiplied by 1e12 at loading for display => has do to be done in another way * Shared kernels: do the "get bad channels" operation in a different way (reading all the files is too slow) * Write shepards.m with new algorithm for nearest neighbors * Use Matlab Coder to compile / optimize some processes |
* !FieldTrip structures: In / Out * BCI2000 Input (via EEGLAB plugin) * Export TF maps to SPM / volumes * EEGLAB import: Selection of conditions in script mode |
Line 179: | Line 187: |
* Version with big fonts for live demos | |
Line 181: | Line 188: |
* Send email to registered users to anounce major improvements * Script tutorials: * Update them to reflect all the recent changes * Script for the time-frequency computation |
|
Line 186: | Line 189: |
* Estimate time to complete each tutorial | * Processes: Describe all the processes |
Line 188: | Line 191: |
* Anatomy: Segmentation with !FreeSurfer | |
Line 191: | Line 193: |
* Coordinate sytems: How to convert between the different coordinates systems in scripts | |
Line 194: | Line 195: |
* Scouts: Atlases of Tzourio-Mazoyer and Brodman * Processes: Describe all the processes * Processes: How to write your own processes (user folder for processes) * Processes: Processing RAW files |
|
Line 199: | Line 196: |
* Temporary folder | |
Line 200: | Line 198: |
* MNE sample dataset | |
Line 202: | Line 199: |
* MRI segmentation with !FreeSurfer => David Wheland * How to make and compress a movie (Brainstorm + !VirtualDub + XVid) * Display the "What's new" page after downloading new version of brainstorm * Ask users to send their channel files, align on Colin, distribute |
== Current bugs == * 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) == Geeky programming details == * 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) |
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
- Full EEG/epilepsy online tutorial
- 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
- Allow multiple RAW windows(columns display)
- Screen setups: save configuration of windows
- 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"
- Documentation: Add definition of bad segments
- Set the amplitude scale for the time series
2DLayout: Doesn't work when changing page => need refresh of GlobalData.Preferences.TopoLayoutOptions.TimeWindow
EEG reference/storage:
- Intracranial electrodes: Define and display in the MRI viewer + 3D figures
- Re-reference recordings
- 2D topo plots: Display bipolar montages with electrodes positions = average of the two electrodes
- Bad locations: indicate with [NaN NaN NaN] instead of [0 0 0]
- Location of bipolar electrodes: Save the 2 positions in the channel file, with weights=[1 -1]
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
Connectivity
- Inter-trial coherence
- Figures: interaction with sensor selection / scouts selections
- Display NxN as 1xN
- PLV: Add a time integration
- Correlation with delays?
- Work on progress bars
- Circle plot: Display Neuromag sensors
Processes
- Time-frequency:
- Frequency bands: extended syntax (ex: [2 3 4], 10:5:90, ...)
- Scouts values for timefreq on surfaces
- 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
- When processing multiple files: waitbar is not behaving well
- Improve interface for SSP on imported recordings
- Define SSP from one time point
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
- Spectrum figure: use bst_selection
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
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 (contact Zainab Fatima):
- Read CTF SAM .svl
- Create new file type in the database
- 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?
Anatomy
- Scouts:
- Mix constrained/unconstrained/volume sources, using the "Source model" atlas
- 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
- 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
- Advanced tutorials:
- EEG (How to import an EEG cap)
Current bugs
- 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)
Geeky programming details
- 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)