3D Shape Analysis for Connectomics

3D Shape Analysis for Connectomics

Connectomics has been described by Profs. Lichtman, Pfister, and Shavit as the field concerned with mapping networks of interconnected cells inside a nervous system, at the synaptic level. The field has seen a surge in the amount of available data with advances in electron microscopy producing detail rich 3D models of cellular organelle like neurons. During Spring 2019, I had the fortune to work with such exciting data at Prof. Hanspeter Pfister's lab, the Visual Computing Group at Harvard. I worked with nano-scale resolution 3D models of mitochondria and pyramidal neurons and attempted to make a tiny contribution towards shape analysis of these organelle. I'm grateful to my advisors Donglai Wei and Zudi Lin for giving me the opportuntiy to work on this project, and for their guidance and patience throughout the project.

Understanding the variance in shapes of an organelle like mitochondria may be seen as the first step towards relating shape to biological function. In this project we attempted to answer questions such as: How many types of organelle shapes exist? Can all these shapes be decomposed into simpler primitives? Towards that end, we have developed a simple shape descriptor called the "Reduced Graph", which provides an intuitive picture of complex 3D shapes, and at the same time preserves adequate topological details to be used for shape sensitive tasks such as shape matching and retrieval. we also proposed a nomenclature scheme for our Reduced Graph descriptor. We have submitted our work to a peer reviewed conference and are awaiting feedback. My contribution to the project has been designing and implementing an algorithm for extracting the Reduced Graph from a shape.

This website is organized into three sections:

  1. This page describes the problem at hand and provides an overview of the proposed solution.
  2. My experiments with existing shape descriptors are described here: [Different Shape Descriptors]
  3. Our proposed Reduced Graph descriptor, for whose extraction I designed and implemented an algorithm, is described here: [Reduced Graph Descriptor]

Dataset

Connectomics is now done at nano-scale resolution, i.e. we can tell apart two points of an organelle which are just a few nanometers (nm) apart. For this project I dealt with typical resolutions of 60nm. One implication of the availability of such high resolution data is that we can now observe never-seen-before details of organelle. Below are two examples of such rich 3D models from the dataset I worked with.

Left: A mitochondria with complex features such as loops. Right: A neuron with a dense dendrite structure.

Given these types of complex 3D shapes, our big picture goal is to come up with a shape descriptor which captures the essential features of a shape (such as connectivity between two points on the shape, length of paths between two points, loops etc.) and discards non-essential features (such as spines on a neuron's dendrites, loops with a tiny circumference in mitochondria).

Problem Statement and Proposed Solution

The aim is twofolds:

The figure below shows the pipeline developed by our team, in which we start with a 3D voxelized model of an organelle ("3D Shape" in the figure), extract its shape descriptor (the "Reduced Graph"), and propose a nomenclature scheme for naming such Reduced Graphs (the "Topological Nomenclature"). My contribution to this pipeline has been designing and implementing the Graph Reduction algorithm.

Pipeline for extracting Topological Nomenclature of a 3D Organelle Model (Source: Our paper.)

Before settling on the Reduced Graph shape descriptor, I tried some existing approaches, which can be seen here: [Different Shape Descriptors]. The details of the Graph Reduction algorithm can be found here: [Reduced Graph Descriptor].