1 #!/bin/csh -f

     2 # RCSid: $Id: phisto.csh,v 3.5 2008/08/25 04:50:32 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 -d -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 -d -b >> $tf

    16 		if ( $status ) exit 1

    17 	end

    18 endif

    19 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 	| histo $Lmin $Lmax 100

    23 quit:

    24 rm -f $tf