Histogram = gd.ImportGeoVol.getHistogram()

# get list of tuples describing

# the histogram and assign it

# to the variable Histogram

file = open('Histogram.csv', 'w')

# open output file for writing (

# create new file with the given name,

# if file does not exist)

file.write('Value,Count\n')

# write titles for columns in csv file

for i in Histogram:

# loop over all tuples i of Histogram

 file.write(f'{i[0]},{i[1]}\n')

# write values and counts

# into the csv file

file.close()

# close the csv file

nx = gd.ImportGeoVol.getNewImageDimensions(0)

# get number of voxels in

# X-direction (direction 0)

gd.msgBox(f"In X-direction there are {nx} voxels.")

# show message box

OTSU = gd.ImportGeoVol.getOtsuThreshold()

# get threshold and assign

# it to variable OTSU

gd.msgBox(f"OTSU threshold is {OTSU}")

# show message box

OTSU = gd.ImportGeoVol.getMultiOtsuThreshold()

# get threshold list and

# assign it to variable OTSU

gd.msgBox(f"The OTSU thresholds are {OTSU}")

# show message box

vl_old = gd.ImportGeoVol.getVoxelLength()

# get voxel length of

# currently loaded gray value

# image and assign it to

# variable vl_old

gd.ImportGeoVol.setVoxelLength(1e-6)

# set voxel length of

# currently loaded gray

# value image to 1µm

vl_new = gd.ImportGeoVol.getVoxelLength()

# assign new voxel length

# to variable vl_new

gd.msgBox(f"The voxel length was changed from {vl_old} to {vl_new}")

# show message box

Rot = gd.ImportGeoVol.getRotationSuggestion(False)

# get rotation suggestion

# for suggest plane and

# assign it to variable

# rotation.

# Threshold is omitted -

# default to automatic

# threshold

Rotation_args = gd.getBuiltinDefaults("ImageProcessing:Rotation")

# get Built-in Defaults

# for the Python dictionary

# of the GeoDict command

# Rotation and assign the

# dictionary to variable

# Rotation_args

Rotation_args['Phi']   = Rot[0]

# assign the rotation values

# to the corresponding keys

# in the rotation dictionary

Rotation_args['Theta'] = Rot[1]

Rotation_args['Psi']   = Rot[2]

gd.runCmd("ImageProcessing:Rotation", Rotation_args)

# rotate the gray value image,

# version is omitted –

# default to latest