--- ray/src/px/phisto.csh 1998/03/12 15:46:53 3.1 +++ ray/src/px/phisto.csh 2003/02/22 02:07:27 3.3 @@ -1,5 +1,5 @@ #!/bin/csh -f -# SCCSid "$SunId$ SGI" +# RCSid: $Id: phisto.csh,v 3.3 2003/02/22 02:07:27 greg Exp $ # # Compute foveal histogram for picture set # @@ -16,8 +16,8 @@ else if ( $status ) exit 1 end endif -set Lmin=`total -l $tf.dat | rcalc -e 'L=$1*179;$1=if(L-1e-7,log10(L),-7)'` -set Lmax=`total -u $tf.dat | rcalc -e '$1=log10($1*179)'` +set Lmin=`total -l $tf.dat | rcalc -e 'L=$1*179;$1=if(L-1e-7,log10(L)-.01,-7)'` +set Lmax=`total -u $tf.dat | rcalc -e '$1=log10($1*179)+.01'` rcalc -e 'L=$1*179;cond=L-1e-7;$1=log10(L)' $tf.dat \ | histo $Lmin $Lmax 100 quit: