Scalpel: Extracting neurons from calcium imaging data

Ashley Petersen, Noah Simon, Daniela Witten

Research output: Contribution to journalArticlepeer-review

16 Scopus citations


In the past few years, new technologies in the field of neuroscience have made it possible to simultaneously image activity in large populations of neurons at cellular resolution in behaving animals. In mid-2016, a huge repository of this so-called “calcium imaging” data was made publicly available. The availability of this large-scale data resource opens the door to a host of scientific questions for which new statistical methods must be developed. In this paper we consider the first step in the analysis of calcium imaging data—namely, identifying the neurons in a calcium imaging video. We propose a dictionary learning approach for this task. First, we perform image segmentation to develop a dictionary containing a huge number of candidate neurons. Next, we refine the dictionary using clustering. Finally, we apply the dictionary to select neurons and estimate their corresponding activity over time, using a sparse group lasso optimization problem. We assess performance on simulated calcium imaging data and apply our proposal to three calcium imaging data sets. Our proposed approach is implemented in the R package scalpel, which is available on CRAN.

Original languageEnglish (US)
Pages (from-to)2430-2456
Number of pages27
JournalAnnals of Applied Statistics
Issue number4
StatePublished - Dec 2018

Bibliographical note

Publisher Copyright:
© 2018, Institute of Mathematical Statistics. All rights reserved.


  • Calcium imaging
  • Cell sorting
  • Clustering
  • Dictionary learning
  • Neuron identification
  • Segmentation
  • Sparse group lasso


Dive into the research topics of 'Scalpel: Extracting neurons from calcium imaging data'. Together they form a unique fingerprint.

Cite this