#!/bin/sh

[ -r .conf ] || exit 125

UFRAW='eval ufraw-batch --interpolation=ahd --out-type=tiff --out-depth=16 --overwrite --zip'
DCRAW='dcraw -v -r 2.008866 1.183000 1.540229 1.183000 -S 62976 -H 0 -q 3 -o 0 -n 80 -m 4 -4 -T'
#DCRAW='dcraw -v -w -H 0 -q 3 -n 300 -4 -T'
PMATIC='eval panomatic -n ${CPUS} --fullscale --surfscore 50 --kdtreesteps 100 --kdtreeseconddist 0.3 --ransacdist 100 --minmatches 16'
#PMATIC='eval panomatic -n ${CPUS} --surfscore 10 --kdtreesteps 100 --kdtreeseconddist 0.3 --ransacdist 100 --minmatches 16'
REMAPLDR='eval nona -v -t ${CPUS} -m TIFF_m -r ldr -p UINT${DEPTH} -z DEFLATE -o tmp-'
LDRSFX='tif'
#REMAPHDR='eval nona -v -t ${CPUS} -m EXR_m -r hdr -p FLOAT -z DEFLATE -o tmp-'
#HDRSFX='exr'
REMAPHDR='eval nona -v -t ${CPUS} -m TIFF_m -r hdr -p FLOAT -z DEFLATE -o tmp-'
HDRSFX='tif'
FULLA='eval fulla -v -t ${CPUS} -s'
CCM='color_correct -E -c -p.03 -D64 -a-4096:4096+1024 -k24:40/32 -g28:36/32'
BLEND='eval enblend -v --compression=DEFLATE -m ${MEM} -b 16384 --fine-mask'
#BLEND='eval enfuse -v --compression=DEFLATE -m ${MEM} -b 16384 --wExposure=0.01  --wSaturation=0.01  --wContrast=0.99 --HardMask -w'
EFUSE='eval enfuse -v --compression=DEFLATE -m ${MEM} -b 16384'
FFUSE='eval enfuse -v --compression=DEFLATE -m ${MEM} -b 16384 --wExposure=0.01  --wSaturation=0.01  --wContrast=0.99 --HardMask'
#MERGE='hugin_hdrmerge -v -m khan -i 8'
#MERGE='hugin_hdrmerge -v -m avg -c'
MERGE='enfuse -v --compression=DEFLATE -m ${MEM} -b 16384'
#CONV='convert -limit memory 512 -limit map 256'
CONV='convert'
#INPNT='greycstoration -alpha 0.9 -sigma 0.2 -p 5 -a 1 -dt 20 -iter 5 -bits 16 -m mask.png -inpaint'
#INPNT='greycstoration -alpha 0.9 -sigma 0.2 -iter 512 -bits 16 -inpaint'
INPNT='greycstoration -iter 1024 -bits 16 -inpaint'

RGB='-fx (s<=0.03928)?s/12.92:((s+0.055)/1.055)^2.4'
sRGB='-fx (s<=0.00304)?s*12.92:s^0.416666667*1.055-0.055'

PFSMAX=75013120
PFSHOST=192.168.0.55

LANG=C; export LANG

. .conf

[ -n "${DEPTH}" ] || DEPTH="${2}"
[ -n "${DEPTH}" ] || DEPTH=16

[ -z "${KEEPHDR}" ] || KEEPHDR="tee ${DST}.pfs"
[ -n "${KEEPHDR}" ] || KEEPHDR=cat

[ -n "${CPUS}" ] || CPUS=$(sysctl -n kern.smp.cpus)
[ -n "${CPUS}" ] || CPUS=1

[ -n "${MEM}" ] || MEM="$(($(ulimit -d)/1024))"
[ -n "${MEM}" ] || MEM=512
[ ${MEM} -lt 1024 ] || MEM=1024

