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
#	User	Rev	Content
1	gwlarson	3.1	#!/bin/csh -f
2	greg	3.3	# RCSid: $Id$
3	gwlarson	3.1	#
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	gwlarson	3.2	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	gwlarson	3.1	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