Model fitting, acces at the coefficient value

[Cedric Hemom]

Hi,

I use my SSM to the model fitting and I want to know the coefficient assigned by the gp.posterior method.

I can't figure out how I can access the coefficient values of the modes in my gp.

Thanks for your help.

Sincerely,

Cedric.

[Marcel Luethi]

HI Cedric

The mode (singular) of your posterior has the coefficients (0, ... 0). The posterior is a Gaussian distribution and therefore has only one mode.

In general, you can access the coefficient of any sample of your posterior, by using the coefficients method of the GP as follows:

val posterior: PointDistributionModel[_3D, TriangleMesh] = ???
val sample: TriangleMesh[_3D] = posterior.sample()
val coeffs = posterior.coefficients(sample)

Best regards

Marcel

--
You received this message because you are subscribed to the Google Groups "scalismo" group.
To unsubscribe from this group and stop receiving emails from it, send an email to scalismo+u...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/scalismo/c83c6ca8-a878-4632-aeab-9b24f2c2821an%40googlegroups.com.

[Cedric Hemom]

Ah ok,

I would like to use the modes of my SSM to fit the other meshes. I don't know if such a solution is possible in scalismo. 

My goal is to know the coefficients on the modes assigned for each mesh to use them as ground truth.

Thank you for your answer.

Sincerely,

Cedric.

[Marcel Luethi]

Dear Cedric

Maybe I am not exactly understanding your question. Do you mean by modes the point of the distribution with maximum probability or are you referring to the modes of variation?

I assumed the first in my last answer.

If you mean the latter, you can use the posterior distribution that you obtain from the Gaussian process regression in place of your initial model. It is a Gaussian process model  itself and so can be used as a drop in replacement in any place where a model is required.

Best regards,

Marcel

To view this discussion on the web visit https://groups.google.com/d/msgid/scalismo/1052b438-bff2-46a4-b90e-8ab4b3ea1a94n%40googlegroups.com.

[Cedric Hemom]

I was referring to the modes of variation. The problem is that once I have built my SSM from the mapped data, I would like to find the coefficients on the modes of variation to reconstruct the shapes of my dataset. 

I don't know if this optimization is present in scalismo. Unfortunately I have to keep my initial model and I cannot use each time the new model obtained from the gaussian process.

Thanks for your help.

Sincerely,

Cédric.

[Marcel Luethi]

Dear Cedric,

Yes, this is possible. As the data is in correspondence, you can use the coefficients method on the dataset:

val model : PointDistributionModel[_3D, TriangleMesh] = ???
val mesh : TriangleMesh[_3D] = ???

val coeffs = model.coefficients(mesh)
model.instance(coeffs)  // should reproduce the mesh 

Best regards,

Marcel

To view this discussion on the web visit https://groups.google.com/d/msgid/scalismo/0c398175-6266-479e-8598-b0e122198985n%40googlegroups.com.

[Cedric Hemom]

Ok , Thank you very much.

I will try this method.

Sincerely,

Cedric.

[Cedric Hemom]

Is it possible to change the mean of the an SSM from scalismo?

I need to change the mean shape by applying just a rigid transformation but keep the original deformation field.

[Marcel Luethi]

Hi Cedric,

Yes this is possible by changing the underlying reference mesh. You can do it as follows:

val pdm : PointDistributionModel[_3D, TriangleMesh] = ???
val newReferenceMesh : TriangleMesh[_3D] = ???
val pdmWithNewMean = pdm.newReference(newReferenceMesh, TriangleMeshInterpolator3D())

This will not only change the reference mesh, but also adjust the mean such that the represented distribution remains the same.

Best regards, 

Marcel

To view this discussion on the web visit https://groups.google.com/d/msgid/scalismo/be7fc475-904f-41bb-a787-8700e1cf79c0n%40googlegroups.com.

[Cedric Hemom]

It works but it is not perfect. I tried using the changeReference function and applying a rigid transformation and I can't change the mean shape. 

I don't know if there is a function that can directly change just the mean shape because the newReference function gives me approximate results.

[Andreas Morel-Forster]

Hi Cedric

I am not sure, but maybe you are looking for the pdm.transform( ... ) function. This will move your PDM in space to a different location while keeping the statistical shape space the same. You should then observe the same mean shape at a different location in space.

The method newReference instead only changes the reference, i.e. the domain on which the PDM is defined. This method could be used if you want to remesh your reference and transfer the statistical shape space to your new remeshed reference/domain.

Best, Andreas

To view this discussion on the web visit https://groups.google.com/d/msgid/scalismo/a686f6c7-0622-4c98-b1e0-752755dd78a2n%40googlegroups.com.

[Cedric Hemom]

Thanks for the explanation of the newReference method.

I misread the documentation and did not see that the transform method was defined for the PDM class.
Sorry, for this waste of time and thank you very much always for your quick help.

Best, 

Cedric.