prep () {
    for e in ${EXP}
    do
        #eval "SEQ=\${SEQ${e}}"
        eval "SRC=\${SRC${e}}"
        eval "SFX=\${SFX${e}}"
        eval "UCF=\${UCF${e}}"
        N=0
        for n in ${SRC}
        do
            o=$(printf '%02d%s' ${N} ${SFX})
            N=$((${N}+1))
            if [ "${n}" = SSKIP ]
            then
                if [ ${o} = 03-e02 -o ${o} = 04-e02 ]
                then
                    eval "NSFX=\${SFX$((${e}+1))}"
                    ln -s $(printf '%02d%s.tiff' ${N} ${NSFX}) ${o}.tiff
                fi
            fi
            [ -z "${n}" -o "${n}" = SKIP -o "${n}" = SSKIP ] && continue
            [ -r ${o}.tiff ] && continue
            [ -r ${n}.RAW.bz2 ] && bzcat ${n}.RAW.bz2 >${n}.RAW
            [ -r ${n}.RAW.gz ] && zcat ${n}.RAW.gz >${n}.RAW
            if [ -r ${n}.RAW ]
            then
                if [ -z "${LIN}" ]
                then
                    conf=${n}.ufraw
                    [ -z "${UCONF}" ] || conf=${UCONF}.ufraw
                    [ -z "${UCF}" ] || conf=${UCF}.ufraw
                    [ -r ${conf} ] || exit 125
                    ${UFRAW} --conf=${conf} --output=${n}.tiff ${n}.RAW
                else
                    ${DCRAW} ${n}.RAW
                fi
                fl=$(dcraw -i -v ${n}.RAW|awk '/^Focal length: /{print $3}')
                tca="$( ( cd ../../tools/tca && echo "disp(tca('tca.db','fz50',${fl},1,0))"|octave -q|cat|head -1 ) )"
                #echo "Focal length: ${fl}; TCA correction: ${tca}"
                ${FULLA} ${tca} -o ${n}.tiff ${n}.tiff
                eval "$(identify -format 'w=%w; h=%h' ${n}.tiff)"
                ${CONV} ${n}.tiff \
                    -depth 16 \
                    -crop $((${w}-24))x$((${h}-12))+12+6 +repage \
                    -orient undefined \
                    -unsharp 0x0.3 \
                    -compress zip ${o}.tiff
                rm -f ${n}.tiff
                if [ -n "${LIN}" -a "${1}" = prep -a "${DEPTH}" = 8 ]
                then
                    if [ -n "${3}" ]
                    then
                        case ${e} in
                        0) m='1.5625,50';;
                        1) m='6.25,50';;
                        2) m='25,50';;
                        3) m='100,50';;
                        4) m='400,50';;
                        5) m='1600,50';;
                        6) m='6400,50';;
                        *) m='100,100';;
                        esac
                    fi
                    echo "${N}/${e}: ${m}"
                    [ -z "${m}" ] || m="-modulate ${m}"
                    ${CONV} ${o}.tiff \
                        ${m} ${sRGB} \
                        -compress zip ${o}.tiff
                fi
                [ "${DEPTH}" = 16 ] \
                || ${CONV} ${o}.tiff -gamma 1.2 -depth 8 -compress zip ${o}.tiff
                exiftool -overwrite_original -P -TagsFromFile ${n}.RAW ${o}.tiff
                [ -r ${n}.RAW.bz2 -o -r ${n}.RAW.gz ] && rm -f ${n}.RAW
                [ -r ${n}.JPG ] \
                && exiftool -overwrite_original -P -TagsFromFile ${n}.JPG ${o}.tiff
            elif [ -r ${n}.JPG ]
            then
                # sRGB => RGB
                [ -z "${LIN}" ] || rgb="${RGB}"
                ${CONV} ${n}.JPG \
                    -unsharp 0x0.3 \
                    ${rgb} \
                    -depth ${DEPTH} \
                    -compress zip ${o}.tiff
                exiftool -overwrite_original -P -TagsFromFile ${n}.JPG ${o}.tiff
            else
                exit 126
            fi
            [ -r mask-${o}.png -a \( "${1}" != prep -o "${DEPTH}" != 8 \) ] \
            && ${CONV} ${o}.tiff -matte \
                mask-${o}.png -compose multiply -composite \
                mask-${o}.png -compose copyopacity -composite \
                -compress zip ${o}.tiff
        done
    done
}

