Hello,
I have measured some points on some piles of dirt that I need to calculate the volumes of. I have got some points by the foot of the piles and some on top of them to define their geometries. Now is the question how to use this data to get the volumes? I have managed to get a TIN surface from the TIN Generator but I cannot calculate the volume unless it is a solid, right? Any ideas..?
Best answer by takashi
Hi @tobiasp, I looked at the TINVolumeCalculator. This transformer seems to always require a feature having a boundary area geometry from the Boundary_INPUT. If you already have a TIN surface and need to calculate its volume,
The transformer should work with the setting above, but in my quick test, the resulting volume was wrong. The reason is that there are these fatal bugs in the implementation :-0
I therefore would not recommend you to use the TINVolumeCalculator transformer. @trentatsafe, please check carefully the implementation of the TINVolumeCalculator.
This is a workaround using the BRepSolidBoundaryCreator.Beta (from FME Hub) to create a solid below (minimum z and above) the input TIN surface.
Alternatively, this small Python script can be used to directly compute the volume below the input Mesh surface (minimum z and above).
# PythonCaller Script Example
# Assume the input feature has a Mesh geometry (typically TIN).
import fmeobjects
def calculateMeshVolume(feature):
mesh = feature.getGeometry()
if isinstance(mesh, fmeobjects.FMEMesh):
zmin = mesh.boundingCube()[0][2]
pool = mesh.getVertices()
volume = 0.0
for part in mesh:
a, h = 0.0, 0.0
vertices = [pool[i] for i in part.getVertexIndices()]
for i, (x0, y0, z0) in enumerate(vertices[:-1]):
x1, y1, _ = vertices[i+1]
a += (x1 - x0) * (y0 + y1)
h += (z0 - zmin)
volume += abs(a) * 0.5 * (h / (len(vertices) - 1))
feature.setAttribute('_volume', volume)
+5Have you tried the TINVolumeCalculator from FME Hub?
https://hub.safe.com/transformers/tinvolumecalculator
I'm not sure you can use this directly, but, if you edit it, it might give you some pointers on how to solve your particular solution.
+6Hello @tobiasp
To extrapolate off what @courtney_m stated.
The TINVolumeCalculator should provide you with the tool to produce the volume. The transformer itself will calculate the volume of each individual triangle in the TIN down to a reference elevation. It then sums all the triangles together to get the total.
You also have the option of using a boundary polygon which can be used in the Digital Terrain Model generation. This would allow you to use the same data(pile of dirt) and model changes over time.
Hello @tobiasp
To extrapolate off what @courtney_m stated.
The TINVolumeCalculator should provide you with the tool to produce the volume. The transformer itself will calculate the volume of each individual triangle in the TIN down to a reference elevation. It then sums all the triangles together to get the total.
You also have the option of using a boundary polygon which can be used in the Digital Terrain Model generation. This would allow you to use the same data(pile of dirt) and model changes over time.
Hello @tobiasp
To extrapolate off what @courtney_m stated.
The TINVolumeCalculator should provide you with the tool to produce the volume. The transformer itself will calculate the volume of each individual triangle in the TIN down to a reference elevation. It then sums all the triangles together to get the total.
You also have the option of using a boundary polygon which can be used in the Digital Terrain Model generation. This would allow you to use the same data(pile of dirt) and model changes over time.
+19Hi @tobiasp, I looked at the TINVolumeCalculator. This transformer seems to always require a feature having a boundary area geometry from the Boundary_INPUT. If you already have a TIN surface and need to calculate its volume,
The transformer should work with the setting above, but in my quick test, the resulting volume was wrong. The reason is that there are these fatal bugs in the implementation :-0
I therefore would not recommend you to use the TINVolumeCalculator transformer. @trentatsafe, please check carefully the implementation of the TINVolumeCalculator.
This is a workaround using the BRepSolidBoundaryCreator.Beta (from FME Hub) to create a solid below (minimum z and above) the input TIN surface.
Alternatively, this small Python script can be used to directly compute the volume below the input Mesh surface (minimum z and above).
# PythonCaller Script Example
# Assume the input feature has a Mesh geometry (typically TIN).
import fmeobjects
def calculateMeshVolume(feature):
mesh = feature.getGeometry()
if isinstance(mesh, fmeobjects.FMEMesh):
zmin = mesh.boundingCube()[0][2]
pool = mesh.getVertices()
volume = 0.0
for part in mesh:
a, h = 0.0, 0.0
vertices = [pool[i] for i in part.getVertexIndices()]
for i, (x0, y0, z0) in enumerate(vertices[:-1]):
x1, y1, _ = vertices[i+1]
a += (x1 - x0) * (y0 + y1)
h += (z0 - zmin)
volume += abs(a) * 0.5 * (h / (len(vertices) - 1))
feature.setAttribute('_volume', volume)
Hi @tobiasp, I looked at the TINVolumeCalculator. This transformer seems to always require a feature having a boundary area geometry from the Boundary_INPUT. If you already have a TIN surface and need to calculate its volume,
The transformer should work with the setting above, but in my quick test, the resulting volume was wrong. The reason is that there are these fatal bugs in the implementation :-0
I therefore would not recommend you to use the TINVolumeCalculator transformer. @trentatsafe, please check carefully the implementation of the TINVolumeCalculator.
This is a workaround using the BRepSolidBoundaryCreator.Beta (from FME Hub) to create a solid below (minimum z and above) the input TIN surface.
Alternatively, this small Python script can be used to directly compute the volume below the input Mesh surface (minimum z and above).
# PythonCaller Script Example
# Assume the input feature has a Mesh geometry (typically TIN).
import fmeobjects
def calculateMeshVolume(feature):
mesh = feature.getGeometry()
if isinstance(mesh, fmeobjects.FMEMesh):
zmin = mesh.boundingCube()[0][2]
pool = mesh.getVertices()
volume = 0.0
for part in mesh:
a, h = 0.0, 0.0
vertices = [pool[i] for i in part.getVertexIndices()]
for i, (x0, y0, z0) in enumerate(vertices[:-1]):
x1, y1, _ = vertices[i+1]
a += (x1 - x0) * (y0 + y1)
h += (z0 - zmin)
volume += abs(a) * 0.5 * (h / (len(vertices) - 1))
feature.setAttribute('_volume', volume)
+19Hi @tobiasp, I looked at the TINVolumeCalculator. This transformer seems to always require a feature having a boundary area geometry from the Boundary_INPUT. If you already have a TIN surface and need to calculate its volume,
The transformer should work with the setting above, but in my quick test, the resulting volume was wrong. The reason is that there are these fatal bugs in the implementation :-0
I therefore would not recommend you to use the TINVolumeCalculator transformer. @trentatsafe, please check carefully the implementation of the TINVolumeCalculator.
This is a workaround using the BRepSolidBoundaryCreator.Beta (from FME Hub) to create a solid below (minimum z and above) the input TIN surface.
Alternatively, this small Python script can be used to directly compute the volume below the input Mesh surface (minimum z and above).
# PythonCaller Script Example
# Assume the input feature has a Mesh geometry (typically TIN).
import fmeobjects
def calculateMeshVolume(feature):
mesh = feature.getGeometry()
if isinstance(mesh, fmeobjects.FMEMesh):
zmin = mesh.boundingCube()[0][2]
pool = mesh.getVertices()
volume = 0.0
for part in mesh:
a, h = 0.0, 0.0
vertices = [pool[i] for i in part.getVertexIndices()]
for i, (x0, y0, z0) in enumerate(vertices[:-1]):
x1, y1, _ = vertices[i+1]
a += (x1 - x0) * (y0 + y1)
h += (z0 - zmin)
volume += abs(a) * 0.5 * (h / (len(vertices) - 1))
feature.setAttribute('_volume', volume)
If you can choose only base points on the boundary of the TIN, (the volume of the TIN created from all the points - the volume of the TIN created from only the base points) would be your desired volume, I think.
Enter your E-mail address. We'll send you an e-mail with instructions to reset your password.
