ClassesClassesClassesClasses | | | | Operators

fit_circle_contour_xldfit_circle_contour_xldFitCircleContourXldfit_circle_contour_xldFitCircleContourXldFitCircleContourXld (Operator)

Name

fit_circle_contour_xldfit_circle_contour_xldFitCircleContourXldfit_circle_contour_xldFitCircleContourXldFitCircleContourXld — Approximate XLD contours by circles.

Signature

fit_circle_contour_xld(Contours : : Algorithm, MaxNumPoints, MaxClosureDist, ClippingEndPoints, Iterations, ClippingFactor : Row, Column, Radius, StartPhi, EndPhi, PointOrder)

Herror fit_circle_contour_xld(const Hobject Contours, const char* Algorithm, const Hlong MaxNumPoints, double MaxClosureDist, const Hlong ClippingEndPoints, const Hlong Iterations, double ClippingFactor, double* Row, double* Column, double* Radius, double* StartPhi, double* EndPhi, char* PointOrder)

Herror T_fit_circle_contour_xld(const Hobject Contours, const Htuple Algorithm, const Htuple MaxNumPoints, const Htuple MaxClosureDist, const Htuple ClippingEndPoints, const Htuple Iterations, const Htuple ClippingFactor, Htuple* Row, Htuple* Column, Htuple* Radius, Htuple* StartPhi, Htuple* EndPhi, Htuple* PointOrder)

Herror fit_circle_contour_xld(Hobject Contours, const HTuple& Algorithm, const HTuple& MaxNumPoints, const HTuple& MaxClosureDist, const HTuple& ClippingEndPoints, const HTuple& Iterations, const HTuple& ClippingFactor, double* Row, double* Column, double* Radius, double* StartPhi, double* EndPhi, char* PointOrder)

Herror fit_circle_contour_xld(Hobject Contours, const HTuple& Algorithm, const HTuple& MaxNumPoints, const HTuple& MaxClosureDist, const HTuple& ClippingEndPoints, const HTuple& Iterations, const HTuple& ClippingFactor, HTuple* Row, HTuple* Column, HTuple* Radius, HTuple* StartPhi, HTuple* EndPhi, HTuple* PointOrder)

double HXLDCont::FitCircleContourXld(const HTuple& Algorithm, const HTuple& MaxNumPoints, const HTuple& MaxClosureDist, const HTuple& ClippingEndPoints, const HTuple& Iterations, const HTuple& ClippingFactor, double* Column, double* Radius, double* StartPhi, double* EndPhi, char* PointOrder) const

HTuple HXLDContArray::FitCircleContourXld(const HTuple& Algorithm, const HTuple& MaxNumPoints, const HTuple& MaxClosureDist, const HTuple& ClippingEndPoints, const HTuple& Iterations, const HTuple& ClippingFactor, HTuple* Column, HTuple* Radius, HTuple* StartPhi, HTuple* EndPhi, HTuple* PointOrder) const

void FitCircleContourXld(const HObject& Contours, const HTuple& Algorithm, const HTuple& MaxNumPoints, const HTuple& MaxClosureDist, const HTuple& ClippingEndPoints, const HTuple& Iterations, const HTuple& ClippingFactor, HTuple* Row, HTuple* Column, HTuple* Radius, HTuple* StartPhi, HTuple* EndPhi, HTuple* PointOrder)

void HXLDCont::FitCircleContourXld(const HString& Algorithm, Hlong MaxNumPoints, double MaxClosureDist, Hlong ClippingEndPoints, Hlong Iterations, double ClippingFactor, HTuple* Row, HTuple* Column, HTuple* Radius, HTuple* StartPhi, HTuple* EndPhi, HTuple* PointOrder) const

void HXLDCont::FitCircleContourXld(const HString& Algorithm, Hlong MaxNumPoints, double MaxClosureDist, Hlong ClippingEndPoints, Hlong Iterations, double ClippingFactor, double* Row, double* Column, double* Radius, double* StartPhi, double* EndPhi, HString* PointOrder) const

void HXLDCont::FitCircleContourXld(const char* Algorithm, Hlong MaxNumPoints, double MaxClosureDist, Hlong ClippingEndPoints, Hlong Iterations, double ClippingFactor, double* Row, double* Column, double* Radius, double* StartPhi, double* EndPhi, HString* PointOrder) const