auto () {
    autoexp () {
        x=${1}; shift
        e=${1}; shift
        [ -r ${DST}${e}${x}-ao.pto ] && continue
        [ ${#} -lt 2 ] && continue
        S=''
        for a in ${*}
        do
            s="$(printf '%s%s.tiff' ${a} ${e})"
            [ -r "${s}" ] && S="${S} ${s}"
        done
        printf 'e=%s; S=%s\n' ${e} "${S}"
        ${PMATIC} \
            --linearmatch --linearmatchlen 20 \
            --sieve1width 64 --sieve1height 48 --sieve1size 8 \
            --sieve2width 16 --sieve2height 16 --sieve2size 4 \
            -o ${DST}${e}${x}-pm.pto ${S}
        autooptimiser -a -l -s -o ${DST}${e}${x}-ao.pto ${DST}${e}${x}-pm.pto
    }
    autoexpex () {
        e="${1}"; shift
        [ -r ${DST}${e}-ao.pto ] && continue
        if [ -z "${ROWS}" ]
        then
            autoexp '' "${e}" "${@}"
            return ${?}
        fi
        eval "${ROWS}"
        for re in ${REXP}
        do
            eval "RSFX=\${RSFX${re}}"
            [ "${e}" = "${RSFX}" ] && break
        done
        if [ ${?} -eq 0 ]
        then
            eval "RROWS=\${ROWS${re}}"
            for r in ${RROWS}
            do
                eval "RSFX=\${RSFX${re}${r}}"
                eval "RSRC=\${RSRC${re}${r}}"
                autoexp "${RSFX}" "${e}" ${RSRC}
            done
            ptomerge ${DST}${e}[^-]*-ao.pto ${DST}${e}-ao.pto
            ptosort -i n ${DST}${e}-ao.pto ${DST}${e}-ao.pto
            autooptimiser -a -l -s -o ${DST}${e}-ao.pto ${DST}${e}-ao.pto
            return ${?}
        fi
        autoexp '' "${e}" "${@}"
    }
    autosect () {
        a=${1}; shift
        a="$(printf '%02d' ${a})"
        [ -r ${DST}-${a}-ao.pto ] && continue
        echo "${DST}-${a}-ao.pto: ${#}"
        [ ${#} -lt 2 ] && continue
        S=''
        for e in ${*}
        do
            s="$(printf '%s%s.tiff' ${a} ${e})"
            [ -r "${s}" ] && S="${S} ${s}"
        done
        printf 'a=%s; S=%s\n' ${a} "${S}"
        ${PMATIC} \
            --sieve1width 64 --sieve1height 48 --sieve1size 8 \
            --sieve2width 32 --sieve2height 24 --sieve2size 4 \
            -o ${DST}-${a}-pm.pto ${S}
        autooptimiser -a -o ${DST}-${a}-ao.pto ${DST}-${a}-pm.pto
    }
    for e in ${EXP}
    do
        eval "SRC=\${SRC${e}}"
        eval "SFX=\${SFX${e}}"
        A=''
        N=0
        for n in ${SRC}
        do
            [ -z "${n}" -o "${n}" = SKIP ] \
            || A="${A} $(printf '%02d' ${N})"
            N=$((${N}+1))
        done
        autoexpex "${SFX}" ${A}
    done
    for a in ${SECT}
    do
        eval "AEXP=\${EXP${a}}"
        E=''
        for e in ${AEXP}
        do
            eval "E=\"${E} \${SFX${e}}\""
        done
        autosect ${a} ${E}
    done
}

cluster () {
    clexp () {
        e=${1}
        [ -r ${DST}${e}-ao.pto ] || return
        [ -z "${2}" ] || D="-D${2}"
        [ -z "${3}" ] || d="-d${3}"
        yrp=y,r,p/0-
        [ -z "${4}" ] || yrp="y,r,p/0-$((${4}-1)),$((${4}+1))-:r,p/${4}"
        perl ../../tools/ptsed -O0xff \
            -s -Vb,v/0:a,c,d,e/0-:${yrp} \
            -o ${DST}${e}-ao.txt ${DST}${e}-ao.pto
        PToptimizer ${DST}${e}-ao.txt
        perl ../../tools/ptscluster -O0xff \
            -x48 -y36 ${D} ${d} \
            -o ${DST}${e}-ao.txt ${DST}${e}-ao.txt \
        && perl ../../tools/ptsed -O0xff \
            -s -r \
            -o ${DST}${e}-ao.txt ${DST}${e}-ao.txt
        PToptimizer ${DST}${e}-ao.txt
        perl ../../tools/ptscluster -O0xff \
            -x24 -y18 ${D} ${d} \
            -o ${DST}${e}-ao.txt ${DST}${e}-ao.txt \
        && perl ../../tools/ptsed -O0xff \
            -s -r \
            -o ${DST}${e}-ao.txt ${DST}${e}-ao.txt
        PToptimizer ${DST}${e}-ao.txt
        perl ../../tools/ptscluster -O0xff \
            ${D} ${d} \
            -o ${DST}${e}-ao.txt ${DST}${e}-ao.txt \
        && perl ../../tools/ptsed -O0xff \
            -s -r \
            -o ${DST}${e}-ao.txt ${DST}${e}-ao.txt
        PToptimizer ${DST}${e}-ao.txt
        perl ../../tools/ptsed -O0xff \
            -i \
            -o ${DST}${e}.pto ${DST}${e}-ao.txt
    }
    clsect () {
        a=${1}
        a="$(printf '%02d' ${a})"
        [ -r ${DST}-${a}-ao.pto ] || return
        perl ../../tools/ptsed -O0xff \
            -Ib,v/1-:=0:a,c,d,e/0-:0 \
            -Vb/0:y,r,p/1- \
            -s -o ${DST}-${a}-ao.txt ${DST}-${a}-ao.pto
        PToptimizer ${DST}-${a}-ao.txt
        perl ../../tools/ptscluster -O0xff \
            -r0x21 -x48 -y36 -d12 -D24 \
            -o ${DST}-${a}-ao.txt ${DST}-${a}-ao.txt \
        && perl ../../tools/ptsed -O0xff \
            -s -r \
            -o ${DST}-${a}-ao.txt ${DST}-${a}-ao.txt
        PToptimizer ${DST}-${a}-ao.txt
        perl ../../tools/ptscluster -O0xff \
            -r0x21 -x48 -y36 -d8 -D16 \
            -o ${DST}-${a}-ao.txt ${DST}-${a}-ao.txt \
        && perl ../../tools/ptsed -O0xff \
            -s -r \
            -o ${DST}-${a}-ao.txt ${DST}-${a}-ao.txt
        PToptimizer ${DST}-${a}-ao.txt
        perl ../../tools/ptscluster -O0xff \
            -r0x21 -x48 -y36 -D10 \
            -o ${DST}-${a}-ao.txt ${DST}-${a}-ao.txt \
        && perl ../../tools/ptsed -O0xff \
            -s -r \
            -Vb/0:a,c,d,e/1-:y,r,p/1- \
            -o ${DST}-${a}-ao.txt ${DST}-${a}-ao.txt
        PToptimizer ${DST}-${a}-ao.txt
        perl ../../tools/ptscluster -O0xff \
            -r0x21 -x24 -y18 -D7 \
            -o ${DST}-${a}-ao.txt ${DST}-${a}-ao.txt \
        && perl ../../tools/ptsed -O0xff \
            -s -r \
            -o ${DST}-${a}-ao.txt ${DST}-${a}-ao.txt
        PToptimizer ${DST}-${a}-ao.txt
        perl ../../tools/ptscluster -O0xff \
            -r0x21 -D5 \
            -o ${DST}-${a}-ao.txt ${DST}-${a}-ao.txt \
        && perl ../../tools/ptsed -O0xff \
            -s -r \
            -o ${DST}-${a}-ao.txt ${DST}-${a}-ao.txt
        #PToptimizer ${DST}-${a}-ao.txt
        #perl ../../tools/ptscluster -O0xff \
        #    -r0x21 -D3 \
        #    -o ${DST}-${a}-ao.txt ${DST}-${a}-ao.txt \
        #&& perl ../../tools/ptsed -O0xff \
        #    -s -r \
        #    -o ${DST}-${a}-ao.txt ${DST}-${a}-ao.txt
        PToptimizer ${DST}-${a}-ao.txt
        perl ../../tools/ptsed -O0xff \
            -i \
            -o ${DST}-${a}.pto ${DST}-${a}-ao.txt
    }
    for e in ${EXP}
    do
        eval "SFX=\${SFX${e}}"
        clexp "${SFX}" "${DMAX}" "${DMIN}" "${ANC}"
    done
    for a in ${SECT}
    do
        clsect ${a}
    done
}

mkldrmf () {
    RM='rm -f'
    [ -z "${3}" ] || RM='echo'
    f=$(awk '/^p /{if(match($0,/ w[0-9]+/)){w=substr($0,RSTART+2,RLENGTH-2);match($0,/ h[0-9]+/);h=substr($0,RSTART+2,RLENGTH-2);print w"x"h}nextfile}' ${DST}.pto)
    [ -z "${f}" ] || f="-f ${f}"
    c=$(awk '/^p /{if(match($0,/ S[0-9]+,[0-9]+,[0-9]+,[0-9]+/)){split(substr($0,RSTART+2,RLENGTH-2),a,/,/);print(a[2]-a[1])"x"(a[4]-a[3])"+"a[1]"+"a[3]}nextfile}' ${DST}.pto)
    [ -z "${c}" ] || f="-f ${c}"
    w=''
    [ "${3}" = a ] || w='-w'
    cat <<EOT
DST?=${DST}

CPUS?=${CPUS}
MEM?=${MEM}
DEPTH?=${DEPTH}

CONV?=${CONV}
REMAP?=${REMAPLDR}
BLEND?=${BLEND} ${w} ${f}
EFUSE?=${EFUSE}

RM?=${RM}
MV?=mv

all: \${DST}.tiff

EOT
    SRC=''
    BLIN=''
    for a in ${SECT}
    do
        FUIN=''
        FOUT=''
        DEP=''
        t=''
        eval "AEXP=\${EXP${a}}"
        nexp=$(echo "${AEXP}" | wc -w)
        for e in ${AEXP}
        do
            DEP="${FOUTDEP} ${t}"
            eval "$(echo ${e} | awk 'BEGIN{FS=":"}{print "e="$1"; aa="$2}')"
            eval "SFX=\${SFX${e}}"
            [ -n "${aa}" ] || aa=${a}
            n="$(printf '%02d%s.tiff' ${aa} ${SFX})"
            SRC="${SRC} ${n}"
            eval "$(awk 'BEGIN{i=0}/^i /{if(match($0," n\"?"n"\"?")){print "I="i;if(match($0,/ Eev[^ ]+/))print "E="substr($0,RSTART+4,RLENGTH-4);nextfile}++i}' n=${n} ${DST}.pto)"
            [ -n "${I}" ] || return 1
            [ -n "${E}" ] || E=0
            E="-e ${E}"
            [ ${nexp} -gt 1 ] || E=''
            t=$(printf 'tmp-%04d.%s' ${I} ${LDRSFX})
            FUIN="${FUIN} ${t}"
            [ ${a} -eq ${aa} ] || continue
            FOUT=${t}
            DEP="$(echo "${DEP}"|sed -e 's/^ *\(.*\) *$/\1/')"
            cat <<EOT
${t}: \${DST}.pto ${n} ${DEP}
	@\${REMAP} -i ${I} ${E} \${DST}.pto
EOT
        done
        FUIN="$(echo "${FUIN}"|sed -e 's/^ *\(.*\) *$/\1/')"
        [ "${FUIN}" = "${FOUT}" ] || \
        cat <<EOT
	#@env CONV='\${CONV}' RM='\${RM}' sh ../../tools/stackmask ${FUIN}
	@\${EFUSE} -o .${FOUT} ${FUIN}
	@\${RM} ${FUIN}
	@\${MV} .${FOUT} ${FOUT}
EOT
        BLIN="${BLIN} ${FOUT}"
    done
    SRC="$(echo "${SRC}"|sed -e 's/^ *\(.*\) *$/\1/')"
    BLIN="$(echo "${BLIN}"|sed -e 's/^ *\(.*\) *$/\1/')"
    cat <<EOT

SRC=${SRC}
BLIN=${BLIN}

tmp.tiff: \${BLIN}
	@\${RM} \${SRC}
	@\${BLEND} -o tmp.tiff \${BLIN}
	@\${RM} \${BLIN}
\${DST}.tiff: tmp.tiff
	@pfsin tmp.tiff \\
	 | pfstmo_mantiuk06 -v -f 0.6 -s 1 \\
	 | pfsgamma -m 0.8 -g 2.2 \\
	 | pfsoutimgmagick \${DST}.tiff
	@\${RM} tmp.tiff
	@\${CONV} \${DST}.tiff \\
		+matte +repage \\
		-unsharp 0x0.3 \\
		-compress zip \${DST}.tiff
EOT
}

mkhdrmf () {
    RM='rm -f'
    [ -z "${3}" ] || RM='echo'
    f=$(awk '/^p /{if(match($0,/ w[0-9]+/)){w=substr($0,RSTART+2,RLENGTH-2);match($0,/ h[0-9]+/);h=substr($0,RSTART+2,RLENGTH-2);print w"x"h}nextfile}' ${DST}.pto)
    [ -z "${f}" ] || f="-f ${f}"
    c=$(awk '/^p /{if(match($0,/ S[0-9]+,[0-9]+,[0-9]+,[0-9]+/)){split(substr($0,RSTART+2,RLENGTH-2),a,/,/);print(a[2]-a[1])"x"(a[4]-a[3])"+"a[1]"+"a[3]}nextfile}' ${DST}.pto)
    [ -z "${c}" ] || f="-f ${c}"
    w=''
    [ "${3}" = a ] || w='-w'
    cat <<EOT
DST?=${DST}

CPUS?=${CPUS}
MEM?=${MEM}
DEPTH?=${DEPTH}

CONV?=${CONV}
REMAP?=${REMAPHDR}
BLEND?=${BLEND} ${w} ${f}
MERGE?=${MERGE}

RM?=${RM}
MV?=mv

all: \${DST}.tiff

EOT
    SRC=''
    BLIN=''
    for a in ${SECT}
    do
        FUIN=''
        MASK=''
        FOUT=''
        DEP=''
        t=''
        eval "AEXP=\${EXP${a}}"
        for e in ${AEXP}
        do
            DEP="${FOUTDEP} ${t}"
            eval "$(echo ${e} | awk 'BEGIN{FS=":"}{print "e="$1"; aa="$2}')"
            eval "SFX=\${SFX${e}}"
            [ -n "${aa}" ] || aa=${a}
            n="$(printf '%02d%s.tiff' ${aa} ${SFX})"
            SRC="${SRC} ${n}"
            eval "$(awk 'BEGIN{i=0}/^i /{if(match($0," n\"?"n"\"?")){print "I="i;if(match($0,/ Eev[^ ]+/))print "E="substr($0,RSTART+4,RLENGTH-4);nextfile}++i}' n=${n} ${DST}.pto)"
            [ -n "${I}" ] || return 1
            [ -n "${E}" ] || E=0
            t=$(printf 'tmp-%04d.%s' ${I} ${HDRSFX})
            m=$(printf 'tmp-%04d_gray.pgm' ${I})
            FUIN="${FUIN} ${t}"
            MASK="${MASK} ${m}"
            [ ${a} -eq ${aa} ] || continue
            FOUT=${t}
            DEP="$(echo "${DEP}"|sed -e 's/^ *\(.*\) *$/\1/')"
            cat <<EOT
