#!/bin/sh

[ -r .conf ] || exit 125

#UFRAW='eval ufraw-batch --interpolation=ahd --out-type=tiff --out-depth=16 --overwrite --zip --wavelet-denoising-threshold=300'
DCRAW='dcraw -v -r 1.110000 1.000000 2.895000 1.000000 -S 62976 -H 0 -q 3 -n 150 -m 4 -4 -T'
#DCRAW='dcraw -v -w -H 0 -q 3 -n 300 -4 -T'
NONA='eval nona -v -t ${CPUS} -m TIFF_m -r ldr -p UINT${DEPTH} -z DEFLATE -o tmp-'
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 -otmp.tiff'
ENBL='enblend -z -m 896 --fine-mask'
#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'

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

prep () {
    for e in ${EXP}
    do
        eval "SEQ=\${SEQ${e}}"
        eval "SRC=\${SRC${e}}"
        eval "SFX=\${SFX${e}}"
        N=0
        for n in ${SRC}
        do
            [ -z "${n}" -o "${n}" = SKEEP ] && continue
            o=$(printf '%02d%s' ${N} ${SFX})
            N=$((${N}+1))
            [ -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 "${DCRAW}" ]
                then
                    conf=${n}.ufraw
                    [ -z "${CONF}" ] || conf=${CONF}.ufraw
                    [ -r ${conf} ] || exit 125
                    ${UFRAW} --conf=${conf} --output=${n}.tiff ${n}.RAW
                else
                    ${DCRAW} ${n}.RAW
                fi
                exiftool -overwrite_original -P -TagsFromFile ${n}.RAW ${n}.tiff
                fl=$(dcraw -i -v ${n}.RAW|awk '/^Focal length: /{print $3}')
                [ -r ${n}.RAW.bz2 -o -r ${n}.RAW.gz ] && rm -f ${n}.RAW
                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 "${DCRAW}" -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} -gamma 2.2 \
                        -compress zip ${o}.tiff
                fi
                [ "${DEPTH}" = 16 ] \
                || ${CONV} ${o}.tiff -depth 8 -compress zip ${o}.tiff
            elif [ -r ${n}.JPG ]
            then
                ${CONV} ${n}.JPG \
                    -depth ${DEPTH} \
                    -orient undefined \
                    -unsharp 0x0.3 \
                    -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 -channel A -fx 'v.r' \
                -compress zip ${o}.tiff
        done
    done
}

auto () {
    SIZE=2048
    E=''
    S=''
    for e in ${EXP}
    do
        eval "SRC=\${SRC${e}}"
        eval "SFX=\${SFX${e}}"
        E="${E} ${SFX}"
        A=''
        N=0
        for n in ${SRC}
        do
            A="${A} $(printf '%02d' ${N})"
            [ -z "${n}" -o "${n}" = SKEEP ] && continue
            S="${S} $(printf '%02d%s.tiff' ${N} ${SFX})"
            N=$((${N}+1))
        done
    done
    #printf 'E=%s\nA=%s\nS=%s\n' "${E}" "${A}" "${S}"; return 0
    N=0
    for n in ${S}
    do
        [ -r ${n}.key.gz ] && continue
        mono ../../tools/autopano-sift-2.4/bin/generatekeys-sd.exe \
            ${n} ${n}.key.gz ${SIZE} &
        N=$((${N}+1))
        [ $((${N}%${CPUS})) -eq 0 ] && wait
    done
    N=0
    for a in ${A}
    do
        [ -r ${DST}-ap${a}.pto ] && continue
        S=''
        for e in ${E}
        do
            s="$(printf '%s%s.tiff.key.gz' ${a} ${e})"
            [ -r "${s}" ] && S="${S} ${s}"
        done
        mono ../../tools/autopano-sift-2.4/bin/autopano.exe \
            --ransac on --maxmatches 200 \
            --refine --keep-unrefinable off \
            ${DST}-ap${a}.pto ${S} 2>&1 \
            | grep -v '^peak-search at level ' &
        N=$((${N}+1))
        [ $((${N}%${CPUS})) -eq 0 ] && wait
    done
    N=0
    for e in ${E}
    do
        [ -r ${DST}-ap${e}.pto ] && continue
        S=''
        for a in ${A}
        do
            s="$(printf '%s%s.tiff.key.gz' ${a} ${e})"
            [ -r "${s}" ] && S="${S} ${s}"
        done
        mono ../../tools/autopano-sift-2.4/bin/autopano.exe \
            --ransac on --maxmatches 200 \
            --refine --keep-unrefinable off \
            ${DST}-ap${e}.pto ${S} 2>&1 \
            | grep -v '^peak-search at level ' &
        N=$((${N}+1))
        [ $((${N}%${CPUS})) -eq 0 ] && wait
    done
    wait
    rm *.key.gz
}

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 "SEQ=\${SEQ${e}}"
        NIN=''
        EIN=''
        for s in ${SEQ}
        do
            NIN="${NIN} -i ${s}"
            EIN=$(printf "${EIN} tmp-%04d.tif" ${s})
        done
        eval "SFX=\${SFX${e}}"
        ${NONA} ${NIN} ${DST}.pto \
        && ${ENBL} ${f} -o tmp${SFX}.tiff ${EIN} \
        && rm ${EIN} \
        && ${CONV} tmp${SFX}.tiff \
            +matte \
            -unsharp 0x0.25 \
            -compress zip ${DST}${SFX}.tiff \
        && rm tmp${SFX}.tiff
    done

    [ -n "${3}" ] \
    || for e in ${EXP}
    do
        eval "SFX=\${SFX${e}}"
        rm [0-9][0-9]${SFX}.tiff
    done
}

