- This project was part of the Reylab spikesorting pipeline.
The following is a processing pipeline for analyzing neural recordings using MATLAB and Python. Here are the requirements for each:
- MATLAB version 9.6 (R2019a) or later
- Signal Processing Toolbox
- DSP System Toolbox
- Parallel Computing Toolbox
- MATLAB Parallel Server
- Polyspace Bug Finder
- NPMK for Blackrock recordings
- Neuroshare for Ripple recordings
You can install the required Python packages using the following command:
pip install -r requirements.txt
-
Parsing the Recording
- For Ripple recordings, use the
parse_ripple.mscript. - For Blackrock recordings, use the
parse_NSx.mscript. - Note: The Neuroshare library is required for Ripple recordings.
- For Ripple recordings, use the
-
Analyzing Power Spectrum and Calculating Notches
- Use the
new_check_lfp_power_NSX.mscript to generate figures for analyzing the power spectrum and calculating notches. - If you prefer to use Python for this step, there is a function implemented in Python called
filter_freq_peaksthat performs notch filtering. If you used thenew_check_lfp_power_NSX.mscript, you can setload_mat_notches=Falsein the Python pipeline.
- Use the
-
Configuration Setup
- Before starting the pipeline, make sure to set up the
config.Yamlfile in theconfigfolder. - Update the paths in the configuration file to match the paths on your computer.
- Before starting the pipeline, make sure to set up the
-
Execution
- Open the
main.ipynbnotebook and follow the cells to execute the processing pipeline. - The notebook contains comments specifically for new Python users and provides details on how to use the
spikeinterfacelibrary.
- Open the
Notes We need to change the bandpass filter from butterworth to