interface IG22MasksHolder<A> {
  replace Mask with IG22Mask<A>.
  
  WidthAndHeight maskSize();
  L<Image2B> maskImages();
  //LPair<Image2B, A> maskImagesWithLabels();
  L<? extends Mask> masks();
  G22GhostImage ghost();
  double certainty();
  default int nMasks() { ret l(masks()); }
  
  default S renderCountAndCertainty() {
    ret main nMasks(nMasks()) + " with c=" + formatDouble3(certainty());
  }
  
  default bool canSplit() { ret nMasks() > 1; }
  
  default L<IG22MasksHolder<A>> subHolders() { null; }
  default void transformSubHolders(IF1<IG22MasksHolder<A>> f) {}
  
  default PStackComputable<Void> findSimilarMasks(
    Image2B maskImage, ProbabilisticList<Mask> outList) {
    ret new FindSimilarMasks_BruteForce<A>(this, maskImage, outList);
  }

  srecord noeq FindSimilarMasks_BruteForce<A>(IG22MasksHolder<A> masksHolder,
    Image2B maskImage,
    ProbabilisticList<Mask> outList) extends PStackComputableWithStep<Void> {
    Mask currentMask;
    
    toString { ret formatFunctionCall(shortClassName(this), masksHolder, "step=" + step); }
    
    void step(IPStack stack) {
      if (step == 0) {
        ++step;
        stack.optionsOrReturn(this, map(masksHolder.masks(), mask ->
          self -> self.currentMask = mask));
      } else {
        long diff = binaryImagePixelDifference_sameSize(
          maskImage,
          currentMask.image());
          
        // TODO: This formula will drop images which are
        // complete inverts of each other since then p is 0
        // and ProbabilisticList discards this
        double p = 1-doubleRatio(diff, area(currentMask.image()));
        printVars("Found mask", +diff, +p, +currentMask);
        outList.at(p, currentMask);
        stack.ret();
      }
    }
  }
}