# Orthogonal Slices
# =================

# View three orthogonal slices from a mesh.
# Use the `pyvista.DataObjectFilters.slice_orthogonal()` filter to create these
# slices simultaneously.
import pyvista as pv
from pyvista import examples

mesh = examples.download_embryo()
mesh.bounds

# Create three slices. Easily control their locations with the `x`, `y`,
# and `z` arguments.
slices = mesh.slice_orthogonal(x=100, z=75)

cpos = pv.CameraPosition(
    position=(540.9115516905358, -617.1912234499737, 180.5084853429126),
    focal_point=(128.31920055083387, 126.4977720785509, 111.77682599082095),
    viewup=(-0.1065160140819035, 0.032750075477590124, 0.9937714884722322),
)
dargs = dict(cmap='gist_ncar_r')

pl = pv.Plotter()
pl.add_mesh(slices, **dargs)
pl.show_grid()
pl.show(cpos=cpos)

pl = pv.Plotter(shape=(2, 2))
# XYZ - show 3D scene first
pl.subplot(1, 1)
pl.add_mesh(slices, **dargs)
pl.show_grid()
pl.camera_position = cpos
# XY
pl.subplot(0, 0)
pl.add_mesh(slices, **dargs)
pl.show_grid()
pl.camera_position = 'xy'
pl.enable_parallel_projection()
# ZY
pl.subplot(0, 1)
pl.add_mesh(slices, **dargs)
pl.show_grid()
pl.camera_position = 'zy'
pl.enable_parallel_projection()
# XZ
pl.subplot(1, 0)
pl.add_mesh(slices, **dargs)
pl.show_grid()
pl.camera_position = 'xz'
pl.enable_parallel_projection()

pl.show()

# ----------------------------------------------------------------------
# Generated by sphinx-examples-as-code https://github.com/pyvista/sphinx-examples-as-code