Whole cell tracking and movement reconstruction through an optimal - - PowerPoint PPT Presentation

whole cell tracking and movement reconstruction through
SMART_READER_LITE
LIVE PREVIEW

Whole cell tracking and movement reconstruction through an optimal - - PowerPoint PPT Presentation

Whole cell tracking and movement reconstruction through an optimal control problem Feng Wei Yang with Anotida Madzvamuse and Chandrasekhar Venkataraman University of Sussex F.W.Yang@sussex.ac.uk 8 June 2015 at IBiDi 1 / 24 Outline Whole


slide-1
SLIDE 1

Whole cell tracking and movement reconstruction through an optimal control problem

Feng Wei Yang with Anotida Madzvamuse and Chandrasekhar Venkataraman

University of Sussex F.W.Yang@sussex.ac.uk

8 June 2015 at IBiDi

1 / 24

slide-2
SLIDE 2

Outline

Whole cell tracking through optimal control problem

The mathematical model Toy models for proving concepts Real world application 3-D simulation

IBiDi data

Our automatic segmentation algorithm Particle tracking results

2 / 24

slide-3
SLIDE 3

Objectives

To track the morphology of cells and reconstruct their movements:

  • V. Peschetola et al. Cytoskeleton, 2013

3 / 24

slide-4
SLIDE 4

What is our signature

Particle tracking: Particle tracking: Particle tracking: The morphology of cells are not considered Manually tracking is slow Automatic tracking algorithms are often flawed

Segmentation is suboptimal for real data Tracking through patten recognition is challenging

Pure geometric math models: Pure geometric math models: Pure geometric math models: Resolution of the data matters Typically no cell-setting is considered It is a complicated procedure to

  • btain the results

Computational power and advanced numerical methods have to be included for 3-D real-life cell tracking

4 / 24

slide-5
SLIDE 5

Our optimal control model

The volume conserved mean curvature flow:

  • V

V V (x x x, t) = (−σH(x x x, t) + η(x x x, t) + λV (t))v v v(x x x, t) on Γ(t), t ∈ (0, T], Γ(0) = Γ0.

The phase-field approximation of the above equation - Allen-Cahn:

     ∂tφ(x x x, t) = △φ(x x x, t) − 1

ǫ2 G ′(φ(x

x x, t)) − 1

ǫ(cGη(x

x x, t) − λ(t)) in Ω × (0, T], ∇φ · ν ν νΩ = 0 on ∂Ω × (0, T], φ(·, 0) = φ0 in Ω.

5 / 24

slide-6
SLIDE 6

Our optimal control model cont.

The objective functional:

J(φ, η) = 1 2

(φ(x x x, T) − φobs(x x x))2 dx x x + θ 2 T

η(x x x, t)2dx x xdt,

and now we solve the minimisation problem:

minηJ(φ, η), with J given above.

6 / 24

slide-7
SLIDE 7

Obtaining solutions

We are using one of the most efficient solution methods, combining most advanced adaptive techniques. Meanwhile, parallelism is employed and the computation has been carried out on large computer cluster with multiple number of computational cores.

7 / 24

slide-8
SLIDE 8

Toy model for proving concepts

A circle becomes two ellipses. Initial data Desired data

8 / 24

slide-9
SLIDE 9

Toy model for proving concepts

9 / 24

slide-10
SLIDE 10

Toy model for proving concepts cont.

A circle becomes 4 children circles. Initial data Desired data

10 / 24

slide-11
SLIDE 11

Toy model for proving concepts cont.

11 / 24

slide-12
SLIDE 12

Toy model for proving concepts cont.

12 / 24

slide-13
SLIDE 13

Real world application

Two segmented cell images. Initial data Desired data

13 / 24

slide-14
SLIDE 14

Real world application cont.

14 / 24

slide-15
SLIDE 15

Real world application cont.

15 / 24

slide-16
SLIDE 16

3-D simulation

(a) (b) (c) (d)

16 / 24

slide-17
SLIDE 17

3-D simulation cont.

17 / 24

slide-18
SLIDE 18

IBiDi data

18 / 24

slide-19
SLIDE 19

IBiDi data cont.

19 / 24

slide-20
SLIDE 20

Issues with basic segmentation

20 / 24

slide-21
SLIDE 21

Our simple solution to segmentation

21 / 24

slide-22
SLIDE 22

Particle tracking

  • N. Stuurman. Mtrack2 at valelab.ucsf.edu, 2009

22 / 24

slide-23
SLIDE 23

Particle tracking cont.

23 / 24

slide-24
SLIDE 24

The end

Thank you.

24 / 24