Parallel Discrete Convolutions on Adaptive Particle Representations of Images

Research output: Contribution to journalResearch articleContributedpeer-review

Contributors

  • Joel Jonsson - , Chair of Scientific Computing for Systems Biology, Max Planck Institute of Molecular Cell Biology and Genetics, Center for Systems Biology Dresden (CSBD), Center for Scalable Data Analytics and Artificial Intelligence (ScaDS.AI) Dresden/Leipzig (Author)
  • Bevan L. Cheeseman - , Chair of Scientific Computing for Systems Biology, Max Planck Institute of Molecular Cell Biology and Genetics, Center for Systems Biology Dresden (CSBD), Oxford Nanoimaging (ONI) (Author)
  • Suryanarayana Maddu - , Chair of Scientific Computing for Systems Biology, Max Planck Institute of Molecular Cell Biology and Genetics, Center for Systems Biology Dresden (CSBD), Center for Scalable Data Analytics and Artificial Intelligence (ScaDS.AI) Dresden/Leipzig, Flatiron Institute, Harvard University (Author)
  • Krzysztof Gonciarz - , TUD Dresden University of Technology, Max Planck Institute of Molecular Cell Biology and Genetics, Center for Systems Biology Dresden (CSBD) (Author)
  • Ivo F. Sbalzarini - , Chair of Scientific Computing for Systems Biology, Clusters of Excellence PoL: Physics of Life, Max Planck Institute of Molecular Cell Biology and Genetics, Center for Systems Biology Dresden (CSBD), Center for Scalable Data Analytics and Artificial Intelligence (ScaDS.AI) Dresden/Leipzig (Author)

Abstract

We present data structures and algorithms for native implementations of discrete convolution operators over Adaptive Particle Representations (APR) of images on parallel computer architectures. The APR is a content-adaptive image representation that locally adapts the sampling resolution to the image signal. It has been developed as an alternative to pixel representations for large, sparse images as they typically occur in fluorescence microscopy. It has been shown to reduce the memory and runtime costs of storing, visualizing, and processing such images. This, however, requires that image processing natively operates on APRs, without intermediately reverting to pixels. Designing efficient and scalable APR-native image processing primitives, however, is complicated by the APR's irregular memory structure. Here, we provide the algorithmic building blocks required to efficiently and natively process APR images using a wide range of algorithms that can be formulated in terms of discrete convolutions. We show that APR convolution naturally leads to scale-adaptive algorithms that efficiently parallelize on multi-core CPU and GPU architectures. We quantify the speedups in comparison to pixel-based algorithms and convolutions on evenly sampled data. We achieve pixel-equivalent throughputs of up to 1TB/s on a single Nvidia GeForce RTX 2080 gaming GPU, requiring up to two orders of magnitude less memory than a pixel-based implementation.

Details

Original languageEnglish
Article number9796006
Pages (from-to)4197-4212
Number of pages16
JournalIEEE Transactions on Image Processing
Volume31
Publication statusPublished - 1 Jan 2022
Peer-reviewedYes

External IDs

Scopus 85132705028
ORCID /0000-0003-4414-4340/work/142252173

Keywords

Keywords

  • Image reconstruction, Image resolution, Image processing, Microscopy, Convolution, Data structures, Signal resolution

Library keywords