from guf import GUF

# import GUF library

guf_file = GUF("StokesResult/FlowField_z.vap")

# access GUF file

# FlowField_z.vap

guf_header = guf_file.getHeader()

# assign the header

# stringmap to the

# variable guf_header

print(guf_header)

# print the complete

# header to the

# GeoDict console

imagenumber = guf_header["NumberOfImages"]

# assign the number of

# images to the

# variable imagenumber

gd.msgBox(f"The file contains {imagenumber} images.")

# show message dialog

from guf import GUF

# import GUF library

guf_file = GUF("StokesResult/FlowField_z.vap")

# access GUF file

# FlowField_z.vap

guf_image = guf_file.getImage("Velocity")

# assign the numpy array

# corresponding to the

# image Velocity to the

# variable guf_image

guf_field = guf_image["VelocityX"]

# assign the numpy array

# corresponding to the

# flow field VelocityX to

# the variable guf_field

gd.msgBox(f"The velocity at position (50,50,50) in the Velocity X field is {guf_field[50][50][50]}.")

# show a message dialog

# for the velocity at

# position (50,50,50)

from guf import GUF

# import GUF library

guf_file = GUF("FilterLifeTime/Batch00001/TrackerFinalParticles.gpp")

# access FilterDict result

# file TrackerFinalParticles.gpp

guf_array = guf_file.getArray("ParticlePositions")

# assign the numpy array

# containing the particle

# positions to the variable

# guf_array

id_5  = guf_array["ID"][5]

# assign fifth element in

# the column ID to the

# variable id_5

# (count starts with 0)

pos_5 = guf_array["Position X"][5]

# assign fifth element in

# the column Position X

# to the variable pos_5

# (count starts with 0)

time_5 = guf_array["Time"][5]

# assign fifth element

# in the column Time to

# the variable time_5

# (count starts with 0)

gd.msgBox(f"The particle with ID {id_5} has the X-position {pos_5} at time {time_5}.")

# show message dialog

guf_row = guf_array[5]

# assign the numpy array

# containing the sixth

# entry of all columns

#  to the variable guf_row

gd.msgBox(f"The particle with ID {guf_row[0]} has the X-position {guf_row[2]} at time {guf_row[8]}.")

# show the same message

# dialog as before

from guf import GUF

# import GUF library

guf_file = GUF("FiberGeo/Structure.gdt")

# access GUF file Structure.gdt

# in the folder FiberGeo

guf_map = guf_file.getMap("GADStats")

# assign the stringmap of the

# GAD statistics for all 2D

# slices to the variable guf_map

print(guf_map)

# print the stringmap to the

# GeoDict console

objectscount_Z = guf_map["PerSliceObjectCountsZ"]

# assign the string containing

# statistics for the Z-slices to

# the variable objectscount_Z

objectscount_Z = objectscount_Z.split(',')

# split the string by commas,

# and obtain a list

count_Z_slice_13 = objectscount_Z[13]

# assign the 14th entry in the

# list (index 13 as counting

# starts with 0) to the variable

# count_Z_slice

gd.msgBox(f" In Z-slice 14 there are {count_Z_slice_13} objects.")

# show a message dialog.