src/px/phisto.csh

#!/bin/csh -f
# RCSid: $Id$
#
# Compute foveal histogram for picture set
#
set tf=/usr/tmp/ph$$
onintr quit
if ( $#argv == 0 ) then
        pfilt -1 -x 128 -y 128 -p 1 \
                        | pvalue -o -h -H -d -b > $tf.dat
else
        rm -f $tf.dat
        foreach i ( $* )
                pfilt -1 -x 128 -y 128 -p 1 $i \
                                | pvalue -o -h -H -d -b >> $tf.dat
                if ( $status ) exit 1
        end
endif
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:
rm -f $tf.dat
Revision:	3.3
Committed:	Sat Feb 22 02:07:27 2003 UTC (21 years, 2 months ago) by greg
Content type:	application/x-csh
Branch:	MAIN
CVS Tags:	rad3R5, rad3R6, rad3R6P1
Changes since 3.2:	+1 -1 lines
Log Message:	Changes and check-in for 3.5 release Includes new source files and modifications not recorded for many years See ray/doc/notes/ReleaseNotes for notes between 3.1 and 3.5 release
#	Content
1	#!/bin/csh -f
2	# RCSid: $Id$
3	#
4	# Compute foveal histogram for picture set
5	#
6	set tf=/usr/tmp/ph$$
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
11	else
12	rm -f $tf.dat
13	foreach i ( $* )
14	pfilt -1 -x 128 -y 128 -p 1 $i \
15	\| pvalue -o -h -H -d -b >> $tf.dat
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
23	quit:
24	rm -f $tf.dat