void HOperatorSetX.FitCircleContourXld(
[in] IHUntypedObjectX* Contours, [in] VARIANT Algorithm, [in] VARIANT MaxNumPoints, [in] VARIANT MaxClosureDist, [in] VARIANT ClippingEndPoints, [in] VARIANT Iterations, [in] VARIANT ClippingFactor, [out] VARIANT* Row, [out] VARIANT* Column, [out] VARIANT* Radius, [out] VARIANT* StartPhi, [out] VARIANT* EndPhi, [out] VARIANT* PointOrder)

VARIANT HXLDContX.FitCircleContourXld(
[in] BSTR Algorithm, [in] Hlong MaxNumPoints, [in] double MaxClosureDist, [in] Hlong ClippingEndPoints, [in] Hlong Iterations, [in] double ClippingFactor, [out] VARIANT* Column, [out] VARIANT* Radius, [out] VARIANT* StartPhi, [out] VARIANT* EndPhi, [out] VARIANT* PointOrder)

static void HOperatorSet.FitCircleContourXld(HObject contours, HTuple algorithm, HTuple maxNumPoints, HTuple maxClosureDist, HTuple clippingEndPoints, HTuple iterations, HTuple clippingFactor, out HTuple row, out HTuple column, out HTuple radius, out HTuple startPhi, out HTuple endPhi, out HTuple pointOrder)

void HXLDCont.FitCircleContourXld(string algorithm, int maxNumPoints, double maxClosureDist, int clippingEndPoints, int iterations, double clippingFactor, out HTuple row, out HTuple column, out HTuple radius, out HTuple startPhi, out HTuple endPhi, out HTuple pointOrder)

void HXLDCont.FitCircleContourXld(string algorithm, int maxNumPoints, double maxClosureDist, int clippingEndPoints, int iterations, double clippingFactor, out double row, out double column, out double radius, out double startPhi, out double endPhi, out string pointOrder)

Description

fit_circle_contour_xldfit_circle_contour_xldFitCircleContourXldfit_circle_contour_xldFitCircleContourXldFitCircleContourXld approximates the XLD contours ContoursContoursContoursContoursContourscontours by circles. It does not perform a segmentation of the input contours. Thus, one has to make sure that each contour corresponds to one and only one circle. The operator returns for each contour the center (RowRowRowRowRowrow, ColumnColumnColumnColumnColumncolumn), and the RadiusRadiusRadiusRadiusRadiusradius.

The algorithm used for the fitting of circles can be selected via AlgorithmAlgorithmAlgorithmAlgorithmAlgorithmalgorithm:

'algebraic'"algebraic""algebraic""algebraic""algebraic""algebraic"

This approach minimizes the algebraic distance between the contour points and the resulting circle.

'ahuber'"ahuber""ahuber""ahuber""ahuber""ahuber"

Similar to 'algebraic'. Here the contour points are weighted to decrease the impact of outliers based on the approach of Huber (see below).

'atukey'"atukey""atukey""atukey""atukey""atukey"

Similar to 'algebraic'. Here the contour points are weighted to decrease the impact of outliers based on the approach of Tukey (see below).

'geometric'"geometric""geometric""geometric""geometric""geometric"

This approach minimizes the geometric distance between the contour points and the resulting circle. The distance measure is statistically optimal, but takes more computational time. This option is recommended if the contour points are considerably destorted by noise.

'geohuber'"geohuber""geohuber""geohuber""geohuber""geohuber"

Similar to 'geometric'. Here the contour points are weighted to decrease the impact of outliers based on the approach of Huber (see below).

'geotukey'"geotukey""geotukey""geotukey""geotukey""geotukey"

Similar to 'geometric'. Here the contour points are weighted to decrease the impact of outliers based on the approach of Tukey (see below).

For '*huber' and '*tukey' a robust error statistics is used to estimate the standard deviation of the distances from the contour points without outliers from the approximating circle. The parameter ClippingFactorClippingFactorClippingFactorClippingFactorClippingFactorclippingFactor (a scaling factor for the standard deviation) controls the amount of outliers: the smaller the value chosen for ClippingFactorClippingFactorClippingFactorClippingFactorClippingFactorclippingFactor the more outliers are detected. In the Tukey algorithm outliers are removed, whereas in the Huber algorithm outliers are only damped, or more precisely they are weighted linearly. Without any robust weighting the squares of the distances are taken as error values in the optimization, i.e. a least square formulation. In practice, the approach of Tukey is recommended.

