在python中保存3D NetworkX图形,以便以后使用paraview查看

时间:2019-01-10 14:46:04

标签: python graph networkx paraview

我编写了此脚本,该脚本使用python中的NetworkX绘制了随机3D图形。该脚本的输出是一个3D图形,我可以在其中围绕图形结构旋转相机。

import networkx as nx
from mpl_toolkits.mplot3d import Axes3D
import matplotlib.pyplot as plt
import numpy as np
import random
import pickle


def gen_random_3d_graph(n_nodes, radius):
    pos = {i: (random.uniform(0, 1), random.uniform(0, 1), random.uniform(0, 1)) for i in range(n_nodes)}
    graph = nx.random_geometric_graph(n_nodes, radius, pos=pos)
    return graph


def plot_3d_network(graph, angle):
    pos = nx.get_node_attributes(graph, 'pos')

    with plt.style.context("bmh"):
        fig = plt.figure(figsize=(10, 7))
        ax = Axes3D(fig)
        for key, value in pos.items():
            xi = value[0]
            yi = value[1]
            zi = value[2]

            ax.scatter(xi, yi, zi, edgecolor='b', alpha=0.9)
            for i, j in enumerate(graph.edges()):
                x = np.array((pos[j[0]][0], pos[j[1]][0]))
                y = np.array((pos[j[0]][1], pos[j[1]][1]))
                z = np.array((pos[j[0]][2], pos[j[1]][2]))
                ax.plot(x, y, z, c='black', alpha=0.9)
    ax.view_init(30, angle)
    pickle.dump(fig, open('FigureObject.fig.pickle', 'wb'))
    plt.show()


if __name__ == '__main__':

    graph01 = gen_random_3d_graph(15, 0.6)
    plot_3d_network(graph01, 0)

我想保存此图以便以后使用paraview进行查看。我尝试了pickle,但没有成功。 反正可以在paraview中查看3D图形吗?

1 个答案:

答案 0 :(得分:1)

NetworkX和Paraview之间没有通用的数据类型。

通过将vtk导入python环境并使用vtkLines创建图形,您可以非常轻松地编写VTK文件(请查看此示例:https://lorensen.github.io/VTKExamples/site/Python/GeometricObjects/ColoredLines/

或者您可以从NetworkX导出到JSON并在Paraview中编写Python可编程源以读取自定义数据结构(请参见这些示例https://www.paraview.org/Wiki/Python_Programmable_Filter