bilinear: Bilinear interpolator#

Bilinear interpolator, peak finder, line-profile for images

class BilinearImage#

Bilinear interpolator for images … or any data on a regular grid

__call__()#

Function f((y, x)) where f is a continuous function made from the image and (y,x)=(row, column) is the pixel coordinates in natural C-order

Parameters:

x – 2-tuple of float (row, column)

Returns:

Interpolated signal from the image

coarse_local_maxi(self, Py_ssize_t x) Py_ssize_t#

Return the nearest local maximum … without sub-pixel refinement

Parameters:

idx – start index (=row*width+column)

Returns:

local maximum index

local_maxi(self, coord)#

Return the nearest local maximum … with sub-pixel refinement

Nearest maximum search:

steepest ascent

Sub-pixel refinement:

Second order Taylor expansion of the function; At the maximum, the first derivative is null delta = x-i = -Inverse[Hessian].gradient if Hessian is singular or |delta|>1: use a center of mass.

Parameters:

coord – 2-tuple of scalar (row, column)

Returns:

2-tuple of float with the nearest local maximum

map_coordinates(self, coordinates)#

Map coordinates of the array on the image

Parameters:

coordinates – 2-tuple of array of the same size (row_array, column_array)

Returns:

array of values at given coordinates

opp_f(self, coord)#

Function -f((y,x)) for peak finding via minimizer.

Gives large number outside the boundaries to return into the image

Parameters:

x – 2-tuple of float in natural C order, i.e (row, column)

Returns:

Negative interpolated signal from the image

profile_line(self, src, dst, int linewidth=1, method=u'mean')#

Return the mean or sum of intensity profile of an image measured along a scan line.

Parameters:
  • src (2-tuple of numeric scalar) – The start point of the scan line.

  • dst (2-tuple of numeric scalar) – The end point of the scan line. The destination point is included in the profile, in contrast to standard numpy indexing.

  • linewidth (int) – Width of the scanline (unit image pixel).

  • method (str) – ‘mean’ or ‘sum’ depending if we want to compute the mean intensity along the line or the sum.

Returns:

The intensity profile along the scan line. The length of the profile is the ceil of the computed length of the scan line.

Return type:

1d array

Inspired from skimage