Creating a Statistical Shape Model

we can generate a statistical shape model from the algorithms.shapeModelling.CylindricalPointModel method

import sys
import pymethods as pm
import pyvista as pv

model = pm.algorithms.shapeModelling.CylindricalPointModel(
r'../../Datasets\shape_modelling_LADs', structured_shape=(3, 360, 720))
model.build(scale=True)

we start by specifying the folder where the shape vtks are stored and the desired structured shape (dims, azimuthal, longitudinal). We then build the model by calling the build method. When we set scaling to true we ensure that the statistical models have the same mean inlet diameter as the original areteries. To obtain the explained variance we can call the explained_variance method of the model

print("the fraction of explained variance is:", model.explained_variance())

to generate a new shape we supply the modes to the models generateShape method

p = pv.BackgroundPlotter()
shape = model.generateShape([0, -1, 0])
mesh = pv.StructuredGrid(*shape)
p.add_mesh(mesh)
p.show()

p = pv.BackgroundPlotter()
shape = model.generateShape([3, 2.5, -2, -2, 0.1, 2,1,3,1,3,1,2])
mesh = pv.StructuredGrid(*shape)
p.add_mesh(mesh)
p.show()