Autofocusing with the help of the empirical Haar transform Przemysaw - - PowerPoint PPT Presentation

autofocusing with the help of the empirical haar transform
SMART_READER_LITE
LIVE PREVIEW

Autofocusing with the help of the empirical Haar transform Przemysaw - - PowerPoint PPT Presentation

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions Autofocusing with the help of the empirical Haar transform Przemysaw Sliwi nski and Krzysztof Berezowski Institute of Computer Engineering, Control


slide-1
SLIDE 1

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Autofocusing with the help

  • f the empirical Haar transform

Przemysław ´ Sliwi´ nski and Krzysztof Berezowski

Institute of Computer Engineering, Control and Robotics Wrocław University of Technology, POLAND

WASC 2012, Clermont-Ferrand, April 5-6th, 2012

slide-2
SLIDE 2

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Introduction

Presentation schedule

Motivations and inspirations

slide-3
SLIDE 3

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Introduction

Presentation schedule

Motivations and inspirations Model and formal assumptions

slide-4
SLIDE 4

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Introduction

Presentation schedule

Motivations and inspirations Model and formal assumptions Generic algorithm and its properties

slide-5
SLIDE 5

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Introduction

Presentation schedule

Motivations and inspirations Model and formal assumptions Generic algorithm and its properties AF criteria

slide-6
SLIDE 6

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Introduction

Presentation schedule

Motivations and inspirations Model and formal assumptions Generic algorithm and its properties AF criteria Unbalanced Haar Transform and Single-Photon AF

slide-7
SLIDE 7

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Introduction

Presentation schedule

Motivations and inspirations Model and formal assumptions Generic algorithm and its properties AF criteria Unbalanced Haar Transform and Single-Photon AF Experimental results and conclusions

slide-8
SLIDE 8

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Introduction

Motivations and inspirations

Problem A proper and reliable focusing algorithm is a conditio sine qua non

  • f a ’good image’. Not only from an aesthetic vantage point, but

also in automated applications. We exploit a plethora of the ’off-the-shelf’ theoretical results developed in various disciplines:

slide-9
SLIDE 9

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Introduction

Motivations and inspirations

Problem A proper and reliable focusing algorithm is a conditio sine qua non

  • f a ’good image’. Not only from an aesthetic vantage point, but

also in automated applications. We exploit a plethora of the ’off-the-shelf’ theoretical results developed in various disciplines:

signal and image processing, image analysis, harmonic analysis, control theory, or

slide-10
SLIDE 10

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Introduction

Motivations and inspirations

Problem A proper and reliable focusing algorithm is a conditio sine qua non

  • f a ’good image’. Not only from an aesthetic vantage point, but

also in automated applications. We exploit a plethora of the ’off-the-shelf’ theoretical results developed in various disciplines:

signal and image processing, image analysis, harmonic analysis, control theory, or information theory, probability theory and mathematical statistics, as well.

slide-11
SLIDE 11

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Introduction

Alternatives

Stereo-vision Solution Our algorithm works with standard matrix sensors & standard

  • ptics, and employs standard transforms and routines. . .
slide-12
SLIDE 12

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Introduction

Alternatives

Stereo-vision

two sensors

Solution Our algorithm works with standard matrix sensors & standard

  • ptics, and employs standard transforms and routines. . .
slide-13
SLIDE 13

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Introduction

Alternatives

Stereo-vision

two sensors two lenses, etc.

Solution Our algorithm works with standard matrix sensors & standard

  • ptics, and employs standard transforms and routines. . .
slide-14
SLIDE 14

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Introduction

Alternatives

Stereo-vision

two sensors two lenses, etc.

Light-field cameras Solution Our algorithm works with standard matrix sensors & standard

  • ptics, and employs standard transforms and routines. . .
slide-15
SLIDE 15

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Introduction

Alternatives

Stereo-vision

two sensors two lenses, etc.

Light-field cameras

lack resolution/dynamic range

Solution Our algorithm works with standard matrix sensors & standard

  • ptics, and employs standard transforms and routines. . .
