SurfaceProjectionY

Measure that computes the surface of the projection along the Y axis.

This measure returns the surface of the projection along the Y axis for each 3d shapes. The surface is computed on the XZ plane.

Here is an example of the surface projection Y measurement :

Measure the surface of the projection along the Y axis

Measure synthesis :

Measure Type	Measure Unit Type	Parameter Type	Result Type	Shape Requirements
Geometry 3d	None	None	Value (ipsdk::ipUInt64)	Row Intersections

See Shape measurement for additional information on these pictograms

Example of Python code :

Generic example in 3d case :

import PyIPSDK
import PyIPSDK.IPSDKIPLShapeAnalysis as shapeanalysis

    # Create the infoset
    inMeasureInfoSet3d = PyIPSDK.createMeasureInfoSet3d()
    PyIPSDK.createMeasureInfo(inMeasureInfoSet3d, "SurfaceProjectionYMsr")
    
    #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("SurfaceProjectionYMsr")
    
    # 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

    //ImagePtr pInGreyImg3d = loadTiffImageFile(inputGreyImgPath);

    //Shape3dCollPtr pShape3dColl = boost::make_shared<Shape3dColl>();
    //IPSDK_REQUIRE(readFromXmlFile(inputShape3dCollPath, *pShape3dColl) == true);

    ImageGeometryPtr pImageGeometry = geometry3d(eImageBufferType::eIBT_Label16, 100, 100, 50);
    boost::shared_ptr<MemoryImage> pImg(boost::make_shared<MemoryImage>());
    pImg->init(*pImageGeometry);

    FormulaInOptImages inOptImages1;
    inOptImages1._pInOptImg1 = pImg;
    ImagePtr pOutBinImg = formula3dImg("if (x>5&&x<40&&y>10&&y<30&&z>2&&z<15, True, if (x>16&&x<36&&y>40&&y<70&&z>5&&z<25, True, if (x>1&&x<51&&y>50&&y<60&&z>5&&z<25, True, if (sqrt((x-30-cx)*(x-30-cx)+(y-20-cy)*(y-20-cy)+(z-cz)*(z-cz))<=20, True, False))))", inOptImages1);

    ImagePtr pInLabelImg = connectedComponent3dImg(pOutBinImg);

    Shape3dCollPtr pShape3dColl = labelShapeExtraction3d(pInLabelImg);

    // define a measure info set
    MeasureInfoSetPtr pMeasureInfoSet = MeasureInfoSet::create3dInstance();
    createMeasureInfo(pMeasureInfoSet, "SurfaceProjectionYMsr");
    
    // compute measure on shape 3d collection
    MeasureSetPtr pOutMeasureSet = shapeAnalysis3d(pInLabelImg, pShape3dColl, pMeasureInfoSet);

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