[{"data":1,"prerenderedAt":-1},["ShallowReactive",2],{"skill:01M7EPK220KGEYRM6MYNETMP2A":3,"skill-same:01M7EPK220KGEYRM6MYNETMP2A":96},{"skillNo":4,"name":5,"nameZh":6,"description":7,"descriptionZh":6,"repositoryFullName":8,"ownerLogin":9,"ownerAvatarUrl":10,"ownerAccountNo":11,"ownerDisplayName":12,"repositoryNo":13,"starsCount":14,"licenseSpdxId":15,"licenseFileUrl":16,"sourceHtmlUrl":17,"categories":18,"keywords":19,"publishedAt":20,"linkOnly":21,"markdownContent":22,"markdownBody":23,"frontMatter":24,"refName":25,"commitSha":26,"skillDirectory":27,"forksCount":28,"repositoryHtmlUrl":29,"repositoryHomepageUrl":6,"ownerHtmlUrl":30,"ownerBlogUrl":31,"ownerDescription":6,"updatedAt":32,"install":33,"resources":46},"01M7EPK220KGEYRM6MYNETMP2A","neuropixels-analysis",null,"Neuropixels neural recording analysis. Load SpikeGLX\u002FOpenEphys data, preprocess, motion correction, Kilosort4 spike sorting, quality metrics, Allen\u002FIBL curation, AI-assisted visual analysis, for Neuropixels 1.0\u002F2.0 extracellular electrophysiology. Use when working with neural recordings, spike sorting, extracellular electrophysiology, or when the user mentions Neuropixels, SpikeGLX, Open Ephys, Kilosort, quality metrics, or unit curation.","TanNhatCMS\u002Fjarvis","TanNhatCMS","https:\u002F\u002Favatars.githubusercontent.com\u002Fu\u002F113796420?v=4","01M753MV17KG5CD3152G3WQEQY","KietNT","01M7EKEDYVYQ6PWVJJMZ911VMX",0,"MIT","https:\u002F\u002Fgithub.com\u002FTanNhatCMS\u002Fjarvis\u002Fblob\u002Fmain\u002Fskills\u002Fams-imp-neuropixels-analysis\u002FLICENSE.txt","https:\u002F\u002Fgithub.com\u002FTanNhatCMS\u002Fjarvis\u002Fblob\u002Fmain\u002Fskills\u002Fams-imp-neuropixels-analysis\u002FSKILL.md",[],[],"2026-11-15 14:19:17",false,"---\nname: neuropixels-analysis\ndescription: \"Neuropixels neural recording analysis. Load SpikeGLX\u002FOpenEphys data, preprocess, motion correction, Kilosort4 spike sorting, quality metrics, Allen\u002FIBL curation, AI-assisted visual analysis, for Neuropixels 1.0\u002F2.0 extracellular electrophysiology. Use when working with neural recordings, spike sorting, extracellular electrophysiology, or when the user mentions Neuropixels, SpikeGLX, Open Ephys, Kilosort, quality metrics, or unit curation.\"\n---\n\n# Neuropixels Data Analysis\n\n## Overview\n\nComprehensive toolkit for analyzing Neuropixels high-density neural recordings using current best practices from SpikeInterface, Allen Institute, and International Brain Laboratory (IBL). Supports the full workflow from raw data to publication-ready curated units.\n\n## When to Use This Skill\n\nThis skill should be used when:\n- Working with Neuropixels recordings (.ap.bin, .lf.bin, .meta files)\n- Loading data from SpikeGLX, Open Ephys, or NWB formats\n- Preprocessing neural recordings (filtering, CAR, bad channel detection)\n- Detecting and correcting motion\u002Fdrift in recordings\n- Running spike sorting (Kilosort4, SpykingCircus2, Mountainsort5)\n- Computing quality metrics (SNR, ISI violations, presence ratio)\n- Curating units using Allen\u002FIBL criteria\n- Creating visualizations of neural data\n- Exporting results to Phy or NWB\n\n## Supported Hardware & Formats\n\n| Probe | Electrodes | Channels | Notes |\n|-------|-----------|----------|-------|\n| Neuropixels 1.0 | 960 | 384 | Requires phase_shift correction |\n| Neuropixels 2.0 (single) | 1280 | 384 | Denser geometry |\n| Neuropixels 2.0 (4-shank) | 5120 | 384 | Multi-region recording |\n\n| Format | Extension | Reader |\n|--------|-----------|--------|\n| SpikeGLX | `.ap.bin`, `.lf.bin`, `.meta` | `si.read_spikeglx()` |\n| Open Ephys | `.continuous`, `.oebin` | `si.read_openephys()` |\n| NWB | `.nwb` | `si.read_nwb()` |\n\n## Quick Start\n\n### Basic Import and Setup\n\n```python\nimport spikeinterface.full as si\nimport neuropixels_analysis as npa\n\n# Configure parallel processing\njob_kwargs = dict(n_jobs=-1, chunk_duration='1s', progress_bar=True)\n```\n\n### Loading Data\n\n```python\n# SpikeGLX (most common)\nrecording = si.read_spikeglx('\u002Fpath\u002Fto\u002Fdata', stream_id='imec0.ap')\n\n# Open Ephys (common for many labs)\nrecording = si.read_openephys('\u002Fpath\u002Fto\u002FRecord_Node_101\u002F')\n\n# Check available streams\nstreams, ids = si.get_neo_streams('spikeglx', '\u002Fpath\u002Fto\u002Fdata')\nprint(streams)  # ['imec0.ap', 'imec0.lf', 'nidq']\n\n# For testing with subset of data\nrecording = recording.frame_slice(0, int(60 * recording.get_sampling_frequency()))\n```\n\n### Complete Pipeline (One Command)\n\n```python\n# Run full analysis pipeline\nresults = npa.run_pipeline(\n    recording,\n    output_dir='output\u002F',\n    sorter='kilosort4',\n    curation_method='allen',\n)\n\n# Access results\nsorting = results['sorting']\nmetrics = results['metrics']\nlabels = results['labels']\n```\n\n## Standard Analysis Workflow\n\n### 1. Preprocessing\n\n```python\n# Recommended preprocessing chain\nrec = si.highpass_filter(recording, freq_min=400)\nrec = si.phase_shift(rec)  # Required for Neuropixels 1.0\nbad_ids, _ = si.detect_bad_channels(rec)\nrec = rec.remove_channels(bad_ids)\nrec = si.common_reference(rec, operator='median')\n\n# Or use our wrapper\nrec = npa.preprocess(recording)\n```\n\n### 2. Check and Correct Drift\n\n```python\n# Check for drift (always do this!)\nmotion_info = npa.estimate_motion(rec, preset='kilosort_like')\nnpa.plot_drift(rec, motion_info, output='drift_map.png')\n\n# Apply correction if needed\nif motion_info['motion'].max() > 10:  # microns\n    rec = npa.correct_motion(rec, preset='nonrigid_accurate')\n```\n\n### 3. Spike Sorting\n\n```python\n# Kilosort4 (recommended, requires GPU)\nsorting = si.run_sorter('kilosort4', rec, folder='ks4_output')\n\n# CPU alternatives\nsorting = si.run_sorter('tridesclous2', rec, folder='tdc2_output')\nsorting = si.run_sorter('spykingcircus2', rec, folder='sc2_output')\nsorting = si.run_sorter('mountainsort5', rec, folder='ms5_output')\n\n# Check available sorters\nprint(si.installed_sorters())\n```\n\n### 4. Postprocessing\n\n```python\n# Create analyzer and compute all extensions\nanalyzer = si.create_sorting_analyzer(sorting, rec, sparse=True)\n\nanalyzer.compute('random_spikes', max_spikes_per_unit=500)\nanalyzer.compute('waveforms', ms_before=1.0, ms_after=2.0)\nanalyzer.compute('templates', operators=['average', 'std'])\nanalyzer.compute('spike_amplitudes')\nanalyzer.compute('correlograms', window_ms=50.0, bin_ms=1.0)\nanalyzer.compute('unit_locations', method='monopolar_triangulation')\nanalyzer.compute('quality_metrics')\n\nmetrics = analyzer.get_extension('quality_metrics').get_data()\n```\n\n### 5. Curation\n\n```python\n# Allen Institute criteria (conservative)\ngood_units = metrics.query(\"\"\"\n    presence_ratio > 0.9 and\n    isi_violations_ratio \u003C 0.5 and\n    amplitude_cutoff \u003C 0.1\n\"\"\").index.tolist()\n\n# Or use automated curation\nlabels = npa.curate(metrics, method='allen')  # 'allen', 'ibl', 'strict'\n```\n\n### 6. AI-Assisted Curation (For Uncertain Units)\n\nWhen using this skill with Claude Code, Claude can directly analyze waveform plots and provide expert curation decisions. For programmatic API access:\n\n```python\nfrom anthropic import Anthropic\n\n# Setup API client\nclient = Anthropic()\n\n# Analyze uncertain units visually\nuncertain = metrics.query('snr > 3 and snr \u003C 8').index.tolist()\n\nfor unit_id in uncertain:\n    result = npa.analyze_unit_visually(analyzer, unit_id, api_client=client)\n    print(f\"Unit {unit_id}: {result['classification']}\")\n    print(f\"  Reasoning: {result['reasoning'][:100]}...\")\n```\n\n**Claude Code Integration**: When running within Claude Code, ask Claude to examine waveform\u002Fcorrelogram plots directly - no API setup required.\n\n### 7. Generate Analysis Report\n\n```python\n# Generate comprehensive HTML report with visualizations\nreport_dir = npa.generate_analysis_report(results, 'output\u002F')\n# Opens report.html with summary stats, figures, and unit table\n\n# Print formatted summary to console\nnpa.print_analysis_summary(results)\n```\n\n### 8. Export Results\n\n```python\n# Export to Phy for manual review\nsi.export_to_phy(analyzer, output_folder='phy_export\u002F',\n                 compute_pc_features=True, compute_amplitudes=True)\n\n# Export to NWB\nfrom spikeinterface.exporters import export_to_nwb\nexport_to_nwb(rec, sorting, 'output.nwb')\n\n# Save quality metrics\nmetrics.to_csv('quality_metrics.csv')\n```\n\n## Common Pitfalls and Best Practices\n\n1. **Always check drift** before spike sorting - drift > 10μm significantly impacts quality\n2. **Use phase_shift** for Neuropixels 1.0 probes (not needed for 2.0)\n3. **Save preprocessed data** to avoid recomputing - use `rec.save(folder='preprocessed\u002F')`\n4. **Use GPU** for Kilosort4 - it's 10-50x faster than CPU alternatives\n5. **Review uncertain units manually** - automated curation is a starting point\n6. **Combine metrics with AI** - use metrics for clear cases, AI for borderline units\n7. **Document your thresholds** - different analyses may need different criteria\n8. **Export to Phy** for critical experiments - human oversight is valuable\n\n## Key Parameters to Adjust\n\n### Preprocessing\n- `freq_min`: Highpass cutoff (300-400 Hz typical)\n- `detect_threshold`: Bad channel detection sensitivity\n\n### Motion Correction\n- `preset`: 'kilosort_like' (fast) or 'nonrigid_accurate' (better for severe drift)\n\n### Spike Sorting (Kilosort4)\n- `batch_size`: Samples per batch (30000 default)\n- `nblocks`: Number of drift blocks (increase for long recordings)\n- `Th_learned`: Detection threshold (lower = more spikes)\n\n### Quality Metrics\n- `snr_threshold`: Signal-to-noise cutoff (3-5 typical)\n- `isi_violations_ratio`: Refractory violations (0.01-0.5)\n- `presence_ratio`: Recording coverage (0.5-0.95)\n\n## Bundled Resources\n\n### scripts\u002Fpreprocess_recording.py\nAutomated preprocessing script:\n```bash\npython scripts\u002Fpreprocess_recording.py \u002Fpath\u002Fto\u002Fdata --output preprocessed\u002F\n```\n\n### scripts\u002Frun_sorting.py\nRun spike sorting:\n```bash\npython scripts\u002Frun_sorting.py preprocessed\u002F --sorter kilosort4 --output sorting\u002F\n```\n\n### scripts\u002Fcompute_metrics.py\nCompute quality metrics and apply curation:\n```bash\npython scripts\u002Fcompute_metrics.py sorting\u002F preprocessed\u002F --output metrics\u002F --curation allen\n```\n\n### scripts\u002Fexport_to_phy.py\nExport to Phy for manual curation:\n```bash\npython scripts\u002Fexport_to_phy.py metrics\u002Fanalyzer --output phy_export\u002F\n```\n\n### assets\u002Fanalysis_template.py\nComplete analysis template. Copy and customize:\n```bash\ncp assets\u002Fanalysis_template.py my_analysis.py\n# Edit parameters and run\npython my_analysis.py\n```\n\n### reference\u002Fstandard_workflow.md\nDetailed step-by-step workflow with explanations for each stage.\n\n### reference\u002Fapi_reference.md\nQuick function reference organized by module.\n\n### reference\u002Fplotting_guide.md\nComprehensive visualization guide for publication-quality figures.\n\n## Detailed Reference Guides\n\n| Topic | Reference |\n|-------|-----------|\n| Full workflow | [reference\u002Fstandard_workflow.md](reference\u002Fstandard_workflow.md) |\n| API reference | [reference\u002Fapi_reference.md](reference\u002Fapi_reference.md) |\n| Plotting guide | [reference\u002Fplotting_guide.md](reference\u002Fplotting_guide.md) |\n| Preprocessing | [PREPROCESSING.md](PREPROCESSING.md) |\n| Spike sorting | [SPIKE_SORTING.md](SPIKE_SORTING.md) |\n| Motion correction | [MOTION_CORRECTION.md](MOTION_CORRECTION.md) |\n| Quality metrics | [QUALITY_METRICS.md](QUALITY_METRICS.md) |\n| Automated curation | [AUTOMATED_CURATION.md](AUTOMATED_CURATION.md) |\n| AI-assisted curation | [AI_CURATION.md](AI_CURATION.md) |\n| Waveform analysis | [ANALYSIS.md](ANALYSIS.md) |\n\n## Installation\n\n```bash\n# Core packages\npip install spikeinterface[full] probeinterface neo\n\n# Spike sorters\npip install kilosort          # Kilosort4 (GPU required)\npip install spykingcircus     # SpykingCircus2 (CPU)\npip install mountainsort5     # Mountainsort5 (CPU)\n\n# Our toolkit\npip install neuropixels-analysis\n\n# Optional: AI curation\npip install anthropic\n\n# Optional: IBL tools\npip install ibl-neuropixel ibllib\n```\n\n## Project Structure\n\n```\nproject\u002F\n├── raw_data\u002F\n│   └── recording_g0\u002F\n│       └── recording_g0_imec0\u002F\n│           ├── recording_g0_t0.imec0.ap.bin\n│           └── recording_g0_t0.imec0.ap.meta\n├── preprocessed\u002F           # Saved preprocessed recording\n├── motion\u002F                 # Motion estimation results\n├── sorting_output\u002F         # Spike sorter output\n├── analyzer\u002F               # SortingAnalyzer (waveforms, metrics)\n├── phy_export\u002F             # For manual curation\n├── ai_curation\u002F            # AI analysis reports\n└── results\u002F\n    ├── quality_metrics.csv\n    ├── curation_labels.json\n    └── output.nwb\n```\n\n## Additional Resources\n\n- **SpikeInterface Docs**: https:\u002F\u002Fspikeinterface.readthedocs.io\u002F\n- **Neuropixels Tutorial**: https:\u002F\u002Fspikeinterface.readthedocs.io\u002Fen\u002Fstable\u002Fhow_to\u002Fanalyze_neuropixels.html\n- **Kilosort4 GitHub**: https:\u002F\u002Fgithub.com\u002FMouseLand\u002FKilosort\n- **IBL Neuropixel Tools**: https:\u002F\u002Fgithub.com\u002Fint-brain-lab\u002Fibl-neuropixel\n- **Allen Institute ecephys**: https:\u002F\u002Fgithub.com\u002FAllenInstitute\u002Fecephys_spike_sorting\n- **Bombcell (Automated QC)**: https:\u002F\u002Fgithub.com\u002FJulie-Fabre\u002Fbombcell\n- **SpikeAgent (AI Curation)**: https:\u002F\u002Fgithub.com\u002FSpikeAgent\u002FSpikeAgent\n","# Neuropixels Data Analysis\n\n## Overview\n\nComprehensive toolkit for analyzing Neuropixels high-density neural recordings using current best practices from SpikeInterface, Allen Institute, and International Brain Laboratory (IBL). Supports the full workflow from raw data to publication-ready curated units.\n\n## When to Use This Skill\n\nThis skill should be used when:\n- Working with Neuropixels recordings (.ap.bin, .lf.bin, .meta files)\n- Loading data from SpikeGLX, Open Ephys, or NWB formats\n- Preprocessing neural recordings (filtering, CAR, bad channel detection)\n- Detecting and correcting motion\u002Fdrift in recordings\n- Running spike sorting (Kilosort4, SpykingCircus2, Mountainsort5)\n- Computing quality metrics (SNR, ISI violations, presence ratio)\n- Curating units using Allen\u002FIBL criteria\n- Creating visualizations of neural data\n- Exporting results to Phy or NWB\n\n## Supported Hardware & Formats\n\n| Probe | Electrodes | Channels | Notes |\n|-------|-----------|----------|-------|\n| Neuropixels 1.0 | 960 | 384 | Requires phase_shift correction |\n| Neuropixels 2.0 (single) | 1280 | 384 | Denser geometry |\n| Neuropixels 2.0 (4-shank) | 5120 | 384 | Multi-region recording |\n\n| Format | Extension | Reader |\n|--------|-----------|--------|\n| SpikeGLX | `.ap.bin`, `.lf.bin`, `.meta` | `si.read_spikeglx()` |\n| Open Ephys | `.continuous`, `.oebin` | `si.read_openephys()` |\n| NWB | `.nwb` | `si.read_nwb()` |\n\n## Quick Start\n\n### Basic Import and Setup\n\n```python\nimport spikeinterface.full as si\nimport neuropixels_analysis as npa\n\n# Configure parallel processing\njob_kwargs = dict(n_jobs=-1, chunk_duration='1s', progress_bar=True)\n```\n\n### Loading Data\n\n```python\n# SpikeGLX (most common)\nrecording = si.read_spikeglx('\u002Fpath\u002Fto\u002Fdata', stream_id='imec0.ap')\n\n# Open Ephys (common for many labs)\nrecording = si.read_openephys('\u002Fpath\u002Fto\u002FRecord_Node_101\u002F')\n\n# Check available streams\nstreams, ids = si.get_neo_streams('spikeglx', '\u002Fpath\u002Fto\u002Fdata')\nprint(streams)  # ['imec0.ap', 'imec0.lf', 'nidq']\n\n# For testing with subset of data\nrecording = recording.frame_slice(0, int(60 * recording.get_sampling_frequency()))\n```\n\n### Complete Pipeline (One Command)\n\n```python\n# Run full analysis pipeline\nresults = npa.run_pipeline(\n    recording,\n    output_dir='output\u002F',\n    sorter='kilosort4',\n    curation_method='allen',\n)\n\n# Access results\nsorting = results['sorting']\nmetrics = results['metrics']\nlabels = results['labels']\n```\n\n## Standard Analysis Workflow\n\n### 1. Preprocessing\n\n```python\n# Recommended preprocessing chain\nrec = si.highpass_filter(recording, freq_min=400)\nrec = si.phase_shift(rec)  # Required for Neuropixels 1.0\nbad_ids, _ = si.detect_bad_channels(rec)\nrec = rec.remove_channels(bad_ids)\nrec = si.common_reference(rec, operator='median')\n\n# Or use our wrapper\nrec = npa.preprocess(recording)\n```\n\n### 2. Check and Correct Drift\n\n```python\n# Check for drift (always do this!)\nmotion_info = npa.estimate_motion(rec, preset='kilosort_like')\nnpa.plot_drift(rec, motion_info, output='drift_map.png')\n\n# Apply correction if needed\nif motion_info['motion'].max() > 10:  # microns\n    rec = npa.correct_motion(rec, preset='nonrigid_accurate')\n```\n\n### 3. Spike Sorting\n\n```python\n# Kilosort4 (recommended, requires GPU)\nsorting = si.run_sorter('kilosort4', rec, folder='ks4_output')\n\n# CPU alternatives\nsorting = si.run_sorter('tridesclous2', rec, folder='tdc2_output')\nsorting = si.run_sorter('spykingcircus2', rec, folder='sc2_output')\nsorting = si.run_sorter('mountainsort5', rec, folder='ms5_output')\n\n# Check available sorters\nprint(si.installed_sorters())\n```\n\n### 4. Postprocessing\n\n```python\n# Create analyzer and compute all extensions\nanalyzer = si.create_sorting_analyzer(sorting, rec, sparse=True)\n\nanalyzer.compute('random_spikes', max_spikes_per_unit=500)\nanalyzer.compute('waveforms', ms_before=1.0, ms_after=2.0)\nanalyzer.compute('templates', operators=['average', 'std'])\nanalyzer.compute('spike_amplitudes')\nanalyzer.compute('correlograms', window_ms=50.0, bin_ms=1.0)\nanalyzer.compute('unit_locations', method='monopolar_triangulation')\nanalyzer.compute('quality_metrics')\n\nmetrics = analyzer.get_extension('quality_metrics').get_data()\n```\n\n### 5. Curation\n\n```python\n# Allen Institute criteria (conservative)\ngood_units = metrics.query(\"\"\"\n    presence_ratio > 0.9 and\n    isi_violations_ratio \u003C 0.5 and\n    amplitude_cutoff \u003C 0.1\n\"\"\").index.tolist()\n\n# Or use automated curation\nlabels = npa.curate(metrics, method='allen')  # 'allen', 'ibl', 'strict'\n```\n\n### 6. AI-Assisted Curation (For Uncertain Units)\n\nWhen using this skill with Claude Code, Claude can directly analyze waveform plots and provide expert curation decisions. For programmatic API access:\n\n```python\nfrom anthropic import Anthropic\n\n# Setup API client\nclient = Anthropic()\n\n# Analyze uncertain units visually\nuncertain = metrics.query('snr > 3 and snr \u003C 8').index.tolist()\n\nfor unit_id in uncertain:\n    result = npa.analyze_unit_visually(analyzer, unit_id, api_client=client)\n    print(f\"Unit {unit_id}: {result['classification']}\")\n    print(f\"  Reasoning: {result['reasoning'][:100]}...\")\n```\n\n**Claude Code Integration**: When running within Claude Code, ask Claude to examine waveform\u002Fcorrelogram plots directly - no API setup required.\n\n### 7. Generate Analysis Report\n\n```python\n# Generate comprehensive HTML report with visualizations\nreport_dir = npa.generate_analysis_report(results, 'output\u002F')\n# Opens report.html with summary stats, figures, and unit table\n\n# Print formatted summary to console\nnpa.print_analysis_summary(results)\n```\n\n### 8. Export Results\n\n```python\n# Export to Phy for manual review\nsi.export_to_phy(analyzer, output_folder='phy_export\u002F',\n                 compute_pc_features=True, compute_amplitudes=True)\n\n# Export to NWB\nfrom spikeinterface.exporters import export_to_nwb\nexport_to_nwb(rec, sorting, 'output.nwb')\n\n# Save quality metrics\nmetrics.to_csv('quality_metrics.csv')\n```\n\n## Common Pitfalls and Best Practices\n\n1. **Always check drift** before spike sorting - drift > 10μm significantly impacts quality\n2. **Use phase_shift** for Neuropixels 1.0 probes (not needed for 2.0)\n3. **Save preprocessed data** to avoid recomputing - use `rec.save(folder='preprocessed\u002F')`\n4. **Use GPU** for Kilosort4 - it's 10-50x faster than CPU alternatives\n5. **Review uncertain units manually** - automated curation is a starting point\n6. **Combine metrics with AI** - use metrics for clear cases, AI for borderline units\n7. **Document your thresholds** - different analyses may need different criteria\n8. **Export to Phy** for critical experiments - human oversight is valuable\n\n## Key Parameters to Adjust\n\n### Preprocessing\n- `freq_min`: Highpass cutoff (300-400 Hz typical)\n- `detect_threshold`: Bad channel detection sensitivity\n\n### Motion Correction\n- `preset`: 'kilosort_like' (fast) or 'nonrigid_accurate' (better for severe drift)\n\n### Spike Sorting (Kilosort4)\n- `batch_size`: Samples per batch (30000 default)\n- `nblocks`: Number of drift blocks (increase for long recordings)\n- `Th_learned`: Detection threshold (lower = more spikes)\n\n### Quality Metrics\n- `snr_threshold`: Signal-to-noise cutoff (3-5 typical)\n- `isi_violations_ratio`: Refractory violations (0.01-0.5)\n- `presence_ratio`: Recording coverage (0.5-0.95)\n\n## Bundled Resources\n\n### scripts\u002Fpreprocess_recording.py\nAutomated preprocessing script:\n```bash\npython scripts\u002Fpreprocess_recording.py \u002Fpath\u002Fto\u002Fdata --output preprocessed\u002F\n```\n\n### scripts\u002Frun_sorting.py\nRun spike sorting:\n```bash\npython scripts\u002Frun_sorting.py preprocessed\u002F --sorter kilosort4 --output sorting\u002F\n```\n\n### scripts\u002Fcompute_metrics.py\nCompute quality metrics and apply curation:\n```bash\npython scripts\u002Fcompute_metrics.py sorting\u002F preprocessed\u002F --output metrics\u002F --curation allen\n```\n\n### scripts\u002Fexport_to_phy.py\nExport to Phy for manual curation:\n```bash\npython scripts\u002Fexport_to_phy.py metrics\u002Fanalyzer --output phy_export\u002F\n```\n\n### assets\u002Fanalysis_template.py\nComplete analysis template. Copy and customize:\n```bash\ncp assets\u002Fanalysis_template.py my_analysis.py\n# Edit parameters and run\npython my_analysis.py\n```\n\n### reference\u002Fstandard_workflow.md\nDetailed step-by-step workflow with explanations for each stage.\n\n### reference\u002Fapi_reference.md\nQuick function reference organized by module.\n\n### reference\u002Fplotting_guide.md\nComprehensive visualization guide for publication-quality figures.\n\n## Detailed Reference Guides\n\n| Topic | Reference |\n|-------|-----------|\n| Full workflow | [reference\u002Fstandard_workflow.md](reference\u002Fstandard_workflow.md) |\n| API reference | [reference\u002Fapi_reference.md](reference\u002Fapi_reference.md) |\n| Plotting guide | [reference\u002Fplotting_guide.md](reference\u002Fplotting_guide.md) |\n| Preprocessing | [PREPROCESSING.md](PREPROCESSING.md) |\n| Spike sorting | [SPIKE_SORTING.md](SPIKE_SORTING.md) |\n| Motion correction | [MOTION_CORRECTION.md](MOTION_CORRECTION.md) |\n| Quality metrics | [QUALITY_METRICS.md](QUALITY_METRICS.md) |\n| Automated curation | [AUTOMATED_CURATION.md](AUTOMATED_CURATION.md) |\n| AI-assisted curation | [AI_CURATION.md](AI_CURATION.md) |\n| Waveform analysis | [ANALYSIS.md](ANALYSIS.md) |\n\n## Installation\n\n```bash\n# Core packages\npip install spikeinterface[full] probeinterface neo\n\n# Spike sorters\npip install kilosort          # Kilosort4 (GPU required)\npip install spykingcircus     # SpykingCircus2 (CPU)\npip install mountainsort5     # Mountainsort5 (CPU)\n\n# Our toolkit\npip install neuropixels-analysis\n\n# Optional: AI curation\npip install anthropic\n\n# Optional: IBL tools\npip install ibl-neuropixel ibllib\n```\n\n## Project Structure\n\n```\nproject\u002F\n├── raw_data\u002F\n│   └── recording_g0\u002F\n│       └── recording_g0_imec0\u002F\n│           ├── recording_g0_t0.imec0.ap.bin\n│           └── recording_g0_t0.imec0.ap.meta\n├── preprocessed\u002F           # Saved preprocessed recording\n├── motion\u002F                 # Motion estimation results\n├── sorting_output\u002F         # Spike sorter output\n├── analyzer\u002F               # SortingAnalyzer (waveforms, metrics)\n├── phy_export\u002F             # For manual curation\n├── ai_curation\u002F            # AI analysis reports\n└── results\u002F\n    ├── quality_metrics.csv\n    ├── curation_labels.json\n    └── output.nwb\n```\n\n## Additional Resources\n\n- **SpikeInterface Docs**: https:\u002F\u002Fspikeinterface.readthedocs.io\u002F\n- **Neuropixels Tutorial**: https:\u002F\u002Fspikeinterface.readthedocs.io\u002Fen\u002Fstable\u002Fhow_to\u002Fanalyze_neuropixels.html\n- **Kilosort4 GitHub**: https:\u002F\u002Fgithub.com\u002FMouseLand\u002FKilosort\n- **IBL Neuropixel Tools**: https:\u002F\u002Fgithub.com\u002Fint-brain-lab\u002Fibl-neuropixel\n- **Allen Institute ecephys**: https:\u002F\u002Fgithub.com\u002FAllenInstitute\u002Fecephys_spike_sorting\n- **Bombcell (Automated QC)**: https:\u002F\u002Fgithub.com\u002FJulie-Fabre\u002Fbombcell\n- **SpikeAgent (AI Curation)**: https:\u002F\u002Fgithub.com\u002FSpikeAgent\u002FSpikeAgent\n",{"name":5,"description":7},"main","af388a77c2fe5e0b4fd879cca4300a0ffa9a6579","skills\u002Fams-imp-neuropixels-analysis",1,"https:\u002F\u002Fgithub.com\u002FTanNhatCMS\u002Fjarvis","https:\u002F\u002Fgithub.com\u002FTanNhatCMS","https:\u002F\u002Fmrkiet.name.vn","2026-11-07 14:23:08",{"source":8,"skillName":5,"commands":34,"prompt":44,"zipAvailable":45},[35,38,41],{"tool":36,"command":37},"npx","npx skills add TanNhatCMS\u002Fjarvis --skill neuropixels-analysis",{"tool":39,"command":40},"bunx","bunx skills add TanNhatCMS\u002Fjarvis --skill neuropixels-analysis",{"tool":42,"command":43},"pnpm","pnpm dlx skills add TanNhatCMS\u002Fjarvis --skill neuropixels-analysis","安装技能，技能地址：https:\u002F\u002Fgithub.com\u002FTanNhatCMS\u002Fjarvis\u002Fblob\u002Fmain\u002Fskills\u002Fams-imp-neuropixels-analysis\u002FSKILL.md",true,[47,53,58,64,68,72,76,80,84,88,92],{"path":48,"resourceType":49,"resourceTypeDesc":50,"sizeBytes":51,"sourceHtmlUrl":52},"assets\u002Fanalysis_template.py","SCRIPT","脚本",9014,"https:\u002F\u002Fgithub.com\u002FTanNhatCMS\u002Fjarvis\u002Fblob\u002Fmain\u002Fskills\u002Fams-imp-neuropixels-analysis\u002Fassets\u002Fanalysis_template.py",{"path":54,"resourceType":55,"resourceTypeDesc":56,"sizeBytes":57,"sourceHtmlUrl":16},"LICENSE.txt","LICENSE","许可证",1065,{"path":59,"resourceType":60,"resourceTypeDesc":61,"sizeBytes":62,"sourceHtmlUrl":63},"references\u002Fapi_reference.md","REFERENCE","参考文档",8175,"https:\u002F\u002Fgithub.com\u002FTanNhatCMS\u002Fjarvis\u002Fblob\u002Fmain\u002Fskills\u002Fams-imp-neuropixels-analysis\u002Freferences\u002Fapi_reference.md",{"path":65,"resourceType":60,"resourceTypeDesc":61,"sizeBytes":66,"sourceHtmlUrl":67},"references\u002Fplotting_guide.md",10837,"https:\u002F\u002Fgithub.com\u002FTanNhatCMS\u002Fjarvis\u002Fblob\u002Fmain\u002Fskills\u002Fams-imp-neuropixels-analysis\u002Freferences\u002Fplotting_guide.md",{"path":69,"resourceType":60,"resourceTypeDesc":61,"sizeBytes":70,"sourceHtmlUrl":71},"references\u002Fstandard_workflow.md",8646,"https:\u002F\u002Fgithub.com\u002FTanNhatCMS\u002Fjarvis\u002Fblob\u002Fmain\u002Fskills\u002Fams-imp-neuropixels-analysis\u002Freferences\u002Fstandard_workflow.md",{"path":73,"resourceType":49,"resourceTypeDesc":50,"sizeBytes":74,"sourceHtmlUrl":75},"scripts\u002Fcompute_metrics.py",5226,"https:\u002F\u002Fgithub.com\u002FTanNhatCMS\u002Fjarvis\u002Fblob\u002Fmain\u002Fskills\u002Fams-imp-neuropixels-analysis\u002Fscripts\u002Fcompute_metrics.py",{"path":77,"resourceType":49,"resourceTypeDesc":50,"sizeBytes":78,"sourceHtmlUrl":79},"scripts\u002Fexplore_recording.py",5529,"https:\u002F\u002Fgithub.com\u002FTanNhatCMS\u002Fjarvis\u002Fblob\u002Fmain\u002Fskills\u002Fams-imp-neuropixels-analysis\u002Fscripts\u002Fexplore_recording.py",{"path":81,"resourceType":49,"resourceTypeDesc":50,"sizeBytes":82,"sourceHtmlUrl":83},"scripts\u002Fexport_to_phy.py",2477,"https:\u002F\u002Fgithub.com\u002FTanNhatCMS\u002Fjarvis\u002Fblob\u002Fmain\u002Fskills\u002Fams-imp-neuropixels-analysis\u002Fscripts\u002Fexport_to_phy.py",{"path":85,"resourceType":49,"resourceTypeDesc":50,"sizeBytes":86,"sourceHtmlUrl":87},"scripts\u002Fneuropixels_pipeline.py",13286,"https:\u002F\u002Fgithub.com\u002FTanNhatCMS\u002Fjarvis\u002Fblob\u002Fmain\u002Fskills\u002Fams-imp-neuropixels-analysis\u002Fscripts\u002Fneuropixels_pipeline.py",{"path":89,"resourceType":49,"resourceTypeDesc":50,"sizeBytes":90,"sourceHtmlUrl":91},"scripts\u002Fpreprocess_recording.py",4235,"https:\u002F\u002Fgithub.com\u002FTanNhatCMS\u002Fjarvis\u002Fblob\u002Fmain\u002Fskills\u002Fams-imp-neuropixels-analysis\u002Fscripts\u002Fpreprocess_recording.py",{"path":93,"resourceType":49,"resourceTypeDesc":50,"sizeBytes":94,"sourceHtmlUrl":95},"scripts\u002Frun_sorting.py",2573,"https:\u002F\u002Fgithub.com\u002FTanNhatCMS\u002Fjarvis\u002Fblob\u002Fmain\u002Fskills\u002Fams-imp-neuropixels-analysis\u002Fscripts\u002Frun_sorting.py",[97,104,112,119,126,133],{"skillNo":98,"name":99,"nameZh":6,"description":100,"descriptionZh":6,"repositoryFullName":8,"ownerLogin":9,"ownerAvatarUrl":10,"ownerAccountNo":11,"ownerDisplayName":12,"repositoryNo":13,"starsCount":14,"licenseSpdxId":15,"licenseFileUrl":6,"sourceHtmlUrl":101,"categories":102,"keywords":103,"publishedAt":20},"01M7EWYJR0J8DARY7E59Z6Y9MV","web-scrape","Turn web pages into LLM-ready markdown via the Firecrawl API (renders JS, anti-bot handling). Use when browser_open\u002Ffetch fails on JS-heavy pages, or when you need clean markdown\u002Fstructured extraction from a site. Requires FIRECRAWL_API_KEY.","https:\u002F\u002Fgithub.com\u002FTanNhatCMS\u002Fjarvis\u002Fblob\u002Fmain\u002Fskills\u002Fweb-scrape\u002FSKILL.md",[],[],{"skillNo":105,"name":106,"nameZh":6,"description":107,"descriptionZh":6,"repositoryFullName":8,"ownerLogin":9,"ownerAvatarUrl":10,"ownerAccountNo":11,"ownerDisplayName":12,"repositoryNo":13,"starsCount":14,"licenseSpdxId":15,"licenseFileUrl":108,"sourceHtmlUrl":109,"categories":110,"keywords":111,"publishedAt":20},"01M7EWYBN5P6F2AJG3FR24T9BP","design-taste-frontend","tasteskill v2, the anti-slop frontend skill. Art direction for websites, distilled from a study of 247 curated gallery sites (One Page Love, Awwwards Sites of the Day). Deals researched directions (idea, composition, type, palette, photographs, section blueprints) from a library, enforces a strict no-slop house style with a real-browser inspector, and runs a look-and-fix loop on screenshots. Use this whenever you build, design, redesign, restyle, polish or review a landing page, marketing site, portfolio, product page, one-pager, hero or any public-facing web page, in any stack (HTML\u002FCSS, React, Next.js, Vue, Svelte, Astro, Tailwind) — also when the user says a page looks generic, templated, \"AI-made\", sloppy, or asks to make it look designed, premium, award-level or less like every other site. Not for dashboards or app screens.","https:\u002F\u002Fgithub.com\u002FTanNhatCMS\u002Fjarvis\u002Fblob\u002Fmain\u002Fskills\u002Ftaste-design\u002FLICENSE","https:\u002F\u002Fgithub.com\u002FTanNhatCMS\u002Fjarvis\u002Fblob\u002Fmain\u002Fskills\u002Ftaste-design\u002FSKILL.md",[],[],{"skillNo":113,"name":114,"nameZh":6,"description":115,"descriptionZh":6,"repositoryFullName":8,"ownerLogin":9,"ownerAvatarUrl":10,"ownerAccountNo":11,"ownerDisplayName":12,"repositoryNo":13,"starsCount":14,"licenseSpdxId":15,"licenseFileUrl":6,"sourceHtmlUrl":116,"categories":117,"keywords":118,"publishedAt":20},"01M7EWY8FQFNQ2GB7NZZWGSSDX","test-driven-development","Use when implementing any feature or bugfix, before writing implementation code","https:\u002F\u002Fgithub.com\u002FTanNhatCMS\u002Fjarvis\u002Fblob\u002Fmain\u002Fskills\u002Fsuperpowers-tdd\u002FSKILL.md",[],[],{"skillNo":120,"name":121,"nameZh":6,"description":122,"descriptionZh":6,"repositoryFullName":8,"ownerLogin":9,"ownerAvatarUrl":10,"ownerAccountNo":11,"ownerDisplayName":12,"repositoryNo":13,"starsCount":14,"licenseSpdxId":15,"licenseFileUrl":6,"sourceHtmlUrl":123,"categories":124,"keywords":125,"publishedAt":20},"01M7EWY5SHBFRH1HDHVWM9K7TV","systematic-debugging","Use when encountering any bug, test failure, or unexpected behavior, before proposing fixes","https:\u002F\u002Fgithub.com\u002FTanNhatCMS\u002Fjarvis\u002Fblob\u002Fmain\u002Fskills\u002Fsuperpowers-debugging\u002FSKILL.md",[],[],{"skillNo":127,"name":128,"nameZh":6,"description":129,"descriptionZh":6,"repositoryFullName":8,"ownerLogin":9,"ownerAvatarUrl":10,"ownerAccountNo":11,"ownerDisplayName":12,"repositoryNo":13,"starsCount":14,"licenseSpdxId":15,"licenseFileUrl":6,"sourceHtmlUrl":130,"categories":131,"keywords":132,"publishedAt":20},"01M7EWY3NM0Z0JAKS7HDZ24K66","requesting-code-review","Use when completing tasks, implementing major features, or before merging to verify work meets requirements","https:\u002F\u002Fgithub.com\u002FTanNhatCMS\u002Fjarvis\u002Fblob\u002Fmain\u002Fskills\u002Fsuperpowers-code-review\u002FSKILL.md",[],[],{"skillNo":134,"name":135,"nameZh":6,"description":136,"descriptionZh":6,"repositoryFullName":8,"ownerLogin":9,"ownerAvatarUrl":10,"ownerAccountNo":11,"ownerDisplayName":12,"repositoryNo":13,"starsCount":14,"licenseSpdxId":15,"licenseFileUrl":6,"sourceHtmlUrl":137,"categories":138,"keywords":139,"publishedAt":20},"01M7EWY0WRMWX4M9ZFY0WKC593","brainstorming","You MUST use this before any creative work - creating features, building components, adding functionality, modifying behavior, or planning anything new, in software or out of it (a talk, a business, a renovation).","https:\u002F\u002Fgithub.com\u002FTanNhatCMS\u002Fjarvis\u002Fblob\u002Fmain\u002Fskills\u002Fsuperpowers-brainstorming\u002FSKILL.md",[],[]]