src/px/falsecolor.csh

#!/bin/csh -fe
# RCSid: $Id: falsecolor.csh,v 2.15 2005/11/18 22:50:21 greg Exp $
#
# Create false color image with legend
#
# Added user-definable legend 2004/01/20 Rob Guglielmetti

set td=/tmp/fc$$
onintr quit
set mult=179
set label=Nits
set scale=1000
set decades=0
set redv='def_red(v)'
set grnv='def_grn(v)'
set bluv='def_blu(v)'
set ndivs=8
set picture='-'
set cpict=
set loff=0
set legwidth=100
set legheight=200
while ($#argv > 0)
        switch ($argv[1])
        case -lw:
                shift argv
                set legwidth="$argv[1]"
                breaksw
        case -lh:
                shift argv
                set legheight="$argv[1]"
                breaksw
        case -m:
                shift argv
                set mult="$argv[1]"
                breaksw
        case -s:
                shift argv
                set scale="$argv[1]"
                if ("$scale" =~ [aA]*) set needfile
                breaksw
        case -l:
                shift argv
                set label="$argv[1]"
                breaksw
        case -log:
                shift argv
                set decades=$argv[1]
                breaksw
        case -r:
                shift argv
                set redv="$argv[1]"
                breaksw
        case -g:
                shift argv
                set grnv="$argv[1]"
                breaksw
        case -b:
                shift argv
                set bluv="$argv[1]"
                breaksw
        case -spec:
                shift argv
                set redv='1.6*v-.6'
                set grnv='if(v-.375,1.6-1.6*v,8/3*v)'
                set bluv='1-8/3*v'
                breaksw
        case -i:
                shift argv
                set picture="$argv[1]"
                breaksw
        case -p:
                shift argv
                set cpict="$argv[1]"
                breaksw
        case -ip:
        case -pi:
                shift argv
                set picture="$argv[1]"
                set cpict="$argv[1]"
                breaksw
        case -cl:
                set docont=a
                set loff=12
                breaksw
        case -cb:
                set docont=b
                set loff=13
                breaksw
        case -e:
                set doextrem
                set needfile
                breaksw
        case -n:
                shift argv
                set ndivs="$argv[1]"
                breaksw
        default:
                echo bad option "'$argv[1]'"
                exit 1
        endsw
        shift argv
end
mkdir $td
if ($?needfile && "$picture" == '-') then
        cat > $td/picture
        set picture=$td/picture
endif
if ("$scale" =~ [aA]*) then
        set LogLmax=`phisto $picture | tail -2 | sed -n '1s/    [0-9]*$//p'`
        set scale=`ev "$mult/179*10^$LogLmax"`
endif
cat > $td/pc0.cal <<_EOF_
PI : 3.14159265358979323846 ;
scale : $scale ;
mult : $mult ;
ndivs : $ndivs ;

or(a,b) : if(a,a,b);
EPS : 1e-7;
neq(a,b) : if(a-b-EPS,1,b-a-EPS);
btwn(a,x,b) : if(a-x,-1,b-x);
clip(x) : if(x-1,1,if(x,x,0));
frac(x) : x - floor(x);
boundary(a,b) : neq(floor(ndivs*a+.5),floor(ndivs*b+.5));

interp_arr2(i,x,f):(i+1-x)*f(i)+(x-i)*f(i+1);
interp_arr(x,f):if(x-1,if(f(0)-x,interp_arr2(floor(x),x,f),f(f(0))),f(1));
def_redp(i):select(i,0.18848,0.05468174,
0.00103547,8.311144e-08,7.449763e-06,0.0004390987,0.001367254,
0.003076,0.01376382,0.06170773,0.1739422,0.2881156,0.3299725,
0.3552663,0.372552,0.3921184,0.4363976,0.6102754,0.7757267,
0.9087369,1,1,0.9863);
def_red(x):interp_arr(x/0.0454545+1,def_redp);
def_grnp(i):select(i,0.0009766,2.35501e-05,
0.0008966244,0.0264977,0.1256843,0.2865799,0.4247083,0.4739468,
0.4402732,0.3671876,0.2629843,0.1725325,0.1206819,0.07316644,
0.03761026,0.01612362,0.004773749,6.830967e-06,0.00803605,
0.1008085,0.3106831,0.6447838,0.9707);
def_grn(x):interp_arr(x/0.0454545+1,def_grnp);
def_blup(i):select(i,0.2666,0.3638662,0.4770437,
0.5131397,0.5363797,0.5193677,0.4085123,0.1702815,0.05314236,
0.05194055,0.08564082,0.09881395,0.08324373,0.06072902,
0.0391076,0.02315354,0.01284458,0.005184709,0.001691774,
2.432735e-05,1.212949e-05,0.006659406,0.02539);
def_blu(x):interp_arr(x/0.0454545+1,def_blup);

isconta = if(btwn(0,v,1),or(boundary(vleft,vright),boundary(vabove,vbelow)),-1);
iscontb = if(btwn(0,v,1),btwn(.4,frac(ndivs*v),.6),-1); 

ra = 0;
ga = 0;
ba = 0;

in = 1;

ro = if(in,clip($redv),ra);
go = if(in,clip($grnv),ga);
bo = if(in,clip($bluv),ba);
_EOF_
cat > $td/pc1.cal <<_EOF_
norm : mult/scale/le(1);

v = map(li(1)*norm);

vleft = map(li(1,-1,0)*norm);
vright = map(li(1,1,0)*norm);
vabove = map(li(1,0,1)*norm);
vbelow = map(li(1,0,-1)*norm);

map(x) = x;

ra = ri(nfiles);
ga = gi(nfiles);
ba = bi(nfiles);
_EOF_
set pc0args=(-f $td/pc0.cal)
set pc1args=(-f $td/pc1.cal)
if ($?docont) then
        set pc0args=($pc0args -e "in=iscont$docont")
endif
if ("$cpict" == "") then
        set pc1args=($pc1args -e 'ra=0;ga=0;ba=0')
else if ("$cpict" == "$picture") then
        set cpict=
endif
if ("$decades" != "0") then
        set pc1args=($pc1args -e "map(x)=if(x-10^-$decades,log10(x)/$decades+1,0)")
        set imap="imap(y)=10^((y-1)*$decades)"
else
        set imap="imap(y)=y"
endif
if ( $legwidth > 20 && $legheight > 40 ) then
pcomb $pc0args -e 'v=(y+.5)/yres;vleft=v;vright=v' \
                -e 'vbelow=(y-.5)/yres;vabove=(y+1.5)/yres' \
                -x $legwidth -y $legheight > $td/scol.pic
( echo "$label"; cnt $ndivs \
                | rcalc -e '$1='"($scale)*imap(($ndivs-.5-"'$1'")/$ndivs)" \
                -e "$imap" | sed -e 's/\(\.[0-9][0-9][0-9]\)[0-9]*/\1/' ) \
        | psign -s -.15 -cf 1 1 1 -cb 0 0 0 \
                -h `ev "floor($legheight/$ndivs+.5)"` > $td/slab.pic
else
        set legwidth=0
        set legheight=0
        (echo "" ; echo "-Y 1 +X 1" ; echo "aaa" ) > $td/scol.pic
        cp $td/scol.pic $td/slab.pic
endif
if ( $?doextrem ) then
        pextrem -o $picture > $td/extrema
        set minpos=`sed 2d $td/extrema | rcalc -e '$2=$2;$1=$1+'"$legwidth"`
        set minval=`rcalc -e '$1=($3*.27+$4*.67+$5*.06)*'"$mult" $td/extrema | sed -e 2d -e 's/\(\.[0-9][0-9][0-9]\)[0-9]*/\1/'`
        set maxpos=`sed 1d $td/extrema | rcalc -e '$2=$2;$1=$1+'"$legwidth"`
        set maxval=`rcalc -e '$1=($3*.27+$4*.67+$5*.06)*'"$mult" $td/extrema | sed -e 1d -e 's/\(\.[0-9][0-9][0-9]\)[0-9]*/\1/'`
        psign -s -.15 -a 2 -h 16 $minval > $td/minv.pic
        psign -s -.15 -a 2 -h 16 $maxval > $td/maxv.pic
        pcomb $pc0args $pc1args $picture $cpict \
                | pcompos $td/scol.pic 0 0 \
                        +t .1 "\!pcomb -e 'lo=1-gi(1)' $td/slab.pic" \
                        `ev 2 $loff-1` -t .5 $td/slab.pic 0 $loff \
                  - $legwidth 0 $td/minv.pic $minpos $td/maxv.pic $maxpos
else
        pcomb $pc0args $pc1args $picture $cpict \
                | pcompos $td/scol.pic 0 0 \
                        +t .1 "\!pcomb -e 'lo=1-gi(1)' $td/slab.pic" \
                        `ev 2 $loff-1` -t .5 $td/slab.pic 0 $loff \
                        - $legwidth 0
endif
quit:
rm -rf $td
Revision:	2.16
Committed:	Mon Nov 13 02:41:46 2006 UTC (17 years, 6 months ago) by greg
Content type:	application/x-csh
Branch:	MAIN
Changes since 2.15:	+16 -8 lines
Log Message:	Added -spec and "-s auto" options to falsecolor
#	User	Rev	Content
1	greg	2.8	#!/bin/csh -fe
2	greg	2.16	# RCSid: $Id: falsecolor.csh,v 2.15 2005/11/18 22:50:21 greg Exp $
3	greg	1.1	#
4			# Create false color image with legend
5			#
6	greg	2.9	# Added user-definable legend 2004/01/20 Rob Guglielmetti
7
8	greg	2.12	set td=/tmp/fc$$
9	greg	1.1	onintr quit
10	greg	1.9	set mult=179
11	greg	1.1	set label=Nits
12			set scale=1000
13	greg	1.3	set decades=0
14	greg	2.13	set redv='def_red(v)'
15			set grnv='def_grn(v)'
16			set bluv='def_blu(v)'
17	greg	1.1	set ndivs=8
18			set picture='-'
19			set cpict=
20	greg	2.4	set loff=0
21	greg	1.5	set legwidth=100
22			set legheight=200
23	greg	1.1	while ($#argv > 0)
24			switch ($argv[1])
25	greg	2.9	case -lw:
26			shift argv
27			set legwidth="$argv[1]"
28			breaksw
29			case -lh:
30			shift argv
31			set legheight="$argv[1]"
32			breaksw
33	greg	1.1	case -m:
34			shift argv
35			set mult="$argv[1]"
36			breaksw
37			case -s:
38			shift argv
39			set scale="$argv[1]"
40	greg	2.16	if ("$scale" =~ [aA]*) set needfile
41	greg	1.1	breaksw
42			case -l:
43			shift argv
44			set label="$argv[1]"
45			breaksw
46	greg	1.3	case -log:
47			shift argv
48			set decades=$argv[1]
49			breaksw
50	greg	1.1	case -r:
51			shift argv
52			set redv="$argv[1]"
53			breaksw
54			case -g:
55			shift argv
56			set grnv="$argv[1]"
57			breaksw
58			case -b:
59			shift argv
60			set bluv="$argv[1]"
61			breaksw
62	greg	2.16	case -spec:
63			shift argv
64			set redv='1.6*v-.6'
65			set grnv='if(v-.375,1.6-1.6v,8/3v)'
66			set bluv='1-8/3*v'
67			breaksw
68	greg	1.1	case -i:
69			shift argv
70			set picture="$argv[1]"
71			breaksw
72			case -p:
73			shift argv
74			set cpict="$argv[1]"
75			breaksw
76			case -ip:
77			case -pi:
78	greg	1.2	shift argv
79	greg	1.1	set picture="$argv[1]"
80			set cpict="$argv[1]"
81			breaksw
82			case -cl:
83			set docont=a
84	greg	2.5	set loff=12
85	greg	1.1	breaksw
86			case -cb:
87			set docont=b
88	greg	2.5	set loff=13
89	greg	1.1	breaksw
90	greg	2.3	case -e:
91			set doextrem
92			set needfile
93			breaksw
94	greg	1.1	case -n:
95	greg	1.2	shift argv
96	greg	1.1	set ndivs="$argv[1]"
97	greg	1.2	breaksw
98	greg	1.1	default:
99	greg	2.11	echo bad option "'$argv[1]'"
100	greg	1.1	exit 1
101			endsw
102			shift argv
103			end
104	greg	1.5	mkdir $td
105	greg	2.16	if ($?needfile && "$picture" == '-') then
106			cat > $td/picture
107			set picture=$td/picture
108			endif
109			if ("$scale" =~ [aA]*) then
110			set LogLmax=`phisto $picture \| tail -2 \| sed -n '1s/ [0-9]*$//p'`
111			set scale=`ev "$mult/179*10^$LogLmax"`
112			endif
113	greg	1.5	cat > $td/pc0.cal <<_EOF_
114	greg	1.8	PI : 3.14159265358979323846 ;
115	greg	1.1	scale : $scale ;
116			mult : $mult ;
117			ndivs : $ndivs ;
118
119			or(a,b) : if(a,a,b);
120			EPS : 1e-7;
121			neq(a,b) : if(a-b-EPS,1,b-a-EPS);
122	greg	1.7	btwn(a,x,b) : if(a-x,-1,b-x);
123	greg	1.5	clip(x) : if(x-1,1,if(x,x,0));
124	greg	1.1	frac(x) : x - floor(x);
125	greg	2.2	boundary(a,b) : neq(floor(ndivsa+.5),floor(ndivsb+.5));
126	greg	1.1
127	greg	2.13	interp_arr2(i,x,f):(i+1-x)f(i)+(x-i)f(i+1);
128			interp_arr(x,f):if(x-1,if(f(0)-x,interp_arr2(floor(x),x,f),f(f(0))),f(1));
129			def_redp(i):select(i,0.18848,0.05468174,
130	greg	2.15	0.00103547,8.311144e-08,7.449763e-06,0.0004390987,0.001367254,
131			0.003076,0.01376382,0.06170773,0.1739422,0.2881156,0.3299725,
132			0.3552663,0.372552,0.3921184,0.4363976,0.6102754,0.7757267,
133			0.9087369,1,1,0.9863);
134	greg	2.13	def_red(x):interp_arr(x/0.0454545+1,def_redp);
135	greg	2.14	def_grnp(i):select(i,0.0009766,2.35501e-05,
136	greg	2.15	0.0008966244,0.0264977,0.1256843,0.2865799,0.4247083,0.4739468,
137			0.4402732,0.3671876,0.2629843,0.1725325,0.1206819,0.07316644,
138			0.03761026,0.01612362,0.004773749,6.830967e-06,0.00803605,
139			0.1008085,0.3106831,0.6447838,0.9707);
140	greg	2.13	def_grn(x):interp_arr(x/0.0454545+1,def_grnp);
141			def_blup(i):select(i,0.2666,0.3638662,0.4770437,
142	greg	2.15	0.5131397,0.5363797,0.5193677,0.4085123,0.1702815,0.05314236,
143			0.05194055,0.08564082,0.09881395,0.08324373,0.06072902,
144			0.0391076,0.02315354,0.01284458,0.005184709,0.001691774,
145			2.432735e-05,1.212949e-05,0.006659406,0.02539);
146	greg	2.13	def_blu(x):interp_arr(x/0.0454545+1,def_blup);
147
148	greg	1.6	isconta = if(btwn(0,v,1),or(boundary(vleft,vright),boundary(vabove,vbelow)),-1);
149	greg	1.5	iscontb = if(btwn(0,v,1),btwn(.4,frac(ndivs*v),.6),-1);
150	greg	1.1
151	greg	1.5	ra = 0;
152			ga = 0;
153			ba = 0;
154
155			in = 1;
156	greg	2.13
157			ro = if(in,clip($redv),ra);
158			go = if(in,clip($grnv),ga);
159			bo = if(in,clip($bluv),ba);
160	greg	1.5	_EOF_
161			cat > $td/pc1.cal <<_EOF_
162			norm : mult/scale/le(1);
163
164	greg	1.4	v = map(li(1)*norm);
165	greg	1.5
166	greg	1.4	vleft = map(li(1,-1,0)*norm);
167			vright = map(li(1,1,0)*norm);
168			vabove = map(li(1,0,1)*norm);
169			vbelow = map(li(1,0,-1)*norm);
170	greg	1.1
171	greg	1.5	map(x) = x;
172	greg	1.1
173			ra = ri(nfiles);
174			ga = gi(nfiles);
175			ba = bi(nfiles);
176			_EOF_
177	greg	1.5	set pc0args=(-f $td/pc0.cal)
178			set pc1args=(-f $td/pc1.cal)
179	greg	1.1	if ($?docont) then
180	greg	1.5	set pc0args=($pc0args -e "in=iscont$docont")
181	greg	1.1	endif
182			if ("$cpict" == "") then
183	greg	1.5	set pc1args=($pc1args -e 'ra=0;ga=0;ba=0')
184	greg	1.1	else if ("$cpict" == "$picture") then
185			set cpict=
186			endif
187	greg	2.6	if ("$decades" != "0") then
188	greg	1.5	set pc1args=($pc1args -e "map(x)=if(x-10^-$decades,log10(x)/$decades+1,0)")
189	greg	1.3	set imap="imap(y)=10^((y-1)*$decades)"
190			else
191			set imap="imap(y)=y"
192			endif
193	greg	2.9	if ( $legwidth > 20 && $legheight > 40 ) then
194	greg	1.5	pcomb $pc0args -e 'v=(y+.5)/yres;vleft=v;vright=v' \
195			-e 'vbelow=(y-.5)/yres;vabove=(y+1.5)/yres' \
196			-x $legwidth -y $legheight > $td/scol.pic
197	greg	2.10	( echo "$label"; cnt $ndivs \
198	greg	1.3	\| rcalc -e '$1='"($scale)*imap(($ndivs-.5-"'$1'")/$ndivs)" \
199	greg	2.7	-e "$imap" \| sed -e 's/\(\.[0-9][0-9][0-9]\)[0-9]*/\1/' ) \
200	greg	2.3	\| psign -s -.15 -cf 1 1 1 -cb 0 0 0 \
201			-h `ev "floor($legheight/$ndivs+.5)"` > $td/slab.pic
202	greg	2.9	else
203			set legwidth=0
204			set legheight=0
205			(echo "" ; echo "-Y 1 +X 1" ; echo "aaa" ) > $td/scol.pic
206			cp $td/scol.pic $td/slab.pic
207			endif
208	greg	2.3	if ( $?doextrem ) then
209			pextrem -o $picture > $td/extrema
210			set minpos=`sed 2d $td/extrema \| rcalc -e '$2=$2;$1=$1+'"$legwidth"`
211	greg	2.7	set minval=`rcalc -e '$1=($3.27+$4.67+$5.06)'"$mult" $td/extrema \| sed -e 2d -e 's/\(\.[0-9][0-9][0-9]\)[0-9]*/\1/'`
212	greg	2.3	set maxpos=`sed 1d $td/extrema \| rcalc -e '$2=$2;$1=$1+'"$legwidth"`
213	greg	2.7	set maxval=`rcalc -e '$1=($3.27+$4.67+$5.06)'"$mult" $td/extrema \| sed -e 1d -e 's/\(\.[0-9][0-9][0-9]\)[0-9]*/\1/'`
214	greg	2.3	psign -s -.15 -a 2 -h 16 $minval > $td/minv.pic
215			psign -s -.15 -a 2 -h 16 $maxval > $td/maxv.pic
216			pcomb $pc0args $pc1args $picture $cpict \
217	greg	2.13	\| pcompos $td/scol.pic 0 0 \
218			+t .1 "\!pcomb -e 'lo=1-gi(1)' $td/slab.pic" \
219			`ev 2 $loff-1` -t .5 $td/slab.pic 0 $loff \
220	greg	2.3	- $legwidth 0 $td/minv.pic $minpos $td/maxv.pic $maxpos
221			else
222			pcomb $pc0args $pc1args $picture $cpict \
223	greg	2.13	\| pcompos $td/scol.pic 0 0 \
224			+t .1 "\!pcomb -e 'lo=1-gi(1)' $td/slab.pic" \
225			`ev 2 $loff-1` -t .5 $td/slab.pic 0 $loff \
226			- $legwidth 0
227	greg	2.3	endif
228	greg	1.1	quit:
229	greg	1.5	rm -rf $td