1 |
{ RCSid $Id$ } |
2 |
{ |
3 |
Fit Spherical Harmonic Coefficients for a Light Probe |
4 |
|
5 |
2/16/2005 G. Ward |
6 |
|
7 |
Constants: |
8 |
xres - horizontal input resolution |
9 |
yres - vertical input resolution |
10 |
map - light probe mapping (1=angular, 2=mirror) |
11 |
Usage: |
12 |
pvalue -h -H -pR -d probe.hdr \ |
13 |
| rcalc -e 'xres:512;yres:512;map:1' \ |
14 |
-f spharm.cal -f fitSH.cal \ |
15 |
| total > coef_R.dat |
16 |
Repeat for green and blue components if desired |
17 |
} |
18 |
map : 1; { just angular for now } |
19 |
|
20 |
mod(n,d) : n - floor(n/d)*d; |
21 |
sq(x) : x*x; |
22 |
|
23 |
x = mod(recno-1, xres); |
24 |
y = yres-1 - floor((recno-1)/xres); |
25 |
h = 1/xres*x - .5; |
26 |
v = 1/yres*y - .5; |
27 |
|
28 |
beta_ang = 2*PI*sqrt(h*h + v*v); |
29 |
gamma_ang = atan2(v, h); |
30 |
sine_beta_ang = sin(beta_ang); |
31 |
Dx_ang = cos(beta_ang); |
32 |
Dy_ang = cos(gamma_ang)*sine_beta_ang; |
33 |
Dz_ang = sin(gamma_ang)*sine_beta_ang; |
34 |
theta_ang = acos(Dz_ang); |
35 |
phi_ang = atan2(Dy_ang, Dx_ang); |
36 |
omega_ang = sq(4*PI/(xres+yres))*if(beta_ang-.01, sine_beta_ang/beta_ang, 1); |
37 |
|
38 |
cond = select(map, PI-1e-4-beta_ang); |
39 |
theta = select(map, theta_ang); |
40 |
phi = select(map, phi_ang); |
41 |
omega = select(map, omega_ang); |
42 |
|
43 |
{ Level 0-2 harmonics } |
44 |
$1 = $1*SphericalHarmonicB(1,theta,phi)*omega; |
45 |
$2 = $1*SphericalHarmonicB(2,theta,phi)*omega; |
46 |
$3 = $1*SphericalHarmonicB(3,theta,phi)*omega; |
47 |
$4 = $1*SphericalHarmonicB(4,theta,phi)*omega; |
48 |
$5 = $1*SphericalHarmonicB(5,theta,phi)*omega; |
49 |
$6 = $1*SphericalHarmonicB(6,theta,phi)*omega; |
50 |
$7 = $1*SphericalHarmonicB(7,theta,phi)*omega; |
51 |
$8 = $1*SphericalHarmonicB(8,theta,phi)*omega; |
52 |
$9 = $1*SphericalHarmonicB(9,theta,phi)*omega; |
53 |
{ Level 3 harmonics } |
54 |
$10 = $1*SphericalHarmonicB(10,theta,phi)*omega; |
55 |
$11 = $1*SphericalHarmonicB(11,theta,phi)*omega; |
56 |
$12 = $1*SphericalHarmonicB(12,theta,phi)*omega; |
57 |
$13 = $1*SphericalHarmonicB(13,theta,phi)*omega; |
58 |
$14 = $1*SphericalHarmonicB(14,theta,phi)*omega; |
59 |
$15 = $1*SphericalHarmonicB(15,theta,phi)*omega; |
60 |
$16 = $1*SphericalHarmonicB(16,theta,phi)*omega; |
61 |
{ Level 4 harmonics } |
62 |
$17 = $1*SphericalHarmonicB(17,theta,phi)*omega; |
63 |
$18 = $1*SphericalHarmonicB(18,theta,phi)*omega; |
64 |
$19 = $1*SphericalHarmonicB(19,theta,phi)*omega; |
65 |
$20 = $1*SphericalHarmonicB(20,theta,phi)*omega; |
66 |
$21 = $1*SphericalHarmonicB(21,theta,phi)*omega; |
67 |
$22 = $1*SphericalHarmonicB(22,theta,phi)*omega; |
68 |
$23 = $1*SphericalHarmonicB(23,theta,phi)*omega; |
69 |
$24 = $1*SphericalHarmonicB(24,theta,phi)*omega; |
70 |
$25 = $1*SphericalHarmonicB(25,theta,phi)*omega; |
71 |
{ Level 5 harmonics } |
72 |
$26 = $1*SphericalHarmonicB(26,theta,phi)*omega; |
73 |
$27 = $1*SphericalHarmonicB(27,theta,phi)*omega; |
74 |
$28 = $1*SphericalHarmonicB(28,theta,phi)*omega; |
75 |
$29 = $1*SphericalHarmonicB(29,theta,phi)*omega; |
76 |
$30 = $1*SphericalHarmonicB(30,theta,phi)*omega; |
77 |
$31 = $1*SphericalHarmonicB(31,theta,phi)*omega; |
78 |
$32 = $1*SphericalHarmonicB(32,theta,phi)*omega; |
79 |
$33 = $1*SphericalHarmonicB(33,theta,phi)*omega; |
80 |
$34 = $1*SphericalHarmonicB(34,theta,phi)*omega; |
81 |
$35 = $1*SphericalHarmonicB(35,theta,phi)*omega; |
82 |
$36 = $1*SphericalHarmonicB(36,theta,phi)*omega; |
83 |
{ Level 6 harmonics } |
84 |
$37 = $1*SphericalHarmonicB(37,theta,phi)*omega; |
85 |
$38 = $1*SphericalHarmonicB(38,theta,phi)*omega; |
86 |
$39 = $1*SphericalHarmonicB(39,theta,phi)*omega; |
87 |
$40 = $1*SphericalHarmonicB(40,theta,phi)*omega; |
88 |
$41 = $1*SphericalHarmonicB(41,theta,phi)*omega; |
89 |
$42 = $1*SphericalHarmonicB(42,theta,phi)*omega; |
90 |
$43 = $1*SphericalHarmonicB(43,theta,phi)*omega; |
91 |
$44 = $1*SphericalHarmonicB(44,theta,phi)*omega; |
92 |
$45 = $1*SphericalHarmonicB(45,theta,phi)*omega; |
93 |
$46 = $1*SphericalHarmonicB(46,theta,phi)*omega; |
94 |
$47 = $1*SphericalHarmonicB(47,theta,phi)*omega; |
95 |
$48 = $1*SphericalHarmonicB(48,theta,phi)*omega; |
96 |
$49 = $1*SphericalHarmonicB(49,theta,phi)*omega; |
97 |
{ Level 7 harmonics } |
98 |
$50 = $1*SphericalHarmonicB(50,theta,phi)*omega; |
99 |
$51 = $1*SphericalHarmonicB(51,theta,phi)*omega; |
100 |
$52 = $1*SphericalHarmonicB(52,theta,phi)*omega; |
101 |
$53 = $1*SphericalHarmonicB(53,theta,phi)*omega; |
102 |
$54 = $1*SphericalHarmonicB(54,theta,phi)*omega; |
103 |
$55 = $1*SphericalHarmonicB(55,theta,phi)*omega; |
104 |
$56 = $1*SphericalHarmonicB(56,theta,phi)*omega; |
105 |
$57 = $1*SphericalHarmonicB(57,theta,phi)*omega; |
106 |
$58 = $1*SphericalHarmonicB(58,theta,phi)*omega; |
107 |
$59 = $1*SphericalHarmonicB(59,theta,phi)*omega; |
108 |
$60 = $1*SphericalHarmonicB(60,theta,phi)*omega; |
109 |
$61 = $1*SphericalHarmonicB(61,theta,phi)*omega; |
110 |
$62 = $1*SphericalHarmonicB(62,theta,phi)*omega; |
111 |
$63 = $1*SphericalHarmonicB(63,theta,phi)*omega; |
112 |
$64 = $1*SphericalHarmonicB(64,theta,phi)*omega; |
113 |
{ Level 8 harmonics } |
114 |
$65 = $1*SphericalHarmonicB(65,theta,phi)*omega; |
115 |
$66 = $1*SphericalHarmonicB(66,theta,phi)*omega; |
116 |
$67 = $1*SphericalHarmonicB(67,theta,phi)*omega; |
117 |
$68 = $1*SphericalHarmonicB(68,theta,phi)*omega; |
118 |
$69 = $1*SphericalHarmonicB(69,theta,phi)*omega; |
119 |
$70 = $1*SphericalHarmonicB(70,theta,phi)*omega; |
120 |
$71 = $1*SphericalHarmonicB(71,theta,phi)*omega; |
121 |
$72 = $1*SphericalHarmonicB(72,theta,phi)*omega; |
122 |
$73 = $1*SphericalHarmonicB(73,theta,phi)*omega; |
123 |
$74 = $1*SphericalHarmonicB(74,theta,phi)*omega; |
124 |
$75 = $1*SphericalHarmonicB(75,theta,phi)*omega; |
125 |
$76 = $1*SphericalHarmonicB(76,theta,phi)*omega; |
126 |
$77 = $1*SphericalHarmonicB(77,theta,phi)*omega; |
127 |
$78 = $1*SphericalHarmonicB(78,theta,phi)*omega; |
128 |
$79 = $1*SphericalHarmonicB(79,theta,phi)*omega; |
129 |
$80 = $1*SphericalHarmonicB(80,theta,phi)*omega; |
130 |
$81 = $1*SphericalHarmonicB(81,theta,phi)*omega; |