200 likes | 302 Views
Fast spherical hormonic transform. 歐永俊 資工三 b92b02053. Background. Direct computation: O(N^4). Separable!!. BUT. Can be reduced to O(N^3). Background. Outer IDFT:. O(N^3). O(n^2lgn). FFT. Inner SHM:. O(N^3). O(n^2lgn). Fast SHM. Background.
E N D
Fast spherical hormonic transform 歐永俊 資工三 b92b02053
Background Direct computation: O(N^4) Separable!! BUT Can be reduced to O(N^3)
Background Outer IDFT: O(N^3) O(n^2lgn) FFT Inner SHM: O(N^3) O(n^2lgn) Fast SHM
Background The problem of a fast spherical harmonic transform is now reduced to the efficient calculation of these discrete Legendre transforms. for an arbitrary input vector s with kth component [s]k.
STEP1:using the Legendre recurrence The recurrence relation satisfied by degree:L Stands at level L, iterating the recurrence formula forward r steps associated legendre functions, satisfy follows: degree: r,r-1
STEP1:using the Legendre recurrence In matrix form( recurrence relation): Purpose : rewrite projections onto high degree Legendre polynomials as sums of projections onto (shifted) legendre Polynomials of lower degree .
STEP1:using the Legendre recurrence stored at early time degree at most r Higher degree inner product can be compute as inner products of stored data and (precomputed) sampled values of the polynomials A,B( remember it satisfies recurrence relation of asso.legendre) “A hint to divide and conquer but not enough!”
STEP2:smoothing and subsampling(working inthe “cosine transform” domain) What is not enough? Projections of a data vector onto lower degree trigonometric polynomials should be computed more efficiently. How? Can be done by first lowpass filtering(smoothing) and subsampling the data vector(in cosine transform domain).Then the complexity of computing these inner products grows linearly with n.
STEP2:smoothing and subsampling(working inthe “cosine transform” domain) denote the n-dinmensional orthogonal DCT matrix comprised of normalized sampled cosines DCT
STEP2:smoothing and subsampling(working inthe “cosine transform” domain) fact1: orthogonality of C implies fact2: for any trigonometric polynomial Q of degree at most N-1 the cosine coefficients vanish for n>deg(Q) Computation of the low degree Legendre projections accomplished with very few operations
STEP2:smoothing and subsampling(working inthe “cosine transform” domain) define: critically sampled lowpass operator (bandwidth p) define: truncation operator (keeps first p coordinates of given input vector) is to first to compute cosine representation of a vector the effect of of length N, then remove all frequency components beyond p from s (smoothing), and finally, to keep only those samples necessary to represent this smoothed version(subsampling).
STEP2:smoothing and subsampling(working inthe “cosine transform” domain) Q: degree p Note: the effect of applying to is to simply sample Q on coarser grid O(nlgn) by fast DCT/IDCT
STEP3:divide and conquer low degree transform coefficients: high degree transform coefficients: lower part: (a discrete transform of size N/2)
STEP3:divide and conquer higher part: (Page 7) ( ) ( ) with degree at most N/2 (a discrete transform of size N/2)
Summary O(NlgN) page 15 two problems of size N/2 Page 16,17 Overall is
reference 1.A fast spherical harmonics transform algorithm R Suda, M Takami - Math. Comp, 2002 2. D. Healy, S. Moore, and D. Rockmore, “An FFT for the 2-sphere - improvements and variations” Proc. of ICASSP-96 volume 3.pp.1323-1326