${t}: \${DST}.pto ${n} ${DEP}
	@\${REMAP} -i ${I} \${DST}.pto
EOT
        done
        FUIN="$(echo "${FUIN}"|sed -e 's/^ *\(.*\) *$/\1/')"
        MASK="$(echo "${MASK}"|sed -e 's/^ *\(.*\) *$/\1/')"
        [ "${FUIN}" = "${FOUT}" ] || \
        cat <<EOT
	@\${MERGE} -o .${FOUT} ${FUIN}
	@\${RM} ${FUIN}
	@\${MV} .${FOUT} ${FOUT}
EOT
        cat <<EOT
	@\${RM} ${MASK}
EOT
        BLIN="${BLIN} ${FOUT}"
    done
    SRC="$(echo "${SRC}"|sed -e 's/^ *\(.*\) *$/\1/')"
    BLIN="$(echo "${BLIN}"|sed -e 's/^ *\(.*\) *$/\1/')"
    resize=''
    [ "${3}" = a ] \
    || resize='-unsharp 0x0.3 -resize x4096 -unsharp 0x0.75+0.6+0.05'
    [ "${3}" = b ] \
    || resize='-unsharp 0x0.75+0.6+0.05'
    cat <<EOT

SRC=${SRC}
BLIN=${BLIN}

