StructureModelIndex3d

Indicates the relative prevalence of rod and plate structures in three-dimensional binary image.

Access to parameter description

For an introduction:

• section Image Analysis
• section Morphometry

As the fragmentation algorithm, this measurement is based on surface convexity analysis.
This algorithm dilates the input binarized image as shown in the fragmentation index algorithm documentation. The kernelRadius parameter is the number of iterations of dilation performed with a cubic structuring element of size 3.

The Structure Model Index $SMI$ is given by: $$ SMI = 6\times\frac{S_d -S}{S\times \Delta_d}\times\frac{V}{S} $$ where

• $S$ is the object surface area before dilation
• $S_d$ is the object surface area after dilation
• $\Delta_d$ is the radius of dilation structuring element, in calibration unit
• $V$ is the initial, undilated object volume

An ideal plate, cylinder, and sphere have Structure Model Index values of 0, 3 and 4 respectively.

Reference:
T.Hildebrand, P.Ruegsegger. "Quantification of Bone Microarchitecture with the Structure Model Index". Comput Methods Biomech Biomed Engin., vol. 1, pp. 15-23, 1997.

See also

• FragmentationIndex
• AverageObjectThickness3d

// Function prototype.
StructureModelIndexMsr::Ptr
structureModelIndex3d( std::shared_ptr< iolink::ImageView > inputImage,
                       int32_t kernelRadius,
                       StructureModelIndexMsr::Ptr outputMeasurement = NULL );

// Function prototype.
structure_model_index_3d( input_image, kernel_radius = 3, output_measurement = None )

// Function prototype.
public static StructureModelIndexMsr
StructureModelIndex3d( IOLink.ImageView inputImage,
                       Int32 kernelRadius = 3,
                       StructureModelIndexMsr outputMeasurement = null );

// Command constructor.
StructureModelIndex3d();

// Get method of the inputImage parameter.
std::shared_ptr< iolink::ImageView > inputImage() const;
// Set method of the inputImage parameter.
void setInputImage( std::shared_ptr< iolink::ImageView > inputImage );

// Get method of the kernelRadius parameter.
int32_t kernelRadius() const;
// Set method of the kernelRadius parameter.
void setKernelRadius( const int32_t& kernelRadius );

// Get method of the outputMeasurement parameter.
StructureModelIndexMsr::Ptr outputMeasurement() const;

// Method to launch the command.
void execute();

// Property of the inputImage parameter.
StructureModelIndex3d.input_image
// Property of the kernelRadius parameter.
StructureModelIndex3d.kernel_radius
// Property of the outputMeasurement parameter.
StructureModelIndex3d.output_measurement

// Method to launch the command.
execute()

// Command constructor.
StructureModelIndex3d()

// Property of the inputImage parameter.
StructureModelIndex3d.inputImage
// Property of the kernelRadius parameter.
StructureModelIndex3d.kernelRadius
// Property of the outputMeasurement parameter.
StructureModelIndex3d.outputMeasurement

// Method to launch the command.
Execute()

auto foam_sep = readVipImage( std::string( IMAGEDEVDATA_IMAGES_FOLDER ) + "foam_sep.vip" );

StructureModelIndex3d structureModelIndex3dAlgo;
structureModelIndex3dAlgo.setInputImage( foam_sep );
structureModelIndex3dAlgo.setKernelRadius( 3 );
structureModelIndex3dAlgo.execute();

std::cout << "index: " << structureModelIndex3dAlgo.outputMeasurement()->index( ) ;

foam_sep = imagedev.read_vip_image(imagedev_data.get_image_path("foam_sep.vip"))

structure_model_index_3d_algo = imagedev.StructureModelIndex3d()
structure_model_index_3d_algo.input_image = foam_sep
structure_model_index_3d_algo.kernel_radius = 3
structure_model_index_3d_algo.execute()

print( "index: ", str( structure_model_index_3d_algo.output_measurement.index( ) ) )

ImageView foam_sep = Data.ReadVipImage( @"Data/images/foam_sep.vip" );

StructureModelIndex3d structureModelIndex3dAlgo = new StructureModelIndex3d
{
    inputImage = foam_sep,
    kernelRadius = 3
};
structureModelIndex3dAlgo.Execute();

Console.WriteLine( "index: " + structureModelIndex3dAlgo.outputMeasurement.index( ) );

auto foam_sep = readVipImage( std::string( IMAGEDEVDATA_IMAGES_FOLDER ) + "foam_sep.vip" );

auto result = structureModelIndex3d( foam_sep, 3 );

std::cout << "index: " << result->index( ) ;

foam_sep = imagedev.read_vip_image(imagedev_data.get_image_path("foam_sep.vip"))

result = imagedev.structure_model_index_3d( foam_sep, 3 )

print( "index: ", str( result.index( ) )  )

ImageView foam_sep = Data.ReadVipImage( @"Data/images/foam_sep.vip" );

StructureModelIndexMsr result = Processing.StructureModelIndex3d( foam_sep, 3 );

Console.WriteLine(  "index: " + result.index( )  );