slide-16
SLIDE 16

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Introduction

Alternatives

Stereo-vision

two sensors two lenses, etc.

Light-field cameras

lack resolution/dynamic range computational photography devices

Solution Our algorithm works with standard matrix sensors & standard

  • ptics, and employs standard transforms and routines. . .
slide-17
SLIDE 17

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Introduction

Alternatives

Stereo-vision

two sensors two lenses, etc.

Light-field cameras

lack resolution/dynamic range computational photography devices

Femtosecond lasers Solution Our algorithm works with standard matrix sensors & standard

  • ptics, and employs standard transforms and routines. . .
slide-18
SLIDE 18

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Introduction

Alternatives

Stereo-vision

two sensors two lenses, etc.

Light-field cameras

lack resolution/dynamic range computational photography devices

Femtosecond lasers

comparatively slow (like line scanners)

Solution Our algorithm works with standard matrix sensors & standard

  • ptics, and employs standard transforms and routines. . .
slide-19
SLIDE 19

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Introduction

Alternatives

Stereo-vision

two sensors two lenses, etc.

Light-field cameras

lack resolution/dynamic range computational photography devices

Femtosecond lasers

comparatively slow (like line scanners) computational photography devices

Solution Our algorithm works with standard matrix sensors & standard

  • ptics, and employs standard transforms and routines. . .
slide-20
SLIDE 20

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Problem statement

AF model

CAPTURED SCENE CAPTURED SCENE LENS LENS IMAGE SENSOR IMAGE SENSOR FOCUS FUNCTION

CALCULATOR

FOCUS FUNCTION

CALCULATOR

AF CONTROL AF CONTROL

MFD/INF MFD/INF

RANDOM FIELD BLOCK/IMPULSE SAMPLER

Q

R

Q

m n

mn

2 LOW-PASS FILTER f

slide-21
SLIDE 21

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Generic AF algorithm steps

1

