static BufferedImage scaledBufferedImageFromIntegralImage_direct(int w, IntegralImage img) {
  ret scaledBufferedImageFromIntegralImage_direct(img, w);
}

static BufferedImage scaledBufferedImageFromIntegralImage_direct(IntegralImage img, int w) {
  ret scaledBufferedImageFromIntegralImage_direct(img, w, rowsForColsWithAspectRatio(w, img.w, img.h));
}

static BufferedImage scaledBufferedImageFromIntegralImage_direct(IntegralImage img, int w, int h) {
  int w1 = img.w, h1 = img.h;
  BufferedImage out = newBufferedImage(w, h);
  for y to h:
    for x to w: {
      int x1 = x*w1/w, x2 = max(x1+1, (x+1)*w1/w);
      int y1 = y*h1/h, y2 = max(y1+1, (y+1)*h1/h);
      int area = (x2-x1)*(y2-y1);
      int r = integralImage_sumRect(img, x1, y1, x2, y2, 0)/area;
      int g = integralImage_sumRect(img, x1, y1, x2, y2, 1)/area;
      int b = integralImage_sumRect(img, x1, y1, x2, y2, 2)/area;
      out.setRGB(x, y, (r << 16) | (g << 8) | b);
    }
  ret out;
}

Began life as a copy of #1025282