IPSDK  4_1_0_2
IPSDK : Image Processing Software Development Kit
Classes
Skeleton2dDiameterTortuosity

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

Classes

class  ipsdk::imaproc::shape::analysis::Skeleton2dDiameterTortuosityMsr
 Measurement object for measure Skeleton2dDiameterTortuosity. More...
 
class  ipsdk::imaproc::shape::analysis::Skeleton2dDiameterTortuosityMsrInfo
 Information object for measure Skeleton2dDiameterTortuosity. More...
 
class  ipsdk::imaproc::shape::analysis::Skeleton2dDiameterTortuosityMsrParams
 Parameter object for measure Skeleton2dDiameterTortuosity. More...
 

Detailed Description

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 $L$ of the path to the distance between the ends of it $D$ (Arc-chord ratio) :

\[ \tau = \frac{L}{D} \]

See also
https://en.wikipedia.org/wiki/Tortuosity
Note
Tortuosity equals 1 for a straight line and is infinite for a circle (this case cannot append for a diameter path)
skeleton2dDiameterTortuosityMsr1.png
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.

skeleton2dDiameterTortuosityMsr2.png
Author
E.Noirfalise
Date
2016/08/22

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.png
Generic
none.png
None
parameter.png
Skeleton2dDiameterTortuosityMsrParams
Value.png
Value (ipsdk::ipReal64)
BoundaryApproximation.png
Boundary Approximation
See Shape measurement for additional information on these pictograms

Measure Type :

This is a generic measure

This measure can only be used in 2d case

Measure Unit Type:

Measure Skeleton2dDiameterTortuosity is not associated to any unit [ipsdk::shape::analysis::eMsrUnitFormat::eMUF_NoUnit]

Measure Parameter Type :

Measure Skeleton2dDiameterTortuosity is associated to Skeleton2dDiameterTortuosityMsrParams parameters

Measure Result Type :

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

Measure Shape Requirements :

Measure Skeleton2dDiameterTortuosity requires boundary approximation from shape data

Measure Dependencies :

Measure Skeleton2dDiameterTortuosity depends on following measures :

Measure Mode Measure Name Measure Type Measure Parameters
eMVP_2d Skeleton Skeleton2d createSkeleton2dMsrParams(_pMsrParams->getValue<Skeleton2dDiameterTortuosityMsrParams::HolesPolicy>())
eMVP_2d Diameter Skeleton2dDiameter createSkeleton2dDiameterMsrParams(_pMsrParams->getValue<Skeleton2dDiameterTortuosityMsrParams::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, "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());