Heat Equation and its applications in imaging processing and mathematical biology

Transcription

1 Heat Equation and its applications in imaging processing and mathematical biology Yongzhi Xu Department of Mathematics University of Louisville Louisville, KY 40292

2 1. Introduction

3 System of heat equations

4 2. Derivation of heat equation 2.1. Random walks

5

6 2.2. Fick s law

7

8 Then the balance law implies

9 3. Applications of heat equation 1. Heat equation in image processing 2. Heat equation in cancer model and spatial ecological model

10 3.1. Heat equation in image processing Sampling an image: f(x i )

11 Three components of image processing: 1. Image Compression

12 2. Image Denoising

13 3. Image Analysis

14 One common need: Smoothing Smoothing is a necessary part of image formation. An image can be correctly represented by a discrete set of values, the samples, only if it has been previously smoothed.

15 How to smooth a image? Convolution with a bell-shaped function

16

17

18

19 We blur the image! Questions: 1. What does it have to do with heat equation? 2. How will it be helpful for image processing?

20

21

22

23

24

25

26 Answer to question 1:

27 Modified heat equations and applications: Deblur an image by reversing time in the heat equation:

28 Smoothing to detect edges

29 Level curves

30

31

32 3.2: Reaction-diffusion Models and Pattern Formations Let U(x,t) and V(x,t) be the density functions of two chemicals or species which interact or react Alan Turing, The Chemical Basis of Morphogenesis, Phil. Trans. Roy. Soc. (1952). Acknowledgement: Pictures of animal patterns are from Junping Shi s website.

33 Mathematical Biology by James Murray Emeritus Professor University of Washington, Seattle Oxford University, Oxford

34 Why do animals coats have different patterns?

35 Murray s theory Murray suggests that a single mechanism could be responsible for generating all of the common patterns observed. This mechanism is based on a reaction-diffusion system of the morphogen prepatterns, and the subsequent differentiation of the cells to produce melanin simply reflects the spatial patterns of morphogen concentration.

36 Reaction-diffusion systems Domain: rectangle (0, a) X (0,b) Boundary conditions: head and tail (no flux), body side (periodic)

37 Theorem 1 : Snakes always have striped (ring) patterns, but not spotted patterns. Turing-Murray Theory: snake is the example of b/a is large.

38 Snake pictures (stripe patterns)

39 Theorem 2 : There is no animal with striped body and spotted tail, but there is animal with spotted body and striped tail. Turing-Murray theory: The same reaction-diffusion mechanism should be responsible for the patterns on both body and tail. The body is always wider than the tail. If the body is striped, then the tail must also be striped.

40 More examples: Spotted body and striped tail Genet (left), Giraffe (right)

41 Cancer model: Free boundary problem model of DCIS Let the tumor to be within a rigid cylinder occupying a region

42

43 --- the growing boundary of the tumor --- the dimensionless nutrient concentration in the surrounding --- the rate of surrounding transfer per unit length --- the transfer of nutrient from surrounding --- the nutrient consumption rate --- the ratio of the nutrient diffusion time scale To the tumor growth time scale. C<< the nutrient concentration

44 Micropapillary DCIS Cribiform DCIS

45

46 Solid DCIS (1-D model) Cribiform DCIS (Spread evenly) Papillary & MicroPapillary DCIS (Baby tree) Moving, but not growing Stable

47 More sophisticate models may be considered. (I) mixed models

48 (II) segregated models

49 Inverse Problems in cancer modeling Inverse Problem Use one incisional biopsy

50 Inverse Problem use a sequence of needle biopsy

51 Inverse Problem use a sequence of tomography

52 Acknowledgement: Some of the graphs and pictures are copied from the manuscript of Frederic Guichard and Jean-Michel Morel, and from the website of Junping Shi.

53 Thank You!