Many slices around a point at different bearings about pyvista-support

[banesullivan]

The surface from the .slice filter can be oriented in any direction and centered at any location!

EDIT: I've made a number of edits to this because I messed it up the first time 🤦‍♂

Here's an example of how to get many slices at different bearings all centered around a user-chosen location:

import pyvista as pv
from pyvista import examples
import numpy as np

# Load RectilinearGrid of geo-data
mesh = examples.load_channels()

# Create a point to orient slices around
ranges = np.array(mesh.bounds).reshape(-1, 2).ptp(axis=1)
point = np.array(mesh.center) - ranges*0.25

# now generate a few normal vectors to rotate a slice around the z-axis
# use equation for circle since its about the Z-axis
increment = np.pi/6.
# use a container to hold all the slices
slices = pv.MultiBlock() # treat like a dictionary/list
for theta in np.arange(0, np.pi, increment):
    normal = np.array([np.cos(theta), np.sin(theta), 0.0]).dot(np.pi/2.)
    name = 'Bearing: {:.2f}'.format(np.rad2deg(theta))
    slices[name] = mesh.slice(origin=point, normal=normal)
slices

Information	Blocks
MultiBlock Values N Blocks 6 X Bounds 0.000, 250.000 Y Bounds 0.000, 250.000 Z Bounds 0.000, 100.000	Index Name Type 0 Bearing: 0.00 PolyData 1 Bearing: 30.00 PolyData 2 Bearing: 60.00 PolyData 3 Bearing: 90.00 PolyData 4 Bearing: 120.00 PolyData 5 Bearing: 150.00 PolyData

# and now display it!
p = pv.Plotter()
p.add_mesh(slices, cmap='rainbow', n_colors=4,)
p.add_mesh(mesh.outline())
p.show()