The parameter IterationsIterationsIterationsIterationsIterationsiterations specifies the number of iterations for the algorithms 'algebraic', 'ahuber' and 'atukey'. This parameter is ignored for the algorithms 'geometric', 'geohuber' and 'geotukey'. If IterationsIterationsIterationsIterationsIterationsiterations is set to zero, the algorithm does not perform iterative improvements on the fitted circle, but only checks if the initial guess was already close enough depending on the chosen treatment of outliers.

To reduce the computational load, the fitting of circles can be restricted to a subset of the contour points: If a value other than -1 is assigned to MaxNumPointsMaxNumPointsMaxNumPointsMaxNumPointsMaxNumPointsmaxNumPoints, only up to MaxNumPointsMaxNumPointsMaxNumPointsMaxNumPointsMaxNumPointsmaxNumPoints points - uniformly distributed over the contour - are used.

For circular arcs, the points on the circle closest to the start points and end points of the original contours are chosen as start and end points. The corresponding angles refering to the X-axis are returned in StartPhiStartPhiStartPhiStartPhiStartPhistartPhi and EndPhiEndPhiEndPhiEndPhiEndPhiendPhi, see also gen_ellipse_contour_xldgen_ellipse_contour_xldGenEllipseContourXldgen_ellipse_contour_xldGenEllipseContourXldGenEllipseContourXld. Contours, for which the distance between their start points and their end points is less or equal MaxClosureDistMaxClosureDistMaxClosureDistMaxClosureDistMaxClosureDistmaxClosureDist are considered to be closed. Thus, they are approximated by circles instead of circular arcs. Due to artefacts in the pre-processing the start and end points of a contour might be faulty. Therefore, it is possible to exclude ClippingEndPointsClippingEndPointsClippingEndPointsClippingEndPointsClippingEndPointsclippingEndPoints points at the beginning and at the end of a contour from the fitting of circles. However, they are still used for the determination of StartPhiStartPhiStartPhiStartPhiStartPhistartPhi and EndPhiEndPhiEndPhiEndPhiEndPhiendPhi.

The minimun necessary number of contour points for fitting a circle is three. Therefore, it is required that the number of contour points is at least 3 + 2*ClippingEndPointsClippingEndPointsClippingEndPointsClippingEndPointsClippingEndPointsclippingEndPoints .

Parallelization

Parameters

ContoursContoursContoursContoursContourscontours (input_object)  xld_cont(-array) objectHXLDContHXLDContHXLDContHXLDContXHobject

Input contours.

AlgorithmAlgorithmAlgorithmAlgorithmAlgorithmalgorithm (input_control)  string HTupleHTupleHTupleVARIANTHtuple (string) (string) (HString) (char*) (BSTR) (char*)

Algorithm for the fitting of circles.

Default value: 'algebraic' "algebraic" "algebraic" "algebraic" "algebraic" "algebraic"

List of values: 'ahuber'"ahuber""ahuber""ahuber""ahuber""ahuber", 'algebraic'"algebraic""algebraic""algebraic""algebraic""algebraic", 'atukey'"atukey""atukey""atukey""atukey""atukey", 'geohuber'"geohuber""geohuber""geohuber""geohuber""geohuber", 'geometric'"geometric""geometric""geometric""geometric""geometric", 'geotukey'"geotukey""geotukey""geotukey""geotukey""geotukey"

MaxNumPointsMaxNumPointsMaxNumPointsMaxNumPointsMaxNumPointsmaxNumPoints (input_control)  integer HTupleHTupleHTupleVARIANTHtuple (integer) (int / long) (Hlong) (Hlong) (Hlong) (Hlong)

Maximum number of contour points used for the computation (-1 for all points).

Default value: -1

Restriction: MaxNumPoints >= 3

MaxClosureDistMaxClosureDistMaxClosureDistMaxClosureDistMaxClosureDistmaxClosureDist (input_control)  real HTupleHTupleHTupleVARIANTHtuple (real) (double) (double) (double) (double) (double)

Maximum distance between the end points of a contour to be considered as 'closed'.

Default value: 0.0

Restriction: MaxClosureDist >= 0.0

ClippingEndPointsClippingEndPointsClippingEndPointsClippingEndPointsClippingEndPointsclippingEndPoints (input_control)  integer HTupleHTupleHTupleVARIANTHtuple (integer) (int / long) (Hlong) (Hlong) (Hlong) (Hlong)

Number of points at the beginning and at the end of the contours to be ignored for the fitting.

Default value: 0

Restriction: ClippingEndPoints >= 0

IterationsIterationsIterationsIterationsIterationsiterations (input_control)  integer HTupleHTupleHTupleVARIANTHtuple (integer) (int / long) (Hlong) (Hlong) (Hlong) (Hlong)

