# Cuberille Mesh
**URL:** https://discourse.itk.org/t/cuberille-mesh/6130
**Category:** Beginner Questions
**Created:** [August 7, 2023, 11:47pm UTC](https://discourse.itk.org/t/cuberille-mesh/6130 "2023-08-07T23:47:50Z")
**Posts on this page:** 14
**Page:** 1
### Author:  [@CrownMella](https://discourse.itk.org/u/CrownMella)
#### Post date: [August 7, 2023, 11:47pm UTC](https://discourse.itk.org/t/cuberille-mesh/6130/1 "2023-08-07T23:47:50Z")
Hi,
Is there a way to convert the mesh from cuberille image to mesh filter back to an image? I saw something like itk.TriangleMeshToBinaryImageFilter but my python brain doesnt understand how
instance = itk.TriangleMeshToBinaryImageFilter[itk.Mesh[itk.UC,3], itk.Image[itk.UC,3]].New()
works…is it like
image\_from\_mesh = instance(mesh) because then it says
ITK ERROR: TriangleMeshToBinaryImageFilter(0000016AD471C630): Must Set Image Size
---
### Author:  [@dzenanz](https://discourse.itk.org/u/dzenanz)
#### Post date: [August 8, 2023, 2:35pm UTC](https://discourse.itk.org/t/cuberille-mesh/6130/2 "2023-08-08T14:35:40Z")
Mesh does not explicitly contain the resolution of the image from which it originated. And implicitly determining it is hard. So you need to specify the image grid parameters (size, spacing, origin, orientation).
---
### Author:  [@dzenanz](https://discourse.itk.org/u/dzenanz)
#### Post date: [August 8, 2023, 2:36pm UTC](https://discourse.itk.org/t/cuberille-mesh/6130/3 "2023-08-08T14:36:15Z")
If you have the original image, you could take the information from it.
---
### Author:  [@CrownMella](https://discourse.itk.org/u/CrownMella)
#### Post date: [August 8, 2023, 9:53pm UTC](https://discourse.itk.org/t/cuberille-mesh/6130/4 "2023-08-08T21:53:16Z")
I got it to work, thanks for the help!
---
### Author:  [@dzenanz](https://discourse.itk.org/u/dzenanz)
#### Post date: [August 9, 2023, 1:13pm UTC](https://discourse.itk.org/t/cuberille-mesh/6130/5 "2023-08-09T13:13:51Z")
Perhaps share your solution in case someone in the future runs into the same issue?
---
### Author:  [@CrownMella](https://discourse.itk.org/u/CrownMella)
#### Post date: [August 9, 2023, 4:47pm UTC](https://discourse.itk.org/t/cuberille-mesh/6130/6 "2023-08-09T16:47:49Z")
Sure, no problem! I was trying to transform the 3D segmentations in physical space so this was my way of doing it. Please let me know if there are other ways!
```auto
import itk
MeshType = itk.Mesh[itk.UC,3]
ImageType = itk.Image[itk.UC, 3]
PixelType = itk.ctype("unsigned char")
segmentation= itk.imread(path, PixelType)
inv = itk.transformread(invpath)
size = itk.size(segmentation)
origin = segmentation.GetOrigin()
spacing = segmentation.GetSpacing()
direction = segmentation.GetDirection()
mesh = itk.cuberille_image_to_mesh_filter(segmentation)
transformed_mesh = itk.transform_mesh_filter(
mesh,
transform=inv
)
mesh_to_image_filter = itk.TriangleMeshToBinaryImageFilter[MeshType, ImageType].New()
mesh_to_image_filter.SetInput(transformed_mesh)
mesh_to_image_filter.SetSize(size)
mesh_to_image_filter.SetOrigin(origin)
mesh_to_image_filter.SetSpacing(spacing)
mesh_to_image_filter.SetDirection(direction)
mesh_to_image_filter.Update()
tf_image_from_mesh = mesh_to_image_filter.GetOutput()
```
---
### Author:  [@dzenanz](https://discourse.itk.org/u/dzenanz)
#### Post date: [August 9, 2023, 5:41pm UTC](https://discourse.itk.org/t/cuberille-mesh/6130/7 "2023-08-09T17:41:38Z")
If your transform is rigid or affine, you can do this operation much more elegantly via [TransformGeometryImageFilter](https://itk.org/Doxygen/html/classitk_1_1TransformGeometryImageFilter.html), not to mention there is no loss and it is nearly instantaneous compared to cuberille filter workaround.
---
### Author:  [@dzenanz](https://discourse.itk.org/u/dzenanz)
#### Post date: [August 9, 2023, 5:44pm UTC](https://discourse.itk.org/t/cuberille-mesh/6130/8 "2023-08-09T17:44:34Z")
If your transform is not affine, I would still recommend the image resampling approach, without a diversion into mesh. You might want a fancy interpolator such as [Label Gaussian](https://itk.org/Doxygen/html/classitk_1_1LabelImageGaussianInterpolateImageFunction.html) instead of the usual nearest which is the normal choice for label resampling.
---
### Author:  [@CrownMella](https://discourse.itk.org/u/CrownMella)
#### Post date: [August 10, 2023, 4:13pm UTC](https://discourse.itk.org/t/cuberille-mesh/6130/9 "2023-08-10T16:13:37Z")
Hi Dzenanz,
Thanks for these two links! I still have yet to understand itk stuff coming from sitk. As a follow up, I get the concerns for time efficiency of my code but where would the loss (what is loss also) come from?
Also, I’m not sure if this is what you mean but I did try the SITK resampling to my images but when I converted them to mesh files, they looked quite different from what I would get say from applying a transform on the segmentations/meshes via 3DSlicer.
---
### Author:  [@dzenanz](https://discourse.itk.org/u/dzenanz)
#### Post date: [August 10, 2023, 5:16pm UTC](https://discourse.itk.org/t/cuberille-mesh/6130/10 "2023-08-10T17:16:32Z")
Conversion `mesh` ↔ `label image` is not fully reversible in general case, so these conversions should be avoided when possible.
> SITK resampling to my images but when I converted them to mesh files, they looked quite different from what I would get say from applying a transform on the segmentations/meshes via 3DSlicer.
If they were quite different, I suspect you were not doing it right. When resampling an image, direction of transform is opposite (inverted) from the one used to transform points (mesh). Try resampling with an inverse of the transform you were originally trying.
---
### Author:  [@CrownMella](https://discourse.itk.org/u/CrownMella)
#### Post date: [August 10, 2023, 5:46pm UTC](https://discourse.itk.org/t/cuberille-mesh/6130/11 "2023-08-10T17:46:38Z")
Hi dzenanz,
The differences were not drastic but they were noticeable at least (the resampled image turned mesh looked grainier)! I used the transform produced by registering the two images for the resampling. I inverted this transform for the meshes (for 3dslicer I did not need to do this).
---
### Author:  [@dzenanz](https://discourse.itk.org/u/dzenanz)
#### Post date: [August 10, 2023, 6:04pm UTC](https://discourse.itk.org/t/cuberille-mesh/6130/12 "2023-08-10T18:04:07Z")
> [@CrownMella](#):
>
> the resampled image turned mesh looked grainier
It sounds like you might have chosen lower resolution for resampling, or you like the smoothing implicit in cuberille converter. If you resample to an image of higher resolution, and/or use smooth label interpolator, how do the results compare then?
---
### Author:  [@CrownMella](https://discourse.itk.org/u/CrownMella)
#### Post date: [August 10, 2023, 6:12pm UTC](https://discourse.itk.org/t/cuberille-mesh/6130/13 "2023-08-10T18:12:33Z")
Apologies, my terminology is a bit weird. When I say grainier, I mean two labels can get interspersed with one another whereas in the mesh it does not. This was using SITK.resample so I may have to go up a level and do it by learning and using itk stuff to achieve a better result? (like the fancy interpolator you mentioned)
---
### Author:  [@dzenanz](https://discourse.itk.org/u/dzenanz)
#### Post date: [August 10, 2023, 6:41pm UTC](https://discourse.itk.org/t/cuberille-mesh/6130/14 "2023-08-10T18:41:51Z")
If labels get interspersed, that is probably an artifact of label interpolation. Try using `nearest` interpolator instead of the default `linear`.