3,495 research outputs found
Discrete Signal Processing on Graphs: Frequency Analysis
Signals and datasets that arise in physical and engineering applications, as
well as social, genetics, biomolecular, and many other domains, are becoming
increasingly larger and more complex. In contrast to traditional time and image
signals, data in these domains are supported by arbitrary graphs. Signal
processing on graphs extends concepts and techniques from traditional signal
processing to data indexed by generic graphs. This paper studies the concepts
of low and high frequencies on graphs, and low-, high-, and band-pass graph
filters. In traditional signal processing, there concepts are easily defined
because of a natural frequency ordering that has a physical interpretation. For
signals residing on graphs, in general, there is no obvious frequency ordering.
We propose a definition of total variation for graph signals that naturally
leads to a frequency ordering on graphs and defines low-, high-, and band-pass
graph signals and filters. We study the design of graph filters with specified
frequency response, and illustrate our approach with applications to sensor
malfunction detection and data classification
Greedy low-rank algorithm for spatial connectome regression
Recovering brain connectivity from tract tracing data is an important
computational problem in the neurosciences. Mesoscopic connectome
reconstruction was previously formulated as a structured matrix regression
problem (Harris et al., 2016), but existing techniques do not scale to the
whole-brain setting. The corresponding matrix equation is challenging to solve
due to large scale, ill-conditioning, and a general form that lacks a
convergent splitting. We propose a greedy low-rank algorithm for connectome
reconstruction problem in very high dimensions. The algorithm approximates the
solution by a sequence of rank-one updates which exploit the sparse and
positive definite problem structure. This algorithm was described previously
(Kressner and Sirkovi\'c, 2015) but never implemented for this connectome
problem, leading to a number of challenges. We have had to design judicious
stopping criteria and employ efficient solvers for the three main sub-problems
of the algorithm, including an efficient GPU implementation that alleviates the
main bottleneck for large datasets. The performance of the method is evaluated
on three examples: an artificial "toy" dataset and two whole-cortex instances
using data from the Allen Mouse Brain Connectivity Atlas. We find that the
method is significantly faster than previous methods and that moderate ranks
offer good approximation. This speedup allows for the estimation of
increasingly large-scale connectomes across taxa as these data become available
from tracing experiments. The data and code are available online
- …