sclass AI_PixelFindabilityMap1 {

  int trials = 16;

  double gain = 1;

  BufferedImage inputImage;

  BWImage scaledInputImage;

  int inputWidth = 128;

  int coreW = 8, coreH = coreW;

  int outerW = 16, outerH = outerW;

  int leewayX, leewayY;

  double[] heatMap;

  run {

    scaledInputImage = BWImage(scaleToWidth(inputWidth, inputImage));

    leewayX = scaledInputImage.getWidth()-coreW;

    leewayY = scaledInputImage.getHeight()-coreH;

    //heatMap = BWImage(leewayX, leewayY);

    heatMap = new double[leewayX*leewayY];

    for y to leewayY:

      for x to leewayX:

        calcPixel(x, y);

  }

  void calcPixel(int x, int y) {

    Rect full = imageRect(scaledInputImage);

    Rect core = rect(x, y, coreW, coreH);

    Pt center = centerOfRect(core);

    //print(+core);

    BWImage bwCore = bwClip(scaledInputImage, core);

    double sum = 0;

    repeat trials {

      Rect outer = randomRect_enclosing_in(outerW, outerH, core, full);

      IBWImage bwClip = clipIBWImage(scaledInputImage, outer);

      Pt solution = Pt(center.x-outer.x, center.y-outer.y);

      //print("Solution: " + solution);

      Pt guess = translatePt(coreW/2, coreH/2, bwSpiralSearch(bwCore, bwClip));

      //print("Guessed:  " + guess);

      if (guess == null) guess = Pt(outerW/2, outerH/2);

      sum += ptSquaredDistance(solution, guess);

    }

    heatMap[y*leewayX+x] = sum;

  }

  BWImage heatMapImage() {

    //double max = doubleMax(heatMap);

    double max = (outerW-coreW)*(outerH-coreH)/gain*trials;

    BWImage img = new(leewayX, leewayY);

    int i = 0;

    for y to leewayY:

      for x to leewayX:

        img.setPixel(x, y, /*1f-*/(float) doubleRatio(heatMap[i++], max));

    ret img;

  }

  BufferedImage mixedWithHeatMap() {

    ret combineBWImagesToColor(scaledInputImage, Color.red, bwBlackBorder(coreW/2.0, coreH/2.0, heatMapImage()), Color.green);

  }

}