tmp.tiff: \${BLIN}
	@\${RM} \${SRC}
	@\${BLEND} -o tmp.tiff \${BLIN}
	@\${RM} \${BLIN}
\${DST}.tiff: tmp.tiff
	@pfsin tmp.tiff \\
	 | pfsclamp --min 1e-20: --max 1e+20 \\
	 | pfstmo_mantiuk06 -v -f 0.6 -s 1 \\
	 | pfsgamma -m 0.8 -g 2.2 \\
	 | pfsoutimgmagick \${DST}.tiff
	@\${RM} tmp.tiff
	@\${CONV} \${DST}.tiff \\
		+matte +repage \\
		${resize} \\
		-compress zip \${DST}.tiff
EOT
#	 | pfstmo_mantiuk06 -v -f 0.45 -s 0.75 \\
#	 | pfsgamma -m 0.825 -g 2.2 \\
#	 | pfstmo_mantiuk08 -v -d pd=crt -c 0.40 -e 1.00 \\
#	 | pfsgamma -m 0.65 -g 1.00 \\
#	 | pfstmo_reinhard02 -v -k 0.090 \\
#	 | pfsgamma -m 0.825 -g 2.185 \\
#	 | pfstmo_reinhard02 -v -k 0.045 \\
#	 | pfsgamma -m 0.825 -g 2.185 \\
}