mkhdr () {
    #pfsinhdrgen ${DST}.hdrgen \
    #| pfshdrcalibrate -v -r linear -g 8.8 -b 16 -s ${DST}.0.m
    PFSCAL='pfshdrcalibrate -v -r linear -g 8.8'
    [ -r ${DST}.m ] \
    && PFSCAL="pfshdrcalibrate -v -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)
    PFSIN="pfsinimgmagick --linear ${DST}.tiff"
    [ -r ${DST}.hdrgen ] \
    && PFSIN="pfsinhdrgen ${DST}.hdrgen | ${PFSCAL} -b ${DEPTH}"
    PFSEVAL=eval
    [ ${PFSSIZE} -le ${PFSMAX} ] || PFSEVAL="rsh ${PFSHOST}"
#    PFSCHAIN="cd $(pwd) \
#&& ${PFSIN} \
#| pfstmo_mantiuk06 -v -f 0.45 -s 0.7 \
#| pfsgamma -m 0.75 -g 2.2 \
#| pfsoutimgmagick ${DST}.tiff \
#"
#    PFSCHAIN="cd $(pwd) \
#&& ${PFSIN} \
#| pfstmo_mantiuk06 -v -f 0.55 -s 0.7 \
#| pfsgamma -m 0.7 -g 2.2 \
#| pfsoutimgmagick ${DST}.tiff \
#"
    PFSCHAIN="cd $(pwd) \
&& ${PFSIN} \
| pfswb -v -k -b 0.9 \
| pfstmo_reinhard05 -v -l 0.8 -c 0.7 -b -0.33 \
| pfsoutimgmagick ${DST}.tiff \
"
#    [ "${3}" = e -o "${3}" = m ] \
#    || PFSCHAIN="cd $(pwd) \
#&& ${PFSIN} \
#| pfstmo_reinhard02 -v -k 0.0225 -s -p 9 -r 6 -u 17 \
#| pfsgamma -m 0.75 -g 2.2 \
#| pfsoutimgmagick ${DST}.tiff \
#"
#    [ "${3}" = e -o "${3}" = m ] \
#    || PFSCHAIN="cd $(pwd) \
#&& ${PFSIN} \
#| pfswb -v -k -b 0.9 \
#| pfstmo_reinhard05 -v -l 0.7 -c 0.7 -b -1.33 \
#| pfsoutimgmagick ${DST}.tiff \
#"
    [ "${3}" = e -o "${3}" = m ] \
    || PFSCHAIN="cd $(pwd) \
&& ${PFSIN} \
| pfswb -v -k -r 1.05 -b 0.9 \
| pfstmo_reinhard05 -v -l 0.55 -c 0.7 -b -0.67 \
| pfsoutimgmagick ${DST}.tiff \
"

    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}" \
    && for e in ${EXP}
       do
        eval "SFX=\${SFX${e}}"
        rm ${DST}${SFX}.tiff
       done
    #&& greycstoration -restore ${DST}.tiff \
    #-a 0.9 -p 2 -dt 20.0 -iter 8 \
    #-v -bits ${DEPTH} \
    #-o ${DST}.tiff

    if [ "${3}" = e ]
    then
        g=512x416+1408+640
        ${CONV} ${DST}.tiff \
            -crop ${g} \
            -compress zip ${DST}.1.tiff 
        mv ${DST}.tiff ${DST}.0.tiff
        greycstoration -inpaint ${DST}.1.tiff \
            -iter 1024 -bits 16 \
            -m mask-${DST}.1.png \
            -o ${DST}.1.tiff
        ${CONV} ${DST}.0.tiff ${DST}.1.tiff \
            -geometry ${g} \
            -compose src-over -composite \
            -compress zip ${DST}.tiff
        rm ${DST}.0.tiff ${DST}.1.tiff
    elif [ "${3}" = m ]
    then
        g=544x512+1408+800
        ${CONV} ${DST}.tiff \
            -crop ${g} \
            -compress zip ${DST}.1.tiff 
        mv ${DST}.tiff ${DST}.0.tiff
        greycstoration -inpaint ${DST}.1.tiff \
            -iter 1024 -bits 16 \
            -m mask-${DST}.1.png \
            -o ${DST}.1.tiff
        ${CONV} ${DST}.0.tiff ${DST}.1.tiff \
            -geometry ${g} \
            -compose src-over -composite \
            -compress zip ${DST}.tiff
        rm ${DST}.0.tiff ${DST}.1.tiff
    else
        g=544x480+640+608
        ${CONV} ${DST}.tiff \
            -crop ${g} \
            -compress zip ${DST}.1.tiff 
        mv ${DST}.tiff ${DST}.0.tiff
        greycstoration -inpaint ${DST}.1.tiff \
            -iter 1024 -bits 16 \
            -m mask-${DST}.1.png \
            -o ${DST}.1.tiff
        ${CONV} ${DST}.0.tiff ${DST}.1.tiff \
            -geometry ${g} \
            -compose src-over -composite \
            -compress zip ${DST}.tiff
        rm ${DST}.0.tiff ${DST}.1.tiff
    fi
}

