cédric dufour ( lts-ibcm collaboration )

1

Signal Processing Laboratory Swiss Federal Institute of Technology, Lausanne

Cédric Dufour( LTS-IBCM Collaboration )

The ‘The ‘microtubulesmicrotubules’ project’ project

2


The ‘The ‘microtubulesmicrotubules’ project’ project Morphological filteringMorphological filtering Markers extractionMarkers extraction Microtubules extractionMicrotubules extraction ResultsResults Algorithm testingAlgorithm testing Final assessmentFinal assessment

Content

3


Obtain Obtain specific proteins’ density statisticsspecific proteins’ density statistics, related to , related to the neural cell microtubules structures.the neural cell microtubules structures.

Markers

The ‘microtubules’ project The goal

Microtubules

4


The ‘microtubules’ project The various steps

Isolate the markers Isolate the markers maskmask..

Isolate the microtubules Isolate the microtubules mask mask andand skeleton skeleton..

Compute the microtubules’ surface and length.Compute the microtubules’ surface and length.

Compute the markers’ quantity, overall and near the Compute the markers’ quantity, overall and near the

microtubules.microtubules.

5


The ‘microtubules’ project The way it is done

Markers extraction :Markers extraction :

Morphological filtering and local maximum Morphological filtering and local maximum

detection.detection.

Microtubules extraction :Microtubules extraction :

Selective filtering using linear oriented structuring Selective filtering using linear oriented structuring

element’s correlation with thresholded image.element’s correlation with thresholded image.

6


Morphological filtering What is it ?

Morphological filteringMorphological filtering is a filtering method is a filtering method

originated from the theory of originated from the theory of mathematical mathematical

morphologymorphology..

The base of all morphological processing are the The base of all morphological processing are the

‘‘erosionerosion’ and ‘’ and ‘dilationdilation’ morphological functions.’ morphological functions.

7


Morphological filtering Binary ‘erosion’

Bb

B bXXbyByyX

,,

Original set Eroded set

Structuring elmt.

8


Morphological filtering Binary ‘dilation’

Bb

B bXBbXxbxX

,,

Original set Dilated set

Structuring elmt.

9


Morphological filtering Binary ‘opening’

XX BBB

Original set Opened set

Structuring elmt.

10


Morphological filtering Binary ‘closing’

XX BBB

Original set Closed set

Structuring elmt.

11


Replace the Replace the intersectionintersection//unionunion operators with operators with infimuminfimum//supremumsupremum operators operators

yxfxf

yxfxf

ByB

ByB

sup

inf

Morphological filtering Function (or ‘gray scale’) morphology

Original fct.

Dilated fct.

Eroded fct.

Original fct.

Closed fct.

Opened fct.

12


Base imageBase image

Markers extraction Step by step : 1

13


Opening with a 3x3 squareOpening with a 3x3 square


14


Subtraction with original image (Tophat)Subtraction with original image (Tophat)


15


NormalizationNormalization


16


‘‘Loosy’ thresholdLoosy’ threshold


17


Remove small elements (artefacts)Remove small elements (artefacts)


18


Mask base image ( and normalize )Mask base image ( and normalize )


19


Local maximum detection ( in 5x5 disc neighborhood )Local maximum detection ( in 5x5 disc neighborhood )


20


Resulting Resulting markers maskmarkers mask

Markers extraction In the end

21


Base imageBase image

Microtubules extraction Step by step : 1

22


Opening with 21x21 square (Opening with 21x21 square ( background) background)


23


Subtraction with original image (Tophat)Subtraction with original image (Tophat)


24


Markers removalMarkers removal


25


ThresholdThreshold


26


Resulting binary imageResulting binary image


27


Oriented linear element correlationOriented linear element correlation


28


Small elements (artefacts) removalSmall elements (artefacts) removal


29


Repeat last 2 steps for different orientationsRepeat last 2 steps for different orientations


30


Thresholding of filter accumulatorThresholding of filter accumulator


31


Closing with a 3x3 cross to remove irregularitiesClosing with a 3x3 cross to remove irregularities


32


Skeleton by thinning and cntd. points suppressingSkeleton by thinning and cntd. points suppressing


33


Resulting Resulting microtubules mask microtubules mask andand skeleton skeleton

Microtubules extraction In the end

34


Results Image 1 (demo image)

35


Results Image 2 (low density)

36


Results Image 3 (high density)

37


Errors Errors markers: 0% / 4.8% markers: 0% / 4.8%

Algorithm testing Real image (PSNR = 35dB)

38


Errors Errors microtubules: 0.7% / 0.5%; markers: 0% / 0% microtubules: 0.7% / 0.5%; markers: 0% / 0%

Algorithm testing Synthetic image (PSNR = 35dB)

39


We’ve been able to offer the biologists…We’ve been able to offer the biologists…– a fully automatic analysis program,

– running in a powerful and wide-spread environment

(MatLab),

– giving good results, according to the biologists’ needs.

Final assessment In general

40


Microtubules mask extraction :Microtubules mask extraction :

The poor quality of the input images (very low The poor quality of the input images (very low

contrast) leads to a low-efficiency microtubules contrast) leads to a low-efficiency microtubules

mask extraction procedure (mask extraction procedure (the algorithm misses the algorithm misses

the most evanescent structuresthe most evanescent structures).).

Final assessment Problems

41


Microtubules mask extraction :Microtubules mask extraction :

Use tracking algorithm to follow the full microtubule Use tracking algorithm to follow the full microtubule

(pseudo-linear) structure.(pseudo-linear) structure.

N.B.N.B. This is not easy because of the variable This is not easy because of the variable

number of microtubules that may cross in one point number of microtubules that may cross in one point

(resulting in tracking uncertainties)(resulting in tracking uncertainties)

Final assessment Improvements

42


Thank you for your attentionThank you for your attention

cédric dufour ( lts-ibcm collaboration )

Documents