Skip to main content

A class of generalized pseudo-splines

Abstract

In this paper, a class of refinable functions is given by smoothening pseudo-splines in order to get divergence free and curl free wavelets. The regularity and stability of them are discussed. Based on that, the corresponding Riesz wavelets are constructed.

1 Introduction

We denote by and the set of integers and real numbers, respectively. Let L p (R) stand for the classical Lebesgue space

L p (R):= { f , R | f ( t ) | p d t < + }

with the norm f p = ( R | f ( t ) | p d t ) 1 p and L (R) consisting of all Lebesgue measurable and bounded functions on . Similarly, the discrete space p (Z):={{ a n }, n | a n | p <+,nZ} with { a n } p = ( n | a n | p ) 1 p . As usual, given f L 1 (R) L 2 (R), its Fourier transform is defined by

f ˆ (ω):= R f(x) e i x ω dx

on . The Fourier transform of a function in L 2 (R) is understood as the unitary extension. We write h=fg for the convolution h(x)= R f(xt)g(t)dt, defined for any pair of functions f and g such that the integral exists almost everywhere. Clearly, h ˆ (ω)= f ˆ (ω) g ˆ (ω) in the frequency domain, when all the Fourier transforms exist in that formula. Given g L 2 (R), {g(xk),kZ} is called a Riesz basis of its linearly generating space, if for each { λ k } 2 (Z) there exist two positive constants A and B such that

A k Z λ k 2 k Z λ k g ( x k ) 2 2 B k Z λ k 2 .
(1.1)

The numbers A, B are called lower Riesz bound and upper Riesz bound, respectively.

Multiresolution analysis provides a classical method to construct wavelets.

Definition 1 A multiresolution analysis of L 2 (R) means a sequence of closed linear subspaces V j of L 2 (R) which satisfies

  1. (i)

    V j V j + 1 , jZ,

  2. (ii)

    f(x) V j if and only if f(2x) V j + 1 ,

  3. (iii)

    j Z V j ¯ = L 2 (R) and j Z V j ={0},

  4. (iv)

    there exists a function ϕ L 2 (R) such that {ϕ(xk),kZ} forms a Riesz basis of  V 0 .

The function ϕ in Definition 1 is said to be a scaling function, if it satisfies

ϕ(x)= k a k ϕ(2xk)
(1.2)

for some sequence { a k } 2 (Z). Define the Fourier series c ˆ of a sequence { c k } 2 (Z) by

c ˆ (ξ):= k Z c k e i k ξ ,ξR.

Then the refinement equation (1.2) becomes

ϕ ˆ (ξ)= a ˆ (ξ/2) ϕ ˆ (ξ/2),ξR.

The function a ˆ is called the refinement mask of ϕ. The pseudo-spline of Type I was first introduced in [1] to construct tight framelets. The pseudo-spline of Type II was first studied by Dong and Shen in [2]. There have been many developments in the theory of pseudo-splines over the past ten years [3, 4]. Its applications in image denoising and image in-painting are also very extensive [5, 6]. The pseudo-spline is defined by its refinement mask. The refinement mask of a pseudo-spline of Type I with order (m,) is given by

| 1 a ˆ (ξ) | 2 :=| 1 a ˆ m , (ξ) | 2 := cos 2 m (ξ/2) j = 0 ( m + j ) sin 2 j (ξ/2) cos 2 ( j ) (ξ/2)
(1.3)

and the refinement of a pseudo-spline of Type II with order (m,) is given by

a ˆ 2 (ξ):= 2 a ˆ m , (ξ):= cos 2 m (ξ/2) j = 0 ( m + j ) sin 2 j (ξ/2) cos 2 ( j ) (ξ/2).
(1.4)

The mask of Type I is obtained by taking the square root of the mask of Type II using the Fejér-Riesz lemma [7], i.e. a ˆ 2 (ξ)= | 1 a ˆ ( ξ ) | 2 . The corresponding pseudo-spline can be defined in terms of their Fourier transform, i.e.

ϕ ˆ m , k (ξ):= j = 1 k a ˆ m , ( 2 j ξ ) ,k=1,2.

In order to smoothen the pseudo-spline, one can use the convolution method. Take the smoothed pseudo-spline

