Operations on vectors

This example shows how to perform some basic operations slicing and selecting vectors. This is designed to be very flexible and powerful. Many operations such as slicing with a boolean array are supported.

Additionally, lazy operations are supported and can be chained together. These are often faster than their non-lazy counterparts as dask very effectively prunes the computation graph.

import pyxem as pxm
import hyperspy.api as hs

hs.set_log_level("ERROR")

s = pxm.data.tilt_boundary_data()
temp = s.template_match_disk(disk_r=5, subtract_min=False)

vectors = s.get_diffraction_vectors(threshold_abs=0.4, min_distance=5)

# Plotting all the vectors

s.plot()
all_vectors = vectors.to_markers(color="red", sizes=10, alpha=0.5)
s.add_marker(all_vectors)

• 
• 

[                                        ] | 0% Completed | 150.95 us
[                                        ] | 0% Completed | 100.48 ms
[                                        ] | 0% Completed | 200.85 ms
[                                        ] | 0% Completed | 301.22 ms
[########################################] | 100% Completed | 401.62 ms

[                                        ] | 0% Completed | 133.99 us
[                                        ] | 0% Completed | 100.64 ms
[                                        ] | 0% Completed | 200.94 ms
[                                        ] | 0% Completed | 301.26 ms
[                                        ] | 0% Completed | 401.66 ms
[                                        ] | 0% Completed | 502.35 ms
[                                        ] | 0% Completed | 602.66 ms
[                                        ] | 0% Completed | 703.09 ms
[########################################] | 100% Completed | 803.45 ms

[                                        ] | 0% Completed | 140.50 us
[########################################] | 100% Completed | 100.46 ms

[                                        ] | 0% Completed | 251.13 us
[########################################] | 100% Completed | 101.37 ms

slic_vectors = (vectors.ivec[:, vectors.ivec[0] < 10]).to_markers(
    color="green", sizes=5, alpha=0.5
)

s.plot()
s.add_marker([all_vectors, slic_vectors])

• 
• 

[                                        ] | 0% Completed | 142.52 us
[########################################] | 100% Completed | 100.52 ms

[                                        ] | 0% Completed | 138.47 us
[########################################] | 100% Completed | 100.46 ms

[                                        ] | 0% Completed | 151.19 us
[########################################] | 100% Completed | 100.49 ms

[                                        ] | 0% Completed | 139.46 us
[########################################] | 100% Completed | 100.43 ms

slic_vectors = (
    vectors.ivec[:, (vectors.ivec[0] > 0) * (vectors.ivec[0] < 10)]
).to_markers(color="w", sizes=5, alpha=0.5)
s.plot()
s.add_marker([all_vectors, slic_vectors])

• 
• 

[                                        ] | 0% Completed | 144.08 us
[########################################] | 100% Completed | 100.47 ms

[                                        ] | 0% Completed | 138.28 us
[########################################] | 100% Completed | 100.42 ms

[                                        ] | 0% Completed | 130.58 us
[########################################] | 100% Completed | 100.41 ms

[                                        ] | 0% Completed | 148.56 us
[########################################] | 100% Completed | 100.45 ms

[                                        ] | 0% Completed | 145.80 us
[########################################] | 100% Completed | 100.45 ms

[                                        ] | 0% Completed | 138.30 us
[########################################] | 100% Completed | 100.47 ms

vect_magnitudes = (vectors.ivec[0] ** 2 + vectors.ivec[1] ** 2) ** 0.5
slic_vectors = vectors.ivec[:, vect_magnitudes < 20].to_markers(
    color="w", sizes=5, alpha=0.5
)
s.plot()
s.add_marker([all_vectors, slic_vectors])
s.add_marker([all_vectors, slic_vectors])

• 
• 

[                                        ] | 0% Completed | 148.72 us
[########################################] | 100% Completed | 100.55 ms

[                                        ] | 0% Completed | 145.56 us
[########################################] | 100% Completed | 100.47 ms

[                                        ] | 0% Completed | 141.18 us
[########################################] | 100% Completed | 100.47 ms

[                                        ] | 0% Completed | 150.62 us
[########################################] | 100% Completed | 100.47 ms

[                                        ] | 0% Completed | 133.19 us
[########################################] | 100% Completed | 100.42 ms

slic_vectors = (vectors.ivec[:, vectors.ivec["intensity"] < 0.5]).to_markers(
    color="w", sizes=5, alpha=0.5
)
s.plot()
s.add_marker([all_vectors, slic_vectors])
s.add_marker([all_vectors, slic_vectors])

• 
• 

[                                        ] | 0% Completed | 141.81 us
[########################################] | 100% Completed | 100.46 ms

[                                        ] | 0% Completed | 139.28 us
[########################################] | 100% Completed | 100.42 ms

[                                        ] | 0% Completed | 148.53 us
[########################################] | 100% Completed | 100.45 ms

[                                        ] | 0% Completed | 129.68 us
[########################################] | 100% Completed | 100.42 ms