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