two_point.rst

Two Point (Dasch)

Introduction

The "Dasch two-point" deconvolution algorithm is one of several described in the Dasch [1] paper. See also the three_point and onion_peeling descriptions.

How it works

The Abel integral is broken into intervals between the r_j points, and P^\prime(r) is assumed constant between r_j and r_{j+1}.

When to use it

This method is simple and computationally very efficient. The method incorporates no smoothing.

How to use it

To complete the inverse transform of a full image with the two_point method, simply use the :class:`abel.Transform` class:

abel.Transform(myImage, method='two_point').transform

If you would like to access the two_point algorithm directly (to transform a right-side half-image), you can use :func:`abel.dasch.two_point_transform`.

Example

.. plot:: ../examples/example_dasch_methods.py
    :include-source:

or more information on the PyAbel implementation of the two_point algorithm, please see Pull Request #155.

Citation

[1] Dasch, Applied Optics, Vol 31, No 8, March 1992, Pg 1146-1152.