Maximum number of iterations for the robust weighted fitting.

Default value: 3

Restriction: Iterations >= 0

ClippingFactorClippingFactorClippingFactorClippingFactorClippingFactorclippingFactor (input_control)  real HTupleHTupleHTupleVARIANTHtuple (real) (double) (double) (double) (double) (double)

Clipping factor for the elimination of outliers (typical: 1.0 for Huber and 2.0 for Tukey).

Default value: 2.0

List of values: 1.0, 1.5, 2.0, 2.5, 3.0

Restriction: ClippingFactor > 0

RowRowRowRowRowrow (output_control)  circle.center.y(-array) HTupleHTupleHTupleVARIANTHtuple (real) (double) (double) (double) (double) (double)

Row coordinate of the center of the circle.

ColumnColumnColumnColumnColumncolumn (output_control)  circle.center.x(-array) HTupleHTupleHTupleVARIANTHtuple (real) (double) (double) (double) (double) (double)

Column coordinate of the center of the circle.

RadiusRadiusRadiusRadiusRadiusradius (output_control)  circle.radius(-array) HTupleHTupleHTupleVARIANTHtuple (real) (double) (double) (double) (double) (double)

Radius of circle.

StartPhiStartPhiStartPhiStartPhiStartPhistartPhi (output_control)  real(-array) HTupleHTupleHTupleVARIANTHtuple (real) (double) (double) (double) (double) (double)

Angle of the start point [rad].

EndPhiEndPhiEndPhiEndPhiEndPhiendPhi (output_control)  real(-array) HTupleHTupleHTupleVARIANTHtuple (real) (double) (double) (double) (double) (double)

Angle of the end point [rad].

PointOrderPointOrderPointOrderPointOrderPointOrderpointOrder (output_control)  string(-array) HTupleHTupleHTupleVARIANTHtuple (string) (string) (HString) (char*) (BSTR) (char*)

Point order along the boundary.

List of values: 'negative'"negative""negative""negative""negative""negative", 'positive'"positive""positive""positive""positive""positive"

Result

fit_circle_contour_xldfit_circle_contour_xldFitCircleContourXldfit_circle_contour_xldFitCircleContourXldFitCircleContourXld returns 2 (H_MSG_TRUE) if all parameter values are correct, and circles could be fitted to the input contours. If the input is empty the behaviour can be set via set_system('no_object_result',<Result>)set_system("no_object_result",<Result>)SetSystem("no_object_result",<Result>)set_system("no_object_result",<Result>)SetSystem("no_object_result",<Result>)SetSystem("no_object_result",<Result>). If necessary, an exception is raised. If the parameter ClippingFactorClippingFactorClippingFactorClippingFactorClippingFactorclippingFactor is chosen too small, i.e., all points are classified as outliers, the error 3264 is raised.

Possible Predecessors

gen_contours_skeleton_xldgen_contours_skeleton_xldGenContoursSkeletonXldgen_contours_skeleton_xldGenContoursSkeletonXldGenContoursSkeletonXld, lines_gausslines_gaussLinesGausslines_gaussLinesGaussLinesGauss, lines_facetlines_facetLinesFacetlines_facetLinesFacetLinesFacet, edges_sub_pixedges_sub_pixEdgesSubPixedges_sub_pixEdgesSubPixEdgesSubPix, smooth_contours_xldsmooth_contours_xldSmoothContoursXldsmooth_contours_xldSmoothContoursXldSmoothContoursXld

Possible Successors

gen_ellipse_contour_xldgen_ellipse_contour_xldGenEllipseContourXldgen_ellipse_contour_xldGenEllipseContourXldGenEllipseContourXld, disp_circledisp_circleDispCircledisp_circleDispCircleDispCircle, get_points_ellipseget_points_ellipseGetPointsEllipseget_points_ellipseGetPointsEllipseGetPointsEllipse

See also

fit_ellipse_contour_xldfit_ellipse_contour_xldFitEllipseContourXldfit_ellipse_contour_xldFitEllipseContourXldFitEllipseContourXld, fit_line_contour_xldfit_line_contour_xldFitLineContourXldfit_line_contour_xldFitLineContourXldFitLineContourXld, fit_rectangle2_contour_xldfit_rectangle2_contour_xldFitRectangle2ContourXldfit_rectangle2_contour_xldFitRectangle2ContourXldFitRectangle2ContourXld

Module

Foundation


ClassesClassesClassesClasses | | | | Operators