Skeleton2dDiameterTortuositySee full documentation

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

Skeleton 2d diameter tortuosity measure is based on skeleton 2d diameter measure (see [Concept] Skeleton 2d Diameter). This measure allows to compute tortuosity of skeleton diameter path.

Tortuosity is in this case defined as the ratio of the length of the path to the distance between the ends of it (Arc-chord ratio) :

Note

Tortuosity equals 1 for a straight line and is infinite for a circle (this case cannot append for a diameter path)

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 tortuosity are taken into account with respect to associated edge length for global tortuosity computation.

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

Tortuosity 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	None	Skeleton2dDiameterTortuosityMsrParams	Value (ipsdk::ipReal64)	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.createGeometricCalibration2d(1,1,'px'))
    PyIPSDK.createMeasureInfo(inMeasureInfoSet2d, "Skeleton2dDiameterMeanCurvatureMsr", shapeanalysis.createSkeleton2dDiameterMeanCurvatureMsrParams(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("Skeleton2dDiameterMeanCurvatureMsr")
    
    # 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, "Skeleton2dDiameterTortuosityMsr", createSkeleton2dDiameterTortuosityMsrParams(skeleton2dHolesPolicy));
    
    // compute measure on shape 2d collection
    MeasureSetPtr pOutMeasureSet = shapeAnalysis2d(pInGreyImg2d, pShape2dColl, pMeasureInfoSet);

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