Measure allowing to compute the barycenter ponderated by the gray level of each pixel for shape.

The coordinates $\left[ x_b, y_b \right]^T$ of the grey barycenter for a given 2d or 3d shape are computed as follows :

$\left\lbrace \begin{array}{ccc} x_b & = & \frac{1}{nbElt} \left( \sum_{\left\{x, y, z\right\} \in Shape}{Shape(x, y, z) \times x} \right)\\ y_b & = & \frac{1}{nbElt} \left( \sum_{\left\{x, y, z\right\} \in Shape}{Shape(x, y, z) \times y} \right)\\ z_b & = & \frac{1}{nbElt} \left( \sum_{\left\{x, y, z\right\} \in Shape}{Shape(x, y, z) \times z} \right) \end{array} \right.$

Where $ Shape $ is the set of pixels in the shape and $ nbElt $ is the number of elements (i.e. pixels in 2d case or voxels in 3d case) of the shape.

Here is an example of 2d grey barycenter measurement :

Measure allowing to compute the barycenter ponderated by the gray level of each pixel for shape

Measure synthesis :

Measure Type	Measure Unit Type	Parameter Type	Result Type	Shape Requirements
Intensity	None	None	Custom	Row Intersections

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, "GreyBarycenterMsr")
    
    #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("GreyBarycenterMsr")
    
    # retrieve measure values
    outMsrValues = outMsr.getMeasureResult().getColl(0)
    print("First label measurement equal " + str(outMsrValues[1]))

Generic example in 3d case :

import PyIPSDK

import PyIPSDK.IPSDKIPLShapeAnalysis as shapeanalysis

    # Create the infoset
    inMeasureInfoSet3d = PyIPSDK.createMeasureInfoSet3d()
    PyIPSDK.createMeasureInfo(inMeasureInfoSet3d, "GreyBarycenterMsr")
    
    #Perform the analysis
    outMeasureSet = shapeanalysis.labelAnalysis3d(inGreyImg, inLabelImg, inMeasureInfoSet3d)
    
    # save results to csv format
    PyIPSDK.saveCsvMeasureFile(os.path.join(tmpPath,  "shape_analysis_results.csv"), outMeasureSet)
    
    # retrieve measure results
    outMsr = outMeasureSet.getMeasure("GreyBarycenterMsr")
    
    # retrieve measure values
    outMsrValues = outMsr.getMeasureResult().getColl(0)
    print("First label measurement equal " + str(outMsrValues[1]))

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>();
    readFromXmlFile(inputShape2dCollPath, *pShape2dColl);
    // define a measure info set
    MeasureInfoSetPtr pMeasureInfoSet = MeasureInfoSet::create2dInstance();
    createMeasureInfo(pMeasureInfoSet, "GreyBarycenterMsr");
    
    // compute measure on shape 2d collection
    MeasureSetPtr pOutMeasureSet = shapeAnalysis2d(pInGreyImg2d, pShape2dColl, pMeasureInfoSet);
    // retrieve associated results
    const MeasureConstPtr& pGreyBarycenterOutMsr = pOutMeasureSet->getMeasure("GreyBarycenterMsr");
    const GreyBarycenterMsrResults& outResults = static_cast<const GreyBarycenterMsrResults&>(pGreyBarycenterOutMsr->getMeasureResult());