PVGeo: Voxelize model with known block size

Hello everyone,

lately I have been trying to create a voxelized (“block”) model. The data I want to visualize is based on x,y,z data that was grided before, but does not cover a whole cuboid space. When the orginial grid that the data is retrieved from is regularly spaced, I am able to get a nice voxelized model using PVGeo.filters.VoxelizePoints() . However, in my real application case the orginal grid is not cubic, but instead I have different length in different directions. I do know the resolution of the original grid (thus also the voxel size I am aiming for). When using PVGeo.filters.VoxelizePoints() I get a result that doesnt cover the whole space - there are gaps between voxels and maybe also some overlapping.

Maybe I am missing a simple solution (maybe even jsut in pyvista) for this as I know the original resolution. Some help or ideas would be highly appreciated.

Thanks a lot for any help and support and of course for your time!