| 3 |
|
# |
| 4 |
|
# Compute foveal histogram for picture set |
| 5 |
|
# |
| 6 |
< |
set tf=/tmp/ph$$ |
| 6 |
> |
set tf=`mktemp /tmp/phdat.XXXXX` |
| 7 |
|
onintr quit |
| 8 |
|
if ( $#argv == 0 ) then |
| 9 |
|
pfilt -1 -x 128 -y 128 -p 1 \ |
| 10 |
< |
| pvalue -o -h -H -d -b > $tf.dat |
| 10 |
> |
| pvalue -o -h -H -df -b > $tf |
| 11 |
|
else |
| 12 |
< |
rm -f $tf.dat |
| 12 |
> |
rm -f $tf |
| 13 |
|
foreach i ( $* ) |
| 14 |
|
pfilt -1 -x 128 -y 128 -p 1 $i \ |
| 15 |
< |
| pvalue -o -h -H -d -b >> $tf.dat |
| 15 |
> |
| pvalue -o -h -H -df -b >> $tf |
| 16 |
|
if ( $status ) exit 1 |
| 17 |
|
end |
| 18 |
|
endif |
| 19 |
< |
set Lmin=`total -l $tf.dat | rcalc -e 'L=$1*179;$1=if(L-1e-7,log10(L)-.01,-7)'` |
| 20 |
< |
set Lmax=`total -u $tf.dat | rcalc -e '$1=log10($1*179)+.01'` |
| 21 |
< |
rcalc -e 'L=$1*179;cond=L-1e-7;$1=log10(L)' $tf.dat \ |
| 22 |
< |
| histo $Lmin $Lmax 100 |
| 19 |
> |
set Lmin=`total -if -l $tf | rcalc -e 'L=$1*179;$1=if(L-1e-7,log10(L)-.01,-7)'` |
| 20 |
> |
set Lmax=`total -if -u $tf | rcalc -e '$1=log10($1*179)+.01'` |
| 21 |
> |
rcalc -if -e 'L=$1*179;cond=L-1e-7;$1=log10(L)' $tf \ |
| 22 |
> |
| histo $Lmin $Lmax 777 |
| 23 |
|
quit: |
| 24 |
< |
rm -f $tf.dat |
| 24 |
> |
rm -f $tf |
