| 1 |
#!/bin/csh -fe |
| 2 |
# Convert Radiance animation frames to TIFF output |
| 3 |
# |
| 4 |
set histosiz=200 |
| 5 |
set pfwt=0.9 |
| 6 |
set outdir="" |
| 7 |
set pcopts=() |
| 8 |
set tfopts=() |
| 9 |
if (! $#argv) set argv=(DUMMY) |
| 10 |
# Process options for pcond and ra_tiff |
| 11 |
while ("$argv[1]" =~ -*) |
| 12 |
switch ("$argv[1]") |
| 13 |
case -W: |
| 14 |
shift argv |
| 15 |
set pfwt=$argv[1] |
| 16 |
breaksw |
| 17 |
case -H: |
| 18 |
shift argv |
| 19 |
set histof=$argv[1]:q |
| 20 |
breaksw |
| 21 |
case -D: |
| 22 |
shift argv |
| 23 |
if (! -d $argv[1]:q ) then |
| 24 |
echo "Directory $argv[1] does not exist" |
| 25 |
exit 1 |
| 26 |
endif |
| 27 |
set outdir=$argv[1]:q/ |
| 28 |
breaksw |
| 29 |
case -h*: |
| 30 |
case -a*: |
| 31 |
case -v*: |
| 32 |
case -s*: |
| 33 |
case -c*: |
| 34 |
case -l*: |
| 35 |
set pcopts=($pcopts $argv[1]) |
| 36 |
breaksw |
| 37 |
case -u: |
| 38 |
case -d: |
| 39 |
case -f: |
| 40 |
set pcopts=($pcopts $argv[1-2]) |
| 41 |
shift argv |
| 42 |
breaksw |
| 43 |
case -p: |
| 44 |
shift argv |
| 45 |
set pcopts=($pcopts -p $argv[1-6]) |
| 46 |
shift argv; shift argv; shift argv; shift argv; shift argv |
| 47 |
breaksw |
| 48 |
case -z: |
| 49 |
case -b: |
| 50 |
case -w: |
| 51 |
set tfopts=($tfopts $argv[1]) |
| 52 |
breaksw |
| 53 |
case -g: |
| 54 |
shift argv |
| 55 |
set tfopts=($tfopts -g $argv[1]) |
| 56 |
breaksw |
| 57 |
default: |
| 58 |
echo "$0: bad option: $argv[1]" |
| 59 |
exit 1 |
| 60 |
endsw |
| 61 |
shift argv |
| 62 |
end |
| 63 |
if ($#argv < 2) then |
| 64 |
echo Usage: "$0 [-W prev_frame_wt][-H histo][-D dir][pcond opts][ra_tiff opts] frame1 frame2 .." |
| 65 |
exit 1 |
| 66 |
endif |
| 67 |
# Get shrunken image luminances |
| 68 |
set vald=`mktemp -d /tmp/val.XXXXXX` |
| 69 |
foreach inp ($argv:q) |
| 70 |
set datf="$inp:t" |
| 71 |
set datf="$vald/$datf:r.dat" |
| 72 |
pfilt -1 -x 128 -y 128 -p 1 $inp:q \ |
| 73 |
| pvalue -o -h -H -b -df \ |
| 74 |
| rcalc -if1 -e 'L=$1*179;cond=L-1e-7;$1=log10(L)' \ |
| 75 |
> $datf:q |
| 76 |
end |
| 77 |
# Get Min. and Max. log values |
| 78 |
set Lmin=`cat $vald/*.dat | total -l | rcalc -e '$1=$1-.01'` |
| 79 |
set Lmax=`cat $vald/*.dat | total -u | rcalc -e '$1=$1+.01'` |
| 80 |
if ($?histof) then |
| 81 |
if (-r $histof) then |
| 82 |
# Fix min/max and translate histogram |
| 83 |
set Lmin=`sed -n '1p' $histof | rcalc -e 'min(a,b):if(a-b,b,a);$1=min($1,'"$Lmin)"` |
| 84 |
set Lmax=`sed -n '$p' $histof | rcalc -e 'max(a,b):if(a-b,a,b);$1=max($1,'"$Lmax)"` |
| 85 |
tabfunc -i hfunc < $histof > $vald/oldhist.cal |
| 86 |
cnt $histosiz \ |
| 87 |
| rcalc -e "L10=$Lmin+($Lmax-$Lmin)/$histosiz"'*($1+.5)' \ |
| 88 |
-f $vald/oldhist.cal -e '$1=L10;$2=hfunc(L10)' \ |
| 89 |
> $vald/oldhisto.dat |
| 90 |
endif |
| 91 |
endif |
| 92 |
foreach inp ($argv:q) |
| 93 |
set datf="$inp:t" |
| 94 |
set datf="$vald/$datf:r.dat" |
| 95 |
set outp="$inp:t" |
| 96 |
set outp="$outdir$outp:r.tif" |
| 97 |
endif |
| 98 |
histo $Lmin $Lmax $histosiz < $datf > $vald/newhisto.dat |
| 99 |
if (-f $vald/oldhisto.dat) then |
| 100 |
rlam $vald/newhisto.dat $vald/oldhisto.dat \ |
| 101 |
| rcalc -e '$1=$1;$2=$2+$4*'$pfwt \ |
| 102 |
> $vald/histo.dat |
| 103 |
else |
| 104 |
mv $vald/{new,}histo.dat |
| 105 |
endif |
| 106 |
pcond $pcopts -I $inp:q < $vald/histo.dat \ |
| 107 |
| ra_tiff $tfopts - $outp:q |
| 108 |
mv $vald/{,old}histo.dat |
| 109 |
end |
| 110 |
if ($?histof) then |
| 111 |
cp -f $vald/oldhisto.dat $histof |
| 112 |
endif |
| 113 |
rm -rf $vald |
