File: 3D_binary_structure.py

package info (click to toggle)
scipy 1.17.0-1exp2
  • links: PTS, VCS
  • area: main
  • in suites: experimental
  • size: 235,340 kB
  • sloc: cpp: 506,914; python: 357,038; ansic: 215,028; javascript: 89,566; fortran: 19,308; cs: 3,081; f90: 1,150; sh: 860; makefile: 519; pascal: 284; lisp: 134; xml: 56; perl: 51
file content (19 lines) | stat: -rw-r--r-- 454 bytes parent folder | download | duplicates (4) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
 
import matplotlib.pyplot as plt
import scipy.ndimage


def plot_voxels(varray, ax, title):
    ax.view_init(20, 200)
    ax.voxels(varray, edgecolor="k")
    ax.set_title(title, fontsize=30)


fig = plt.figure(figsize=(16, 9))

for i in [1, 2, 3]:
    ax = fig.add_subplot(1, 3, i, projection="3d")
    arrray = scipy.ndimage.generate_binary_structure(3, i)
    plot_voxels(arrray, ax, title=f"rank=3 \n connectivity={i}")

plt.tight_layout()
plt.show()