where χ [ 1 2 , 1 2 ] denotes the characteristic function of interval [ 1 2 , 1 2 ] and nm. This is equivalent to

ϕ ˆ n , m , (ξ)= ϕ ˆ m , (ξ) ( sin ( ξ / 2 ) ξ / 2 ) n m .

Thus the symbol of ϕ n , m , becomes

a ˆ n , m , (ξ)= ϕ ˆ n , m , (2ξ)/ ϕ ˆ n , m , (ξ)= a ˆ m , (ξ) ( cos ( ξ / 2 ) ) n m .

Therefore, we define the smoothed pseudo-spline by its refinement mask for Type I:

| 1 a ˆ n , m , | 2 := cos 2 n (ξ/2) j = 0 ( m + j ) sin 2 j (ξ/2) cos 2 ( j ) (ξ/2)
(1.5)

and for Type II:

ϕ ˆ r , m , 2 := cos r (ξ/2) j = 0 ( m + j ) sin 2 j (ξ/2) cos 2 ( j ) (ξ/2),

where r2m. When r=2m, it is the pseudo-spline. When r2m, it can be considered as an extension of pseudo-spline. Define the translated form of the Type II by

ϕ ˆ r , m , T (ξ):= e i r ξ 2 2 ϕ ˆ r , m , (ξ).

Then we get the differential relation

ϕ r + 1 , m , T (x)= T ϕ r , m , (x) T ϕ r , m , (x1).
(1.6)

This inherits the property of a B-spline.

Remark 1 One may think that smoothing the pseudo-splines by convolving them with B-splines seems unnecessary since one can simply increase m of the original pseudo-splines. However, by increasing m, we cannot get the differential relation (1.6), which is important for the construction of divergence free wavelets and curl free wavelets in the analysis of incompressible turbulent flows [8, 9].

Remark 2 Similar to the definition of (1.5), we can define a smoothed dual pseudo-spline by its refinement mask,

b ˆ n , m , (ξ)= e i ξ / 2 cos 2 n + 1 (ξ/2) j = 0 ( m + 1 / 2 + j ) sin 2 j (ξ/2) cos 2 ( j ) (ξ/2),
(1.7)

as an extension of dual pseudo-splines in [3] and get the corresponding wavelets.

Remark 3 In addition, one can assume nR in (1.5) and (1.7), as a generalization of fractional splines in [10].

2 Some lemmas

This section gives some lemmas that will be used to prove several results of this paper. We start with some results from [2].

Lemma 1 [2]

For given nonnegative integers m, j, ,

2 ( m + ) 1 2 j = 0 ( m + j ) for all m1 and 0m1.

This lemma will be used in Section 4 in order to prove the Riesz basis property of wavelets. Define P m , (y):= j = 0 ( m + j ) y j ( 1 y ) j , R m , (y):= ( 1 y ) m P m , (y) and R r , m , = ( 1 y ) r 2 P m , (y) where y= sin 2 (ξ/2), r, m, are nonnegative integers and r2m. Then one can find that

R m , ( sin 2 ( ξ / 2 ) ) = 2 a ˆ m , (ξ)and R r , m , ( sin 2 ( ξ / 2 ) ) = 2 a ˆ r , m , (ξ)

and the following lemma holds.

Lemma 2 [2]

For nonnegative integers m and with m1, let P m , and R m , be the polynomials defined above. Then

  1. (1)

    P m , (y)= j = 0 ( m 1 + j j ) y j ;

  2. (2)

    R m , (y)=(m+) ( m + 1 ) y ( 1 y ) m 1 .

With the two lemmas in hand, the basic property of the polynomial R r , m , , which will be used in Section 4, is given.

Lemma 3 For nonnegative integers r, m and ,

  1. (1)

    define Q(y):= R r , m , (y)+ R r , m , (1y); then

    min y [ 0 , 1 ] Q(y)=Q ( 1 2 ) = 2 1 r 2 j = 0 ( m + j ) ;
  2. (2)

    define S(y):= R r , m , 2 (y)+ R r , m , 2 (1y); then

    min y [ 0 , 1 ] S(y)=S ( 1 2 ) = 2 1 r 2 ( j = 0 ( m + j ) ) 2 .