final () {
    SFX=${1}
    shift
    #${CONV} ${DST}.tiff \
    #    -level 6.25%,100%,1.2 \
    #    -sigmoidal-contrast 3.5x8.6% \
    #    ${@} \
    #    -modulate 100,112.5 \
    #    -quality 90 ${DST}${SFX}.jpeg
    ${CONV} ${DST}.tiff \
        -strip \
        -gamma 1.1 \
        ${@} \
        -quality 90 ${DST}${SFX}.jpeg
    r=${?}
    return ${r}
}

prep "${@}"
[ "${1}" = prep ] && exit
[ "${1}" = auto ] && { auto "${@}"; exit; }
mkpano "${@}" || exit ${?}
[ "${1}" = pano ] && exit
#mkhdr "${@}" || exit ${?}
#[ "${1}" = hdr ] && exit
#HH='-resize x2816'
#[ "${3}" = e ] || HH=''
${CONV} ${DST}.tiff -resize x3072 -unsharp 0x0.75+0.6+0.05 -compress zip ${DST}.tiff \
&& final -huge \
&& ${CONV} ${DST}.tiff -resize x1536 -unsharp 0x0.75+0.6+0.05 -compress zip ${DST}.tiff \
&& final -big \
&& ${CONV} ${DST}.tiff -resize x768 -unsharp 0x0.75+0.5+0.04 -compress zip ${DST}.tiff \
&& final -middle \
&& ${CONV} ${DST}.tiff -resize x384 -unsharp 0x0.65+0.5+0.04 -compress zip ${DST}.tiff \
&& ${CONV} ${DST}.tiff -resize x256 -unsharp 0x0.65+0.4+0.03 -compress zip ${DST}.tiff \
&& final -small -gamma 1.05 \
&& rm ${DST}.tiff
exit