inpaint () {
    g=${1}
    n=${2}
    ${CONV} ${DST}.tiff \
        -crop ${g} \
        -compress zip ${DST}.${n}.tiff 
    mv ${DST}.tiff ${DST}.0.tiff
    greycstoration -inpaint ${DST}.${n}.tiff \
        -iter 1024 -bits 16 \
        -m mask-${DST}.${n}.png \
        -o ${DST}.${n}.tiff
    ${CONV} ${DST}.0.tiff ${DST}.${n}.tiff \
        -geometry ${g} \
        -compose src-over -composite \
        -compress zip ${DST}.tiff
    rm ${DST}.0.tiff ${DST}.${n}.tiff
}

apply-patch () 
{
    g=${1}; shift
    se=${1}; shift
    de=${1}; shift
    n=${1}; shift
    [ -n ${n} ] || n=1
    rgb=''; srgb=''
    [ -n ${LIN} ] || { rgb="${RGB}"; srgb="${sRGB}"; }
    ${CONV} ${DST}-e${se}.tiff \
        -crop ${g} +repage \
        ${rgb} "${@}" ${srgb} \
        -compress zip ${DST}-e${de}.${n}.tiff 
    [ -r mask-${DST}-e${de}.${n}.png ] \
    && ${CONV} ${DST}-e${de}.${n}.tiff \
        mask-${DST}-e${de}.${n}.png -compose copyopacity -composite \
        -compress zip ${DST}-e${de}.${n}.tiff
    mv ${DST}-e${de}.tiff ${DST}-e${de}.0.tiff
    ${CONV} ${DST}-e${de}.0.tiff ${DST}-e${de}.${n}.tiff \
        -geometry ${g} \
        -compose src-over -composite \
        -compress zip ${DST}-e${de}.tiff
    rm ${DST}-e${de}.0.tiff ${DST}-e${de}.${n}.tiff
}