Compute the focus function (with optional:

slide-22
SLIDE 22

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Generic AF algorithm steps

1

Compute the focus function (with optional:

1

denoising and

slide-23
SLIDE 23

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Generic AF algorithm steps

1

Compute the focus function (with optional:

1

denoising and

2

sensor output linearization).

slide-24
SLIDE 24

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Generic AF algorithm steps

1

Compute the focus function (with optional:

1

denoising and

2

sensor output linearization).

2

Shift the lens accordingly:

slide-25
SLIDE 25

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Generic AF algorithm steps

1

Compute the focus function (with optional:

1

denoising and

2

sensor output linearization).

2

Shift the lens accordingly:

1

determine the direction

slide-26
SLIDE 26

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Generic AF algorithm steps

1

Compute the focus function (with optional:

1

denoising and

2

sensor output linearization).

2

Shift the lens accordingly:

1

determine the direction

2

set the step-size

slide-27
SLIDE 27

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Generic AF algorithm steps

1

Compute the focus function (with optional:

1

denoising and

2

sensor output linearization).

2

Shift the lens accordingly:

1

determine the direction

2

set the step-size

3

Make it reliable in noisy environments!

slide-28
SLIDE 28

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Problem statement

Assumptions

1

The scene is a 2D homogenous second-order stationary process (thus an ergodic (in the wide sense) random field) with unknown distribution and unknown correlation function.

slide-29
SLIDE 29

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Problem statement

Assumptions

1

The scene is a 2D homogenous second-order stationary process (thus an ergodic (in the wide sense) random field) with unknown distribution and unknown correlation function.

2

The lens system is modeled with the help of the first-order

  • ptics laws, that is, the lens is merely a simple centered

moving average filter with an order proportional to the distance of the sensor from the image plane and to the size of the lens aperture.

slide-30
SLIDE 30

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Problem statement

Assumptions

1

The scene is a 2D homogenous second-order stationary process (thus an ergodic (in the wide sense) random field) with unknown distribution and unknown correlation function.

2

The lens system is modeled with the help of the first-order

  • ptics laws, that is, the lens is merely a simple centered

moving average filter with an order proportional to the distance of the sensor from the image plane and to the size of the lens aperture.

3

The image sensor acts as a block sampler, that the lens-produced image is orthogonally projected onto the space

  • f piecewise constant functions.
slide-31
SLIDE 31

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

AF algorithm foundations

The AF algorithm is based on the following lemmas: Lemma Under assumptions 1-3, the variance of the captured image is a unimodal function w.r.t. the order of the lens filter and attains its maximum value for the in-focus image. Lemma The variance estimate is tantamount to the orthogonal expansion

  • f the image acquired by the sensor.
slide-32
SLIDE 32

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

AF algorithm routine

The algorithm seeks the maximum of the (noised) image variance.

slide-33
SLIDE 33

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

AF algorithm routine

The algorithm seeks the maximum of the (noised) image variance. The focus functions is global!

slide-34
SLIDE 34

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

AF algorithm routine

The algorithm seeks the maximum of the (noised) image variance. The focus functions is global! Various (bi-)orthogonal expansions can be used to estimate the variance:

slide-35
SLIDE 35

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

AF algorithm routine

The algorithm seeks the maximum of the (noised) image variance. The focus functions is global! Various (bi-)orthogonal expansions can be used to estimate the variance:

trigonometric (DCT, Fourier, Hartley),

slide-36
SLIDE 36

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

AF algorithm routine

The algorithm seeks the maximum of the (noised) image variance. The focus functions is global! Various (bi-)orthogonal expansions can be used to estimate the variance:

trigonometric (DCT, Fourier, Hartley), Walsh-Hadamard (additions and subtractions only!),

slide-37
SLIDE 37

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

AF algorithm routine

The algorithm seeks the maximum of the (noised) image variance. The focus functions is global! Various (bi-)orthogonal expansions can be used to estimate the variance:

trigonometric (DCT, Fourier, Hartley), Walsh-Hadamard (additions and subtractions only!), wavelet, e.g. orthogonal Haar, Daubechies, or biorthogonal LeGall (5/3) and Cohen-Daubechies-Vial (9/7),

slide-38
SLIDE 38

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

AF algorithm routine

The algorithm seeks the maximum of the (noised) image variance. The focus functions is global! Various (bi-)orthogonal expansions can be used to estimate the variance:

trigonometric (DCT, Fourier, Hartley), Walsh-Hadamard (additions and subtractions only!), wavelet, e.g. orthogonal Haar, Daubechies, or biorthogonal LeGall (5/3) and Cohen-Daubechies-Vial (9/7), polynomial, e.g. Chebyshev, Legendre, Zernike (in general–any ’people’s polynomials’).

slide-39
SLIDE 39

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Possible AF function computation implementations

The following discrete orthogonal series transforms are available in the transform coders: DCT transform, (JPEG),

slide-40
SLIDE 40

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Possible AF function computation implementations

The following discrete orthogonal series transforms are available in the transform coders: DCT transform, (JPEG), Haar wavelet transform (JPEG 2K (Part II)), and

slide-41
SLIDE 41

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Possible AF function computation implementations

The following discrete orthogonal series transforms are available in the transform coders: DCT transform, (JPEG), Haar wavelet transform (JPEG 2K (Part II)), and Walsh-Hadamard transform (JPEG XR).

slide-42
SLIDE 42

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Maximum search algorithms

We can apply standard algorithms to find the function’s maximum in a noisy environment: Golden section-search (GSS) performed on:

slide-43
SLIDE 43

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Maximum search algorithms

We can apply standard algorithms to find the function’s maximum in a noisy environment: Golden section-search (GSS) performed on:

averaged image, or

slide-44
SLIDE 44

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Maximum search algorithms

We can apply standard algorithms to find the function’s maximum in a noisy environment: Golden section-search (GSS) performed on:

averaged image, or smoothed image (e.g. by any de-noising routine).

slide-45
SLIDE 45

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Maximum search algorithms

We can apply standard algorithms to find the function’s maximum in a noisy environment: Golden section-search (GSS) performed on:

averaged image, or smoothed image (e.g. by any de-noising routine).

Stochastic approximation (SA) exploiting:

slide-46
SLIDE 46

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Maximum search algorithms

We can apply standard algorithms to find the function’s maximum in a noisy environment: Golden section-search (GSS) performed on:

averaged image, or smoothed image (e.g. by any de-noising routine).

Stochastic approximation (SA) exploiting:

smoothed image, or

slide-47
SLIDE 47

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Maximum search algorithms

We can apply standard algorithms to find the function’s maximum in a noisy environment: Golden section-search (GSS) performed on:

averaged image, or smoothed image (e.g. by any de-noising routine).

Stochastic approximation (SA) exploiting:

smoothed image, or smoothed focus function (e.g. by using standard kernel convolutions).

slide-48
SLIDE 48

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

AF criteria

1

Unimodality — holds in theory. In practice, unimodality can be lost (aperture control!).

slide-49
SLIDE 49

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

AF criteria

1

Unimodality — holds in theory. In practice, unimodality can be lost (aperture control!).

2

Accuracy — corresponds to the resolution of the sensor.

slide-50
SLIDE 50

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

AF criteria

1

Unimodality — holds in theory. In practice, unimodality can be lost (aperture control!).

2

Accuracy — corresponds to the resolution of the sensor.

3

Reproducibility — a sharp top of the extremum holds in theory.

slide-51
SLIDE 51

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

AF criteria

1

Unimodality — holds in theory. In practice, unimodality can be lost (aperture control!).

2

Accuracy — corresponds to the resolution of the sensor.

3

Reproducibility — a sharp top of the extremum holds in theory.

4

Range — global. The variance of the image does not vanish.

slide-52
SLIDE 52

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

AF criteria

1

Unimodality — holds in theory. In practice, unimodality can be lost (aperture control!).

2

Accuracy — corresponds to the resolution of the sensor.

3

Reproducibility — a sharp top of the extremum holds in theory.

4

Range — global. The variance of the image does not vanish.

5

General applicability — a generic class of processes is admitted (ARMA models, Markov fields, and piecewise-smooth models).

slide-53
SLIDE 53

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

AF criteria

1

Unimodality — holds in theory. In practice, unimodality can be lost (aperture control!).

2

Accuracy — corresponds to the resolution of the sensor.

3

Reproducibility — a sharp top of the extremum holds in theory.

4

Range — global. The variance of the image does not vanish.

5

General applicability — a generic class of processes is admitted (ARMA models, Markov fields, and piecewise-smooth models).

6

Insensitivity to other parameters — particularly robust against the noise.

slide-54
SLIDE 54

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

AF criteria

1

Unimodality — holds in theory. In practice, unimodality can be lost (aperture control!).

2

Accuracy — corresponds to the resolution of the sensor.

3

Reproducibility — a sharp top of the extremum holds in theory.

4

Range — global. The variance of the image does not vanish.

5

General applicability — a generic class of processes is admitted (ARMA models, Markov fields, and piecewise-smooth models).

6

Insensitivity to other parameters — particularly robust against the noise.

7

Video signal compatibility — holds

slide-55
SLIDE 55

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

AF criteria

1

Unimodality — holds in theory. In practice, unimodality can be lost (aperture control!).

2

Accuracy — corresponds to the resolution of the sensor.

3

Reproducibility — a sharp top of the extremum holds in theory.

4

Range — global. The variance of the image does not vanish.

5

General applicability — a generic class of processes is admitted (ARMA models, Markov fields, and piecewise-smooth models).

6

Insensitivity to other parameters — particularly robust against the noise.

7

Video signal compatibility — holds

Block samplers = Foveon X3 sensor.

slide-56
SLIDE 56

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

AF criteria

1

Unimodality — holds in theory. In practice, unimodality can be lost (aperture control!).

2

Accuracy — corresponds to the resolution of the sensor.

3

Reproducibility — a sharp top of the extremum holds in theory.

4

Range — global. The variance of the image does not vanish.

5

General applicability — a generic class of processes is admitted (ARMA models, Markov fields, and piecewise-smooth models).

6

Insensitivity to other parameters — particularly robust against the noise.

7

Video signal compatibility — holds

Block samplers = Foveon X3 sensor. Low-pass filter + Dirac comb-based sampler = Bayer CFA sensors.

slide-57
SLIDE 57

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

AF criteria

1

Unimodality — holds in theory. In practice, unimodality can be lost (aperture control!).

2

Accuracy — corresponds to the resolution of the sensor.

3

Reproducibility — a sharp top of the extremum holds in theory.

4

Range — global. The variance of the image does not vanish.

5

General applicability — a generic class of processes is admitted (ARMA models, Markov fields, and piecewise-smooth models).

6

Insensitivity to other parameters — particularly robust against the noise.

7

Video signal compatibility — holds

Block samplers = Foveon X3 sensor. Low-pass filter + Dirac comb-based sampler = Bayer CFA sensors.

8

Fast implementation — all algorithms exploit ’fast’

  • transforms. . .
slide-58
SLIDE 58

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Experimental results

AF against aperture 15 17 19 21 23 25 27 29 f/1.2 f/1.6 f/2 f/4 f/8 f/16

slide-59
SLIDE 59

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Experimental results

AF against transform coder 5 10 15 20 25 JPG H-DCT DCT CD XR MR-CD

slide-60
SLIDE 60

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Experimental results

AF against image size 5 10 15 20 25 25 50 75 100 125 512 256 128 64 XR/256 XR/512 XR/128 XR/64

slide-61
SLIDE 61

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Single-Photon AF

Problem Can the generic algorithm be adapted to the Single-Photon Imagery? There are several noise sources Ylk = Ilk + Poissonlk + Gaussianlk + PRNUlk +crosstalklk + quantizationlk + . . .

slide-62
SLIDE 62

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Single-Photon AF

Problem Can the generic algorithm be adapted to the Single-Photon Imagery? There are several noise sources Ylk = Ilk + Poissonlk + Gaussianlk + PRNUlk +crosstalklk + quantizationlk + . . . Some of them are pure (random) noises, e.g.:

slide-63
SLIDE 63

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Single-Photon AF

Problem Can the generic algorithm be adapted to the Single-Photon Imagery? There are several noise sources Ylk = Ilk + Poissonlk + Gaussianlk + PRNUlk +crosstalklk + quantizationlk + . . . Some of them are pure (random) noises, e.g.:

Shot noise (of Poisson distribution),

slide-64
SLIDE 64

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Single-Photon AF

Problem Can the generic algorithm be adapted to the Single-Photon Imagery? There are several noise sources Ylk = Ilk + Poissonlk + Gaussianlk + PRNUlk +crosstalklk + quantizationlk + . . . Some of them are pure (random) noises, e.g.:

Shot noise (of Poisson distribution), Thermal noise (of Gaussian distribution).

slide-65
SLIDE 65

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Single-Photon AF

Problem Can the generic algorithm be adapted to the Single-Photon Imagery? There are several noise sources Ylk = Ilk + Poissonlk + Gaussianlk + PRNUlk +crosstalklk + quantizationlk + . . . Some of them are pure (random) noises, e.g.:

Shot noise (of Poisson distribution), Thermal noise (of Gaussian distribution).

Some are random but fixed, e.g. Photo-Response Non-uniformity

slide-66
SLIDE 66

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Single-Photon AF

Problem Can the generic algorithm be adapted to the Single-Photon Imagery? There are several noise sources Ylk = Ilk + Poissonlk + Gaussianlk + PRNUlk +crosstalklk + quantizationlk + . . . Some of them are pure (random) noises, e.g.:

Shot noise (of Poisson distribution), Thermal noise (of Gaussian distribution).

Some are random but fixed, e.g. Photo-Response Non-uniformity For the others, it is just convenient to model them as a

  • noise. . .
slide-67
SLIDE 67

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Unbalanced Haar Transform

In order to model PRNU for each pixel we use the Unbalanced Haar Transform instead of the classic one

slide-68
SLIDE 68

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Unbalanced Haar Transform

In order to model PRNU for each pixel we use the Unbalanced Haar Transform instead of the classic one The basic transform step becomes a little bit more complicated than ˆ αm−1,n =

√ 2 2 ˆ

αm,2n +

√ 2 2 ˆ

αm,2n+1 vs. ¯ αm−1,n =

  • Im,2n

Im−1,n ¯

αm,2n +

  • Im,2n+1

Im−1,n ¯

αm,2n+1, with Im−1,n = Im,2n+1 + Im,2n+1 (where Im,2n+1, Im−1,n are the non-uniformity indices).

slide-69
SLIDE 69

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Unbalanced Haar Transform

Application of UHT has some advantages: Can be plugged-in into the standard AF algorithm.

slide-70
SLIDE 70

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Unbalanced Haar Transform

Application of UHT has some advantages: Can be plugged-in into the standard AF algorithm. Remains fast, i.e. linear with number of pixels.

slide-71
SLIDE 71

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Unbalanced Haar Transform

Application of UHT has some advantages: Can be plugged-in into the standard AF algorithm. Remains fast, i.e. linear with number of pixels. Allows for in situ image denoising.

slide-72
SLIDE 72

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Unbalanced Haar Transform

Application of UHT has some advantages: Can be plugged-in into the standard AF algorithm. Remains fast, i.e. linear with number of pixels. Allows for in situ image denoising. Can be computed in parallel.

slide-73
SLIDE 73

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Unbalanced Haar Transform

Application of UHT has some advantages: Can be plugged-in into the standard AF algorithm. Remains fast, i.e. linear with number of pixels. Allows for in situ image denoising. Can be computed in parallel.

  • But. . . requires computing square roots. . .
slide-74
SLIDE 74

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Single-photon AF

We have now Ylk ∼ Ilk + Poissonlk + Gaussianlk. The single-photon-denoising algorithm has two simple steps:

slide-75
SLIDE 75

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Single-photon AF

We have now Ylk ∼ Ilk + Poissonlk + Gaussianlk. The single-photon-denoising algorithm has two simple steps:

’removal’ of the Gaussian part by UHT transform with a

  • thresholding. Then

Ylk ∼ Ilk + Poissonlk

slide-76
SLIDE 76

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Single-photon AF

We have now Ylk ∼ Ilk + Poissonlk + Gaussianlk. The single-photon-denoising algorithm has two simple steps:

’removal’ of the Gaussian part by UHT transform with a

  • thresholding. Then

Ylk ∼ Ilk + Poissonlk application of the Anscombe transform and repeating the previous step (i.e. the UHT transform with a thresholding). Thus Ylk ∼ Ilk

slide-77
SLIDE 77

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Single-photon AF

We have now Ylk ∼ Ilk + Poissonlk + Gaussianlk. The single-photon-denoising algorithm has two simple steps:

’removal’ of the Gaussian part by UHT transform with a

  • thresholding. Then

Ylk ∼ Ilk + Poissonlk application of the Anscombe transform and repeating the previous step (i.e. the UHT transform with a thresholding). Thus Ylk ∼ Ilk

The rest of the AF algorithm remains unchanged.

slide-78
SLIDE 78

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Final conclusions

The proposed AF algorithm: Is robust against a noise.

slide-79
SLIDE 79

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Final conclusions

The proposed AF algorithm: Is robust against a noise. Works with a standard (cheap) equipment.

slide-80
SLIDE 80

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Final conclusions

The proposed AF algorithm: Is robust against a noise. Works with a standard (cheap) equipment. Can reuse existing IPs.

slide-81
SLIDE 81

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Final conclusions

The proposed AF algorithm: Is robust against a noise. Works with a standard (cheap) equipment. Can reuse existing IPs. Can effectively be implemented (e.g. in situ).

slide-82
SLIDE 82

Introduction Problem statement and algorithm Properties Single-photon AF Conclusions

Example

The demonstration movie can be found at: http://diuna.iiar.pwr.wroc.pl/sliwinski/gss-af.avi