domain = gd.getDomain('2025')

# gets the domain parameters

# as a dictionary

nx = domain['Domain']['NX']

# assigns the number of

# voxels in x-direction

# to the variable nx

vl = domain['Domain']['VoxelLength']

# assigns the voxel length to

# variable vl

material = domain['Domain']'Material']['Name']

# assigns the material name

# of material ID00

# to variable material

print(f'The loaded structure has {nx} voxels in X-direction, a voxel length of {vl} and contains the material {material}.')

# show message in the console

nx, ny, nz = gd.getVolDimensions()

# assigns the number of voxels in

# x-direction to the variable nx,

# and the number of voxels in

# y- and z-direction to ny and nz,

# respectively

vl = gd.getVoxelLength()

# assigns the voxel length to the variable vl

voxelCounts = gd.getVoxelCounts2D(2,1)

# get the number of voxels

# of material ID 1 in

# Z-direction for each slice

# and assign them

# as a list to the variable

# voxelCounts

ZSlice_5 = voxelCounts[4]

# compute the total number

# of voxels in the structure

# and assign it to the

# variable TotalVoxels

gd.msgBox(f"In slice 5 {ZSlice_5} voxels have MaterialID 1.")

# show message dialog

nx , ny, nz = gd.getVolDimensions()

# get the number of

# voxels in all three

# directions and assign

# them to variables

TotalVoxels = nx * ny * nz

# compute the total

# number of voxels

# in the structure

# and assign it to

# the variable

# TotalVoxels

Voxels = gd.getVoxelCounts3D()

# gets list of voxel

# counts for material IDs

ID_1 = Voxels[1]/TotalVoxels * 100

# computes volume

# percentage of material

# ID 1 and assign it to

# variable ID_1

gd.msgBox(f"MaterialID 1 is assigned to {ID_1}% of the structure.")

# show message dialog

# of result

svp = gd.getSolidVolumeFraction() * 100

# compute solid volume percentage

gd.msgBox(f"The SVP is {svp}%.")

# show message dialog of result

Voxels = gd.getSelectedVoxels()

# assign list of

# selected voxels to

# variable Voxels

gd.msgBox(f"The first selected voxel is located at position {Voxels[0]}")

# show message dialog

with open("Structure.csv", "w") as fd:

# open output file for writing

# (create new file with the

# given name, if file does not exist)

# and assign it to fd.

# The file stays open for the

# following indented section.

 Structure = gd.getStructure()

# assign 3D numpy array of

# currently loaded structure to

# variable Structure. data type

# is 8-bit unsigned (uint8)

 nx, ny, nz = gd.getVolDimensions()

# assign structure volume dimensions

# to variables nx, ny and nz

 fd.write(f"{nx},{ny},{nz}\n")

# write dimensions of volume

# in first row

for x in range(nx):

# loop over all x-coordinates

  for y in range(ny):

# loop over all y-coordinates

     row = Structure[x,y,:]

# assign the z-row with

# x-coordinate x and y-coordinate y

# to the variable row

     strList = [f"{voxel_value}" for voxel_value in row]

# transform all entries of the row

# in strings and write

# them in the string list strList

     fd.write(",".join(strList) + "\n")  

# writes all entries of strList

# in the csv file, separated

# by commas, adds a new line

# at the end of the list

import numpy as np

# import Python module numpy

# to create numpy arrays

with open("Structure.csv", "r") as fd:

# open input file for reading

# and assign it to fd.

# The file is closed after the

# last indented line following

 first_row = fd.readline().strip()

# read first row and

# remove newline \n

 first_row_list = first_row.split(",")

# assign list of first_row

# entries splitted by commas

# to variable first_row_list

 nx, ny, nz = int(first_row_list[0]), int(first_row_list[1]), int(first_row_list[2])

# assign the volume dimensions

# contained in the first row

# to variables nx, ny and nz

 voxel_value_list = []

# an empty list is assigned to

# variable voxel_value_list to

# store integer values of all voxels

for line in fd:

# loop over all lines in the

# *.csv file, starting with the

# second row, as the first

# was already read

   line_stripped = line.strip()

# remove whitespace before

# and after line. in this case,

# remove newline at end of line.

   LineList = line_stripped.split(",")

# assign a list of all entries

# from line separated by

# commas to variable LineList

   LineList = [int(x) for x in LineList]

# convert each voxel_value

# string to an integer number

   voxel_value_list += LineList

# append voxel values of this

# row to list

 voxel_values = np.array(voxel_value_list, dtype=np.uint8)

# convert voxel values to numpy

# array. data type needs to

# be 8-bit unsigned (np.uint8)

# for GeoDict structures

 Structure = voxel_values.reshape(nx, ny, nz)

# reshape the 1-dimensional

# array voxel_values to a

# 3D array of given dimensions

# nx x ny x nz

 gd.setStructure(Structure, 1e-6)

# visualize the structure defined

# in the csv file in GeoDict,

# by passing the 3D numpy array

# and assigning voxel length 1µm

Struc_Des_old = gd.getStructureDescription()

# get current structure description

gd.setStructureDescription("New Description")

# changes description to “New Description”

Struc_Des_new = gd.getStructureDescription()

# gets new structure description

gd.msgBox(f'The structure description was changed from "{Struc_Des_old}" to "{Struc_Des_new}".')

# outputs the description change

import stringmap

# imports the Python module stringmap

gdr = stringmap.parseGDR('FiberGeo.gdr')

# assign the result file as a string

# to the variable gdr

GDR_Hash_64 = gdr['Geometry:Hash64']

# assign the ID of the structure

# corresponding to the result file

# to the variable GDR_Hash

Structure_Hash_64 = gd.getStructureHash64()

# assign the ID of the loaded structure

# to the variable Structure_Hash

if int(GDR_Hash_64) == Structure_Hash_64:

# condition: if the IDs are equal, the

# following indented section is executed

 gd.msgBox("The loaded structure belongs to the result file FiberGeo.gdr")  

# show message dialog

else:

# if the IDs are not equal, the following

# indented section is executed

 gd.msgBox("The loaded structure does not belong to the result file FiberGeo.gdr.")  

# show message dialog

filename = gd.getStructureFileName()

# assigns the currently loaded structure’s

# file name to the variable filename

gd.msgBox(f"Currently {filename} is loaded.")

# show message dialog