1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
|
"""
========================================
Surface seeding for tractography
========================================
Surface seeding is a way to generate initial position for tractography
from cortical surfaces position :footcite:p:`StOnge2018`.
"""
from fury.io import load_polydata
from fury.utils import (
get_actor_from_polydata,
get_polydata_triangles,
get_polydata_vertices,
normals_from_v_f,
)
import numpy as np
from dipy.data import get_fnames
from dipy.tracking.mesh import (
random_coordinates_from_surface,
seeds_from_surface_coordinates,
)
from dipy.viz import actor, window
###############################################################################
# Fetch and load a surface
brain_lh = get_fnames(name="fury_surface")
polydata = load_polydata(brain_lh)
###############################################################################
# Extract the triangles and vertices
triangles = get_polydata_triangles(polydata)
vts = get_polydata_vertices(polydata)
###############################################################################
# Display the surface
# ===============================================
#
# First, create an actor from the polydata, to display in the scene
scene = window.Scene()
surface_actor = get_actor_from_polydata(polydata)
scene.add(surface_actor)
scene.set_camera(position=(-500, 0, 0), view_up=(0.0, 0.0, 1))
# Uncomment the line below to show to display the window
# window.show(scene, size=(600, 600), reset_camera=False)
window.record(scene=scene, out_path="surface_seed1.png", size=(600, 600))
###############################################################################
# .. rst-class:: centered small fst-italic fw-semibold
#
# Initial cortical surface
#
#
# Generate a list of seeding positions
# ==============================================================
#
# Choose the number of seed
nb_seeds = 100000
nb_triangles = len(triangles)
###############################################################################
# Get a list of triangles indices and trilinear coordinates for each seed
tri_idx, trilin_co = random_coordinates_from_surface(nb_triangles, nb_seeds)
###############################################################################
# Get the 3d cartesian position from triangles indices and trilinear
# coordinates
seed_pts = seeds_from_surface_coordinates(triangles, vts, tri_idx, trilin_co)
###############################################################################
# Compute normal and get the normal direction for each seeds
normals = normals_from_v_f(vts, triangles)
seed_n = seeds_from_surface_coordinates(triangles, normals, tri_idx, trilin_co)
###############################################################################
# Create dot actor for seeds (blue)
seed_actors = actor.dot(seed_pts, colors=(0, 0, 1), dot_size=4.0)
###############################################################################
# Create line actors for seeds normals (green outside, red inside)
normal_length = 0.5
normal_in = np.tile(seed_pts[:, np.newaxis, :], (1, 2, 1))
normal_out = np.tile(seed_pts[:, np.newaxis, :], (1, 2, 1))
normal_in[:, 0] -= seed_n * normal_length
normal_out[:, 1] += seed_n * normal_length
normal_in_actor = actor.line(normal_in, colors=(1, 0, 0))
normal_out_actor = actor.line(normal_out, colors=(0, 1, 0))
###############################################################################
# Visualise seeds and normals along the surface
scene = window.Scene()
scene.add(surface_actor)
scene.add(seed_actors)
scene.add(normal_in_actor)
scene.add(normal_out_actor)
scene.set_camera(position=(-500, 0, 0), view_up=(0.0, 0.0, 1))
# Uncomment the line below to show to display the window
# window.show(scene, size=(600, 600), reset_camera=False)
window.record(scene=scene, out_path="surface_seed2.png", size=(600, 600))
###############################################################################
# .. rst-class:: centered small fst-italic fw-semibold
#
# Surface seeds with normal orientation
#
#
# References
# ----------
#
# .. footbibliography::
#
|
