Skeleton2dDiameterMeanCurvature

Mean curvature of the longest shortest path between any two skeleton 2d verticies. More...

Classes
class	ipsdk::imaproc::shape::analysis::Skeleton2dDiameterMeanCurvatureMsr
	Measurement object for measure Skeleton2dDiameterMeanCurvature. More...
class	ipsdk::imaproc::shape::analysis::Skeleton2dDiameterMeanCurvatureMsrInfo
	Information object for measure Skeleton2dDiameterMeanCurvature. More...
class	ipsdk::imaproc::shape::analysis::Skeleton2dDiameterMeanCurvatureMsrParams
	Parameter object for measure Skeleton2dDiameterMeanCurvature. More...

Detailed Description

Mean curvature of the longest shortest path between any two skeleton 2d verticies.

Skeleton 2d diameter leangth measure is based on skeleton 2d diameter measure (see [Concept] Skeleton 2d Diameter). This measure allows to compute mean curvature along skeleton diameter path.

Mean curvature is in this case defined for 3 consecutive points as the reciprocal of circle radius passing through these points :

Note

Note that measure results depend on skeleton 2d holes processing policy See [Concept] Skeleton 2d for more informations on usage of ipsdk::shape::segmentation::eSkeleton2dHolesPolicy enumerate value. If holes processing policy is set to eSHP_ExtractedAndAnalyzed, holes skeleton diameter mean curvature are taken into account with respect to associated edge length for global mean curvature computation.

Here is an example of measure results for different skeleton 2d holes processing policy.

See also

https://en.wikipedia.org/wiki/Curvature

Author

E.Noirfalise

Date

2016/08/22

Copyright

2016 Reactiv'IP, all right reserved.

Mean curvature of the longest shortest path between any two skeleton 2d verticies

Measure synthesis :

Measure Type	Measure Unit Type	Parameter Type	Result Type	Shape Requirements
Generic	Custom	Skeleton2dDiameterMeanCurvatureMsrParams	Value (ipsdk::ipReal64)	Boundary Approximation

See Shape measurement for additional information on these pictograms

Measure Type :

This is a generic measure

• it does not directly use shape or image data
• it can use dependencies to other measures to compute a result (see Measure Dependencies section)

This measure can only be used in 2d case

Measure Unit Type:

Measure Skeleton2dDiameterMeanCurvature is associated to a custom unit, please see definition of Skeleton2dDiameterMeanCurvatureInfo::computeCustomMsrUnitStr function [ipsdk::shape::analysis::eMsrUnitFormat::eMUF_CustomUnit]

Measure Parameter Type :

Measure Skeleton2dDiameterMeanCurvature is associated to Skeleton2dDiameterMeanCurvatureMsrParams parameters

Measure Result Type :

Measure Skeleton2dDiameterMeanCurvature is associated to ipsdk::shape::analysis::ValueMeasureResult<ipsdk::ipReal64> results

Measure Shape Requirements :

Measure Skeleton2dDiameterMeanCurvature requires boundary approximation from shape data

Measure Dependencies :

Measure Skeleton2dDiameterMeanCurvature depends on following measures :

Measure Mode	Measure Name	Measure Type	Measure Parameters
eMVP_2d	Skeleton	Skeleton2d	createSkeleton2dMsrParams(_pMsrParams->getValue<Skeleton2dDiameterMeanCurvatureMsrParams::HolesPolicy>())
eMVP_2d	Diameter	Skeleton2dDiameter	createSkeleton2dDiameterMsrParams(_pMsrParams->getValue<Skeleton2dDiameterMeanCurvatureMsrParams::HolesPolicy>())

Note

See Shape Analysis 2d for more information on general shape 2d analysis and measurement usage.

Example of Python code :

Generic example in 2d case :

import PyIPSDK
import PyIPSDK.IPSDKIPLShapeAnalysis as shapeanalysis

    # Create the infoset
    inMeasureInfoSet2d = PyIPSDK.createMeasureInfoSet2d(PyIPSDK.createGeometricCalibration2d(1,1,'px'))
    PyIPSDK.createMeasureInfo(inMeasureInfoSet2d, "Skeleton2dDiameterTortuosityMsr", shapeanalysis.createSkeleton2dDiameterTortuosityMsrParams(PyIPSDK.eSHP_Ignored))
    
    #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("Skeleton2dDiameterTortuosityMsr")
    
    # 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>();
    IPSDK_REQUIRE(readFromXmlFile(inputShape2dCollPath, *pShape2dColl) == true);

    // define a measure info set
    MeasureInfoSetPtr pMeasureInfoSet = MeasureInfoSet::create2dInstance();
    createMeasureInfo(pMeasureInfoSet, "Skeleton2dDiameterMeanCurvatureMsr", createSkeleton2dDiameterMeanCurvatureMsrParams(skeleton2dHolesPolicy));
    
    // compute measure on shape 2d collection
    MeasureSetPtr pOutMeasureSet = shapeAnalysis2d(pInGreyImg2d, pShape2dColl, pMeasureInfoSet);

    // retrieve associated results
    const MeasureConstPtr& pSkeleton2dDiameterMeanCurvatureOutMsr = pOutMeasureSet->getMeasure("Skeleton2dDiameterMeanCurvatureMsr");
    const ipsdk::shape::analysis::ValueMeasureResult<ipsdk::ipReal64>& outResults = static_cast<const ipsdk::shape::analysis::ValueMeasureResult<ipsdk::ipReal64>&>(pSkeleton2dDiameterMeanCurvatureOutMsr->getMeasureResult());