gap> F:=ReadImageAsFilteredPureCubicalComplex("my_coins.png",40);;
gap> M:=FiltrationTerm(F,24);;  #Chosen after viewing degree 0 barcode for F
gap> M:=PureComplexThickened(M);;
gap> M:=PureComplexThickened(M);;
gap> C:=PureComplexComplement(M);;
gap> T:=ThickeningFiltration(C,50);;
gap> P:=PersistentBettiNumbers(T,1);;
gap> BarCodeCompactDisplay(P);
