src/px/phisto.csh

#!/bin/csh -f
# RCSid: $Id: phisto.csh,v 3.4 2005/02/16 05:40:11 greg Exp $
#
# Compute foveal histogram for picture set
#
set tf=`mktemp /tmp/phdat.XXXXX`
onintr quit
if ( $#argv == 0 ) then
        pfilt -1 -x 128 -y 128 -p 1 \
                        | pvalue -o -h -H -d -b > $tf
else
        rm -f $tf
        foreach i ( $* )
                pfilt -1 -x 128 -y 128 -p 1 $i \
                                | pvalue -o -h -H -d -b >> $tf
                if ( $status ) exit 1
        end
endif
set Lmin=`total -l $tf | rcalc -e 'L=$1*179;$1=if(L-1e-7,log10(L)-.01,-7)'`
set Lmax=`total -u $tf | rcalc -e '$1=log10($1*179)+.01'`
rcalc -e 'L=$1*179;cond=L-1e-7;$1=log10(L)' $tf \
        | histo $Lmin $Lmax 100
quit:
rm -f $tf
Revision:	3.5
Committed:	Mon Aug 25 04:50:32 2008 UTC (15 years, 8 months ago) by greg
Content type:	application/x-csh
Branch:	MAIN
CVS Tags:	rad4R2P2, rad4R2, rad4R1, rad4R0, rad4R2P1
Changes since 3.4:	+9 -9 lines
Log Message:	Made temporary files safe from malicious attackers
#	User	Rev	Content
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