gap> F:=ReadImageAsFilteredPureCubicalComplex(file,20);;
gap> P:=PersistentBettiNumbersAlt(F,1);;
gap> BarCodeCompactDisplay(P);
