least square circle associated to shape tops

This measure allows to compute least square circle associated to each shape.

Given a shape, algorithm will compute a least square circle approximation of shape tops :

Estimated circle can be computed using a simple parametric estimation of shape tops or can be computed using a robust estimation technic. This behavior is ruled by measure parameter, see ipsdk::imaproc::attr::RobustEstimationConfig.

For more information on parametric estimation process (including outliers management) please reports to Parametric estimation.

Here is an example of the least square circle computation on a grey scale image:

least square circle associated to shape tops

Measure synthesis :

Measure Type	Measure Unit Type	Parameter Type	Result Type	Shape Requirements
Geometry 2d	None	LeastSquareCircle2dMsrParams	Custom	Boundary Approximation

See Shape measurement for additional information on these pictograms

Example of Python code :

Generic example in 2d case :

import PyIPSDK

import PyIPSDK.IPSDKIPLShapeAnalysis as shapeanalysis

    # Create the infoset
    inMeasureInfoSet2d = PyIPSDK.createMeasureInfoSet2d()
    PyIPSDK.createMeasureInfo(inMeasureInfoSet2d, "LeastSquareCircle2dMsr")
    
    #Perform the analysis
    outMeasureSet = shapeanalysis.labelAnalysis2d(inGreyImg, inLabelImg2d, inMeasureInfoSet2d)
    
    # save results to csv format
    PyIPSDK.saveCsvMeasureFile(os.path.join(tmpPath,  "shape_analysis_results.csv"), outMeasureSet)
    
    # retrieve measure results
    outMsr = outMeasureSet.getMeasure("LeastSquareCircle2dMsr")
    
    # retrieve measure values
    outMsrValues = outMsr.getMeasureResult().getColl(0)
    print("First label estimation radius : " + str(outMsrValues[1].first.getRadius()))
    print("First label estimation center : " + str(outMsrValues[1].first.centerX()) + ", " + str(outMsrValues[1].first.centerY()))

Example of C++ code :

Example informations

Associated library

IPSDKIPLShapeAnalysis

Code Example

    // opening grey level input image
    ImagePtr pInGreyImg2d = loadTiffImageFile(inputGreyImgPath);
    // read entity shape 2d collection used for processing
    Shape2dCollPtr pShape2dColl = boost::make_shared<Shape2dColl>();
    IPSDK_REQUIRE(readFromXmlFile(inputShape2dCollPath, *pShape2dColl) == true);
    // define a measure info set
    MeasureInfoSetPtr pMeasureInfoSet = MeasureInfoSet::create2dInstance();
    createMeasureInfo(pMeasureInfoSet, "LeastSquareCircle2dMsr");
    
    // compute measure on shape 2d collection
    MeasureSetPtr pOutMeasureSet = shapeAnalysis2d(pInGreyImg2d, pShape2dColl, pMeasureInfoSet);
    // retrieve associated results
    const MeasureConstPtr& pLeastSquareCircle2dOutMsr = pOutMeasureSet->getMeasure("LeastSquareCircle2dMsr");
    const LeastSquareCircle2dMsrResults& outResults = static_cast<const LeastSquareCircle2dMsrResults&>(pLeastSquareCircle2dOutMsr->getMeasureResult());