src/px/phisto.csh

#!/bin/csh -f
# RCSid: $Id: phisto.csh,v 3.6 2015/06/29 21:03:45 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 -df -b > $tf
else
        rm -f $tf
        foreach i ( $* )
                pfilt -1 -x 128 -y 128 -p 1 $i \
                                | pvalue -o -h -H -df -b >> $tf
                if ( $status ) exit 1
        end
endif
set Lmin=`total -if -l $tf | rcalc -e 'L=$1*179;$1=if(L-1e-7,log10(L)-.01,-7)'`
set Lmax=`total -if -u $tf | rcalc -e '$1=log10($1*179)+.01'`
rcalc -if -e 'L=$1*179;cond=L-1e-7;$1=log10(L)' $tf \
        | histo $Lmin $Lmax 777
quit:
rm -f $tf
Revision:	3.7
Committed:	Fri Aug 21 05:48:28 2015 UTC (8 years, 8 months ago) by greg
Content type:	application/x-csh
Branch:	MAIN
CVS Tags:	rad5R2, rad5R0, rad5R1, rad5R3
Changes since 3.6:	+2 -2 lines
Log Message:	Fixed issue with undersampled histogram in phisto
#	User	Rev	Content
1	gwlarson	3.1	#!/bin/csh -f
2	greg	3.7	# RCSid: $Id: phisto.csh,v 3.6 2015/06/29 21:03:45 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.6	\| pvalue -o -h -H -df -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.6	\| pvalue -o -h -H -df -b >> $tf
16	gwlarson	3.1	if ( $status ) exit 1
17			end
18			endif
19	greg	3.6	set Lmin=`total -if -l $tf \| rcalc -e 'L=$1*179;$1=if(L-1e-7,log10(L)-.01,-7)'`
20			set Lmax=`total -if -u $tf \| rcalc -e '$1=log10($1*179)+.01'`
21			rcalc -if -e 'L=$1*179;cond=L-1e-7;$1=log10(L)' $tf \
22	greg	3.7	\| histo $Lmin $Lmax 777
23	gwlarson	3.1	quit:
24	greg	3.5	rm -f $tf