mkpano () {
    f=$(awk '/^p /{if(match($0,/w[0-9]+/)){w=substr($0,RSTART+1,RLENGTH-1);match($0,/h[0-9]+/);h=substr($0,RSTART+1,RLENGTH-1);print w"x"h}nextfile}' ${DST}.pto)
    [ -z "${f}" ] || f="-f ${f}"
    c=$(awk '/^p /{if(match($0,/S[0-9]+,[0-9]+,[0-9]+,[0-9]+/)){split(substr($0,RSTART+1,RLENGTH-1),a,/,/);print(a[2]-a[1])"x"(a[4]-a[3])"+"a[1]"+"a[3]}nextfile}' ${DST}.pto)
    [ -z "${c}" ] || f="-f ${c}"
    #[ -z "${c}" ] || c="-crop ${c} +repage"
    for e in ${EXP}
    do
        eval "SFX=\${SFX${e}}"
        [ -r ${DST}${SFX}.tiff ] && continue
        NIN=''
        BLIN=''
        for a in ${SECT}
        do
            eval "AEXP=\${EXP${a}}"
            for ee in ${AEXP}
            do
                eval "$(echo ${ee} | awk 'BEGIN{FS=":"}{print "ee="$1"; aa="$2}')"
                [ ${e} -eq ${ee} ] || continue
                [ -n "${aa}" ] || aa=${a}
                n="$(printf '%02d%s.tiff' ${aa} ${SFX})"
                eval "$(awk 'BEGIN{i=0}/^i /{if(match($0," n\"?"n"\"?")){print "I="i;if(match($0,/ Eev[^ ]+/))print "E="substr($0,RSTART+4,RLENGTH-4);nextfile}++i}' n=${n} ${DST}.pto)"
                [ -n "${I}" ] || continue
                NIN="${NIN} -i ${I}"
                BLIN=$(printf "${BLIN} tmp-%04d.tif" ${I})
            done
        done
        ${REMAPLDR} ${NIN} ${DST}.pto \
        && ${BLEND} ${f} -o tmp${SFX}.tiff ${BLIN} \
        && rm ${BLIN} \
        && ${CONV} tmp${SFX}.tiff \
            +matte +repage \
            -unsharp 0x0.2 \
            -compress zip ${DST}${SFX}.tiff \
        && rm tmp${SFX}.tiff
            #-modulate 100,95 \
    done

    [ -n "${3}" ] \
    || for e in ${EXP}
    do
        eval "SFX=\${SFX${e}}"
        rm [0-9][0-9]${SFX}.tiff
    done
    #apply-patch 3072x864+0+512 04 05 1  -modulate  50.000,100
    #apply-patch 3072x864+0+512 04 06 1  -modulate  25.000,100
    #apply-patch 64x64+1440+1472 05 04 2 -modulate 200.000,100
    #apply-patch 64x64+1440+1472 05 03 2 -modulate 400.000,100
    #apply-patch 3072x864+0+512 04 05 1 -evaluate multiply 0.495 -modulate 100,126
    #apply-patch 3072x864+0+512 04 06 1 -evaluate multiply 0.225 -modulate 100,129
    #apply-patch 64x64+1440+1472 05 04 2 -modulate 100,75 -evaluate multiply 2.20
    #apply-patch 64x64+1440+1472 05 03 2 +contrast -modulate 100,65 -evaluate multiply 3.70

    #inpaint 416x64+1536+1920 1
}

mkhdr () {
    #pfsinhdrgen ${DST}.hdrgen \
    #| pfshdrcalibrate -v -r gamma -g 8 -b 16 -s ${DST}.0.m
    PFSCAL='pfshdrcalibrate -v -x -r gamma -g 8'
    [ -r ${DST}.m ] \
    && PFSCAL="pfshdrcalibrate -v -x -c none -f ${DST}.m"
    PFSCNT=2
    [ -r ${DST}.hdrgen ] \
    && PFSCNT=$(wc -l ${DST}.hdrgen|sed -e 's/ *\([0-9]*\) *.*/\1/')
    INIMG=${DST}.tiff
    [ -r ${DST}.hdrgen ] \
    && INIMG=$(awk '{print $1;nextfile}' ${DST}.hdrgen)
    PFSSIZE=$(identify -format "%w*%h*${PFSCNT}" ${INIMG}|bc)
    PFSEVAL=eval
    [ ${PFSSIZE} -le ${PFSMAX} ] || PFSEVAL="rsh ${PFSHOST}"
    PFSIN="pfsintiff --linear ${DST}.tiff"
    [ -r ${DST}.hdrgen ] \
    && PFSIN="pfsinhdrgen ${DST}.hdrgen | ${PFSCAL} -b ${DEPTH}"
    [ -r ${DST}.pfs.bz2 ] \
    && PFSIN="bzcat ${DST}.pfs.bz2"
    PFSCHAIN="cd $(pwd) && ${PFSIN}"

    cat >&2 <<EOM
Image size is ${PFSSIZE} pxel(s)
Max allowed for local operation is ${PFSMAX} pxel(s)
Pipeline
  ${PFSCHAIN}
will be executed using '${PFSEVAL}' command
EOM

    ${PFSEVAL} "${PFSCHAIN}" \
    | pfswb -v \
    | pfstmo_mantiuk08 -v -d pd=crt -c 0.75 -e 1.00 \
    | pfsgamma -m 0.95 -g 1.00 \
    | pfsoutimgmagick ${DST}.0.tiff
    [ ${?} -ne 0 -o -n "${3}" ] \
    || for e in ${EXP}
       do
        eval "SFX=\${SFX${e}}"
        [ -z "${SFX}" ] || rm ${DST}${SFX}.tiff
       done
    ${CONV} ${DST}.tiff \
        -sigmoidal-contrast 4x82% \
        -matte -channel A \
        ${DST}.0.tiff -fx '0.4+0.55*(v.luminance^3)' \
        ${DST}.0.tiff +swap \
        -compose src-over -composite \
        -compress zip ${DST}.tiff
    rm ${DST}.0.tiff
    #test -r ${DST}.tiff
}

