Replies: 2 comments 4 replies
-
|
Few things are going awry here. The primary issue is that you have a pipeline. The glyph's you are tying to update are computed from the I recommend avoiding the HTH |
Beta Was this translation helpful? Give feedback.
-
How is updating the meshes directly done? That is, what must the calls to |
Beta Was this translation helpful? Give feedback.
-
Please see #4317 (comment) |
Beta Was this translation helpful? Give feedback.
-
Thanks, but this doesn't help at all. I have asked for an "exact" solution, meaning based on the code that @apcode is given. The example you are referring to --even if it solves the problem --I couldn't know-- has nothing to do with that. |
Beta Was this translation helpful? Give feedback.
-
|
Here's two examples, one with and without glyphs: Without Glyphs"""Glyph-free example.
Particles are rendered uniformly.
"""
import pyvista as pv
import numpy as np
velocity_factor = 1e-3
# Initialize particle data
num_particles = 100
# we're using a list here to allow updating in-place without using a global
velocities = [(np.random.rand(num_particles, 3) - 0.5) * velocity_factor]
# Create a PyVista mesh for the particles
particles = pv.PolyData(np.random.rand(num_particles, 3))
# Create a PyVista plotter
pl = pv.Plotter()
actor = pl.add_mesh(
particles,
show_scalar_bar=False,
name="points",
render_points_as_spheres=True,
point_size=10,
)
pl.open_movie("glyph-free.mp4", framerate=60)
def update_points():
# Update positions based on velocities
particles.points += velocities[0]
# Update velocities based on particle interactions (you can replace this
# with your custom logic)
velocities[0] += (np.random.rand(num_particles, 3) - 0.5) * velocity_factor
pl.show(interactive_update=True)
for n in range(500):
update_points()
# pl.render() # unnecessary if we're writing a frame
pl.write_frame()
pl.close()glyph-free.mp4With Glyphs"""Glyph-free example.
Particles are rendered uniformly.
"""
import pyvista as pv
import numpy as np
velocity_factor = 1e-3
size_factor = 0.05
# Initialize particle data
num_particles = 100
# we're using a list here to allow updating in-place without using a global
velocities = [(np.random.rand(num_particles, 3) - 0.5) * velocity_factor]
# Create a PyVista mesh for the particles
particles = pv.PolyData(np.random.rand(num_particles, 3))
init_size = np.random.rand(num_particles) * size_factor + size_factor
final_size = np.random.rand(num_particles) * size_factor + size_factor * 2
particles["size"] = init_size
sphere = pv.Sphere(radius=0.3)
glyphs = particles.glyph(
geom=sphere, scale="size", factor=1, rng=(size_factor, size_factor * 2)
)
# Create a PyVista plotter
pl = pv.Plotter()
# Add particles to the plotter
actor = pl.add_mesh(
glyphs, show_scalar_bar=False, name="points", cmap="bwr", smooth_shading=True
)
pl.open_movie("particles-glyph.mp4", framerate=60)
n_frames = 500
def update_points(frame_n):
# Update positions based on velocities
particles.points += velocities[0]
# allow the particles to get bigger over time
ratio = frame_n / n_frames
new_size = init_size * (1 - ratio) + final_size * ratio
print(new_size[0])
particles["size"] = new_size
# Update velocities based on particle interactions (you can replace this
# with your custom logic)
velocities[0] += (np.random.rand(num_particles, 3) - 0.5) * velocity_factor
# we have to overwrite the glyph dataset since we're changing the underlying mesh
actor.mapper.dataset = particles.glyph(
geom=sphere, scale="size", factor=1, rng=(size_factor, 0.1)
)
pl.show(interactive_update=True)
for n in range(n_frames):
update_points(n)
# pl.render() # not needed when writing a frame
pl.write_frame()
pl.close()particles-glyph.mp4 |
Beta Was this translation helpful? Give feedback.
-
|
Agree with @banesullivan's assessment. Here's the code to do this: sample.mp4from tqdm import tqdm
import pyvista as pv
import numpy as np
import random
import time
pv.set_plot_theme('document')
# Initialize particle data
num_particles = 100
positions = np.random.rand(num_particles, 3)
# Create a PyVista mesh for the particles
poly = pv.PolyData(positions)
poly['velocities'] = (np.random.rand(num_particles, 3) - 0.5) * 0.001
poly['size'] = np.random.rand(num_particles) * 0.05 + 0.05
sphere = pv.Sphere(radius=0.3)
# Create a PyVista plotter
pl = pv.Plotter(off_screen=True)
# Set up a callback function for updating particle positions
def gen_glyphs():
"""Update point positions and return glyphs of those particles."""
poly.points += poly['velocities']
# Update velocities based on particle interactions (insert logic here...)
poly['velocities'] += (np.random.rand(num_particles, 3) - 0.5) * 0.01
poly['size'] += np.random.rand(num_particles) * 0.05 - 0.025
return poly.glyph(geom=sphere, scale='size', factor=1, rng=(0.05, 0.1), orient=False)
# add the mesh, open the movie, and write the first frame
pl.add_mesh(gen_glyphs(), smooth_shading=True, name='points', show_scalar_bar=False)
pl.enable_anti_aliasing('ssaa')
pl.open_movie('/tmp/sample.mp4', quality=10)
pl.write_frame()
# write the remaining frames
for n in tqdm(range(260)):
pl.remove_actor('points') # Remove the previous actor
pl.add_mesh(gen_glyphs(), smooth_shading=True, name='points', show_scalar_bar=False)
pl.write_frame()
pl.close() |
Beta Was this translation helpful? Give feedback.
Uh oh!
There was an error while loading. Please reload this page.
Uh oh!
There was an error while loading. Please reload this page.
-
Describe the bug, what's wrong, and what you expect:
I am trying to animate a point cloud of points that move over time. Trying to use the plotter with points rendered as spheres in a window. I do not want to do this off-screen as a gif or mov.
I am trying to use the plotter update_coordinates with interactive_update=True. As I understand it I should be able to change the points and have the changes rendered on screen with this approach.
plotter.show(window_size=(1024, 1024), interactive_update=True)then
plotter.update_coordinates(points, mesh=glyphs, render=True)and
plotter.update()The script crashes on the call to update_coordinates.
I would like help to get this working, or if there is a preferred way to animate a pyvista plotter render on screen can you advise how to achieve this.
** To Reproduce **
** Result **
** System Information **
Date: Sun Apr 23 13:59:15 2023 BST
Python 3.9.13 | packaged by conda-forge | (main, May 27 2022, 17:00:33)
[Clang 13.0.1 ]
Beta Was this translation helpful? Give feedback.
All reactions