def setup$()$:

global source, destination, ti$,$ e$,$ tf

size$(8000,6000)$

source $=$ loadImage$("$image$.$jpeg"$)$

# The destination image is created as a blank image the same size as the source.

destination $=$ createImage$($source$.$width, source.height, RGB$)$

e $=5$

ti $=0$

tf $=$ adpt$\_$th$($ti$)$

while$($abs$($ti $-$ tf$)>$e$)$:

temp $=$ adpt$\_$th$($tf$)$

ti $=$ tf

tf $=$ temp

print$($tf$)$

img$1=$ binarize$($tf$)$

img$2=$ crop$($img$1)$

image$($source$,0,0)$

image$($img$1,0,$ source.height$+10)$

image$($img$2,0,2*$source$.$height$+20)$

def adpt$\_$th$($t$)$:

e $=5$

ti $=120$

r $=[20,40,250,70,90,120,170,190,220,5]$

g$1=[]$

g$2=[]$

for i in range$($len$($r$))$:

if r$[$i$]<=$ ti:

g$1.$append$($r$[$i$])$

else:

g$2.$append$($r$[$i$])$

sum$1=$ sum$($g$1)$

sum$2=$ sum$($g$2)$

l$1=$ len$($g$1)$

l$2=$ len$($g$2)$

a$1=$ sum$1/$ l$1$ if l$1!=0$ else $0$

a$2=$ sum$2/$ l$2$ if l$2!=0$ else $0$

tf $=($a$1+$ a$2)/2$ if $($l$1!=0$ and l$2!=0)$ else $0$

return g$1,$ g$2,$ sum$1,$ sum$2,$ l$1,$ l$2,$ a$1,$ a$2,$ tf

# Example usage

g$1,$ g$2,$ sum$1,$ sum$2,$ l$1,$ l$2,$ a$1,$ a$2,$ tf $=$ adpt$\_$th$(0)$

print$("$g$1$:$",$ g$1)$

print$("$g$2$:$",$ g$2)$

print$("$Sum of elements in g$1$:$",$ sum$1)$

print$("$Sum of elements in g$2$:$",$ sum$2)$

print$("$Length of g$1$:$",$ l$1)$

print$("$Length of g$2$:$",$ l$2)$

print$("$Average of elements in g$1$:$",$ a$1)$

print$("$Average of elements in g$2$:$",$ a$2)$

print$("$Threshold:$",$ tf$)$

def binarize$($t$)$:

source.loadPixels$()$

destination.loadPixels$()$

for x in range$($source$.$width$)$:

for y in range$($source$.$height$)$:

loc $=$ x $+$ y$*$source$.$width

r $=$ brightness$($source$.$pixels$[$loc$])$

#print$($r$)$

# Test the brightness against the threshold

if $($r $>$ tf$)$:

destination.pixels$[$loc$]=$ color$(255)$

else:

destination.pixels$[$loc$]=$ color$(0)$

destination.updatePixels$()$

return$($destination$)$

def crop$($destination$)$:

xl $=0$

xr $=0$

yt $=0$

yb $=0$

destination.loadPixels$()$

found $=$ False

for x in range$($destination$.$width$)$:

for y in range$($destination$.$height$)$:

if $($found $==$ False$)$:

loc $=$ loc $=$ x $+$ y$*$source$.$width

r $=$ brightness$($destination$.$pixels$[$loc$])$

#print$($r$)$

if r $==0$:

xl $=$ x

found $=$ True

found $=$ False

for x in reverse$($range$($destination$.$width$))$:

for y in range$($destination$.$height$)$:

if $($found $==$ False$)$:

loc $=$ loc $=$ x $+$ y$*$source$.$width

r $=$ brightness$($destination$.$pixels$[$loc$])$

#print$($r$)$

if r $==0$:

xr $=$ x

found $=$ True

found $=$ False

for y in range$($destination$.$height$)$:

for x in range$($destination$.$width$)$:

if $($found $==$ False$)$:

loc $=$ loc $=$ x $+$ y$*$source$.$width

r $=$ brightness$($destination$.$pixels$[$loc$])$

#print$($r$)$

if r $==0$:

yt $=$ y

found $=$ True

found $=$ False

for y in reverse$($range$($destination$.$height$))$:

for x in range$($destination$.$width$)$:

if $($found $==$ False$)$:

loc $=$ loc $=$ x $+$ y$*$source$.$width

r $=$ brightness$($destination$.$pixels$[$loc$])$

#print$($r$)$

if r $==0$:

yb $=$ y

found $=$ True

print$($xl$,$ xr$,$ yt$,$ yb$)$

i $=$ destination.get$($xl$,$ yt$,$ xr$-$xl$,$ yb$-$yt$)$

return$($i$)$

Alternate code for this in python