qtvr () {
    ${CONV} ${DST}.tiff \
        -recolor '1.25 -0.1 -0.15  -0.1 1.25 -0.15  -0.05 -0.05 1.1' \
        -sigmoidal-contrast 3x32% \
        ${@} \
        -compress zip ${DST}-qtvr.tiff
    erect2qtvr --erect=${DST}-qtvr.tiff \
        --quality=85 --face=2048 \
        --width=1280 --height=960 \
        --outfile=${DST}.mov
    rm ${DST}-qtvr.tiff
}

final-a () {
    SFX=${1}
    shift
    h=$(identify -format '%h' ${DST}.tiff)
    y=$((${h}/2))
    r=$((${h}/64))
    s=$((${r}/3))
    echo "h=${h}; y=${y}; blur ${r}x${s}"
    ${CONV} \
        \( ${DST}.tiff \
           -recolor '1.25 -0.1 -0.15  -0.1 1.25 -0.15  -0.05 -0.05 1.1' \
           -sigmoidal-contrast 3x32% \
        \) \
        \(  -colorspace gray -virtual-pixel black -size ${h}x${h} xc:black \
            -fill white -box white -draw "circle ${y},${y} ${y},1" \
            -blur ${r}x${s} \
        \) \
        -compose multiply -composite \
        ${@} \
        -quality 90 ${DST}${SFX}.jpeg
    r=${?}
    return ${r}
}

final () {
    SFX=${1}
    shift
    ${CONV} ${DST}.tiff \
        -negate -gamma 1.2 -negate -modulate 115 -gamma 1.4 \
        -sigmoidal-contrast 2.25x5% \
        -recolor '1.15 -0.0825 -0.0725  -0.05 1.125 -0.075  -0.0 -0.0 1.0' \
        ${@} \
        -quality 90 ${DST}${SFX}.jpeg
    r=${?}
    return ${r}
}

prep "${@}"
[ "${1}" = prep ] && exit
[ "${1}" = auto ] && { auto "${@}"; exit; }
[ "${1}" = cluster ] && { cluster "${@}"; exit; }
[ "${1}" = mkldrmf ] && { mkldrmf "${@}"; exit; }
[ "${1}" = mkhdrmf ] && { mkhdrmf "${@}"; exit; }
mkpano "${@}" || exit ${?}
#{ mkldrmf "${@}" >Makefile && make && { [ -n "${3}" ] || rm Makefile; } } || exit ${?}
#{ mkhdrmf "${@}" >Makefile && make && { [ -n "${3}" ] || rm Makefile; } } || exit ${?}
#inpaint 64x64+1760+1376 1
[ "${1}" = pano ] && exit
#mkhdr "${@}" || exit ${?}
#[ "${1}" = hdr ] && exit
#[ "${3}" = a -o "${3}" = b ] || qtvr 
#resize='-resize x3072'
#[ "${3}" = a ] && resize=''
final=final
#[ "${3}" = a ] && final='final-a'
echo ${final}
${CONV} ${DST}.tiff ${resize} -unsharp 0x0.75+0.6+0.05 -compress zip ${DST}.tiff \
&& ${final} -huge \
&& ${CONV} ${DST}.tiff -resize x1536 -unsharp 0x0.75+0.6+0.04 -compress zip ${DST}.tiff \
&& ${final} -big \
&& ${CONV} ${DST}.tiff -resize x768 -unsharp 0x0.65+0.5+0.04 -compress zip ${DST}.tiff \
&& ${final} -middle \
&& ${CONV} ${DST}.tiff -resize x384 -unsharp 0x0.55+0.5+0.04 -compress zip ${DST}.tiff \
&& ${CONV} ${DST}.tiff -resize x256 -unsharp 0x0.45+0.4+0.03 -compress zip ${DST}.tiff \
&& ${final} -small -gamma 1.10 \
&& rm ${DST}.tiff
exit