3D Analysis application examples

The Aphelion™ software product is the only commercial product able to process real volumic 3D images, and display these images using a voxel representation.

To better evaluate the power of the 3D Image Processing and 3D Image Display modules, let us consider an application:

Extract volume fraction of zircone grains in aluminum

Related video

Below is an application that was developed by an Aphelion™ users in the field of Material Science. The image is courtesy of ESRF Grenoble, INSA de Lyon, GEMPPM, and Ecole des Mines de Paris.

The goal of the application is to extract volume fraction of zircone grains in aluminum, to compute grain sizing and determine neighbor distribution, working on the 3D volume.

In the past, most of the 3D analysis were performed using the set of 2D sections of a 3D volume. Nowadays, thanks the processing power of computers, and the quality of sensors, it is possible to deal directly with 3D images.

Below is the description of a very innovative technique based on 3D Morphology and 3D Image Understanding to compute grain sizing and neighbor distribution, two analysis which can only be performed on the 3D data.

  • The image is acquired using a X-Ray micro-tomograph with an advanced sensor.
    Loading images/3D_01.gif

    Original Image

  • The computation of the mean value of the binary image gives the zircone volumic fraction.
  • Since the contrast between the two phases is good, a simple threshold is applied to extract the zircone phase.
  • Since zircone particles are almost spherical, they all appear in the image as a stack of spheres touching each other. The use of the ClustersSplitConvex operator will help to segment all the spheres which are actually convex particles. The operator is based on the 3D implementation of the watershed, a morphological operator.
    Loading images/3D_03.gif

    3D label displayed as a volume

    Loading images/3D_04.gif

    3D label displayed as 3 orthogonal sections

  • To avoid biased measurements, a BorderKill operator is applied to remove all particles intersecting the edge of the image. A Miles-Lantuéjoul correction could also be applied to actually get unbiased measurements taking into account the size of the volume, and the operators involved in the process.
  • The binary image is now converted into a 3D ObjectSet based on the 26-Connectivity. Note that cubic and cubic centered face grids are supported in Aphelion 3D. All spheres are now perfectly individualized. A set of measurements based on the shape are computed, such as the sphericity, and the intercept numbers in the main directions of the grid. The size distribution of the particles is computed and displayed as an histogram, as shown in the chart below.
    Loading images/3D_05.png

    Size distribution

  • The last step of the application involves the computation of the number of neighbors for each zircone spherical grain. This number can only be accessed from the 3D data, using the notion of 3D ObjectSets. The classical method is based on the extraction of each sphere in the volume. Each sphere is dilated, and then intersected with the binary image. A geodesic reconstruction is then performed, and the number of components is computed. It has to be done for each grain, and it takes several minutes to run. We are proposing to use another technique, much faster, and based on the Aphelion ObjectSets.

With the help of ObjectSets, the computation is no longer performed on pixels, but on 3D objects. Since objects are already individualized in an ObjectSet, and available as rasters, we perform a dilation on the Objects, with the condition that grains remain individualized even if they overlap after the dilation. The final result gives the number of overlaps for each grain, and is displayed in the standard Aphelion grid as a new attribute. This computation is very quickly performed since no pixel information is required, and it really proves Aphelion has superior capability than any other software when dealing with 3D images.

The following chart gives the neighbor distribution for the Zircone image, and the grid shows the value of the neighbor attribute. As in the 2D version of Aphelion, message passing is available between the grid and the chart.

Loading images/3D_06.png

Neighbor distribution


ReVAnalyzer: Software product to quantify the volume of Retinal Pigment Epithelium Detachment in macular area

Download the flyer Related poster Related video tutorial Download

ReV Analyzer (Retinal Volume Analyzer) is a software product that takes as input OCT images of a patient, detects and analyzes the retinal Pigment Epithelium Detachment (PED) in that series of images, generates a 3D volume in the 3D space, and provides information on the PED evolution over different visits. It also works on other retina-choroidal areas in a semi-automatic way.

ReVAnalyzer can be used for patients with vascularized PEDs in the scope of exudative Age Related Macula Degeneration, when the PED is visible in multi modal images. Typically, an ophthalmologic clinical examination is performed on a patient during multiple visits, and images of fluorescein and ICG angiography are captured, along with an SD OCT examination in Enhanced Depth Imaging (EDI) mode. ReVAnalyzer is fully compatible with the Spectralis HRA+OCT system manufactured by Heidelberg Engineering GmbH, Heidelberg, Germany.

A retinal volume of 61 lines, centered on the fovea is captured using the Spectralis system. The volume is then exported from the Heidelberg system, and automatically imported into the ReVAnalyzer software product.

Once the volume is opened in the software, the PED area is detected by the operator. Two detection modes are available:

  • A semi-automatic mode that lets the user outline the detachment in each section where it is visible, using a User Assisted tool that automatically fits to the boundary of the detachment.
  • An automatic mode is also available that detects the detachment in the central section, once the user has manually specified that section, and then propagates through all the sections where the detachment is present.

Outline detachment using the User Assisted tool that automatically fits to the boundary of the detachment

Eraser tool used to adjust the detection

In both modes, the user is capable of visually validating the detection in each section, and using a Pen/Eraser tool to adjust the detection.

Once the detection is completed, the software automatically reconstructs and displays the PED in 3D. A full set of display parameters is available to let the user adjust the display in terms of color and orientation, display a section along with the 3D volume, zoom in and out, and move the PED in all directions. The computed volume is provided in pixels and calibrated unit (cubic millimeter).

3D volume and section display

ReVAnalyzer can track the evolution of a detachment over multiple visits for a given patient, and provide a report that can be printed and saved.

ReVAnalyzer has been specified and fully validated by the team of Dr. Souied, Hôpital Intercommunal de Créteil, France. A semi-automated analysis has been performed by three independent readers, twice. Inter and intra reader coefficients were computed based on the Bland Altman statistical method. The analysis was performed on 20 eyes of 20 patients. The value of the coefficient was above 0.99, and the location and size of the PED did not have any impact on the results.

ReVAnalyzer is an excellent tool to accurately measure the PED volume, monitor its evolution over time, and evaluate its response to therapeutic treatments. It can be used by researchers in the field of ophthalmology, big pharmaceutical companies dealing with clinical studies, as well as private ophthalmologists performing patient monitoring through multiple visits.

Video tutorial