Get 20M+ Full-Text Papers For Less Than $1.50/day. Start a 14-Day Trial for You or Your Team.

Learn More →

Angular multiselectivity with spherical wavelets

Angular multiselectivity with spherical wavelets We construct spherical wavelets based on approximate identities that are directional, i.e. not rotation-invariant, and have an adaptive angular selectivity. The problem of how to nd a proper representation of distinct kinds of details of real images, ranging from highly directional to fully isotropic ones, was quite intensively studied for the case of signals over the Euclidean space. However, the present paper is the rst attempt to deal with this task in the case of spherical signals. A multiselectivity scheme, similar to that proposed for R -functions, is presented. Key words and phrases: spherical wavelets, Poisson kernel, angular selectiv- ity, directional wavelet frames, angular multiresolution 2010 Mathematics Subject Classi cation Number: 42C40, 42C15, 94A12, 65T60 1 Introduction Recently, directional wavelets based on approximate identities have been intro- duced [4, 5, 7, 8, 9]. Unfortunately, wavelet families constructed so far have a very low angular resolution, e.g., the rst directional derivative of the Poisson kernel is a wavelet with separable spherical variables (#; ') of the form [1] [1] (x) = (#) cos '; where x = (cos #; sin # cos '; sin # sin ') 2 S; # 2 [0; ]; ' 2 [0; 2): (Unless it leads to misunderstandings, we identify x 2 S with its spherical coordinates (#; ').) It is depicted on the left-hand side of Figure 1 and it can be seen that its ability to detect features with close orientations is not satisfactory. In [7] second order wavelets are used to analyse the land cover data of the world, where the analyzed objects are in principle curves. The wavelets are of the form [2] [2] [2] 2 (x) = (#) + (#) cos '; ;1 ;2 West Pomeranian University of Technology in Szczecin, School of Mathematics, al. Pias- t ow 17, 70{310 Szczecin, Poland arXiv:1804.03044v1 [math.CA] 6 Apr 2018 2 Figure 1: The rst directional derivative of the Poisson kernel on S vs. a wavelet with a bigger angular selectivity { behavior around the north pole see [9, Appendix A], i.e. again functions with a low angular resolution. The aim of the present paper is to construct wavelets with a variable angular selectivity such that their ability to detect directional features can be adapted to the ana- lyzed signal. An example of such a wavelet is presented on the right-hand-side of Figure 1. It is a function with separable variables and Figure 3 shows the dependence of the function of the longitudinal variable on a parameter that can be chosen almost freely. The investigation is inspired by the results obtained by Antoine and Jacques in [2, 1, 11] for the case of wavelets over R . The angular selectivity of spherical wavelets was analysed by [3], but in this case wavelets based on the group theory were involved (the di erence between the two constructions is discussed in [8]). The angular resolving power was de ned in [3, Subsection 9.2.4.2]. Nonetheless, we do not intend to adapt it to our purposes, since we nd the parameter introduced in Section 3 to be a better indicator of the angular selectivity in the case of wavelets constructed in the present paper. Moreover, one that can be a priori chosen. The paper is organized as follows. After an introduction of necessary notions in Section 2, we present a construction of two wavelets with a steerable angular selectivity in Section 3. An angular multiselectivity analysis is then discussed in Section 4. 2 Preliminaries 3 3 2 By S we denote the unit two-dimensional sphere in R , S = f 2 R : jj = 1g, with the measure d invariant under the rotation group and such that d = 4. # 2 [0; ] is the polar (colatitudinal) coordinate, and ' 2 [0; 2) { the azimuthal (longitudinal) coordinate. The scalar product of spherical functions , is de ned by h; i = (x) (x) d(x); S 3 and L (S ) is the set of functions such that k k 2 := h ; i < 1: L (S) Note the di erence to the notation usually used in research concerning n- dimensional spheres, where the integral in the de nition of the scalar product is divided by  , the Lebesgue measure of S , see e.g. [8]. The spherical harmonics Y , l 2 N , k = l;l + 1; : : : ; l, are given by 2l + 1 (l k)! k k k ik' Y (#; ') = (1) P (cos #) e ; l l 4 (l + k)! Y = Y l l for 0  k  l, where P are the (associated) Legendre polynomials (functions), 1 d 0 2 l P (t) = P (t) = (t 1) (Legendre polynomials) l l 2 l! dt k k 2 k=2 P (t) = (1) (1 t ) P (t) for 0 < k  l (associated Legendre functions). dt 2 2 The spherical harmonics build an orthonormal basis for L (S ). For an L (S ){ function , the series 1 l X X k k Y ; Y l l l=1 k=l is the Fourier series of in terms of the spherical harmonics. The constants k k = Y ; are called the Fourier coecients of . l l T denotes the one-dimensional torus T = [0; 2). Functions on T will be iden- ti ed with their 2-periodic extension on R, such that C(T) denotes the class of continuous functions on T that can be continuously extended. An L (T)- function can be represented as the Fourier series ikt f (t)  f e ; k=1 where the Fourier coecients f , k 2 Z, are given by ikt f = f (t) e dt: Similarly, the Fourier transform of an L (R) function F is de ned by i!t F (!) = F (t) e dt;  2 R; and it can be inverted by i!t F (t) = F (!) e d: 1 4 The functions we consider in the present paper are both integrable and con- tinuous such that integrals and series in the Fourier and the inverse Fourier transform are (pointwise) convergent. An L (R)-function F can be periodized to an L (T)-function f by f (t) = f (t + 2j); t 2 T; j=1 and their Fourier transforms are linked by the following Poisson summation formula. Theorem 2.1 Let F 2 L (R). Then its periodization f exists for almost every t 2 T, it is integrable and kfk 1  kFk 1 . Further, L (T) L (R) f = F (k): A family of vectors fg ;  2 Ig  H in a Hilbert space H indexed by a measure space I with a positive measure  is called a frame with weight  if the mapping 7! g is weakly measurable, i.e.  7! hg ; ui is measurable and 2 2 2 Akuk  jhg ; uij d(x)  Bkuk for all u 2 H and 0 < A  B. If fg ;  2 Ig is a frame, then the mapping u 7! fhg ; ui ;  2 Ig is invertible. The wavelet de nition we use originates from [10]. In the case of the two- dimensional sphere and the weight function  7! , it is the following one. 2 2 De nition 2.2 The family f g  L (S ) is called a wavelet (family) of 2R order m if it satis es A(2l + 1)  ( )  B(2l + 1) (1) k=l for some positive constants A and B independent of l 2 N , l > m, and ( ) = 0 (2) k=l for l 2 N , l  m. 2 2 The wavelet transform of an L (S )-function is given by W f (; g) = (g x) f (x) d (x); g 2 SO(3); 4 n S 5 Figure 2: Grid construction and it is invertible by methods known from the frame theory. Moreover, only a discrete set of the wavelet coecients W f (; g) is necessary to reconstruct the analyzed signal. The rst step is a scale discretization. According to [10, Theorem 3.1] there exist constants a and X such that for any sequence R = ( ) with   a and 1 <  = < X the family f (g );  2 R; g 2 j j2N 0 0 j j+1  j 0 j 2 2 SO(3)g is a frame for L (S ). In order to characterize discrete sets of rotations, we need the following de nition (compare [10, De nition 4.1]). De nition 2.3  is a grid of type ( ;  ) if it is a discrete measurable set of 2 1 2 1 rotations in SO(3) = S  S , constructed in the following way. Let P = 2 2 fO : = 1; : : : ; K g be a measurable partition of S into simply connected 2 2 sets such that the diameter of each set (measured in the geodesic distance) is 2 2 not larger than  . Choose from each of the sets O an arbitrary point x = 2 2 2 2 1 (# ; ' ). Now, for a fxed , let P ( ) = fO : = 1; : : : ; K ( )g 2 1 2 1 1 2 2 2 ( ; ) 2 1 be a measurable partition of S into K ( ) simply connected sets of a diameter 1 2 not larger than  . Choose from each of the sets O , = 1; : : : ; K ( ), 1 1 1 2 ( ; ) 2 1 an arbitrary point ' . Then,  is the set of rotations given by ( ; ) 2 1 1 2 2 g = g (' ) g (# ) g (' ); ( ; ) 1 1 2 2 1 ( ; ) 2 1 2 2 where g ( ) is the rotation in the plane ( ;  ) with the rotation angle . j j j+1 Loosely speaking, one has a tight enough discrete set of points fx g 2 S , and for each point a tight enough discrete set of rotations around it.  is a set of SO(3)-rotations being a composition of a tranlastion to a point x and a rotation around this point, as illustrated in Figure 2. The next statement re ects the content of [10, Theorem 4.2] for the case of the two-dimensional sphere. 6 Figure 3: f for  = 1; 2; 4; 8; 16, ' 2 [=2; 3=2] Theorem 2.4 Let be a C {wavelet family with the property that f ; j 2 ;x N ; x 2 S g is a semi{continuous frame. Then, for each j 2 N there exist 0 0 j j numbers  ,  such that 2 1 1 j f (g ); j 2 N ; g j j 2  g j j 0 j ( ; ) ( ; ) 2 1 2 1 j j 2 n j is a frame for L (S ), provided that  is a grid of type ( ;  ). 2 1 3 Construction of the wavelet In this section, we shall present a wavelet family with separated variables and a steerable angular selectivity. The construction is based on the Poisson kernel, 1 1 r 1 p (x) = = (2l + 1)r P (cos #); (3) 2 3=2 4 (1 2r cos # + r ) 4 l=0 where r = e , since our aim is to obtain a wavelet with a simple representation as a function of the spherical variables. On the other hand, in order to prove the wavelet property, some knowledge about the Fourier coecients is necessary. The idea is to choose a periodized di erence of Gaussians dependent on a pa- rameter  2 [1;1) as the function of the longitudinal variable ', 2 2 2 2 () 2 2 F () = e e ;  2 R; f (') = F (' + 2j); ' 2 [0; 2): j2Z Figure 3 presents the function f for several values of  . The challenge is to nd an appropriate function of the colatitudinal variable # such that the inequalities (1) are satis ed. The result of our investigation is that the following families: @ @ (x) = (#; ') = ! (#) f ('); ! (#) =  sin # r  r  [p (x)] ; @r @r 5 2 (x) =  (#; ') =  (#) f (');  (#) =  sin # r  [p (x)] ; @r 7 r = e , are wavelets. In order to prove the wavelet condition (1), we need an estimation of the Fourier coecients (with respect to the spherical harmonics) of and , and since they are functions with separated variables, we rst consider the Fourier series expansion (in exponential functions) of f . Since 2 2 ! ! i! b 2 2 2 2 F (!) = e e ; by the Poisson summation formula the Fourier coecients of f are given by 0; for even k; p 2 (f ) = k k 2 2 e for odd k: Thus, X 2 1 2 ik' f (') = e e : (4) k22Z+1 In order to estimate the Fourier coecients of and , we estimate the mag- nitude of the coecients of sin # 2 l ! (#) = (2l + 1)l r P (cos #) (5) l=1 2 4 2 4 2 2 5 r r(10 19r + r ) (3 14r + 5r ) cos # r(9 r ) cos # sin # = ; 2 7=2 4(1 2r cos # + r ) (6) and sin # (#) = (2l + 1)l(l 1)r P (cos #) (7) l=2 2 2 4 2 2 2 5 r 5 23r + 2r + 4r(7 + r ) cos # (15 + r ) cos # sin # = ; (8) 2 7=2 4(1 2r cos # + r ) r = e , developed into series with respect to the associated Legendre polyno- mials of an odd order k. Lemma 3.1 For r 2 (0; 1), ! (#) = P (cos #) and l=1 (#) = P (cos #) l=2 8 with 3 l+5 3 2 l+3 (l + 2)(l + 3)(l + 4)(l + 5) r (l + 2)(l + 3) (5l + 21l 8) r (2l + 3)(2l + 5)(2l + 7)(2l + 9) (2l 1)(2l + 3)(2l + 5)(2l + 9) 2 5 4 3 2 l+1 2(l + 1) (10l + 65l + 83l 134l 186l + 36) r (2l 3)(2l 1)(2l + 3)(2l + 5)(2l + 7) 2 4 3 2 l1 2(l 1) (5l 2l 41l + 14l + 60) r (9) (2l 5)(2l 3)(2l + 3)(2l + 5) 2 4 3 2 l3 (l 3) (l 2)(l 1)(20l 100l + 39l + 239l 48) r (2l 7)(2l 5)(2l 3)(2l 1)(1 + 2l)(2l + 3) 2 l5 (l 5) (l 4)(l 3)(l 2)(l 1) r (2l 7)(2l 5)(2l 3)(2l 1) and 2 2 l+5 2 2 2 l+3 (l + 2)(l + 3)(l + 4) (l + 5) r (l + 2) (l + 3) (5l + 21l 8) r (2l + 3)(2l + 5)(2l + 7)(2l + 9) (2l 1)(2l + 3)(2l + 5)(2l + 9) 5 4 3 2 l+1 2l(l + 1)(10l + 65l + 83l 134l 186l + 36) r (2l 3)(2l 1)(2l + 3)(2l + 5)(2l + 7) 4 3 2 l1 2(l 2)(l 1)(5l 2l 41l + 14l + 60) r (2l 5)(2l 3)(2l + 3)(2l + 5) 4 3 2 l3 (l 4)(l 3)(l 2)(l 1)(20l 100l + 39l + 239l 48) r (2l 7)(2l 5)(2l 3)(2l 1)(1 + 2l)(2l + 3) l5 (l 6)(l 5)(l 4)(l 3)(l 2)(l 1) r (2l 7)(2l 5)(2l 3)(2l 1) (10) for l  9. The coecients , l = 1; 2; : : : ; 8, and , l = 2; 3; : : : ; 8, do not l l vanish. Proof. By the repeated application of [6, formulae (8.733.4) and (8.735.5)]: k k P (cos #) P (cos #) k1 l1 l+1 sin #P (cos #) = ; 2l + 1 (l k)(l k + 1) (l + k)(l + k + 1) k+1 k k sin #P (cos #) = P (cos #) P (cos #) l l+1 l1 2l + 1 2l + 1 9 we obtain (l 3)(l 2)(l 1)l 5 0 1 (2l + 1) sin #P (cos #) = P (cos #) l l5 (2l 7)(2l 5)(2l 3)(2l 1) (l 1)l(5l 9l 26) P (cos #) l3 (2l 7)(2l 3)(2l 1)(2l + 3) 5 4 3 2 2(10l + 15l 77l 93l + 121l + 60) + P (cos #) l1 (2l 5)(2l 3)(2l + 1)(2l + 3)(2l + 5) 4 3 2 2(5l + 18l 17l 54l + 36) P (cos #) l+1 (2l 3)(2l 1)(2l + 5)(2l + 7) 4 3 2 (l + 1)(l + 2)(20l + 140l + 219l 67l 60) + P (cos #) l+3 (2l 1)(2l + 1)(2l + 3)(2l + 5)(2l + 7)(2l + 9) (l + 1)(l + 2)(l + 3)(l + 4) P (cos #) l+5 (2l + 3)(2l + 5)(2l + 7)(2l + 9) for l  4. A substitution of this formula to (5) resp. (7) and index shifts yield coecients of ! resp.  as in (9) resp. (10) for l  9. The coecients , , l l l = 1; 2; : : : ; 8, are computed in the same manner. The di erence with respect to the case l  9 is caused by vanishing polynomials P . 1 1 b b Lemma 3.2 The Fourier coecients , l = 1; 2; : : : , and  , l = 2; 3; : : : , l l satisfy Z Z 1 1 2 2 d d 1 1 b b (2l + 1)A ; (2l + 1)A   (11) l l 0 0 for some constants A , A independent of l. Proof. It follows from Lemma 3.1 and formula (4) that X X 1 2 k 1 ik' (x) = ! (#) f (') = P (cos #) e e l=1 k22Z+1 X 2 2 k 1 ik' = e P (cos #) e : l;k Since 4  l(l + 1) 1 i' 1 P (cos #) e = Y (#; ') l l 2l + 1 and the spherical harmonics are orthogonal to one another, r r 2 1 4  l(l + 1)  l(l + 1) 1 b 2 2 2 2 =  e  =  e : l l 4  2l + 1  2(2l + 1) 10 Similarly, l(l + 1) 1 b 2 =  e : 2(2l + 1) For l  9, the integrals Z Z 1 1 2 2 1 d 1 d 1 1 b b and l l 2l + 1  2l + 1 0 0 are positive rational functions of l. It can be veri ed that Z Z 1 1 2 2 d and  d l l 0 0 are continuous for l  9 and that their limits for l ! 1 both equal . Consequently, Z Z 1 1 2 2 1 d 1 d 9 1 9 1 b b A := min and A := min l  l l2[9;1) 2l + 1  l2[9;1) 2l + 1 0 0 both exist and are positive. Similarly, for l = 1; 2; : : : ; 8 the numbers 1 d l 1 A := 2l + 1 and for l = 2; 3; : : : ; 8 the numbers 1 d l 1 A := 2l + 1 1 2 8 9 are positive. Inequalities (11) are thus satis ed with A = minfA ; A ; : : : ; A ; A g 2 3 8 9 and A = minfA ; A ; : : : ; A ; A g. Lemma 3.3 Let k 2 N. Then 3r (l + k)! ! (#)P (cos #) sin # d#  p and (12) k(l k)! 2 1 r 3r (l + k)! (#)P (cos #) sin # d#  p (13) k(l k)! 1 r k k Proof. Denote the left-hand-sides of (12) and (13) by !e and e , respectively. l l They are bounded from above by 6 k r sin #jp (r; #)j P (cos #) k l !e  d# and 2 7=2 4 (1 2r cos # + r ) 2  k sin #jp (r; #)j P (cos #) e  d#; 2 7=2 4 (1 2r cos # + r ) 0 11 where 2 4 2 4 2 2 p (r; #) = r(10 19r + r ) (3 14r + 5r ) cos # r(9 r ) cos # and 2 4 2 2 2 p (r; #) = 5 23r + 2r + 4r(7 + r ) cos # (15 + r ) cos #; see (6) and (8). The H older inequality together with k 2 [P (cos #)] (l + k)! d# = sin # k(l k)! (see [6, formula (7.122.1)], note that 2 in the denominator disappears because of the extended integration interval) yield !e  sup sin #jp (r; #)j r2(0;1) #2[0;] " # 1=2 Z  Z 2 1=2 12  k P (cos #) sin # d#  d# 2 7 (1 2r cos # + r ) sin # 0 0 1=2 r 231 (l + k)! 3r (l + k)! 21   p : 4 1024(1 r ) k(l k)! k(l k)! 2 (1 r ) Estimation (13) is obtained in a similar way. It is to be respected that the function (r; #) 7! sin #jp (r; #)j is bounded from above by 41. Theorem 3.4 The function family f g is a wavelet of order 0, and the func- tion family f g is a wavelet of order 1 according to De nition 2.2. Proof. The lower bound inequality in (1) is satis ed by Lemma 3.2 for l  1 in the case of the wavelet f g, respectively for l  2 in the case of the wavelet f g. In order to obtain an estimation from above, note that for odd k's Z Z 2l + 1 (ljkj)! jkj k ik' =  ! (#) P (cos #) sin # d# f (') e d' l l 4 (l +jkj)! 0 0 s s 2l + 1 (ljkj)! 3r (l +jkj)! 2 2 p  e 4 (l +jkj)! jkj(ljkj)! 2 1 r 3r 2l + 1 = e 2jkj(1 r ) and, analogously, 6r 2l + 1 e : 2jkj(1 r ) 12 Thus, since Z Z 1 2 2 1 4 2 e   e  6 d = and d = ; 2 2 1 e 24 1 e 24 0 0 we obtain for each l 2 2 k k X X 2 X 2 d 9 (2l + 1) e 2(2l + 1) e 2 2 2 24  jkj  jkj k=l k22Z+1 k22Z+1 and 2 2 k k X X 2 X 2 d 36( 6)(2l + 1) e 3(2l + 1) e 2 2 2 24  jkj  jkj k=l k22Z+1 k22Z+1 Further, since t 7! is a decreasing function for t 2 (0;1), the estimation 2 2 k Z t 2 2 e 1 e dt jkj 2 t k22N +1 is valid, and since Z t Z Z Z p 1 1 1 1 e e du 2 dt = du  + e du = ln  + ; 1 1 t u u 2 1 0 the upper bound inequality in (1) holds with p   p 2  3 B = ln  + and B = ln  + 2 2 2  2 for the families f g and f g, respectively. Condition (2) follows from the representations (5) and (7). 4 Multiselectivity scheme Our goal is to de ne an angular multiselectivity analysis similar to that pre- sented in [2, 1, 11]. The main di erence between our work and that of Antoine and Jacques is that in the R -case, an exact reconstruction is possible with the use of biorthogonal wavelets and scaling functions, whereas the spherical wavelet transform is invertible by frame methods. A multiselectivity analysis is de ned a priori for the discrete case. The con- tinuous wavelet transform is performed with respect to a xed wavelet, and it 13 is certainly possible to choose a function with a proper resolution power, but the choice must be based on some global estimations. In the case of a wavelet transform for a discrete set of scales and positions, the frame property remains unaltered if a wavelet is chosen for each of the parameters separately. The an- gular selectivity can be adapted to the features detected in each point of the analyzed signal by varying the parameter  . However, the discretization den- sity of the sphere S must be chosen in such a way that the frame property is preserved in the case that a wavelet with the highest angular resolution is used. More exactly, according to the proof of [10, Theorem 4.2], the error E made by the discretization of the set of wavelet coecients is bounded by jE j  4E + E ; 1 2 with E = ck k kr k   kfk ; J = 1; 2; J  1  1 j j where is the wavelet on the scale  , and  is an upper bound for the diam- j J eter of the discretization sets on the scale  , and f denotes the analyzed signal. 2 j1 2 Since in general jr j increases with increasing  , and jE j  2 kfk is j 2 required, an upper bound e for  must be set, and  computed with respect to . Suppose the discretization of scales and S -positions is performed. For each scale  and each position x , the 'best frame' for f , i.e. the best angular selectivity level is to be chosen (see [1, Section 5]). The best match between the image and the wavelet can be quanti ed by D E (g ); f ( ; x ) = arg max max ; g2SO(3) k k 2[1;e] j 1 2 g eb=x where eb denotes the north pole of the sphere. Alternatively, one can a priori choose a discrete (and bounded from above) set T of possible  's, and replace arg max 2[1;e] in the above formula by arg max: 2T Note that according to the proof of [10, Theorem 4.2], the discretization error E is an upper bound for a weighted sum of the discretization errors e ( ) made in each grid point 2 S that can be estimated by 1  2 e ( )  c supj j  supjr j kfk   ; j j 2 1 1 2 2 (absolute values are missing in [10, formula (12)]). The value of a single e ( ) j j does not change if one replaces by ,  =  ( ; x ), and  by  ( ), j 2 j  1 1 j 2 14 chosen appropriately. (Since supj j is independent of  , the diameter  ( ) depends only on supjr j. It is small for wavelets with a big angular resolu- tion, and it is big for wavelets with a small angular resolution.) Consequently, the discrete set of wavelets constructed according to the above description is a 2 2 frame for the whole L (S ) (the frame condition is satis ed independently of f ), but its choice is adapted to a concrete function. References [1] J.-P. Antoine and L. Jacques, Angular multiselectivity analysis of images, in wavelet applications in signal and image processing X, Proc. SPIE 5207 (2003), 196{207. [2] J.-P. Antoine and L. Jacques, The 2-D wavelet transform in image process- ing: two novel applications, in J. Govaerts, M. N. Hounkonnou, and A. Z. Msezane (eds.), Contemporary problems in mathematical physics, Proceed- ings of the Third International Workshop, Cotonou, Republic of Benin, 1 { 7 Nov. 2003, World Scienti c, Singapore 2004, pp. 414{435. [3] J.-P. Antoine, R. Murenzi, P. Vandergheynst, and S. T. Ali, Two- dimensional wavelets and their relatives, Cambridge University Press, Cam- bridge 2008. [4] S. Bernstein, Spherical singular integrals, monogenic kernels and wavelets on the three{dimensional sphere, Adv. Appl. Cli ord Algebr. 19 (2009), no. 2, 173{189. [5] S. Ebert, S. Bernstein, P. Cerejeiras, and U. K ahler, Nonzonal wavelets N th on S , 18 International Conference on the Application of Computer Sci- ence and Mathematics in Architecture and Civil Engineering, Weimar 2009. [6] I.S. Gradshteyn and I.M. Ryzhik, Table of integrals, series, and products, Elsevier/Academic Press, Amsterdam, 2007. [7] M. Hayn and M. Holschneider, Directional spherical multipole wavelets, J. Math. Phys. 50 (2009), no. 7, 073512, 11 pp. [8] I. Iglewska{Nowak, Continuous wavelet transforms on n{dimensional spheres, Appl. Comput. Harmon. Anal. 39 (2015), no. 2, 248{276. [9] I. Iglewska{Nowak, Directional wavelets on n-dimensional spheres, Appl. Comput. Harmon. Anal. (2016), dx.doi.org/10/1016/j.acha.2016.04.008. [10] I. Iglewska{Nowak, Frames of directional wavelets on n{dimensional spheres, Appl. Comput. Harmon. Anal. 43 (2017), 148{161. [11] L. Jacques and J.-P. Antoine, Multiselective pyramidal decomposition of images: wavelets with adaptive angular selectivity, Int. J. Wavelets Mul- tiresolut. Inf. Process. 5 (2007), no. 5, 785{814. http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png Mathematical Physics arXiv (Cornell University)

Angular multiselectivity with spherical wavelets

Iglewska-Nowak, Ilona
Mathematical Physics , Volume 2018 (1804) – Apr 6, 2018
Download PDF
14 pages

Loading next page...
 
/lp/arxiv-cornell-university/angular-multiselectivity-with-spherical-wavelets-SVHG0lwojO

References (12)

ISSN
1063-5203
eISSN
ARCH-3337
DOI
10.1016/j.acha.2017.11.002
Publisher site
See Article on Publisher Site

Abstract

We construct spherical wavelets based on approximate identities that are directional, i.e. not rotation-invariant, and have an adaptive angular selectivity. The problem of how to nd a proper representation of distinct kinds of details of real images, ranging from highly directional to fully isotropic ones, was quite intensively studied for the case of signals over the Euclidean space. However, the present paper is the rst attempt to deal with this task in the case of spherical signals. A multiselectivity scheme, similar to that proposed for R -functions, is presented. Key words and phrases: spherical wavelets, Poisson kernel, angular selectiv- ity, directional wavelet frames, angular multiresolution 2010 Mathematics Subject Classi cation Number: 42C40, 42C15, 94A12, 65T60 1 Introduction Recently, directional wavelets based on approximate identities have been intro- duced [4, 5, 7, 8, 9]. Unfortunately, wavelet families constructed so far have a very low angular resolution, e.g., the rst directional derivative of the Poisson kernel is a wavelet with separable spherical variables (#; ') of the form [1] [1] (x) = (#) cos '; where x = (cos #; sin # cos '; sin # sin ') 2 S; # 2 [0; ]; ' 2 [0; 2): (Unless it leads to misunderstandings, we identify x 2 S with its spherical coordinates (#; ').) It is depicted on the left-hand side of Figure 1 and it can be seen that its ability to detect features with close orientations is not satisfactory. In [7] second order wavelets are used to analyse the land cover data of the world, where the analyzed objects are in principle curves. The wavelets are of the form [2] [2] [2] 2 (x) = (#) + (#) cos '; ;1 ;2 West Pomeranian University of Technology in Szczecin, School of Mathematics, al. Pias- t ow 17, 70{310 Szczecin, Poland arXiv:1804.03044v1 [math.CA] 6 Apr 2018 2 Figure 1: The rst directional derivative of the Poisson kernel on S vs. a wavelet with a bigger angular selectivity { behavior around the north pole see [9, Appendix A], i.e. again functions with a low angular resolution. The aim of the present paper is to construct wavelets with a variable angular selectivity such that their ability to detect directional features can be adapted to the ana- lyzed signal. An example of such a wavelet is presented on the right-hand-side of Figure 1. It is a function with separable variables and Figure 3 shows the dependence of the function of the longitudinal variable on a parameter that can be chosen almost freely. The investigation is inspired by the results obtained by Antoine and Jacques in [2, 1, 11] for the case of wavelets over R . The angular selectivity of spherical wavelets was analysed by [3], but in this case wavelets based on the group theory were involved (the di erence between the two constructions is discussed in [8]). The angular resolving power was de ned in [3, Subsection 9.2.4.2]. Nonetheless, we do not intend to adapt it to our purposes, since we nd the parameter introduced in Section 3 to be a better indicator of the angular selectivity in the case of wavelets constructed in the present paper. Moreover, one that can be a priori chosen. The paper is organized as follows. After an introduction of necessary notions in Section 2, we present a construction of two wavelets with a steerable angular selectivity in Section 3. An angular multiselectivity analysis is then discussed in Section 4. 2 Preliminaries 3 3 2 By S we denote the unit two-dimensional sphere in R , S = f 2 R : jj = 1g, with the measure d invariant under the rotation group and such that d = 4. # 2 [0; ] is the polar (colatitudinal) coordinate, and ' 2 [0; 2) { the azimuthal (longitudinal) coordinate. The scalar product of spherical functions , is de ned by h; i = (x) (x) d(x); S 3 and L (S ) is the set of functions such that k k 2 := h ; i < 1: L (S) Note the di erence to the notation usually used in research concerning n- dimensional spheres, where the integral in the de nition of the scalar product is divided by  , the Lebesgue measure of S , see e.g. [8]. The spherical harmonics Y , l 2 N , k = l;l + 1; : : : ; l, are given by 2l + 1 (l k)! k k k ik' Y (#; ') = (1) P (cos #) e ; l l 4 (l + k)! Y = Y l l for 0  k  l, where P are the (associated) Legendre polynomials (functions), 1 d 0 2 l P (t) = P (t) = (t 1) (Legendre polynomials) l l 2 l! dt k k 2 k=2 P (t) = (1) (1 t ) P (t) for 0 < k  l (associated Legendre functions). dt 2 2 The spherical harmonics build an orthonormal basis for L (S ). For an L (S ){ function , the series 1 l X X k k Y ; Y l l l=1 k=l is the Fourier series of in terms of the spherical harmonics. The constants k k = Y ; are called the Fourier coecients of . l l T denotes the one-dimensional torus T = [0; 2). Functions on T will be iden- ti ed with their 2-periodic extension on R, such that C(T) denotes the class of continuous functions on T that can be continuously extended. An L (T)- function can be represented as the Fourier series ikt f (t)  f e ; k=1 where the Fourier coecients f , k 2 Z, are given by ikt f = f (t) e dt: Similarly, the Fourier transform of an L (R) function F is de ned by i!t F (!) = F (t) e dt;  2 R; and it can be inverted by i!t F (t) = F (!) e d: 1 4 The functions we consider in the present paper are both integrable and con- tinuous such that integrals and series in the Fourier and the inverse Fourier transform are (pointwise) convergent. An L (R)-function F can be periodized to an L (T)-function f by f (t) = f (t + 2j); t 2 T; j=1 and their Fourier transforms are linked by the following Poisson summation formula. Theorem 2.1 Let F 2 L (R). Then its periodization f exists for almost every t 2 T, it is integrable and kfk 1  kFk 1 . Further, L (T) L (R) f = F (k): A family of vectors fg ;  2 Ig  H in a Hilbert space H indexed by a measure space I with a positive measure  is called a frame with weight  if the mapping 7! g is weakly measurable, i.e.  7! hg ; ui is measurable and 2 2 2 Akuk  jhg ; uij d(x)  Bkuk for all u 2 H and 0 < A  B. If fg ;  2 Ig is a frame, then the mapping u 7! fhg ; ui ;  2 Ig is invertible. The wavelet de nition we use originates from [10]. In the case of the two- dimensional sphere and the weight function  7! , it is the following one. 2 2 De nition 2.2 The family f g  L (S ) is called a wavelet (family) of 2R order m if it satis es A(2l + 1)  ( )  B(2l + 1) (1) k=l for some positive constants A and B independent of l 2 N , l > m, and ( ) = 0 (2) k=l for l 2 N , l  m. 2 2 The wavelet transform of an L (S )-function is given by W f (; g) = (g x) f (x) d (x); g 2 SO(3); 4 n S 5 Figure 2: Grid construction and it is invertible by methods known from the frame theory. Moreover, only a discrete set of the wavelet coecients W f (; g) is necessary to reconstruct the analyzed signal. The rst step is a scale discretization. According to [10, Theorem 3.1] there exist constants a and X such that for any sequence R = ( ) with   a and 1 <  = < X the family f (g );  2 R; g 2 j j2N 0 0 j j+1  j 0 j 2 2 SO(3)g is a frame for L (S ). In order to characterize discrete sets of rotations, we need the following de nition (compare [10, De nition 4.1]). De nition 2.3  is a grid of type ( ;  ) if it is a discrete measurable set of 2 1 2 1 rotations in SO(3) = S  S , constructed in the following way. Let P = 2 2 fO : = 1; : : : ; K g be a measurable partition of S into simply connected 2 2 sets such that the diameter of each set (measured in the geodesic distance) is 2 2 not larger than  . Choose from each of the sets O an arbitrary point x = 2 2 2 2 1 (# ; ' ). Now, for a fxed , let P ( ) = fO : = 1; : : : ; K ( )g 2 1 2 1 1 2 2 2 ( ; ) 2 1 be a measurable partition of S into K ( ) simply connected sets of a diameter 1 2 not larger than  . Choose from each of the sets O , = 1; : : : ; K ( ), 1 1 1 2 ( ; ) 2 1 an arbitrary point ' . Then,  is the set of rotations given by ( ; ) 2 1 1 2 2 g = g (' ) g (# ) g (' ); ( ; ) 1 1 2 2 1 ( ; ) 2 1 2 2 where g ( ) is the rotation in the plane ( ;  ) with the rotation angle . j j j+1 Loosely speaking, one has a tight enough discrete set of points fx g 2 S , and for each point a tight enough discrete set of rotations around it.  is a set of SO(3)-rotations being a composition of a tranlastion to a point x and a rotation around this point, as illustrated in Figure 2. The next statement re ects the content of [10, Theorem 4.2] for the case of the two-dimensional sphere. 6 Figure 3: f for  = 1; 2; 4; 8; 16, ' 2 [=2; 3=2] Theorem 2.4 Let be a C {wavelet family with the property that f ; j 2 ;x N ; x 2 S g is a semi{continuous frame. Then, for each j 2 N there exist 0 0 j j numbers  ,  such that 2 1 1 j f (g ); j 2 N ; g j j 2  g j j 0 j ( ; ) ( ; ) 2 1 2 1 j j 2 n j is a frame for L (S ), provided that  is a grid of type ( ;  ). 2 1 3 Construction of the wavelet In this section, we shall present a wavelet family with separated variables and a steerable angular selectivity. The construction is based on the Poisson kernel, 1 1 r 1 p (x) = = (2l + 1)r P (cos #); (3) 2 3=2 4 (1 2r cos # + r ) 4 l=0 where r = e , since our aim is to obtain a wavelet with a simple representation as a function of the spherical variables. On the other hand, in order to prove the wavelet property, some knowledge about the Fourier coecients is necessary. The idea is to choose a periodized di erence of Gaussians dependent on a pa- rameter  2 [1;1) as the function of the longitudinal variable ', 2 2 2 2 () 2 2 F () = e e ;  2 R; f (') = F (' + 2j); ' 2 [0; 2): j2Z Figure 3 presents the function f for several values of  . The challenge is to nd an appropriate function of the colatitudinal variable # such that the inequalities (1) are satis ed. The result of our investigation is that the following families: @ @ (x) = (#; ') = ! (#) f ('); ! (#) =  sin # r  r  [p (x)] ; @r @r 5 2 (x) =  (#; ') =  (#) f (');  (#) =  sin # r  [p (x)] ; @r 7 r = e , are wavelets. In order to prove the wavelet condition (1), we need an estimation of the Fourier coecients (with respect to the spherical harmonics) of and , and since they are functions with separated variables, we rst consider the Fourier series expansion (in exponential functions) of f . Since 2 2 ! ! i! b 2 2 2 2 F (!) = e e ; by the Poisson summation formula the Fourier coecients of f are given by 0; for even k; p 2 (f ) = k k 2 2 e for odd k: Thus, X 2 1 2 ik' f (') = e e : (4) k22Z+1 In order to estimate the Fourier coecients of and , we estimate the mag- nitude of the coecients of sin # 2 l ! (#) = (2l + 1)l r P (cos #) (5) l=1 2 4 2 4 2 2 5 r r(10 19r + r ) (3 14r + 5r ) cos # r(9 r ) cos # sin # = ; 2 7=2 4(1 2r cos # + r ) (6) and sin # (#) = (2l + 1)l(l 1)r P (cos #) (7) l=2 2 2 4 2 2 2 5 r 5 23r + 2r + 4r(7 + r ) cos # (15 + r ) cos # sin # = ; (8) 2 7=2 4(1 2r cos # + r ) r = e , developed into series with respect to the associated Legendre polyno- mials of an odd order k. Lemma 3.1 For r 2 (0; 1), ! (#) = P (cos #) and l=1 (#) = P (cos #) l=2 8 with 3 l+5 3 2 l+3 (l + 2)(l + 3)(l + 4)(l + 5) r (l + 2)(l + 3) (5l + 21l 8) r (2l + 3)(2l + 5)(2l + 7)(2l + 9) (2l 1)(2l + 3)(2l + 5)(2l + 9) 2 5 4 3 2 l+1 2(l + 1) (10l + 65l + 83l 134l 186l + 36) r (2l 3)(2l 1)(2l + 3)(2l + 5)(2l + 7) 2 4 3 2 l1 2(l 1) (5l 2l 41l + 14l + 60) r (9) (2l 5)(2l 3)(2l + 3)(2l + 5) 2 4 3 2 l3 (l 3) (l 2)(l 1)(20l 100l + 39l + 239l 48) r (2l 7)(2l 5)(2l 3)(2l 1)(1 + 2l)(2l + 3) 2 l5 (l 5) (l 4)(l 3)(l 2)(l 1) r (2l 7)(2l 5)(2l 3)(2l 1) and 2 2 l+5 2 2 2 l+3 (l + 2)(l + 3)(l + 4) (l + 5) r (l + 2) (l + 3) (5l + 21l 8) r (2l + 3)(2l + 5)(2l + 7)(2l + 9) (2l 1)(2l + 3)(2l + 5)(2l + 9) 5 4 3 2 l+1 2l(l + 1)(10l + 65l + 83l 134l 186l + 36) r (2l 3)(2l 1)(2l + 3)(2l + 5)(2l + 7) 4 3 2 l1 2(l 2)(l 1)(5l 2l 41l + 14l + 60) r (2l 5)(2l 3)(2l + 3)(2l + 5) 4 3 2 l3 (l 4)(l 3)(l 2)(l 1)(20l 100l + 39l + 239l 48) r (2l 7)(2l 5)(2l 3)(2l 1)(1 + 2l)(2l + 3) l5 (l 6)(l 5)(l 4)(l 3)(l 2)(l 1) r (2l 7)(2l 5)(2l 3)(2l 1) (10) for l  9. The coecients , l = 1; 2; : : : ; 8, and , l = 2; 3; : : : ; 8, do not l l vanish. Proof. By the repeated application of [6, formulae (8.733.4) and (8.735.5)]: k k P (cos #) P (cos #) k1 l1 l+1 sin #P (cos #) = ; 2l + 1 (l k)(l k + 1) (l + k)(l + k + 1) k+1 k k sin #P (cos #) = P (cos #) P (cos #) l l+1 l1 2l + 1 2l + 1 9 we obtain (l 3)(l 2)(l 1)l 5 0 1 (2l + 1) sin #P (cos #) = P (cos #) l l5 (2l 7)(2l 5)(2l 3)(2l 1) (l 1)l(5l 9l 26) P (cos #) l3 (2l 7)(2l 3)(2l 1)(2l + 3) 5 4 3 2 2(10l + 15l 77l 93l + 121l + 60) + P (cos #) l1 (2l 5)(2l 3)(2l + 1)(2l + 3)(2l + 5) 4 3 2 2(5l + 18l 17l 54l + 36) P (cos #) l+1 (2l 3)(2l 1)(2l + 5)(2l + 7) 4 3 2 (l + 1)(l + 2)(20l + 140l + 219l 67l 60) + P (cos #) l+3 (2l 1)(2l + 1)(2l + 3)(2l + 5)(2l + 7)(2l + 9) (l + 1)(l + 2)(l + 3)(l + 4) P (cos #) l+5 (2l + 3)(2l + 5)(2l + 7)(2l + 9) for l  4. A substitution of this formula to (5) resp. (7) and index shifts yield coecients of ! resp.  as in (9) resp. (10) for l  9. The coecients , , l l l = 1; 2; : : : ; 8, are computed in the same manner. The di erence with respect to the case l  9 is caused by vanishing polynomials P . 1 1 b b Lemma 3.2 The Fourier coecients , l = 1; 2; : : : , and  , l = 2; 3; : : : , l l satisfy Z Z 1 1 2 2 d d 1 1 b b (2l + 1)A ; (2l + 1)A   (11) l l 0 0 for some constants A , A independent of l. Proof. It follows from Lemma 3.1 and formula (4) that X X 1 2 k 1 ik' (x) = ! (#) f (') = P (cos #) e e l=1 k22Z+1 X 2 2 k 1 ik' = e P (cos #) e : l;k Since 4  l(l + 1) 1 i' 1 P (cos #) e = Y (#; ') l l 2l + 1 and the spherical harmonics are orthogonal to one another, r r 2 1 4  l(l + 1)  l(l + 1) 1 b 2 2 2 2 =  e  =  e : l l 4  2l + 1  2(2l + 1) 10 Similarly, l(l + 1) 1 b 2 =  e : 2(2l + 1) For l  9, the integrals Z Z 1 1 2 2 1 d 1 d 1 1 b b and l l 2l + 1  2l + 1 0 0 are positive rational functions of l. It can be veri ed that Z Z 1 1 2 2 d and  d l l 0 0 are continuous for l  9 and that their limits for l ! 1 both equal . Consequently, Z Z 1 1 2 2 1 d 1 d 9 1 9 1 b b A := min and A := min l  l l2[9;1) 2l + 1  l2[9;1) 2l + 1 0 0 both exist and are positive. Similarly, for l = 1; 2; : : : ; 8 the numbers 1 d l 1 A := 2l + 1 and for l = 2; 3; : : : ; 8 the numbers 1 d l 1 A := 2l + 1 1 2 8 9 are positive. Inequalities (11) are thus satis ed with A = minfA ; A ; : : : ; A ; A g 2 3 8 9 and A = minfA ; A ; : : : ; A ; A g. Lemma 3.3 Let k 2 N. Then 3r (l + k)! ! (#)P (cos #) sin # d#  p and (12) k(l k)! 2 1 r 3r (l + k)! (#)P (cos #) sin # d#  p (13) k(l k)! 1 r k k Proof. Denote the left-hand-sides of (12) and (13) by !e and e , respectively. l l They are bounded from above by 6 k r sin #jp (r; #)j P (cos #) k l !e  d# and 2 7=2 4 (1 2r cos # + r ) 2  k sin #jp (r; #)j P (cos #) e  d#; 2 7=2 4 (1 2r cos # + r ) 0 11 where 2 4 2 4 2 2 p (r; #) = r(10 19r + r ) (3 14r + 5r ) cos # r(9 r ) cos # and 2 4 2 2 2 p (r; #) = 5 23r + 2r + 4r(7 + r ) cos # (15 + r ) cos #; see (6) and (8). The H older inequality together with k 2 [P (cos #)] (l + k)! d# = sin # k(l k)! (see [6, formula (7.122.1)], note that 2 in the denominator disappears because of the extended integration interval) yield !e  sup sin #jp (r; #)j r2(0;1) #2[0;] " # 1=2 Z  Z 2 1=2 12  k P (cos #) sin # d#  d# 2 7 (1 2r cos # + r ) sin # 0 0 1=2 r 231 (l + k)! 3r (l + k)! 21   p : 4 1024(1 r ) k(l k)! k(l k)! 2 (1 r ) Estimation (13) is obtained in a similar way. It is to be respected that the function (r; #) 7! sin #jp (r; #)j is bounded from above by 41. Theorem 3.4 The function family f g is a wavelet of order 0, and the func- tion family f g is a wavelet of order 1 according to De nition 2.2. Proof. The lower bound inequality in (1) is satis ed by Lemma 3.2 for l  1 in the case of the wavelet f g, respectively for l  2 in the case of the wavelet f g. In order to obtain an estimation from above, note that for odd k's Z Z 2l + 1 (ljkj)! jkj k ik' =  ! (#) P (cos #) sin # d# f (') e d' l l 4 (l +jkj)! 0 0 s s 2l + 1 (ljkj)! 3r (l +jkj)! 2 2 p  e 4 (l +jkj)! jkj(ljkj)! 2 1 r 3r 2l + 1 = e 2jkj(1 r ) and, analogously, 6r 2l + 1 e : 2jkj(1 r ) 12 Thus, since Z Z 1 2 2 1 4 2 e   e  6 d = and d = ; 2 2 1 e 24 1 e 24 0 0 we obtain for each l 2 2 k k X X 2 X 2 d 9 (2l + 1) e 2(2l + 1) e 2 2 2 24  jkj  jkj k=l k22Z+1 k22Z+1 and 2 2 k k X X 2 X 2 d 36( 6)(2l + 1) e 3(2l + 1) e 2 2 2 24  jkj  jkj k=l k22Z+1 k22Z+1 Further, since t 7! is a decreasing function for t 2 (0;1), the estimation 2 2 k Z t 2 2 e 1 e dt jkj 2 t k22N +1 is valid, and since Z t Z Z Z p 1 1 1 1 e e du 2 dt = du  + e du = ln  + ; 1 1 t u u 2 1 0 the upper bound inequality in (1) holds with p   p 2  3 B = ln  + and B = ln  + 2 2 2  2 for the families f g and f g, respectively. Condition (2) follows from the representations (5) and (7). 4 Multiselectivity scheme Our goal is to de ne an angular multiselectivity analysis similar to that pre- sented in [2, 1, 11]. The main di erence between our work and that of Antoine and Jacques is that in the R -case, an exact reconstruction is possible with the use of biorthogonal wavelets and scaling functions, whereas the spherical wavelet transform is invertible by frame methods. A multiselectivity analysis is de ned a priori for the discrete case. The con- tinuous wavelet transform is performed with respect to a xed wavelet, and it 13 is certainly possible to choose a function with a proper resolution power, but the choice must be based on some global estimations. In the case of a wavelet transform for a discrete set of scales and positions, the frame property remains unaltered if a wavelet is chosen for each of the parameters separately. The an- gular selectivity can be adapted to the features detected in each point of the analyzed signal by varying the parameter  . However, the discretization den- sity of the sphere S must be chosen in such a way that the frame property is preserved in the case that a wavelet with the highest angular resolution is used. More exactly, according to the proof of [10, Theorem 4.2], the error E made by the discretization of the set of wavelet coecients is bounded by jE j  4E + E ; 1 2 with E = ck k kr k   kfk ; J = 1; 2; J  1  1 j j where is the wavelet on the scale  , and  is an upper bound for the diam- j J eter of the discretization sets on the scale  , and f denotes the analyzed signal. 2 j1 2 Since in general jr j increases with increasing  , and jE j  2 kfk is j 2 required, an upper bound e for  must be set, and  computed with respect to . Suppose the discretization of scales and S -positions is performed. For each scale  and each position x , the 'best frame' for f , i.e. the best angular selectivity level is to be chosen (see [1, Section 5]). The best match between the image and the wavelet can be quanti ed by D E (g ); f ( ; x ) = arg max max ; g2SO(3) k k 2[1;e] j 1 2 g eb=x where eb denotes the north pole of the sphere. Alternatively, one can a priori choose a discrete (and bounded from above) set T of possible  's, and replace arg max 2[1;e] in the above formula by arg max: 2T Note that according to the proof of [10, Theorem 4.2], the discretization error E is an upper bound for a weighted sum of the discretization errors e ( ) made in each grid point 2 S that can be estimated by 1  2 e ( )  c supj j  supjr j kfk   ; j j 2 1 1 2 2 (absolute values are missing in [10, formula (12)]). The value of a single e ( ) j j does not change if one replaces by ,  =  ( ; x ), and  by  ( ), j 2 j  1 1 j 2 14 chosen appropriately. (Since supj j is independent of  , the diameter  ( ) depends only on supjr j. It is small for wavelets with a big angular resolu- tion, and it is big for wavelets with a small angular resolution.) Consequently, the discrete set of wavelets constructed according to the above description is a 2 2 frame for the whole L (S ) (the frame condition is satis ed independently of f ), but its choice is adapted to a concrete function. References [1] J.-P. Antoine and L. Jacques, Angular multiselectivity analysis of images, in wavelet applications in signal and image processing X, Proc. SPIE 5207 (2003), 196{207. [2] J.-P. Antoine and L. Jacques, The 2-D wavelet transform in image process- ing: two novel applications, in J. Govaerts, M. N. Hounkonnou, and A. Z. Msezane (eds.), Contemporary problems in mathematical physics, Proceed- ings of the Third International Workshop, Cotonou, Republic of Benin, 1 { 7 Nov. 2003, World Scienti c, Singapore 2004, pp. 414{435. [3] J.-P. Antoine, R. Murenzi, P. Vandergheynst, and S. T. Ali, Two- dimensional wavelets and their relatives, Cambridge University Press, Cam- bridge 2008. [4] S. Bernstein, Spherical singular integrals, monogenic kernels and wavelets on the three{dimensional sphere, Adv. Appl. Cli ord Algebr. 19 (2009), no. 2, 173{189. [5] S. Ebert, S. Bernstein, P. Cerejeiras, and U. K ahler, Nonzonal wavelets N th on S , 18 International Conference on the Application of Computer Sci- ence and Mathematics in Architecture and Civil Engineering, Weimar 2009. [6] I.S. Gradshteyn and I.M. Ryzhik, Table of integrals, series, and products, Elsevier/Academic Press, Amsterdam, 2007. [7] M. Hayn and M. Holschneider, Directional spherical multipole wavelets, J. Math. Phys. 50 (2009), no. 7, 073512, 11 pp. [8] I. Iglewska{Nowak, Continuous wavelet transforms on n{dimensional spheres, Appl. Comput. Harmon. Anal. 39 (2015), no. 2, 248{276. [9] I. Iglewska{Nowak, Directional wavelets on n-dimensional spheres, Appl. Comput. Harmon. Anal. (2016), dx.doi.org/10/1016/j.acha.2016.04.008. [10] I. Iglewska{Nowak, Frames of directional wavelets on n{dimensional spheres, Appl. Comput. Harmon. Anal. 43 (2017), 148{161. [11] L. Jacques and J.-P. Antoine, Multiselective pyramidal decomposition of images: wavelets with adaptive angular selectivity, Int. J. Wavelets Mul- tiresolut. Inf. Process. 5 (2007), no. 5, 785{814.

Journal

Mathematical Physics – arXiv (Cornell University)

Published: Apr 6, 2018

There are no references for this article.