gap> file:=HapFile("image1.3.2.png");;
gap> F:=ReadImageAsFilteredPureCubicalComplex(file,40);
Filtered pure cubical complex of dimension 2.
gap> P:=PersistentBettiNumbers(F,0);;
gap> BarCodeCompactDisplay(P);
