# Plot Over Line
# ==============

# Plot the values of a dataset over a line through that dataset
# using the `plot_over_line()` filter.
import pyvista as pv
from pyvista import examples

# Volumetric Mesh
# ---------------

# First a 3D mesh example to demonstrate
mesh = examples.download_kitchen()

# Make two points to construct the line between
a = [mesh.bounds.x_min, mesh.bounds.y_min, mesh.bounds.z_min]
b = [mesh.bounds.x_max, mesh.bounds.y_max, mesh.bounds.z_max]

# Preview how this line intersects this mesh
line = pv.Line(a, b)

pl = pv.Plotter()
pl.add_mesh(mesh, style='wireframe', color='w')
pl.add_mesh(line, color='b')
pl.show()

# Run the filter and produce a line plot
mesh.plot_over_line(a, b, resolution=100)

# Flat Surface
# ------------

# We could also plot the values of a mesh that lies on a flat surface
mesh = examples.download_st_helens()

# Make two points to construct the line between
a = [mesh.center[0], mesh.bounds.y_min, mesh.bounds.z_max]
b = [mesh.center[0], mesh.bounds.y_max, mesh.bounds.z_max]

# Preview how this line intersects this mesh
line = pv.Line(a, b)

pl = pv.Plotter()
pl.add_mesh(mesh)
pl.add_mesh(line, color='white', line_width=10)
pl.add_point_labels([a, b], ['A', 'B'], font_size=48, point_color='red', text_color='red')
pl.show()

# Run the filter and produce a line plot
mesh.plot_over_line(
    a,
    b,
    resolution=10000,
    title='Elevation Profile',
    ylabel='Height above sea level',
    figsize=(10, 5),
)

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