Operators

# moments_region_2nd (Operator)

## Name

moments_region_2nd — Calculate the geometric moments of regions.

## Signature

moments_region_2nd(Regions : : : M11, M20, M02, Ia, Ib)

## Description

moments_region_2nd calculates the geometric moments M11, M20 and M02 of the input regions in Regions. Furthermore the major and minor axis of the input regions are calculated and returned in Ia and Ib.

The covariance matrix is given by:

The row-dependent moment of 2nd order is returned in M20 and the column-dependent moment of 2nd order is returned in M02. The moment M11 represents the covariance between the row and column coordinates of the region points.

Calculation:

and are the coordinates of the center of gravity of a region R. Then the moments are defined by:

wherein r and c run over all pixels of the region R.

Furthermore the length of the major and minor axes are defined by:

wherein is defined by:

The equation for the major and minor axes can be derived from the definition of the moments by diagonalizing the covariance matrix and reforming the resulting formula.

If more than one region is passed, the results are returned in tuples. The index of a tuple element corresponds to the index of the respective input region.

If an empty region is passed, 0.0 is returned for all parameters, if no other behavior was set (see set_system).

## Parallelization

• Multithreading type: reentrant (runs in parallel with non-exclusive operators).
• Automatically parallelized on tuple level.

## Parameters

Regions (input_object)  region(-array) object

Input regions.

M11 (output_control)  real(-array) (real)

Product of inertia of the axes through the center parallel to the coordinate axes.

M20 (output_control)  real(-array) (real)

Moment of 2nd order (row-dependent).

M02 (output_control)  real(-array) (real)

Moment of 2nd order (column-dependent).

Ia (output_control)  real(-array) (real)

Length of the major axis of the input region.

Ib (output_control)  real(-array) (real)

Length of the minor axis of the input region.

## Complexity

If F is the area of the region the mean runtime complexity is O(sqrt(F)).

## Result

The operator moments_region_2nd returns the value 2 (H_MSG_TRUE) if the input is not empty. The behavior in case of empty input (no input regions available) is set via the operator set_system('no_object_result',<Result>). The behavior in case of empty region (region is the empty set) is set via set_system('empty_region_result',<Result>). If necessary an exception is raised.