Proof Since R r , m , (y)= ( 1 y ) r 2 m R m , (y), its derivative is

R r , m , (y)= ( r 2 m ) ( 1 y ) r 2 m 1 R m , (y)+ ( 1 y ) r 2 m R m , (y).

So the derivative is Q (y)= R r , m , (y)+ R r , m , (1y)=I+II, where

I= ( r 2 m ) y r 2 m 1 R m , (1y) ( r 2 m ) ( 1 y ) r 2 m 1 R m , (y)

and

II= ( 1 y ) r 2 m R m , (y) y r 2 m R m , (1y).

Now, we compute them, respectively. For I, by using (1) of Lemma 2, one has

I= ( r 2 m ) j = 0 ( m 1 + j j ) [ y r 2 1 ( 1 y ) j y j ( 1 y ) r 2 1 ] .

For II, by using (2) of Lemma 2, one has

II=(m+) ( m + 1 ) [ y r 2 1 ( 1 y ) ( 1 y ) r 2 1 y ] .

For j=0,,, since y r 2 1 ( 1 y ) j ( 1 y ) r 2 1 y j for all y[0, 1 2 ] and y r 2 1 ( 1 y ) j ( 1 y ) r 2 1 y j for all y[ 1 2 ,1], one has

Q (y)=I+II { 0 , y [ 0 , 1 2 ] ; 0 , y [ 1 2 , 1 ] .

This means Q(y) reaches its minimum value at the point y=1/2. Furthermore,

Q(1/2)=2 R r , m , (1/2)= 2 1 r 2 P m , (1/2)= 2 1 r 2 j = 0 ( m + j ) .

This completes the proof of (1). For (2) of this lemma, since

S(y)= R r , m , 2 (y)+ R r , m , 2 (1y)= ( 1 y ) r 2 m R m , 2 (y)+ y r 2 m R m , 2 (1y),

we have S (y)=III+IV, where

III=(r2m) [ y r 2 m 1 R m , 2 ( 1 y ) ( 1 y ) r 2 m 1 R m , 2 ( y ) ]

and

IV=2 ( 1 y ) r 2 m R m , (y) R m , (y)2 y r 2 m R m , (1y) R m , (1y).

For III, by (1) of Lemma 2, we have

III = ( r 2 m ) [ y r 1 P m , 2 ( 1 y ) ( 1 y ) r 1 P m , 2 ( y ) ] = ( r 2 m ) ( ( j = 0 ( m 1 + j j ) y 1 2 ( r 1 ) ( 1 y ) j ) 2 ( j = 0 ( m 1 + j j ) ( 1 y ) 1 2 ( r 1 ) y j ) 2 ) = ( r 2 m ) ( j = 0 ( m 1 + j j ) ( y 1 2 ( r 1 ) ( 1 y ) j + ( 1 y ) 1 2 ( r 1 ) y j ) ) × ( j = 0 ( m 1 + j j ) ( y 1 2 ( r 1 ) ( 1 y ) j ( 1 y ) 1 2 ( r 1 ) y j ) ) .

For IV, by (2) of Lemma 2, we have

IV 2 = ( 1 y ) r 2 m P m , ( y ) R m , ( y ) y r 2 m P m , ( 1 y ) R m , ( 1 y ) = R m , ( y ) j = 0 ( m 1 + j j ) y j ( 1 y ) r 2 m R m , ( 1 y ) j = 0 ( m 1 + j j ) ( 1 y ) j y r 2 m = ( m + ) ( m 1 + j j ) ( j = 0 ( m 1 + j j ) ( y r 1 ( 1 y ) + j y + j ( 1 y ) r 1 ) ) .

Since

y 1 2 ( r 1 ) ( 1 y ) j ( 1 y ) 1 2 ( r 1 ) y j { 0 , y [ 0 , 1 2 ] ; 0 , y [ 1 2 , 1 ] ,

and for every j,

y r 1 ( 1 y ) + j y + j ( 1 y ) r 1 { 0 , y [ 0 , 1 2 ] ; 0 , y [ 1 2 , 1 ] ,

we have

S (y)=III+IV { 0 , y [ 0 , 1 2 ] ; 0 , y [ 1 2 , 1 ] .

This means S(y) reaches its minimum at point y=1/2. Furthermore, we have

S(1/2)=2 ( 2 r 2 P m , ( y ) ) 2 = 2 1 r 2 ( j = 0 ( m 1 + j j ) ) 2 .

This completes the lemma. □

3 Regularity and stability of scaling function

In this section, we discuss the regularity and stability of a scaling function generated by the refinement mask of a smoothed pseudo-spline. Let

ϕ ˆ (ξ):= j = 1 a ˆ ( 2 j ξ ) .

Then the decay of | ϕ ˆ | can be characterized by | a ˆ | as stated in the following theorem.

Theorem 1 [2]

Let a ˆ (ξ) be a refinement mask of the refinable function ϕ of the form

| a ˆ (ξ)|= cos n (ξ/2)|L(ξ)|,ξ[π,π].

Suppose that

| L ( ξ ) | | L ( 2 π 3 ) | for  | ξ | 2 π 3 , | L ( ξ ) L ( 2 ξ ) | | L ( 2 π 3 ) | 2 for  2 π 3 | ξ | π .

Then | ϕ ˆ (ξ)|C ( 1 + | ξ | ) n + κ , with κ=log(L(| 2 π 3 )|)/log2, and this decay is optimal.

In order to use this lemma, one needs to consider the polynomial corresponding to L(ξ). In fact, Dong and Shen give an important proposition to estimate it in the following proposition.

Proposition 1 [2]

Let P m , (y) be defined as in Section  2, where m, are nonnegative integers with m1. Then

P m , ( y ) P m , ( 3 4 ) for  y [ 0 , 3 4 ] , P m , ( y ) P m , ( 4 y ( 1 y ) ) ( P m , ( 3 4 ) ) 2 for  y [ 3 4 , 1 ] .

Combing Theorem 1 and Proposition 1, we have the following theorem, which characterizes the regularity of a smoothed pseudo-spline.

Theorem 2 Let 2ϕ be the smoothed pseudo-spline of Type II with order r, m, , then

| 2 ϕ ˆ (ξ)|C ( 1 + | ξ | ) r + κ ,

where κ=log( P m , ( 3 4 ))/log2. Consequently, ϕ 2 C α 2 ϵ where α 2 =rκ1. Furthermore, let 1ϕ be the smoothed pseudo-spline of Type I with order n, m, . Then

| 1 ϕ ˆ (ξ)|C ( 1 + | ξ | ) n + κ 2 .

Consequently, ϕ 1 C α 1 ϵ with α 1 =n κ 2 1.

Proof Notice that |L(ξ)| in Theorem 1 is exactly P m , ( sin 2 ( ξ 2 )) and 4y(1y)= sin 2 (ξ); one can easily prove this theorem by Theorem 1 and Proposition 1. □

This theorem shows ϕ k L 2 (R) for k=1,2. Since r2m, the regularity of ϕ is better than a pseudo-spline but the support is longer. For r=6, m=2, =1 the smoothed pseudo-spline ϕ r , m , 2 is shown in Figure 1.

Figure 1
figure 1

The scaling function ϕ 3 , 2 , 1 2 .

Now, we consider the stability of the smoothed pseudo-spline. When ϕ is compactly supported in L 2 (R), it was shown by Jia and Micchelli [11] that the upper Riesz bound in (1.1) always exists. Furthermore, they assert that the existence of a lower Riesz bound is equivalent to

( ϕ ˆ ( ξ + 2 k π ) ) k Z 0for all ξR,
(3.1)

where 0 denotes the zero sequence in 2 (Z). Since a smoothed pseudo-spline is compactly supported and belongs to L 2 (R) for k=1,2, the stability is equivalent to (3.1).

Theorem 3 Smoothed pseudo-splines are stable.

Proof By the definition of refinement mask, for each fixed r 2 m1 and for any 0m1, cos 2 m (ξ/2) 2 a ˆ r , m , (ξ) holds for all ξR. Therefore, we have

| B ˆ r (ξ)|| ϕ ˆ r , m , (ξ)|,

where B r stands for the B-spline with order r. Since B r is stable, the vector ( B ˆ r ( ξ + 2 k π ) ) k Z 0. Hence ( ϕ ˆ r , m , ( ξ + 2 k π ) ) k Z 0.

For a smoothed pseudo-spline of Type I, since a ˆ 2 n , m , 2 (ξ)= | 1 a ˆ n , m , | 2 = 1 a ˆ n , m , (ξ) 1 a ˆ n , m , (ξ), one has

ϕ ˆ 2 n , m , 2 (ξ)= 1 ϕ ˆ n , m , (ξ) 1 ϕ ˆ n , m , (ξ).

Therefore, the set of zeros of ϕ ˆ n , m , 1 (ξ) is contained in that of ϕ ˆ 2 n , m , 2 (ξ) and this guarantees the stability of ϕ 1 (ξ). □

This theorem shows the stability of a smoothed pseudo-spline. From the definition of a Riesz basis, one can find that the translate of a smoothed pseudo-spline is also linearly independent.

4 Riesz wavelets

Since all smoothed pseudo-splines are compactly supported, refinable, stable in L 2 (R), the sequence of spaces ( V n ) n Z defined via Definition 1 forms an MRA. The corresponding wavelets can be constructed by the classical method. Define

ψ ˆ (ξ)= b ˆ ( ξ 2 ) ϕ ˆ ( ξ 2 ) ,where  b ˆ (ξ)= e i ξ a ˆ ( ξ + π ) ¯

and X(ψ):={ ψ n , k = 2 n / 2 ψ( 2 n k),n,kZ}. Then X(ψ) is a Riesz basis. To prove this, the following theorem is needed.

Theorem 4 [2]

Let a ˆ (ξ) be a finitely supported refinement mask of a refinable function ϕ L 2 (R) with a ˆ (0)=1 and a ˆ (π)=0, such that a ˆ can be factorized into the form

| a ˆ (ξ)|= cos n (ξ/2)|L(ξ)|,ξ[π,π],

where is the Fourier series of a finitely supported sequence with L(π)0. Suppose that

| a ˆ (ξ) | 2 +| a ˆ (ξ+π) | 2 0,ξ[π,π].

Define ψ ˆ (2ξ):= e i ξ a ˆ ( ξ + π ) ¯ ϕ ˆ (ξ) and L ˜ := L ( ξ ) | a ˆ ( ξ ) | 2 + | a ˆ ( ξ + π ) | 2 . Assume that L ( ξ ) L ( R ) < 2 n 1 and L ˜ ( ξ ) L ( R ) < 2 n 1 . Then X(ψ) is a Riesz basis for L 2 (R).

From the above theorem, the key step is to estimate the upper Riesz bound of |L(ξ)| and | L ˜ (ξ)|. Notice that

| 1 a ˆ n , m , (ξ) | 2 = 2 a ˆ 2 n , m , (ξ)= cos 2 n (ξ/2) P m , ( sin 2 ( ξ / 2 ) ) .

One has | 1 L(ξ)|= ( P m , ( sin 2 ( ξ / 2 ) ) ) 1 2 , | L 2 (ξ)|= P m , ( sin 2 (ξ/2)) and

| L ˜ 1 |= ( P m , ( y ) ) 1 2 R 2 n , m , ( y ) + R 2 n , m , ( 1 y ) ,| L ˜ 2 |= P m , ( y ) R r , m , 2 ( y ) + R r , m , 2 ( 1 y ) .

Thus, we have the following theorem.

Theorem 5 Let k ϕ, k=1,2 be the smoothed pseudo-spline of Types I and II with order (r,n,m,). The refinement masks k a are given in (1.3) and (1.4). Define

ψ ˆ k (2ξ):= e i ξ a ˆ k ( ξ + π ) ¯ k ϕ ˆ (ξ).

Then X(ψ) forms a Riesz basis for L 2 (R).

Proof By (1) of Lemma 3, one obtains

1 L ˜ L ( R ) = sup y [ 0 , 1 ] ( P m , ( y ) ) 1 2 R 2 n , m , ( y ) + R 2 n , m , ( 1 y ) ( m + ) 1 2 min y [ 0 , 1 ] ( R 2 n , m , ( y ) + R 2 n , m , ( 1 y ) ) 2 n + 1 ( m + ) 1 2 j = 0 ( m + j ) .

Applying Lemma 1, one obtains 1 L ˜ 2 n 1 < 2 n 1 2 . Similarly, one can get

2 L ˜ L ( R ) 2 r 1 < 2 r 1 2 .

Notice that

| k a ˆ ( ξ ) | 2 + | k a ˆ ( ξ + π ) | 2 1for all ξR.

Hence, | k L(ξ)|| k L ˜ (ξ)| for k=1,2. By using Theorem 4, one gets the desired result. □

By definition, the wavelets are also in L 2 (R) and have the same regularity as the scaling function. Still, the support is longer than for pseudo-spline wavelets. For r=6, m=2, =1, the smoothed pseudo-spline ψ 6 , 2 , 1 2 is shown in Figure 2.

Figure 2
figure 2

The corresponding wavelets.

References

  1. Daubechies I, Han B, Ron A, Shen Z: Framelets: MRA-based constructions of wavelet frames. Appl. Comput. Harmon. Anal. 2003,14(1):1-46. 10.1016/S1063-5203(02)00511-0

    Article  MathSciNet  MATH  Google Scholar 

  2. Dong B, Shen Z: Pseudo-splines, wavelets and framelets. Appl. Comput. Harmon. Anal. 2007,22(1):78-104. 10.1016/j.acha.2006.04.008

    Article  MathSciNet  MATH  Google Scholar 

  3. Dong B, Dyn N, Hormann K: Properties of dual pseudo-splines. Appl. Comput. Harmon. Anal. 2010,29(1):104-110. 10.1016/j.acha.2009.08.010

    Article  MathSciNet  MATH  Google Scholar 

  4. Shen Y, Li S, Mo Q: Complex wavelets and framelets from pseudo splines. J. Fourier Anal. Appl. 2010,16(6):885-900. 10.1007/s00041-009-9095-8

    Article  MathSciNet  MATH  Google Scholar 

  5. Cai J-F, Shen Z: Framelet based deconvolution. J. Comput. Math. 2010,28(3):289-308.

    MathSciNet  MATH  Google Scholar 

  6. Cai J, Chan R, Shen L, Shen Z: Tight frame based method for high-resolution image reconstruction. Ser. Contemp. Appl. Math. CAM 14. In Wavelet Methods in Mathematical Analysis and Engineering. Higher Ed. Press, Beijing; 2010:1-36.

    Chapter  Google Scholar 

  7. Daubechies I CBMS-NSF Regional Conference Series in Applied Mathematics 61. In Ten Lectures on Wavelets. SIAM, Philadelphia; 1992.

    Chapter  Google Scholar 

  8. Deriaz E, Perrier V: Direct numerical simulation of turbulence using divergence-free wavelets. Multiscale Model. Simul. 2008,7(3):1101-1129.

    Article  MathSciNet  MATH  Google Scholar 

  9. Chuang Z-T, Liu Y: Wavelets with differential relation. Acta Math. Sin. Engl. Ser. 2011,27(5):1011-1022.

    MathSciNet  MATH  Google Scholar 

  10. Unser M, Blu T: Fractional splines and wavelets. SIAM Rev. 2000 (electronic),42(1):43-67. (electronic) 10.1137/S0036144598349435

    Article  MathSciNet  MATH  Google Scholar 

  11. Jia RQ, Micchelli CA: Using the refinement equations for the construction of pre-wavelets. II. Powers of two. In Curves and Surfaces. Academic Press, Boston; 1991: (Chamonix-Mont-Blanc, 1990)209-246. (Chamonix-Mont-Blanc, 1990)

    Google Scholar 

Download references

Acknowledgements

This study was partially supported by the Joint Funds of the National Natural Science Foundation of China (Grant No. U1204103). The authors would like to express their thanks to the reviewers for their helpful comments and suggestions.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Zhitao Zhuang.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

All authors contributed to each part of this work equally and read and approved the final manuscript.

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (https://creativecommons.org/licenses/by/4.0), which permits use, duplication, adaptation, distribution, and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhuang, Z., Yang, J. A class of generalized pseudo-splines. J Inequal Appl 2014, 359 (2014). https://doi.org/10.1186/1029-242X-2014-359

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1029-242X-2014-359

Keywords