#! /bin/bash -e
# Last edited on 2025-08-06 20:48:20 by stolfi

???

temp="/tmp/$$"

??? Change to use the 10^{-5} fractile.

# Get the actual ${maxval}:
echo "  finding actual maxval ..." 1>&2
mfile="${temp}-maxval.txt"
cat ${hist_file} \
  | gawk \
      ' BEGIN { maxval = 0; mcount = -1; tcount = 0; }
        /^[ ]*[0-9]/ {
          val = $1; count = $2;
          if ((count > 0) && (val > maxval)) {
            maxval = val; mcount = count
          }
          tcount += count
        }
        END { print maxval, mcount, tcount; }
      ' \
  > ${mfile}
real_maxval=( `cat ${mfile}` )
echo "  actual maxval = ${real_maxval[0]}" 1>&2
echo "  actual maxval count = ${real_maxval[1]}" 1>&2
echo "  num pixels = ${real_maxval[2]}" 1>&2