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, size_t x) → size_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