| 1 |
#!/bin/csh -f
|
| 2 |
# RCSid: $Id: phisto.csh,v 3.6 2015/06/29 21:03:45 greg Exp $
|
| 3 |
#
|
| 4 |
# Compute foveal histogram for picture set
|
| 5 |
#
|
| 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 -df -b > $tf
|
| 11 |
else
|
| 12 |
rm -f $tf
|
| 13 |
foreach i ( $* )
|
| 14 |
pfilt -1 -x 128 -y 128 -p 1 $i \
|
| 15 |
| pvalue -o -h -H -df -b >> $tf
|
| 16 |
if ( $status ) exit 1
|
| 17 |
end
|
| 18 |
endif
|
| 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
|
