4D Density Plot in Python

Question

I am looking to plot some density maps from some grid-like data:

X,Y,Z = np.mgrids[-5:5:50j, -5:5:50j, -5:5:50j]
rho = np.random.rand(50,50,50) #for the sake of argument

I am interested in producing an interpolated density plot as shown below, from Mathematica here, using Python.

Is there any solution in Matplotlib or another plotting suite for this sort of plot?

To be clear, I do not want a scatterplot of coloured points, which is not suitable the plot I am trying to make. I would like a 3D interpolated density plot, as shown below.

Answer 1

Plotly

Plotly Approach from https://plotly.com/python/3d-volume-plots/ uses np.mgrid

import plotly.graph_objects as go
import numpy as np
X, Y, Z = np.mgrid[-8:8:40j, -8:8:40j, -8:8:40j]
values = np.sin(X*Y*Z) / (X*Y*Z)

fig = go.Figure(data=go.Volume(
    x=X.flatten(),
    y=Y.flatten(),
    z=Z.flatten(),
    value=values.flatten(),
    isomin=0.1,
    isomax=0.8,
    opacity=0.1, # needs to be small to see through all surfaces
    surface_count=17, # needs to be a large number for good volume rendering
    ))
fig.show()

Pyvista

Volume Rendering example: https://docs.pyvista.org/examples/02-plot/volume.html#sphx-glr-examples-02-plot-volume-py

3D-interpolation code you might need with pyvista: interpolate 3D volume with numpy and or scipy