gap> F:=ReadImageAsFilteredPureCubicalComplex(file,500);;
gap> Y:=FiltrationTerm(F,64);            
Pure cubical complex of dimension 2.
gap> BettiNumber(Y,0);
20
gap> BettiNumber(Y,1);
14
gap> Display(Y);
