.MCAD 304020000 1 88 300 0
.CMD PLOTFORMAT
0 0 1 0 1 0 0 0 1 
0 0 1 0 1 0 0 0 1 
0 1 0 0 1 1 NO-TRACE-STRING
0 2 1 0 1 1 NO-TRACE-STRING
0 3 2 0 1 1 NO-TRACE-STRING
0 4 3 0 1 1 NO-TRACE-STRING
0 1 4 0 1 1 NO-TRACE-STRING
0 2 5 0 1 1 NO-TRACE-STRING
0 3 6 0 1 1 NO-TRACE-STRING
0 4 0 0 1 1 NO-TRACE-STRING
0 1 1 0 1 1 NO-TRACE-STRING
0 2 2 0 1 1 NO-TRACE-STRING
0 3 3 0 1 1 NO-TRACE-STRING
0 4 4 0 1 1 NO-TRACE-STRING
0 1 5 0 1 1 NO-TRACE-STRING
0 2 6 0 1 1 NO-TRACE-STRING
0 3 0 0 1 1 NO-TRACE-STRING
0 4 1 0 1 1 NO-TRACE-STRING
0 1 1 21 15 0 0 3 
.CMD FORMAT  rd=d ct=10 im=i et=6 zt=8 pr=3 mass length time charge temperature tr=0 vm=1
.CMD SET ORIGIN 0
.CMD SET TOL 0.001000000000000
.CMD SET PRNCOLWIDTH 6
.CMD SET PRNPRECISION 4
.CMD PRINT_SETUP 0.700000 0.500000 0.000000 0.000000 0
.CMD HEADER_FOOTER 1 1 *empty* *empty* *empty* 0 1 *empty* *empty* *empty*
.CMD HEADER_FOOTER_FONT fontID=14 family=Arial points=10 bold=0 italic=0 underline=0 colrid=13
.CMD HEADER_FOOTER_FONT fontID=15 family=Arial points=10 bold=0 italic=0 underline=0 colrid=13
.CMD DEFAULT_TEXT_PARPROPS 0 0 0
.CMD DEFINE_FONTSTYLE_NAME fontID=0 name=Variables
.CMD DEFINE_FONTSTYLE_NAME fontID=1 name=Constants
.CMD DEFINE_FONTSTYLE_NAME fontID=2 name=Text
.CMD DEFINE_FONTSTYLE_NAME fontID=4 name=User^1
.CMD DEFINE_FONTSTYLE_NAME fontID=5 name=User^2
.CMD DEFINE_FONTSTYLE_NAME fontID=6 name=User^3
.CMD DEFINE_FONTSTYLE_NAME fontID=7 name=User^4
.CMD DEFINE_FONTSTYLE_NAME fontID=8 name=User^5
.CMD DEFINE_FONTSTYLE_NAME fontID=9 name=User^6
.CMD DEFINE_FONTSTYLE_NAME fontID=10 name=User^7
.CMD DEFINE_FONTSTYLE fontID=0 family=Times^New^Roman points=10 bold=0 italic=0 underline=0 colrid=-1
.CMD DEFINE_FONTSTYLE fontID=1 family=Times^New^Roman points=10 bold=0 italic=0 underline=0 colrid=-1
.CMD DEFINE_FONTSTYLE fontID=2 family=Arial points=10 bold=0 italic=0 underline=0 colrid=-1
.CMD DEFINE_FONTSTYLE fontID=4 family=Arial points=10 bold=0 italic=0 underline=0 colrid=-1
.CMD DEFINE_FONTSTYLE fontID=5 family=Courier^New points=10 bold=0 italic=0 underline=0 colrid=-1
.CMD DEFINE_FONTSTYLE fontID=6 family=System points=10 bold=0 italic=0 underline=0 colrid=-1
.CMD DEFINE_FONTSTYLE fontID=7 family=Script points=10 bold=0 italic=0 underline=0 colrid=-1
.CMD DEFINE_FONTSTYLE fontID=8 family=Roman points=10 bold=0 italic=0 underline=0 colrid=-1
.CMD DEFINE_FONTSTYLE fontID=9 family=Modern points=10 bold=0 italic=0 underline=0 colrid=-1
.CMD DEFINE_FONTSTYLE fontID=10 family=Times^New^Roman points=10 bold=0 italic=0 underline=0 colrid=-1
.CMD UNITS U=1
.CMD DIMENSIONS_ANALYSIS 0 0
.CMD COLORTAB_ENTRY 0 0 0
.CMD COLORTAB_ENTRY 128 0 0
.CMD COLORTAB_ENTRY 0 128 0
.CMD COLORTAB_ENTRY 128 128 0
.CMD COLORTAB_ENTRY 0 0 128
.CMD COLORTAB_ENTRY 128 0 128
.CMD COLORTAB_ENTRY 0 128 128
.CMD COLORTAB_ENTRY 128 128 128
.CMD COLORTAB_ENTRY 192 192 192
.CMD COLORTAB_ENTRY 255 0 0
.CMD COLORTAB_ENTRY 0 255 0
.CMD COLORTAB_ENTRY 255 255 0
.CMD COLORTAB_ENTRY 0 0 255
.CMD COLORTAB_ENTRY 255 0 255
.CMD COLORTAB_ENTRY 0 255 255
.CMD COLORTAB_ENTRY 255 255 255
.TXT 2 17 255 0 0
Cg a51.625000,51.625000,48
{\rtf\ansi \deff0{\colortbl;\red0\green0\blue128;}{\fonttbl{\f0
\fcharset0\fnil Arial;}}\plain\cf1\fs20 \pard {\fs24\b Frequency 
Domain Steerable Pyramid Filter Design}}
.TXT 3 -16 257 0 0
Cg a82.500000,82.500000,30
{\rtf\ansi \deff0{\colortbl;\red0\green0\blue128;}{\fonttbl{\f0
\fcharset0\fnil Arial;}}\plain\cf1\fs20 \pard First derivative, k = 2 filter}
.TXT 0 33 5 0 0
Cg a48.625000,48.625000,21
{\rtf\ansi \deff0{\colortbl;\red0\green0\blue128;}{\fonttbl{\f0
\fcharset0\fnil Arial;}}\plain\cf1\fs20 \pard {\b Ken Castleman} 5/28/98}
.TXT 0 40 256 0 0
Cg a18.625000,18.625000,12
{\rtf\ansi \deff0{\colortbl;\red0\green0\blue128;}{\fonttbl{\f0
\fcharset0\fnil Arial;}}\plain\cf1\fs20 \pard StPyr_02.MCD}
.EQN 5 -70 6 0 0
{0:N}NAME:16
.EQN 0 8 220 0 0
{0:a}NAME:{0:floor}NAME(({0:N}NAME)/(2))
.EQN 0 15 7 0 0
{0:i}NAME:0;{0:N}NAME
.EQN 0 14 8 0 0
{0:k}NAME:0;{0:N}NAME
.EQN 0 38 96 0 0
{0:j}NAME:\(-1)
.TXT 5 -78 162 0 0
Cg a85.625000,85.625000,67
{\rtf\ansi \deff0{\colortbl;\red0\green0\blue128;}{\fonttbl{\f0
\fcharset0\fnil Arial;}}\plain\cf1\fs20 \pard Define smooth-edged 
lowpass and highpass functions (raised cosine).}
.EQN 3 54 218 0 0
({0:\q}NAME)[({0:i}NAME,{0:k}NAME):{0:angle}NAME({0:i}NAME+0.001-{0:a}NAME,{0:k}NAME-{0:a}NAME)
.EQN 2 -54 22 0 0
{0:LP}NAME({0:x1}NAME,{0:x2}NAME,{0:x}NAME):{0:if}NAME({0:x}NAME<{0:x1}NAME,1,{1:if}NAME({0:x}NAME>{0:x2}NAME,0,\({0}0.5*(1+{0:cos}NAME({0:\p}NAME*(({0:x}NAME-{0:x1}NAME)/({0:x2}NAME-{0:x1}NAME)))))))
.EQN 3 54 219 0 0
({0:\r}NAME)[({0:i}NAME,{0:k}NAME):\((({0:i}NAME-{0:a}NAME))^(2)+(({0:k}NAME-{0:a}NAME))^(2))
.EQN 3 -54 28 0 0
{0:HP}NAME({0:x1}NAME,{0:x2}NAME,{0:x}NAME):{0:if}NAME({0:x}NAME<{0:x1}NAME,0,{1:if}NAME({0:x}NAME>{0:x2}NAME,1,\({0}0.5*(1-{0:cos}NAME({0:\p}NAME*(({0:x}NAME-{0:x1}NAME)/({0:x2}NAME-{0:x1}NAME)))))))
.EQN 1 54 130 0 0
{0:f1}NAME:0*{0:a}NAME
.EQN 0 8 233 0 0
{0:f2}NAME:(5)/(8)*{0:a}NAME
.EQN 0 10 234 0 0
{0:f3}NAME:(7)/(8)*{0:a}NAME
.EQN 0 8 235 0 0
{0:f4}NAME:1.2*{0:a}NAME
.TXT 5 -79 236 0 0
Cg a86.625000,86.625000,275
{\rtf\ansi \deff0{\colortbl;\red0\green0\blue128;}{\fonttbl{\f0
\fcharset0\fnil Arial;}}\plain\cf1\fs20 \pard Define the transfer 
functions of the highpass, two lowpass and two bandpass filters.\par 
The constants f1 and f2 control the steepness of the cutoffs; a is the 
folding frequency.\par The LP(f3,f4,f) filter is not required for 
invertibility, but it makes the bandpass kernels smaller.}
.EQN 9 -1 77 0 0
({0:B1}NAME)[({0:i}NAME,{0:k}NAME):{0:LP}NAME({0:f3}NAME,{0:f4}NAME,({0:\r}NAME)[({0:i}NAME,{0:k}NAME))*{0:HP}NAME({0:f1}NAME,({0:a}NAME)/(2),({0:\r}NAME)[({0:i}NAME,{0:k}NAME))*{0:cos}NAME(({0:\q}NAME)[({0:i}NAME,{0:k}NAME))
.EQN 0 39 78 0 0
({0:B2}NAME)[({0:i}NAME,{0:k}NAME):{0:LP}NAME({0:f3}NAME,{0:f4}NAME,({0:\r}NAME)[({0:i}NAME,{0:k}NAME))*{0:HP}NAME({0:f1}NAME,({0:a}NAME)/(2),({0:\r}NAME)[({0:i}NAME,{0:k}NAME))*{0:cos}NAME(({0:\q}NAME)[({0:i}NAME,{0:k}NAME)-({0:\p}NAME)/(2))
.EQN 6 -39 132 0 0
({0:L1}NAME)[({0:i}NAME,{0:k}NAME):{0:LP}NAME({0:f1}NAME,({0:a}NAME)/(2),({0:\r}NAME)[({0:i}NAME,{0:k}NAME))
.EQN 0 26 133 0 0
({0:L0}NAME)[({0:i}NAME,{0:k}NAME):{0:LP}NAME({0:f2}NAME,{0:a}NAME,({0:\r}NAME)[({0:i}NAME,{0:k}NAME))
.EQN 0 26 144 0 0
({0:H0}NAME)[({0:i}NAME,{0:k}NAME):{0:HP}NAME({0:f2}NAME,{0:a}NAME,({0:\r}NAME)[({0:i}NAME,{0:k}NAME))
.TXT 5 -52 165 0 0
Cg a87.625000,87.625000,61
{\rtf\ansi \deff0{\colortbl;\red0\green0\blue128;}{\fonttbl{\f0
\fcharset0\fnil Arial;}}\plain\cf1\fs20 \pard Compute functions to 
show that the constraints are satisfied.}
.EQN 4 0 68 0 0
({0:M1}NAME)[({0:i}NAME,{0:k}NAME):((({0:B1}NAME)[({0:i}NAME,{0:k}NAME)))^(2)+((({0:B2}NAME)[({0:i}NAME,{0:k}NAME)))^(2)
.EQN 0 26 241 0 0
({0:M2}NAME)[({0:i}NAME,{0:k}NAME):((({0:H0}NAME)[({0:i}NAME,{0:k}NAME)))^(2)+((({0:L0}NAME)[({0:i}NAME,{0:k}NAME)))^(2)
.EQN 0 28 244 0 0
({0:M3}NAME)[({0:i}NAME,{0:k}NAME):(((({0:L1}NAME)[({0:i}NAME,{0:k}NAME)))^(2)+((({0:B1}NAME)[({0:i}NAME,{0:k}NAME)))^(2)+((({0:B2}NAME)[({0:i}NAME,{0:k}NAME)))^(2))
.EQN 5 -54 247 0 0
({0:M4}NAME)[({0:i}NAME,{0:k}NAME):((({0:H0}NAME)[({0:i}NAME,{0:k}NAME)))^(2)+((({0:L0}NAME)[({0:i}NAME,{0:k}NAME)))^(2)*(((({0:L1}NAME)[({0:i}NAME,{0:k}NAME)))^(2)+((({0:B1}NAME)[({0:i}NAME,{0:k}NAME)))^(2)+((({0:B2}NAME)[({0:i}NAME,{0:k}NAME)))^(2))
.TXT 0 54 248 0 0
Cg a63.625000,63.625000,34
{\rtf\ansi \deff0{\colortbl;\red0\green0\blue128;}{\fonttbl{\f0
\fcharset0\fnil Arial;}}\plain\cf1\fs20 \pard M4 = 1 is Simoncelli's 
constraint.}
.TXT 6 -54 166 0 0
Cg a87.625000,87.625000,61
{\rtf\ansi \deff0{\colortbl;\red0\green0\blue128;}{\fonttbl{\f0
\fcharset0\fnil Arial;}}\plain\cf1\fs20 \pard Plot L{\fs16\dn 0}(u,v), H{
\fs16\dn 0}(u,v) and the sum of their squared magnitudes}
.EQN 3 11 273 0 0
{0:L0}NAME{1 4 3 215 45 1 21 20 0 1 1 1 4 -1 1 0 1 1 1 4 -1 1 0 1 0 1 2 -1 1 6 0 16777215 0 90 2 NO-TITLE}{57}
2 5 21 21 0 1 1.5 7
1 5 4 3 0 4 1 0
3 1 2 0.1
.EQN 0 22 274 0 0
{0:H0}NAME{1 4 3 215 45 1 21 20 0 1 1 1 4 -1 1 0 1 1 1 4 -1 1 0 1 0 1 2 -1 1 6 0 16777215 0 90 2 NO-TITLE}{57}
2 5 21 21 0 1 1.5 7
1 5 4 3 0 4 1 0
3 1 2 0.1
.EQN 0 22 275 0 0
{0:M2}NAME{1 4 3 215 45 1 21 20 0 1 1 1 4 -1 1 0 1 1 1 4 -1 1 0 1 0 0 2 0 1 6 0 16777215 0 90 2 NO-TITLE}{57}
2 5 21 21 0 1 1.5 7
1 5 4 3 0 4 1 0
3 1 2 0.1
.TXT 27 -55 167 0 0
Cg a87.625000,87.625000,76
{\rtf\ansi \deff0{\colortbl;\red0\green0\blue128;}{\fonttbl{\f0
\fcharset0\fnil Arial;}}\plain\cf1\fs20 \pard Plot the two steerable 
bandpass filters and their sum of squared magnitudes.}
.EQN 3 11 276 0 0
{0:B1}NAME{1 4 3 215 45 1 21 20 0 1 1 1 4 -1 1 0 1 1 1 4 -1 1 0 1 1 0 4 -1 1 6 0 16777215 0 90 2 NO-TITLE}{57}
2 5 21 21 0 1 1.5 7
1 5 4 3 0 4 1 0
3 1 2 0.1
.EQN 0 22 277 0 0
{0:B2}NAME{1 4 3 215 45 1 21 20 0 1 1 1 4 -1 1 0 1 1 1 4 -1 1 0 1 1 0 4 -1 1 6 0 16777215 0 90 2 NO-TITLE}{57}
2 5 21 21 0 1 1.5 7
1 5 4 3 0 4 1 0
3 1 2 0.1
.EQN 0 22 278 0 0
{0:M1}NAME{1 4 3 215 45 1 21 20 0 1 1 1 4 -1 1 0 1 1 1 4 -1 1 0 1 0 1 2 -1 1 6 0 16777215 0 90 2 NO-TITLE}{57}
2 5 21 21 0 1 1.5 7
1 5 4 3 0 4 1 0
3 1 2 0.1
.TXT 34 -55 254 0 0
Cg a87.625000,87.625000,103
{\rtf\ansi \deff0{\colortbl;\red0\green0\blue128;}{\fonttbl{\f0
\fcharset0\fnil Arial;}}\plain\cf1\fs20 \pard Plot L{\fs16\dn 1}(u,v), 
the sum of the bandpass squared magnitudes, their sum, and the overall 
transfer function.}
.EQN 2 0 252 0 0
{0:L1}NAME{1 4 3 215 45 1 21 20 0 1 1 1 4 -1 1 0 1 1 1 4 -1 1 0 1 0 1 2 -1 1 6 0 16777215 0 90 2 NO-TITLE}{57}
2 5 21 21 0 1 1.5 7
1 5 4 3 0 4 1 0
3 1 2 0.1
.EQN 0 22 253 0 0
{0:M1}NAME{1 4 3 215 45 1 21 20 0 1 1 1 4 -1 1 0 1 1 1 4 -1 1 0 1 0 1 2 -1 1 6 0 16777215 0 90 2 NO-TITLE}{57}
2 5 21 21 0 1 1.5 7
1 5 4 3 0 4 1 0
3 1 2 0.1
.EQN 0 22 251 0 0
{0:M3}NAME{1 4 3 215 45 1 21 20 0 1 1 1 4 -1 1 0 1 1 1 4 -1 1 0 1 0 0 2 0 1 6 0 16777215 0 90 2 NO-TITLE}{57}
2 5 21 21 0 1 1.5 7
1 5 4 3 0 4 1 0
3 1 2 0.1
.EQN 0 22 250 0 0
{0:M4}NAME{1 4 3 215 45 1 21 20 0 1 1 1 4 -1 1 0 1 1 1 4 -1 1 0 1 0 0 2 0 1 6 0 16777215 0 90 2 NO-TITLE}{57}
2 5 21 21 0 1 1.5 7
1 5 4 3 0 4 1 0
3 1 2 0.1
.TXT 26 -66 291 0 0
Cg a87.625000,87.625000,86
{\rtf\ansi \deff0{\colortbl;\red0\green0\blue128;}{\fonttbl{\f0
\fcharset0\fnil Arial;}}\plain\cf1\fs20 \pard {\b Demonstrate steerability} 
- generate and synthesize 60 degree filters and compare them.}
.EQN 5 1 294 0 0
({0:B3}NAME)[({0:i}NAME,{0:k}NAME):{0:LP}NAME({0:f3}NAME,{0:f4}NAME,({0:\r}NAME)[({0:i}NAME,{0:k}NAME))*{0:HP}NAME({0:f1}NAME,({0:a}NAME)/(2),({0:\r}NAME)[({0:i}NAME,{0:k}NAME))*{0:cos}NAME(({0:\q}NAME)[({0:i}NAME,{0:k}NAME)-({0:\p}NAME)/(3))
.EQN 0 38 295 0 0
{0:B4}NAME:(1)/(2)*{0:B1}NAME+(\(3))/(2)*{0:B2}NAME
.EQN 4 -39 283 0 0
{0:B3}NAME{1 4 3 215 45 1 21 20 0 1 1 1 4 -1 1 0 1 1 1 4 -1 1 0 1 1 0 4 -1 1 6 0 16777215 0 90 2 NO-TITLE}{57}
2 5 21 21 0 1 1.5 7
1 5 4 3 0 4 1 0
3 1 2 0.1
.EQN 0 22 285 0 0
{0:B4}NAME{1 4 3 215 45 1 21 20 0 1 1 1 4 -1 1 0 1 1 1 4 -1 1 0 1 1 0 4 -1 1 6 0 16777215 0 90 2 NO-TITLE}{57}
2 5 21 21 0 1 1.5 7
1 5 4 3 0 4 1 0
3 1 2 0.1
.EQN 0 22 289 0 0
{0:B4}NAME-{0:B3}NAME{1 4 3 215 45 1 21 20 0 1 1 1 4 -1 1 0 1 1 1 4 -1 1 0 1 1 0 4 -1 1 6 0 16777215 0 90 2 NO-TITLE}{57}
2 5 21 21 0 1 1.5 7
1 5 4 3 0 4 1 0
3 1 2 0.1
.TXT 28 -43 158 0 0
Cg a85.500000,85.500000,55
{\rtf\ansi \deff0{\colortbl;\red0\green0\blue128;}{\fonttbl{\f0
\fcharset0\fnil Arial;}}\plain\cf1\fs20 \pard Compute the convolution 
kernels using the centered DFT.}
.EQN 0 49 299 0 0
({0:W}NAME)[({0:i}NAME,{0:k}NAME):(1)/({0:N}NAME+1)*{0:exp}NAME(-{0:j}NAME*2*{0:\p}NAME*({0:i}NAME-{0:a}NAME)*({0:k}NAME-{0:a}NAME)/({0:N}NAME+1))
.EQN 6 -50 135 0 0
{0:b1}NAME:{0:W}NAME*({0:j}NAME*{0:B1}NAME)*{0:W}NAME
.EQN 0 16 136 0 0
{0:b2}NAME:{0:W}NAME*({0:j}NAME*{0:B2}NAME)*{0:W}NAME
.EQN 0 18 139 0 0
{0:l0}NAME:{0:W}NAME*{0:L0}NAME*{0:W}NAME
.EQN 0 13 153 0 0
{0:l1}NAME:{0:W}NAME*{0:L1}NAME*{0:W}NAME
.EQN 0 14 138 0 0
{0:h0}NAME:{0:W}NAME*{0:H0}NAME*{0:W}NAME
.TXT 4 -61 170 0 0
Cg a87.625000,87.625000,47
{\rtf\ansi \deff0{\colortbl;\red0\green0\blue128;}{\fonttbl{\f0
\fcharset0\fnil Arial;}}\plain\cf1\fs20 \pard Plot the two bandpass 
filter impulse responses.}
.EQN 3 8 261 0 0
{0:b1}NAME{1 4 3 215 45 1 34 29 0 1 1 1 4 -1 1 0 1 1 1 4 -1 1 0 1 0 1 2 -1 1 6 0 16777215 0 90 2 NO-TITLE}{57}
2 5 21 21 0 1 1.5 7
1 5 4 3 0 4 1 0
3 1 2 0.1
.EQN 0 35 262 0 0
{0:b2}NAME{1 4 3 215 45 1 34 29 0 1 1 1 4 -1 1 0 1 1 1 4 -1 1 0 1 0 1 2 -1 1 6 0 16777215 0 90 2 NO-TITLE}{57}
2 5 21 21 0 1 1.5 7
1 5 4 3 0 4 1 0
3 1 2 0.1
.TXT 53 -42 171 0 0
Cg a86.625000,86.625000,63
{\rtf\ansi \deff0{\colortbl;\red0\green0\blue128;}{\fonttbl{\f0
\fcharset0\fnil Arial;}}\plain\cf1\fs20 \pard Plot the highpass and 
the two lowpass filter impulse responses.}
.EQN 3 -1 267 0 0
{0:h0}NAME{1 4 3 215 45 1 28 26 0 1 1 1 4 -1 1 0 1 1 1 4 -1 1 0 1 0 1 2 -1 1 6 0 16777215 0 90 2 NO-TITLE}{57}
2 5 21 21 0 1 1.5 7
1 5 4 3 0 4 1 0
3 1 2 0.1
.EQN 0 29 271 0 0
{0:l0}NAME{1 4 3 215 45 1 28 26 0 1 1 1 4 -1 1 0 1 1 1 4 -1 1 0 1 0 1 2 -1 1 6 0 16777215 0 90 2 NO-TITLE}{57}
2 5 21 21 0 1 1.5 7
1 5 4 3 0 4 1 0
3 1 2 0.1
.EQN 0 29 272 0 0
{0:l1}NAME{1 4 3 215 45 1 28 26 0 1 1 1 4 -1 1 0 1 1 1 4 -1 1 0 1 0 1 2 -1 1 6 0 16777215 0 90 2 NO-TITLE}{57}
2 5 21 21 0 1 1.5 7
1 5 4 3 0 4 1 0
3 1 2 0.1
.EQN 34 17 210 0 0
{0:i}NAME${0:k}NAME$({0:l0}NAME)[({0:i}NAME,{0:k}NAME)={0}?_n_u_l_l_
.EQN 7 0 209 0 0
{0:i}NAME${0:k}NAME$({0:l1}NAME)[({0:i}NAME,{0:k}NAME)={0}?_n_u_l_l_
.EQN 7 -73 186 0 0
1000*{0:b2}NAME={150112}?_n_u_l_l_
.EQN 0 73 211 0 0
{0:i}NAME${0:k}NAME$({0:h0}NAME)[({0:i}NAME,{0:k}NAME)={0}?_n_u_l_l_
.EQN 7 0 212 0 0
{0:i}NAME${0:k}NAME$({0:b1}NAME)[({0:i}NAME,{0:k}NAME)={0}?_n_u_l_l_
.EQN 7 0 221 0 0
{0:i}NAME${0:k}NAME$({0:b2}NAME)[({0:i}NAME,{0:k}NAME)={0}?_n_u_l_l_
.EQN 8 -1 201 0 0
{0:d}NAME:({2,2}ö{0:a}NAMEö{0:N}NAME+1ö{0:a}NAMEö{0:N}NAME+1)
.EQN 6 0 202 0 0
{0:d}NAME={0}?_n_u_l_l_
.TXT 9 -73 222 0 0
Cg a86.625000,86.625000,71
{\rtf\ansi \deff0{\colortbl;\red0\green0\blue128;}{\fonttbl{\f0
\fcharset0\fnil Arial;}}\plain\cf1\fs20 \pard Round off to four digits 
and scale the values for writing kernel files.}
.EQN 0 65 208 0 0
{0:b}NAME:1
.EQN 7 -66 187 0 0
({0:K2_L0}NAME)[({0:i}NAME,{0:k}NAME):({0:floor}NAME(10000*{0:Re}NAME(({0:l0}NAME)[({0:i}NAME,{0:k}NAME))+0.5))/({0:b}NAME)
.EQN 0 29 188 0 0
({0:K2_L1}NAME)[({0:i}NAME,{0:k}NAME):({0:floor}NAME(10000*{0:Re}NAME(({0:l1}NAME)[({0:i}NAME,{0:k}NAME))+0.5))/({0:b}NAME)
.EQN 0 29 190 0 0
({0:K2_H0}NAME)[({0:i}NAME,{0:k}NAME):({0:floor}NAME(10000*{0:Re}NAME(({0:h0}NAME)[({0:i}NAME,{0:k}NAME))+0.5))/({0:b}NAME)
.EQN 6 -58 199 0 0
({0:K2_B1}NAME)[({0:i}NAME,{0:k}NAME):({0:floor}NAME(10000*{0:Re}NAME(({0:b1}NAME)[({0:i}NAME,{0:k}NAME))+0.5))/({0:b}NAME)
.EQN 0 29 200 0 0
({0:K2_B2}NAME)[({0:i}NAME,{0:k}NAME):({0:floor}NAME(10000*{0:Re}NAME(({0:b2}NAME)[({0:i}NAME,{0:k}NAME))+0.5))/({0:b}NAME)
.EQN 6 -29 216 0 0
{0:i}NAME${0:k}NAME$({0:K2_L0}NAME)[({0:i}NAME,{0:k}NAME)={0}?_n_u_l_l_
.EQN 0 19 215 0 0
{0:i}NAME${0:k}NAME$({0:K2_L1}NAME)[({0:i}NAME,{0:k}NAME)={0}?_n_u_l_l_
.EQN 0 19 217 0 0
{0:i}NAME${0:k}NAME$({0:K2_H0}NAME)[({0:i}NAME,{0:k}NAME)={0}?_n_u_l_l_
.EQN 0 18 213 0 0
{0:i}NAME${0:k}NAME$({0:K2_B1}NAME)[({0:i}NAME,{0:k}NAME)={0}?_n_u_l_l_
.EQN 0 17 214 0 0
{0:i}NAME${0:k}NAME$({0:K2_B2}NAME)[({0:i}NAME,{0:k}NAME)={0}?_n_u_l_l_
.TXT 9 -73 223 0 0
Cg a85.625000,85.625000,103
{\rtf\ansi \deff0{\colortbl;\red0\green0\blue128;}{\fonttbl{\f0
\fcharset0\fnil Arial;}}\plain\cf1\fs20 \pard Write kernel files for 
input to the WiT 2-D convolution operator (Use b = 10,000 for unscaled 
kernels).}
.EQN 5 3 185 0 0
{0:WRITEPRN}NAME({0:K2_L0.arr}NAME):{0:d}NAME
.EQN 0 27 203 0 0
{0:APPENDPRN}NAME({0:K2_L0.arr}NAME):{0:K2_L0}NAME
.EQN 3 -27 189 0 0
{0:WRITEPRN}NAME({0:K2_L1.arr}NAME):{0:d}NAME
.EQN 0 27 204 0 0
{0:APPENDPRN}NAME({0:K2_L1.arr}NAME):{0:K2_L1}NAME
.EQN 3 -27 191 0 0
{0:WRITEPRN}NAME({0:K2_H0.arr}NAME):{0:d}NAME
.EQN 0 27 205 0 0
{0:APPENDPRN}NAME({0:K2_H0.arr}NAME):{0:K2_H0}NAME
.EQN 3 -27 197 0 0
{0:WRITEPRN}NAME({0:K2_B1.arr}NAME):{0:d}NAME
.EQN 0 27 206 0 0
{0:APPENDPRN}NAME({0:K2_B1.arr}NAME):{0:K2_B1}NAME
.EQN 3 -27 198 0 0
{0:WRITEPRN}NAME({0:K2_B2.arr}NAME):{0:d}NAME
.EQN 0 27 207 0 0
{0:APPENDPRN}NAME({0:K2_B2.arr}NAME):{0:K2_B2}NAME
.TXT 6 -29 224 0 0
Cg a84.375000,84.375000,533
{\rtf\ansi \deff0{\colortbl;\red0\green0\blue128;}{\fonttbl{\f0
\fcharset0\fnil Arial;}}\plain\cf1\fs20 \pard 1.  E. P. Simoncelli, 
and W. T. Freeman, "The Steerable Pyramid: A Flexible Architecture for 
Multi-Scale Derivative Computation," {\i Proc. ICIP-95}: 444-447, 
1995.\par 2. P. J. Burt, and E. H. Adelson, "The Laplacian Pyramid as 
a Compact Image Code," \par {\i IEEE Trans}. {\b C-31}:532-540, 1983.
\par 3. E. P. Simoncelli, W. T. Freeman, E. H. Adelson, and D. J. 
Heeger, "Shiftable Multiscale Transforms," \par {\i I}{\i EEE Trans}. {
\b IT-38}(2):587-607, 1992.\par 4. W. T. Freeman and E. H. Adelson, 
"The Design and Use of Steerable Filters," {\i IEEE Trans} {\b PAMI-13}
(9):891-906, 1991.\par }
.TXT 21 35 298 0 0
Cg a55.625000,55.625000,9
{\rtf\ansi \deff0{\colortbl;\red0\green0\blue128;}{\fonttbl{\f0
\fcharset0\fnil Arial;}}\plain\cf1\fs20 \pard == END ==}