<figure><img src="image_1.png" /></figure>

HAL
open science

# Fractional Wavelet-Based Generative Scattering Networks

Jiasong Wu, Xiang Qiu, Jing Zhang, Fuzhi Wu, Youyong Kong, Guanyu

Yang, Lotfi Senhadji, Huazhong Shu

► To cite this version:

Jiasong Wu, Xiang Qiu, Jing Zhang, Fuzhi Wu, Youyong Kong, et al. Fractional Wavelet-Based Generative Scattering Networks. Frontiers in Neurorobotics, 2021, 15, pp.752752. 10.3389/fnbot.2021.752752. hal-03467433

HAL Id: hal-03467433

https://hal.science/hal-03467433

Submitted on 2 Jun 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Fractional Wavelet-Based Generative Scattering Networks

Jiasong Wu<sup>1,2,3,4*</sup>, Xiang Qiu<sup>1,4</sup>, Jing Zhang<sup>1,4</sup>, Fuzhi Wu<sup>1,4</sup>, Youyong Kong<sup>1,2,4</sup>, Guanyu Yang<sup>1,2,4</sup>, Lotfi Senhadji<sup>3,4</sup> and Huazhong Shu<sup>1,2,4</sup>

¹ Laboratory of Image Science and Technology, Key Laboratory of Computer Network and Information Integration, Southeast University, Ministry of Education, Nanjing, China, ² Jiangsu Provincial Joint International Research Laboratory of Medical Information Processing, School of Computer Science and Engineering, Southeast University, Nanjing, China, ³ Univ Rennes, INSERM, LTSI-UMR 1099, Rennes, France, ⁴ Centre de Recherche en Information Biomédicale Sino-Français (CRIBs), Univ Rennes, INSERM, Rennes, France

Generative adversarial networks and variational autoencoders (VAEs) provide impressive image generation from Gaussian white noise, but both are difficult to train, since they need a generator (or encoder) and a discriminator (or decoder) to be trained simultaneously, which can easily lead to unstable training. To solve or alleviate these synchronous training problems of generative adversarial networks (GANs) and VAEs, researchers recently proposed generative scattering networks (GSNs), which use wavelet scattering networks (ScatNets) as the encoder to obtain features (or ScatNet embeddings) and convolutional neural networks (CNNs) as the decoder to generate an image. The advantage of GSNs is that the parameters of ScatNets do not need to be learned, while the disadvantage of GSNs is that their ability to obtain representations of ScatNets is slightly weaker than that of CNNs. In addition, the dimensionality reduction method of principal component analysis (PCA) can easily lead to overfitting in the training of GSNs and, therefore, affect the quality of generated images in the testing process. To further improve the quality of generated images while keeping the advantages of GSNs, this study proposes generative fractional scattering networks (GFRSNs), which use more expressive fractional wavelet scattering networks (FrScatNets), instead of ScatNets as the encoder to obtain features (or FrScatNet embeddings) and use similar CNNs of GSNs as the decoder to generate an image. Additionally, this study develops a new dimensionality reduction method named feature-map fusion (FMF) instead of performing PCA to better retain the information of FrScatNets,; it also discusses the effect of image fusion on the quality of the generated image. The experimental results obtained on the CIFAR-10 and CelebA datasets show that the proposed GFRSNs can lead to better generated images than the original GSNs on testing datasets. The experimental results of the proposed GFRSNs with deep convolutional GAN (DCGAN), progressive GAN (PGAN), and CycleGAN are also given.

**Keywords:** generative model, fractional wavelet scattering network, image generation, image fusion, feature-map fusion

OPEN ACCESS

**Edited by:**
Dongming Zhou,
Yunnan University, China

**Reviewed by:**
Jianhua Yin,
Shandong University, China
Gao Chen,
Dongguan University of
Technology, China
Jinyu Cai,
Fuzhou University, China

*Correspondence:
Jiasong Wu
jswu@seu.edu.cn

**Received:** 03 August 2021
**Accepted:** 20 September 2021
**Published:** 26 October 2021

**Citation:**
Wu J, Qiu X, Zhang J, Wu F, Kong Y, Yang G, Senhadji L and Shu H (2021) *Fractional Wavelet-Based Generative Scattering Networks*. Front. Neurorobot. 15:752752. doi: 10.3389/fnbot.2021.752752
# INTRODUCTION

Generative models have recently attracted the attention of many researchers, and they are widely used in image synthesis, image restoration, image inpainting, image reconstruction, and other applications. Many generative models have been proposed in the literature. They can be roughly classified into two types (Goodfellow et al., 2014): explicit density and implicit density models.

Among explicit density generative models, variational auto-encoders (VAEs) (Kingma and Welling, 2014) and their variants (Rezende et al., 2014; Salimans et al., 2015; Gregor et al., 2018) are most likely the most commonly used models, since they have useful latent representation, which can be used in inference queries. Kingma and Welling (2014) were the first to propose VAEs, which train an encoder and decoder simultaneously and can perform efficient inference and learning in directed probabilistic models and in the presence of continuous latent variables with intractable posterior distributions. Salimans et al. (2015) bridged the gap between Markov chain Monte Carlo (MCMC) and VAEs, and incorporated one or more steps of MCMC into variational approximation. Sohn et al. (2015) proposed a conditional VAE (CVAE), which joins existing label information in training to generate corresponding category data. Rezende and Mohamed (2015) introduced a new approach for specifying flexible, arbitrarily complex, and scalable approximate posterior distributions and made a clear improvement in the performance and applicability of variational inference. Sønderby et al. (2016) presented a ladder variational autoencoder, which uses a process similar to a ladder network and recursively corrects the generation distribution based on a data-independent approximate likelihood. Higgins et al. (2017) presented a β-VAE, which is a modification of a variational autoencoder (VAE), with special emphasis on discovering disentangled latent factors. Oord et al. (2017) proposed a simple yet powerful generative model that learns discrete representations and allowed the model to circumvent issues of posterior collapse. Gregor et al. (2018) proposed temporal difference VAE (TD-VAE), which is a generative sequence model that learns representations containing explicit beliefs about states in several steps into the future Razavi et al. (2019) proposed vector quantized variational autoencoder (VQ-VAE), which augments with powerful priors over latent codes and is able to generate samples with a quality that rival those of state-of-the-art GANs on multifaceted datasets, such as ImageNet. Simonovsky and Komodakis (2018) proposed Graph VAE, sidesteps the hurdles of linearization of discrete structures by outputting a probabilistic fully connected graph of a predefined maximum size directly at once. For more references on VAEs, see Blei et al. (2017).

Among implicit density generative models, generative adversarial networks (GANs) (Goodfellow et al., 2014) and their variants (Chen et al., 2016; Radford et al., 2016) are probably the most commonly used models, since they provide better generated images than other generative models. Goodfellow et al. (2014) were the first to propose GANs, which estimate generative models *via* an adversarial process, where a generative model G

and a discriminative model D are trained simultaneously without the need for Markov chains or unrolled approximate inference networks during either training or the generation of samples. However, the application of GANs to real-world computer vision problems still encounters at least three significant challenges (Wang et al., 2021): (1) high-quality image generation; (2) diverse image generation; and (3) stable training. Therefore, many variants of GANs have been proposed to handle the three challenges. The variants can be roughly classified into two groups (Wang et al., 2021): architecture variant GANs and loss variant GANs.

In terms of architecture variant GANs, for example, Radford et al. (2016) proposed deep convolutional GAN (DCGAN), which uses a convolutional neural network (CNN) as the discriminator D and deploys a deconvolutional neural network architecture for G; the spatial upsampling ability of the deconvolution operation enables the generation of images with higher resolution compared with the original GANs. Mirza and Osindero (2014) proposed conditional GAN (CGAN), which imposes a condition of additional information, such as a class label, to control the process of data generation in a supervised manner. Chen et al. (2016) presented InfoGAN, which decomposes an input noise vector into a standard incompressible latent vector and another latent variable to capture salient semantic features of real samples. Karras et al. (2018) presented progressive GAN (PGAN) for generative high-resolution images using the idea of progressive neural networks (Rusu et al., 2017), which does not suffer from forgetting and is able to deploy prior knowledge *via* lateral connections to previously learned features. Karras et al. (2020a,b) proposed StyleGAN, which leads to an automatically learned, unsupervised separation of high-level attributes and stochastic variation in generated images and, thus, enables intuitive, scale-specific control of the synthesis. More recently, Hudson and Zitnick (2021) introduced the Generative Adversarial Transformer (GANformer), which is a generalization of the StyleGAN and a simple yet effective generalization of the vanilla transformer, for a visual synthesis task.

In terms of loss-variant GANs, for example, Arjovsky et al. (2017) proposed Wasserstein GAN (WGAN), which uses the Wasserstein distance as the loss measure for optimization instead of Kullback-Leibler divergence. Gulrajani et al. (2017) proposed an improved method for training the discriminator for a WGAN, by penalizing the norm of discriminator gradients with respect to data samples during training rather than performing parameter clipping. Nowozin et al. (2016) proposed an alternative cost, which is a function of the f-divergence, for updating the generator, which is less likely to saturate at the beginning of training. Zhu et al. (2017) proposed CycleGAN for the task of image-to-image translation. Qi (2020) presented loss-sensitive GAN (LS-GAN), which trains the generator to produce realistic samples by minimizing the designated margins between real and generated samples. Miyato et al. (2018) proposed spectral normalization GAN (SN-GAN), which uses a weight normalization technique to train the discriminator more stably. Brock et al. (2019) proposed BigGAN, which uses hinge loss instead of Jensen-Shannon divergence and a large-scale dataset to train the generator to produce more realistic samples.
Although GANs and VAEs are great generative models, they raise many questions. A significant disadvantage of VAEs is that the resulting generative models produce blurred images compared with GANs, since the quality of VAEs crucially relies on the expressiveness of their inference models. A significant disadvantage of GANs is that the training process is very difficult and may lead to unstable training and model collapse. To design a network that can maintain the characteristics of high-quality generated images of GANs as much as possible while reducing the training difficulty of GANs, Angles and Mallat (2018) proposed generative scattering networks (GSNs), which use wavelet scattering networks (ScatNets) (Bruna and Mallat, 2013) as the encoder to obtain features (or ScatNet embeddings) and the deconvolutional neural network of DCGAN (Radford et al., 2016) as the decoder to generate an image. The advantage of GSNs is that there is no need to learn the parameters of ScatNets; therefore, the difficulty of training is reduced when compared with DCGAN, while the disadvantage of GSNs is that generated images can lose details, which affects the quality of the generated images. After careful inspection, we determined that the sources of relatively low-quality generated images of GSNs include at least two aspects: (1) the expression ability of ScatNets is slightly weaker than that of CNNs used in DCGAN; (2) applying PCA (Abdi and Williams, 2010) to reduce the dimension of the feature map of ScatNets in the encoder part of GSNs leads to an overfitting problem in the testing process of GSNs. This finding leads to the central question of our study:

Can we change the feature extraction method of ScatNets to a more powerful one that still does not need learning? Can we develop a more suitable dimensionality reduction method to solve the overfitting problem in the testing process of GSNs?

In an attempt to solve the above questions, in this study, we propose generative fractional scattering networks (GFRSNs), which can be seen as an extension of GSNs. The contributions of this article are as follows:

1) We use, for more expressiveness, fractional wavelet scattering networks (FrScatNets) (Liu et al., 2019) instead of ScatNets (Bruna and Mallat, 2013) to extract features of images, and we use image fusion (Liu et al., 2016; Yang et al., 2017) in GFRSNs to effectively improve the visual quality of the generated images.

2) We propose a new dimensionality reduction method named feature-map fusion (FMF), which is more suitable for reducing the feature dimension of FrScatNets than PCA, since the FMF method greatly alleviates the overfitting problem on the testing datasets using GFRSNs.

3) The image generated by the proposed GFRSN on the test set is better than that produced by the original GSNs.

The remainder of this article is organized as follows: In section Generative Scattering Networks (GSNs), wavelet scattering networks and the architectural components of GSNs are briefly introduced. The main architectural components of GFRSNs, which include fractional wavelet scattering networks, the FMF dimensionality reduction method and generative networks, and an image fusion method are introduced in section Generative Fractional Scattering Networks (GFRSNs). The performance of

<figure><img src="image_3.png" /><figcaption>FIGURE 1 | Structure of generative scattering networks (GSNs) with principal component analysis (PCA) dimensionality reduction method.</figcaption></figure>

GFRSNs is analyzed and compared with that of the original GSNs in section Numerical Experiments. The conclusions and further discussion are presented in section Conclusions.

## GENERATIVE SCATTERING NETWORKS (GSNS)

In this section, we first briefly recall the generative scattering networks (GSNs) (Angles and Mallat, 2018), whose structure is shown in Figure 1.

The input Mth-order tensor $\mathcal{X} \in \mathbb{R}^{N_1 \times N_2 \times \cdots \times N_K}$, where $\mathbb{R}$ denotes the real domain and each $N_i, i = 1, 2, 3, \cdots, K$, addresses the i-mode of $\mathcal{X}$, and is first fed into the feature extraction part of the encoder to obtain the ScatNet features $y \in \mathbb{R}^{M_1 \times M_2 \times \cdots \times M_L}$. The next part of the encoder aims to map the features to a Gaussian latent variable $z \in \mathbb{R}^U$, which is accomplished by whitening and projection to a lower-dimensional space. Inspired by Zou and Lerman (2019), we refer to this process as Gaussianization. Decoder G can be seen as a generator and is trained by minimizing the reconstruction loss between the output $\tilde{\mathcal{X}} \in \mathbb{R}^{N_1 \times N_2 \times \cdots \times N_K}$ and input $\mathcal{X}$. In other words, the generator calculation is regarded as the inverse problem of the scattering transform. The main components of GSNs include ScatNets, Gaussianization with PCA, and the generative network G. These components are recalled as follows.

### Wavelet Scattering Networks (ScatNets)

Let the complex bandpass filter $\psi_{\lambda}$ be constructed by scaling and rotating a filter $\psi$, respectively, by $2^j$ and $\delta$, as follows (Bruna and Mallat, 2013):

$$ \psi_{\lambda}(t) = 2^{2j} \psi(2^j \delta^{-1} t), \quad \lambda = 2^j \delta, \qquad (1) $$

with $0 \le j \le J - 1$, and $\delta = k\pi/K$, $k = 0, 1, \dots, K - 1$.

The wavelet-modulus coefficients of $x$ are given by:

$$ U[\lambda] x = |x * \psi_{\lambda}|. \qquad (2) $$

The scattering propagator $U[p]$ is defined by cascading wavelet-modulus operators

$$ \begin{aligned} U[p] x &= U[\lambda_m] \cdots U[\lambda_2] U[\lambda_1] x \\ &= \left| \left| x * \psi_{\lambda_1} \right| * \psi_{\lambda_2} \right| \cdots * \psi_{\lambda_m} \right|, \quad (3) \end{aligned} $$
where $p = (\lambda_1, \lambda_2, \dots, \lambda_m)$ are the frequency-decreasing paths; in other words, $|\lambda_k| \ge |\lambda_{k+1}|$, $k = 1, 2, \dots, m-1$. Note that $U[\emptyset]x = x$, and $\emptyset$ expresses the empty set.

The scattering operator $S_J$ performs spatial averaging on a domain whose width is proportional to $2^J$:

$$ S[p]x = U[p]x * \phi_J = U[\lambda_m] \cdots U[\lambda_2] U[\lambda_1] x * \phi_J \\ = \left| \left| x * \psi_{\lambda_1} \right| * \psi_{\lambda_2} \cdots * \psi_{\lambda_m} \right| * \phi_J. \quad (4) $$

The network nodes of layer m correspond to the set $P^m$ of all paths $p = (\lambda_1, \lambda_2, \dots, \lambda_m)$ of length m. This m-th layer stores the propagated signals $\{U[p]x\}_{p \in P^m}$ and outputs the scattering coefficients $\{S[p]x\}_{p \in P^m}$. The output is obtained by cascading the scattering coefficients of every layer.

Note that $x$ in (2) can be one-dimensional data $x \in \mathbb{R}^{N_1}$, two-dimensional data $X \in \mathbb{R}^{N_1 \times N_2}$, and third-order tensor $X \in \mathbb{R}^{N_1 \times N_2 \times N_3}$, which can be seen as $N_3$ two-dimensional data $X \in \mathbb{R}^{N_1 \times N_2}$, and ScatNet addresses with these $X_s$ one by one and then superimposes the results as output features. According to Mallat (2012), if we feed the input $X \in \mathbb{R}^{N_1 \times N_2 \times N_3}$ into ScatNet, then we can obtain ScatNet features (or ScatNet embeddings) as follows:

$$ y = S[p]X \in \mathbb{R}^{N_3 \times (1+L^2J(J-1)/2) \times (N_1/2^J) \times (N_2/2^J)}, \quad (5) $$

where $N_3$ is the number of input sample channels, and $N_1$ and $N_2$ are the width and height of the input sample, respectively. $N_1/2^J$ and $N_2/2^J$ are the width and height of the output features. $J$ is a scale factor, and $L$ is the number of rotation angles. Note that the number of feature maps in the first, second, and third layers is 1, $L$, and $L^2J(J-1)/2$, respectively.

## Gaussianization With PCA

As shown in Figure 1, the last step of the encoder maps the transformed features in such a way that we can sample from the Gaussian distribution to generate new images, as required by the generator. Specifically, let $\{y\}_{t=1}^T$ be the output features of the ScatNet embedding, and let $y$ be the representing matrix of $\{y\}_{t=1}^T$, while $z$ is the latent variable of the generator. As advocated in Angles and Mallat (2018), $z$ can be interpreted as an address, with a dimension $d$ lower than that in the input image. Hence, to get a lower-dimensional embedding of the output features, one can perform principal component analysis (PCA) (Abdi and Williams, 2010) to project the features of the scattering transform to a lower-dimensional space.

Next, to whiten them, we choose $u = \frac{1}{T} \sum_{t=1}^T y$, $\Sigma = \frac{1}{T} \sum_{t=1}^T (y - u)(y - u)^*$, and the whitening map $A = \Sigma^{-1/2} (Id - u)$.

Hence, the resulting embedding of the encoder is

$$ z = \sum^{-1/2} (y - u). \quad (6) $$

After the above process, the whitened sample is uncorrelated, and their distribution will be close to a normal one with identity covariance (Angles and Mallat, 2018), which is exactly what we want to feed to the generator.

## Generator Networks in GSNs

The generative network G of GSNs is a neural one, which is similar to the generator of DCGAN (Radford et al., 2016), which inverts the whitened scattering embedding on training samples. The generator network G includes a fully connected layer (FC), batch normalization layer (BN) (Ioffe and Szegedy, 2015), bilinear upsampling (Upsample) layer, and convolutional layer (Conv2d) with a kernel size of $7 \times 7$. Except for the last layer, which uses the tanh activation function, the others use the default ReLU (Nair and Hinton, 2010) activation function.

Generative scattering networks with PCA as the dimensionality reductional method choose the $L_1$-norm loss function and solve the following optimization problem (Zhao et al., 2017):

$$ g_1 = \min \text{Loss}_{L_1}(\mathcal{X}, \tilde{\mathcal{X}}) = \min \frac{1}{N} \sum_{i=1}^{N} |\mathcal{X}^{(i)} - \tilde{\mathcal{X}}^{(i)}|, \quad (7) $$

where $\mathcal{X}$ represents the input data, $\tilde{\mathcal{X}}$ represents the generative data, $\mathcal{X}^{(i)}$ represents the i-th input sample, and $\tilde{\mathcal{X}}^{(i)}$ represents the i-th generative sample:

$$ \tilde{\mathcal{X}} = G(\text{PCA}(S[p]\mathcal{X})), \quad (8) $$

where $S[p]\mathcal{X}$ denotes the feature extraction process with ScatNets, and PCA(.) represents that the feature dimensionality reduction method is PCA. $G(.)$ represents the generative network G. The optimization problems in (7) are then solved with the Adam optimizer (Kingma and Ba, 2015) using the default hyperparameters.

## GENERATIVE FRACTIONAL SCATTERING NETWORKS (GFRSNS)

In this section, we introduce the proposed generative fractional scattering networks (GFRSNs), whose structure is shown in Figure 2.

The input $\mathcal{X} \in \mathbb{R}^{N_1 \times N_2 \times \cdots \times N_K}$ is first fed into the fractional wavelet scattering networks (FrScatNets) to obtain FrScatNet features (or FrScatNet embeddings) $y_\alpha \in \mathbb{R}^{M_1 \times M_2 \times \cdots \times M_L}$, whose dimensions are then reduced by the proposed feature-map fusion (FMF) method to obtain an implicit tensor $z_\alpha \in$

<figure><img src="image_4.png" /><figcaption>FIGURE 2 | Structure of generative fractional scattering networks (GFRSNs).</figcaption></figure>
<figure><img src="image_5.png" /><figcaption>FIGURE 3 | Fractional wavelet scattering network and the feature-map fusion dimensional reduction method.</figcaption></figure>

$\mathbb{R}^{O_1 \times O_2 \times \cdots \times O_K}$, which is then fed into the generator G to obtain the generated output tensor $\tilde{\chi}_{\alpha} \in \mathbb{R}^{N_1 \times N_2 \times \cdots \times N_K}$. In other words, the generative network G is seen as the inverse problem of FrScatNets. The main components of GFRSNs include FrScatNets, Gaussianization with feature-map fusion dimensionality reduction method, and the generative network G. In the following, these components of GFRSNs are introduced.

## Fractional Wavelet Scattering Networks (FrScatNets)

In this subsection, fractional wavelet scattering networks (FrScatNets) (Liu et al., 2019) are briefly introduced. Similar to (2), the fractional wavelet modulus coefficients of $x$ are given by:

$$U_{\alpha}[\lambda]x = |x\Theta_{\alpha}\psi_{\lambda}|, \qquad (9)$$

where $\Theta_{\alpha}$ is the fractional convolution defined by Shi et al. (2010);

$$x(t) \Theta_{\alpha} \psi_{\lambda}(t) = e^{-\frac{i}{2} t^2 \cot \theta} \left[ \left( x(t) e^{\frac{i}{2} t^2 \cot \theta} \right)^* \psi_{\lambda}(t) \right], \quad (10)$$

where the parameter $\alpha$ is the fractional order and $\theta = \alpha\pi/2$ represents the rotation angle. Note that when $\alpha = 1$, the fractional convolution in (10) reduces to conventional convolution in (2).

The fractional scattering propagator $U_{\alpha}[p]$ is defined by cascading fractional wavelet modulus operators

$$\begin{aligned} U_{\alpha}[p]x &= U_{\alpha}[\lambda_m] \cdots U_{\alpha}[\lambda_2] U_{\alpha}[\lambda_1]x \\ &= \left| \left| x\Theta_{\alpha}\psi_{\lambda_1} \mid \Theta_{\alpha}\psi_{\lambda_2} \mid \cdots \mid \Theta_{\alpha}\psi_{\lambda_m} \right| \right|, \quad (11) \end{aligned}$$

where $p = (\lambda_1, \lambda_2, \dots, \lambda_m)$ are the frequency-decreasing paths; in other words, $|\lambda_k| \ge |\lambda_{k+1}|$, $k = 1, 2, \dots, m-1$. Note that $U_{\alpha}[\emptyset]x = x$, and $\emptyset$ expresses the empty set.

The fractional scattering operator $S_{\alpha}$ performs spatial averaging on a domain whose width is proportional to $2^J$:

$$\begin{aligned} S_{\alpha}[p]x &= U_{\alpha}[p]x\Theta_{\alpha}\phi_{J} = U_{\alpha}[\lambda_{m}] \cdots U_{\alpha}[\lambda_{1}]x\Theta_{\alpha}\phi_{J} \\ &= \left| \left| x\Theta_{\alpha}\psi_{\lambda_{1}} \mid \Theta_{\alpha}\psi_{\lambda_{2}} \mid \cdots \mid \Theta_{\alpha}\psi_{\lambda_{m}} \right| \right| \Theta_{\alpha}\phi_{J}. \quad (12) \end{aligned}$$

The structure of FrScatNets is shown on the left of Figure 3.

The network nodes of the layer m correspond to the set $P^m$ of all paths $p = (\lambda_1, \lambda_2, \dots, \lambda_m)$ of length m. This m-th layer stores the propagated signals $\{U_{\alpha}[p]x\}_{p \in P^m}$ and outputs the fractional scattering coefficients $\{S_{\alpha}[p]x\}_{p \in P^m}$. The output is obtained by cascading the fractional scattering coefficients of every layer. Note that when $\alpha = 1$, the FrScatNets in (12) default to conventional ScatNets in (4), since the fractional convolution in (10) reduces to conventional convolution in (2).

Note that FrScatNets retain the advantages of ScatNets, for example, no need for learning, translation-invariant property, linearized deformations, and certain parameters. Compared with ScatNets, FrScatNet adds a free parameter $\alpha$, which represents fractional order. With $\alpha$ continuously growing from 0 to 2, FrScatNets can show the characteristics of an image from time domain to frequency domain. Thus, FrScatNets provide more fractional domain choices for the feature extraction of input data. Furthermore, for the image generation task in this study, we can obtain as many generated images from FrScatNet embeddings as different fractional orders $\alpha_i$, and then they can be fused to further improve the quality of the generated images.
If we feed the input $\mathcal{X} \in \mathbb{R}^{N_1 \times N_2 \times N_3}$ into the FrScatNet, then we can obtain the features of FrScatNet (or FrScatNet embeddings) as follows:

$$y_{\alpha} = S_{\alpha}[p] X \in \mathbb{R}^{N_3 \times (1+LJ+L^2J(J-1)/2) \times (N_1/2^l) \times (N_2/2^l)}. \quad (13)$$

Note that the size of output features of FrScatNets is the same as that of ScatNets, whose size is shown in (5).

## Gaussianization With FMF

In this subsection, we introduce a new method called FMF to reduce the dimensionality of the features after a fractional scattering transformation. We propose such an algorithm based on the hierarchical tree structure of features extracted by the fractional scattering transform to replace PCA to map the features to a low-dimensional space. More specifically, since the output features of different layers from the fractional scattering transform have a hierarchical structure, which is not considered in the PCA algorithm, we need a dimensionality reduction method that can make full use of this hierarchical information. The number of feature maps in the first, second, and third layers is 1, $LJ$, and $L^2J(J-1)/2$, respectively. Obviously, the third layer has the largest number of feature maps. Therefore, we fuse only the feature maps from the third layer of the fractional scattering transform to significantly reduce the data dimension. The fusion method is very simple: we obtain a new feature map by simply taking the average of every $L(J-1)/2$ feature map, which obtains $LJ$ feature maps after applying the FMF method to the output of the third layer of FrScatNets. The dotted box on Figure 3 illustrates the proposed FMF method.

Therefore, an input tensor $\mathcal{X} \in \mathbb{R}^{N_1 \times N_2 \times N_3}$ is fed into the FrScatNets to obtain FrScatNet features $y_{\alpha}$ in (13), which are then processed by the FMF method, obtaining an implicit tensor

$$Z_{\alpha} = \text{FMF}(y_{\alpha}) \in \mathbb{R}^{N_3 \times (1+LJ+LJ) \times (N_1/2^l) \times (N_2/2^l)}, \quad (14)$$

whose size is significantly smaller than the size shown in (13) without using the FMF method. Note that FMF(.) means performing the FMF method.

The obtained implicit tensor $Z_{\alpha}$ is then input to the generator network G, described below, to obtain the generated image.

## Generative Networks in GFRSNs

The generative network G of GFRSNs is also a deconvolutional neural network that has a generator similar to that of DCGAN (Radford et al., 2016), which inverts fractional scattering embeddings on training samples. The generative network G of GFRSNs also includes a fully convolutional layer (Fully Conv) (Long et al., 2015) and several convolution blocks that consist of bilinear upsampling (UP), two convolutional layers (Conv) with a 3 × 3 kernel size, batch normalization, and ReLU (the activation function of the last convolution layer is tanh). GFRSNs also choose the $L_1$-norm loss function and solve the following optimization problem:

$$g_2 = \min \text{Loss}_{L_1}(\mathcal{X}, \tilde{\mathcal{X}}_{\alpha}) = \min \frac{1}{N} \sum_{i=1}^{N} \left| \mathcal{X}^{(i)} - \tilde{\mathcal{X}}_{\alpha}^{(i)} \right|, \quad (15)$$

where $\tilde{\mathcal{X}}_{\alpha}$ represents the generative data and $\tilde{\mathcal{X}}_{\alpha}^{(i)}$ represents the i-th generative sample, and

$$\tilde{\mathcal{X}}_{\alpha} = G(\text{FMF}(S_{\alpha}[p]\mathcal{X})), \quad (16)$$

where $S_{\alpha}[p]\mathcal{X}$ denotes the feature extraction process with FrScatNets, FMF(.) represents the dimensionality reduction process, and $G(.)$ represents the generative network.

The optimization problem in (15) is then solved with the Adam optimizer (Kingma and Ba, 2015).

## Image Fusion

In contrast to GSNs, the proposed generative fractional scattering networks (GFRSNs) embed the input using FrScatNets, which allows for deriving many embeddings, since FrScatNets have an additional fractional order $\alpha$; therefore, we can embed the input in different fractional order domains. These FrScatNet embeddings may extract many different but complementary features from the input. We can effectively use these embeddings to generate many images and further improve the quality of the synthesized images using fusion methods. In this study, as shown at the bottom of Figure 2, we use a simple image fusion method that is weighted average. As examples, we simply use the following:

$$\tilde{\mathcal{X}}_{\alpha_1, \alpha_2} = \lambda \tilde{\mathcal{X}}_{\alpha_1} + (1 - \lambda) \tilde{\mathcal{X}}_{\alpha_2}, \quad (17)$$

where $\lambda$ is the balanced parameter, which is set here to 0.5.

## NUMERICAL EXPERIMENTS

In this section, we evaluate the quality of the generated images by the proposed GFRSNs by means of several experiments. The quality of the generated images is evaluated with two criteria: peak signal to noise ratio (PSNR) (Wang et al., 2003) and structural similarity (SSIM) (Wang et al., 2004).

We performed experiments on two datasets that have different levels of variability: CIFAR-10 (Krizhevsky, 2009) and CelebA (Liu et al., 2015). The CIFAR-10 dataset includes 50,000 training images and 10,000 testing images, whose sizes are 32 × 32 × 3. In all the experiments on the CIFAR-10 dataset, after image grayscale preprocessing, the number of rotation angles $L$ is set to 8, and the fractional scattering averaging scale is set to $2^l = 2^3 = 8$, which means that we linearize translations and deformations of up to 8 pixels. Therefore, the size of the output features from FrScatNets according to Equation (13) is 1 × 217 × 4 × 4, which is then, after the FMF method according to Equation (14), reduced to 1 × 49 × 4 × 4 (the size of implicit tensor $Z_{\alpha}$). In addition, the CelebA dataset contains thousands of images, and we choose 65,536 training images and 16,384 testing images, whose sizes are 128 × 128 × 3. In all the experiments on the CelebA dataset, after image grayscale preprocessing, the number of rotation angle $L$ is set to 8, and the fractional scattering averaging scale is set to $2^l = 2^4 = 16$, which means that we linearize translations and deformations of up to 16 pixels. Thus, the size of the output features from FrScatNets according to (13) is 1 × 417 × 8 × 8, which is then, after FMF method according
<table><caption>**TABLE 1** | Core parameters of FrScatNet with and without feature dimensionality reduction.</caption><thead><tr><th rowspan="2">Parameter</th><th rowspan="2">Descriptions</th><th colspan="2">Dataset</th></tr><tr><th>CIFAR-10</th><th>CelebA</th></tr></thead><tbody><tr><td>N<sub>1</sub> × N<sub>2</sub> ×<br>N<sub>3</sub></td><td>The size of input image</td><td>32 × 32 × 1</td><td>128 × 128 × 1</td></tr><tr><td>J</td><td>The fractional scattering averaging scale</td><td>3</td><td>4</td></tr><tr><td>L</td><td>The number of rotation angle</td><td>8</td><td>8</td></tr><tr><td>N<sub>3</sub> × (1 + LJ +<br>L<sup>2</sup>JJ<sup>-1</sup>) × N<sub>1</sub> / 2J × N<sub>2</sub> / 2J</td><td>The shape of FrScatNets features y<sub>α</sub></td><td>1 × 217 × 4 × 4</td><td>1 × 417 × 8 × 8</td></tr><tr><td>N<sub>3</sub> ×<br>(1 + 2 × LJ) ×<br>N<sub>1</sub> / 2J × N<sub>2</sub> / 2J</td><td>The shape of implicit tensor Z<sub>α</sub> with FMF</td><td>1 × 49 × 4 × 4</td><td>1 × 65 × 8 × 8</td></tr></tbody></table>

to Equation (14), reduced to 1 × 65 × 8 × 8 (the size of implicit tensor $Z_α$). **Table 1** shows the core parameters of FrScatNet and its settings on the CIFAR-10 and CelebA datasets.

In the following, we first compare the visual quality of the generated images with different feature dimensionality reduction methods in the framework of GFRSNs. Then, we compare the visual quality of the generated images with FrScatNets. Finally, we compare the visual quality of the fused images and unfused images. The following experiments are implemented using PyTorch on a PC machine, which sets up an Ubuntu 16.04 operating system and has an Intel (R) Core(TM) i7-8700K CPU with a speed of 3.7 GHz and 32 GB RAM, and has two NVIDIA GeForce GTX1080-Ti GPUs.

## Image Generative Results With Different Dimensionality Reduction Methods

In this subsection, we compare the results on the quality of generative images with two different dimensionality reduction methods: the PCA method and the proposed FMF method. We set the fractional orders to be $α_1 = α_2 = 1$, and use conventional ScatNets to extract features from the input $\mathcal{X}$ for simplicity.

For the PCA-based GFRSNs, the flow chart is shown in **Figure 1**. For the CIFAR-10 dataset, the size of the implicit vector $\mathbf{z}$ is 49 × 4 × 4 = 784, and for the CelebA dataset, the size of the implicit vector $\mathbf{z}$ is 65 × 8 × 8 = 4,160. We use the PyTorch code of generative scattering networks<fnref n="1" /> provided by Tomás Angles. The PSNR and SSIM on the CIFAR-10 and CelebA datasets are shown in the second columns of **Tables 2, 3**, respectively.

As shown in the two tables, the scores of PSNR (Train PSNR) and SSIM (Train SSIM), both in the training dataset, are very good for the PCA-based GFRSNs; however, their corresponding values (test PSNR and test SSIM) in the testing dataset are slightly low. This phenomenon indicates that an overfitting problem has occurred using the PCA-based GFRSNs. We argue the reason behind this phenomenon could be that the output feature of

<table><caption>**TABLE 2** | Peak signal to noise ratio (PSNR) and structural similarity (SSIM) scores of training and testing images from FrScatNets with fractional orders $α_1 = α_2 = 1$ on the CIFAR-10 dataset.</caption><thead><tr><th></th><th>PCA</th><th>Feature-Map Fusion</th><th>Increased (%)</th></tr></thead><tbody><tr><td>Train PSNR</td><td><strong>23.08</strong></td><td>20.1500</td><td>-12.69</td></tr><tr><td>Test PSNR</td><td>17.92</td><td><strong>18.1000</strong></td><td><strong>1.00</strong></td></tr><tr><td>Train SSIM</td><td><strong>0.9428</strong></td><td>0.8859</td><td>-6.08</td></tr><tr><td>Test SSIM</td><td>0.8206</td><td><strong>0.8352</strong></td><td><strong>1.78</strong></td></tr></tbody></table>

Increased means the percentages of relative improvements of FMF over principal component analysis (PCA), the better results are shown in bold.

<table><caption>**TABLE 3** | PSNR and SSIM scores of training and testing images from FrScatNets with fractional orders $α_1 = α_2 = 1$ on the CelebA dataset.</caption><thead><tr><th></th><th>PCA</th><th>Feature-Map Fusion</th><th>Increased (%)</th></tr></thead><tbody><tr><td>Train PSNR</td><td><strong>23.8124</strong></td><td>22.7526</td><td>-4.45</td></tr><tr><td>Test PSNR</td><td><strong>20.5312</strong></td><td>19.7688</td><td>-3.71</td></tr><tr><td>Train SSIM</td><td><strong>0.9529</strong></td><td>0.944</td><td>-0.93</td></tr><tr><td>Test SSIM</td><td><strong>0.9104</strong></td><td>0.8993</td><td>-1.22</td></tr></tbody></table>

Increased means the percentage of relative improvements of FMF over PCA, the better results are shown in bold.

FrScatNets $y_α$ in (16) is a 4th-order tensor, which is performed by PCA to obtain an implicit vector $\mathbf{z}$. This process loses correlations between various dimensions of the data. Therefore, we consider using FMF as the dimensionality reduction method to maintain the structures of the input data better.

For the proposed FMF-based GFRSNs, the flow chart is shown in **Figure 2**. The size of the implicit tensor $Z_{α_i}$ is 1 × 49 × 4 × 4 on CIFAR-10, and for the CelebA dataset, the size of implicit tensor $Z_{α_i}$ is 1 × 65 × 8 × 8. The PSNR and SSIM on the CIFAR-10 and CelebA datasets are shown in the third columns of **Tables 2, 3**, respectively. As can be seen from the two tables, train PSNR and train SSIM of the FMF-based GFRSNs are slightly worse than those of the PCA-based GFRSNs on the CIFAR-10 and CelebA datasets; however, the test PSNR and test SSIM of the proposed FMF-based GFRSNs are better than those of the PCA-based GFRSNs. For example, Test PSNR and Test SSIM have relatively increased by 1 and 1.8%, respectively, when compared with the PCA-based GFRSNs, on the CIFAR-10 dataset. However, with regard to the CelebA dataset, Test PSNR and Test SSIM have decreased by 3.71 and 1.22%, respectively, when compared with the PCA-based GFRSNs. Nevertheless, the experimental results still show that the overfitting problem on the testing datasets can be alleviated with the FMF dimensionality reduction method.

Although the performance of the proposed FMF method on the CIFAR-10 dataset is better than that of PCA and has a similar generation ability on the CelebA dataset, more importantly, FMF has better generalization performance under the framework of GFRSNs. In other words, our generative model will not overfit on the test set. However, in order to better reflect the role of fractional scattering transformation and, hence, abolish the influence of FMF, we still use the PCA method in the following two experiments.

<footnote n="1">¹https://github.com/tomas-angles/generative-scattering-networks</footnote>
<table><caption>**TABLE 4** | Results with FrScatNets on the CIFAR-10 dataset.</caption><thead><tr><th>(α<sub>1</sub>, α<sub>2</sub>)</th><th>Fusion</th><th>Test PSNR</th><th>Increased (%)</th><th>Test SSIM</th><th>Increased (%)</th></tr></thead><tbody><tr><td>Base line with<br>(α<sub>1</sub>, α<sub>2</sub>) =<br>(1.00,1.00)</td><td>No</td><td>18.1000</td><td>0</td><td>0.8352</td><td>0</td></tr><tr><td>(0.10,1.00)</td><td>No</td><td>13.9738</td><td>-22.80</td><td>0.5442</td><td>-34.84</td></tr><tr><td></td><td>Yes</td><td>16.9597</td><td>-6.30</td><td>0.6974</td><td>-16.50</td></tr><tr><td>(0.40,1.00)</td><td>No</td><td><b>18.8280</b></td><td>4.02</td><td><b>0.8514</b></td><td><b>1.94</b></td></tr><tr><td></td><td>Yes</td><td><b>18.9869</b></td><td>4.90</td><td><b>0.8970</b></td><td><b>7.40</b></td></tr><tr><td>(0.70,1.00)</td><td>No</td><td>18.6614</td><td>3.10</td><td>0.8469</td><td>1.40</td></tr><tr><td></td><td>Yes</td><td>18.8421</td><td>4.10</td><td>0.8887</td><td>6.40</td></tr><tr><td>(1.30,1.00)</td><td>No</td><td>18.6169</td><td>2.86</td><td>0.8462</td><td>1.32</td></tr><tr><td></td><td>Yes</td><td>18.8059</td><td>3.90</td><td>0.8870</td><td>6.20</td></tr><tr><td>(1.60,1.00)</td><td>No</td><td><b>18.8209</b></td><td><b>3.98</b></td><td><b>0.8517</b></td><td><b>1.98</b></td></tr><tr><td></td><td>Yes</td><td><b>18.9688</b></td><td><b>4.80</b></td><td><b>0.8987</b></td><td><b>7.60</b></td></tr><tr><td>(1.90,1.00)</td><td>No</td><td>14.0110</td><td>-22.59</td><td>0.5474</td><td>-34.46</td></tr><tr><td></td><td>Yes</td><td>16.9959</td><td>-6.10</td><td>0.7041</td><td>-15.70</td></tr><tr><td>(1.00,0.10)</td><td>No</td><td>14.0099</td><td>-22.60</td><td>0.5498</td><td>-34.17</td></tr><tr><td></td><td>Yes</td><td>16.9054</td><td>-6.60</td><td>0.6941</td><td>-16.90</td></tr><tr><td>(1.00,0.40)</td><td>No</td><td>18.9351</td><td>4.61</td><td>0.8550</td><td>2.37</td></tr><tr><td></td><td>Yes</td><td>18.9507</td><td>4.70</td><td>0.8978</td><td>7.50</td></tr><tr><td>(1.00,0.70)</td><td>No</td><td>18.7289</td><td>3.47</td><td>0.8499</td><td>1.76</td></tr><tr><td></td><td>Yes</td><td>18.7335</td><td>3.50</td><td>0.8753</td><td>4.80</td></tr><tr><td>(1.00,1.30)</td><td>No</td><td>18.6947</td><td>3.29</td><td>0.8434</td><td>0.98</td></tr><tr><td></td><td>Yes</td><td>18.5887</td><td>2.70</td><td>0.8753</td><td>4.80</td></tr><tr><td>(1.00,1.60)</td><td>No</td><td>18.9056</td><td>4.45</td><td>0.8545</td><td>2.31</td></tr><tr><td></td><td>Yes</td><td>18.9507</td><td>4.70</td><td>0.8987</td><td>7.60</td></tr><tr><td>(1.00,1.90)</td><td>No</td><td>14.0487</td><td>-22.38</td><td>0.5520</td><td>-33.91</td></tr><tr><td></td><td>Yes</td><td>16.9778</td><td>-6.20</td><td>0.7074</td><td>-15.30</td></tr><tr><td>Fusing<br>(0.40,1.00)<br>and<br>(1.60,1.00)</td><td>Yes</td><td><b>19.1589</b></td><td><b>5.85</b></td><td><b>0.8927</b></td><td><b>6.89</b></td></tr></tbody></table>

Some better results are shown in bold. In the second column, "No" means un-fused image, and "Yes" means fused image. We also show the percentages of relative improvements on Test PSNR and Test SSIM of FrScatNets of various fractional orders (α₁, α₂) over the conventional ScatNets (the first row), respectively.

<table><caption>**TABLE 5** | Results with FrScatNets on CelebA dataset.</caption><thead><tr><th>(α<sub>1</sub>, α<sub>2</sub>)</th><th>Fusion</th><th>Test PSNR</th><th>Increased (%)</th><th>Test SSIM</th><th>Increased (%)</th></tr></thead><tbody><tr><td>Base line with<br>(α<sub>1</sub>, α<sub>2</sub>) =<br>(1.00,1.00)</td><td>No</td><td>21.1668</td><td>0</td><td>0.9221</td><td>0</td></tr><tr><td>(0.10,1.00)</td><td>No</td><td>18.3728</td><td>-13.2</td><td>0.7709</td><td>-16.4</td></tr><tr><td></td><td>Yes</td><td>21.0186</td><td>-0.7</td><td>0.9156</td><td>-0.7</td></tr><tr><td>(0.40,1.00)</td><td>No</td><td>21.4631</td><td><b>1.1</b></td><td><b>0.9350</b></td><td><b>3.3</b></td></tr><tr><td></td><td>Yes</td><td>22.2040</td><td><b>5.3</b></td><td><b>0.9608</b></td><td><b>6.5</b></td></tr><tr><td>(0.70,1.00)</td><td>No</td><td>21.3996</td><td>1.1</td><td>0.9525</td><td>2.3</td></tr><tr><td></td><td>Yes</td><td>22.3098</td><td>5.4</td><td>0.9820</td><td>6.2</td></tr><tr><td>(1.30,1.00)</td><td>No</td><td>21.3785</td><td>1</td><td>0.9433</td><td>2.4</td></tr><tr><td></td><td>Yes</td><td>22.3098</td><td>5.4</td><td>0.9793</td><td>6.2</td></tr><tr><td>(1.60,1.00)</td><td>No</td><td><b>21.4631</b></td><td><b>1.4</b></td><td><b>0.9571</b></td><td><b>3.8</b></td></tr><tr><td></td><td>Yes</td><td><b>22.3310</b></td><td><b>5.5</b></td><td><b>0.9839</b></td><td><b>6.7</b></td></tr><tr><td>(1.90,1.00)</td><td>No</td><td>18.6268</td><td>-12</td><td>0.7866</td><td>-14.7</td></tr><tr><td></td><td>Yes</td><td>21.1456</td><td>-0.1</td><td>0.9219</td><td>-0.02</td></tr><tr><td>(1.00,0.10)</td><td>No</td><td>18.2458</td><td>-13.8</td><td>0.7561</td><td>-18</td></tr><tr><td></td><td>Yes</td><td>20.9551</td><td>-1</td><td>0.9092</td><td>-1.4</td></tr><tr><td>(1.00,0.40)</td><td>No</td><td>21.5055</td><td>1.6</td><td>0.9405</td><td>2</td></tr><tr><td></td><td>Yes</td><td>22.3098</td><td>5.4</td><td>0.9756</td><td>5.8</td></tr><tr><td>(1.00,0.70)</td><td>No</td><td>21.2515</td><td>0.4</td><td>0.9249</td><td>0.3</td></tr><tr><td></td><td>Yes</td><td>22.2251</td><td>5</td><td>0.9700</td><td>5.2</td></tr><tr><td>(1.00,1.30)</td><td>No</td><td>21.2303</td><td>0.3</td><td>0.9267</td><td>0.5</td></tr><tr><td></td><td>Yes</td><td>22.2251</td><td>5</td><td>0.9581</td><td>3.9</td></tr><tr><td>(1.00,1.60)</td><td>No</td><td>21.4843</td><td>1.5</td><td>0.9433</td><td>2.3</td></tr><tr><td></td><td>Yes</td><td>22.3098</td><td>5.4</td><td>0.9765</td><td>5.9</td></tr><tr><td>(1.00,1.90)</td><td>No</td><td>18.7115</td><td>-11.6</td><td>0.7912</td><td>-14.2</td></tr><tr><td></td><td>Yes</td><td>21.1732</td><td>0.03</td><td>0.9212</td><td>-0.1</td></tr><tr><td>Fusing<br>(0.40,1.00)<br>and<br>(1.60,1.00)</td><td>Yes</td><td><b>22.0770</b></td><td><b>4.3</b></td><td><b>0.9802</b></td><td><b>6.3</b></td></tr></tbody></table>

Some better results are shown in bold. In the second column, "No" means un-fused image, and "Yes" means fused image. We also show the percentages of relative improvements on Test PSNR and Test SSIM of FrScatNets of various fractional orders (α₁, α₂) over the conventional ScatNets (the first row), respectively.

## Image Generative Results With Different Fractional Order α

In this subsection, we explore the impact of fractional order α on the quality of the generated image using the framework of GFRSNs shown in Figure 1. The other parameter settings of FrScatNets are shown in Table 1. We choose the L₁ loss function in (15) and train the generator with the Adam optimizer using the default hyperparameters.

In this subsection, we use a two-dimensional fractional Morlet wavelet to construct the FrScatNets. For the two-dimensional fractional wavelet, two fractional orders, α₁ and α₂, are needed to determine the rotational angle. The angle is defined as $θ = απ/2$, ranging from 0 to $π$; thus, the fractional orders α₁ and α₂ change from 0 to 2. To save computation time, we fix one order as 1 and the other order changes within the range 0–2 for computing the fractional scattering coefficients. The chosen values are 0.1, 0.4, 0.7, 1, 1.3, 1.6, and 1.9. The above parameter settings are

same as those in Liu et al. (2016). Note that FrScatNets reduce to conventional ScatNets when $α₁ = α₂ = 1$. The PSNR and SSIM of the generated images from FrScatNets on the CIFAR-10 and CelebA datasets are shown in Tables 4, 5.

Generally, as shown in Table 4, best results are not obtained using FrScatNets with (α₁, α₂) = (1, 1), which means that FrScatNets with some fractional order choice of (α₁, α₂) obtain better embeddings than the conventional ScatNets. For example, both the PSNR and SSIM results are very good the FrScatNets with (α₁, α₂) = (0.4, 1.00) were used and whose Test PSNR and Test SSIM increased by 4.2 and 1.9%, respectively, compared with those of the ScatNets.

For the CelebA dataset, as shown in Table 5, both the PSNR and SSIM scores in the test set are also very good when FrScatNets with (α₁, α₂) = (1.6, 1) are used. Indeed, Test PSNR and Test SSIM increased by 1.4 and 3.8%, respectively, compared with those of the ScatNets.
<figure><img src="image_6.png" /><figcaption>**FIGURE 4** | Generative images on CIFAR-10 dataset using FrScatNet embeddings. (A) Original training images; (B) generative training images using FrScatNets with $(\alpha_1, \alpha_2) = (0.4, 1)$; (C) generative training images using FrScatNets with $(\alpha_1, \alpha_2) = (1, 1)$; (D) fused training image using FrScatNets with $(\alpha_1, \alpha_2) = (0.4, 1)$ and $(\alpha_1, \alpha_2) = (1, 1)$; (E) original testing images; (F) generative testing images using FrScatNets with $(\alpha_1, \alpha_2) = (0.4, 1)$; (G) generative testing images using FrScatNets with $(\alpha_1, \alpha_2) = (1, 1)$; (H) fused testing image using FrScatNets with $(\alpha_1, \alpha_2) = (0.4, 1)$ and $(\alpha_1, \alpha_2) = (1, 1)$.</figcaption></figure>

The generative images on the CIFAR-10 dataset using FrScatNets with $(\alpha_1, \alpha_2) = (0.4, 1)$ and $(\alpha_1, \alpha_2) = (1, 1)$ are shown in **Figure 4**. The generative images on the CelebA dataset using FrScatNets with $(\alpha_1, \alpha_2) = (1.6, 1)$ and $(\alpha_1, \alpha_2) = (1, 1)$ are shown in https://mmlab.ie.cuhk.edu.hk/projects/CelebA.html.

## Image Generative Results With Image Fusion

In this subsection, we explore the impact of image fusion on the quality of the generated images using the framework of GFRSNs shown in **Figure 2**.

Since conventional GSNs are a good baseline for the framework of GFRSNs with different fractional orders $(\alpha_1, \alpha_2)$, as an example, we consider the case in which the generative images from FrScatNets with different fractional orders $(\alpha_1, \alpha_2)$, where $\alpha_1$ and $\alpha_2$ are not simultaneously equal to 1.00, are fused with the generative images from conventional ScatNets, in other words, FrScatNets with fractional orders $(\alpha_1, \alpha_2) = (1, 1)$. Since the fractional parameters can have multiple choices, naturally, we hope to explore the effect of image fusion under different fractional parameters. All the fused images are achieved using the average method shown in Equation (17), and we choose $\lambda = 0.5$. The PSNR and SSIM results of

fused images on the CIFAR-10 dataset are shown in **Table 4**, and those on the CelebA dataset are shown in **Table 5**. Note that the results are shown in the row where the “Fusion or not?” column is “Yes” in **Tables 4, 5**. As can be seen from the two tables, the results of PSNR and SSIM for the fused images are generally better than those for the unfused images from FrScatNets with different fractional orders $(\alpha_1, \alpha_2)$, where $\alpha_1$ and $\alpha_2$ are not 1 at the same time. For example, when the generative images from FrScatNets with $(\alpha_1, \alpha_2) = (0.4, 1)$ are fused with the generative images from ScatNets, the Test PSNR and Test SSIM are increased from 18.828 and 0.8514 to 18.9869 and.897, respectively, on the CIFAR-10 dataset. The results are also better than those of ScatNet-based GFRSRNs, whose Test PSNR and Test SSIM are 18.1 and 0.8352, respectively. When the generative images from FrScatNets with $(\alpha_1, \alpha_2) = (1.6, 1)$ are fused with the generative images from ScatNets, the test PSNR and test SSIM are increased from 21.4632 and 0.9571 to 22.337 and 0.9839, respectively, on the CelebA dataset. The results are also better than those of ScatNet-based GFRSRNs, whose test PSNR and test SSIM are 21.1668 and 0.944, respectively. The fused images on the CIFAR-10 dataset are shown in **Figures 4D,H** and those on the CelebA dataset are shown in https://mmlab.ie.cuhk.edu.hk/projects/CelebA.html, respectively.
We also consider the generative images from FrScatNets with fractional orders (0.4, 1) and (1.6, 1), and the results are shown in the last row of Tables 4, 5, respectively. As can be seen from the two tables, the test PSNR and test SSIM are better than the fusion results of fractional orders (1.6, 1) and (1, 1) on both the CIFAR-10 and CelebA datasets.

## The Deformation Property of the Proposed GFRSNs

In this section, we evaluate the deformation property of the proposed GFRSNs as generally done in GANs. Specifically, given two images $x_1$ and $x_2$, we modify $\beta$ to get the interpolated images:

$$x_{\beta} = G((1 - \beta)z_1 + \beta z_2), \quad \text{for } z_1 = \Phi(x_1) \text{ and } z_2 = \Phi(x_2), \qquad (18)$$

where $\Phi(\cdot)$ denotes the fixed embedding, that is, the fractional scattering transform and Gaussianization process. The results are shown in https://mmlab.ie.cuhk.edu.hk/projects/CelebA.html.

As Angles and Mallat (2018) point out, the Lipschitz continuity to deformations of the scattering network resulting in the continuous deformation from one image to another image. https://mmlab.ie.cuhk.edu.hk/projects/CelebA.html shows that the proposed GFRSNs improve the capability to extract information while maintaining the deformation properties when compared with GSNs. On the other hand, we reproduce the morphing properties of GANs without learning a discriminator.

Besides, we should note that the generated images have strong similarities with those in the training set and, thus, lead to some unrealistic results; this is partially due to the autoencoder architecture of our model. Although under the autoencoder architecture, regarding the generative model as an inverse problem of FrScatNets, can eliminate our need to train an encoder or a discriminator, however, within this supervised paradigm, the generalization ability of the model may be limited to some extent. Therefore, when we try to recover images from unknown images, the results of the model will generate images that are similar to ones in the training set.

## Comparison Results With GANs

In this section, we compared the results of the proposed GFRSNs with GANs on the CelebA dataset.

### Comparison Results With DCGAN and PGAN

We compare the visual results of the proposed GFRSNs with those of the DCGAN (Radford et al., 2016) and progressive GAN (PGAN)<fnref n="2" /> (Karras et al., 2018), as shown in https://mmlab.ie.cuhk.edu.hk/projects/CelebA.html, from which we can see that DCGAN produces a certain degree of distortion. On the contrary, the proposed GFRSNs and PGAN do not show this kind of problem. PGAN generates more image details than the proposed GFRSNs, and we think that the reasons are:

(1) The proposed GFRSNs still belong to the autoencoder architecture, which is generally inferior to that of the

GANs in terms of image generation quality. However, the autoencoder has its own merits; for example, it can obtain an image code (or a latent vector), which is very helpful for downstream tasks such as image classification. In contrast, the GANs cannot generate this latent vector.

(2) The proposed GFRSNs use learning-free FrScatNets instead of CNNs in the encoder stage, which significantly reduces the parameters (for example, reducing the parameters by half compared with DCGAN). However, it also has a certain impact on image generation performance.

(3) The proposed GFRSNs can maintain the structure of the face but show smoothed results to a certain extent. The reason for this is, maybe, the choice of $L_1$ loss.

(4) PGAN uses a more advanced low-resolution to high-resolution generation paradigm, which is more effective than the generator used in GFRSNs.

Note that we choose DCGAN as one of the compared methods, since we use the same generator architecture as the DCGAN. The reason we choose PGAN rather than the more recent BigGAN (Brock et al., 2019) as the other compared method is that the two models achieved similar results without additional class information.

## Comparison Results With CycleGAN

We compare the objective evaluation criteria (PSNR and SSIM) with CycleGAN<fnref n="3" /> (Zhu et al., 2017) on the CelebA dataset. Note that SSIM and PSNR are not suitable for evaluating the quality of GANs, since GANs, generally, generate images directly from Gaussian white noise. That is, we do not have real images corresponding to the generated images, but real images are needed to calculate the PSNR and SSIM scores.

The reason we choose CycleGAN as the compared method is that it can be seen as a special kind of autoencoder model and, hence, can be used to calculate the PSNR and SSIM scores. The structure of CycleGAN is shown in https://mmlab.ie.cuhk.edu.hk/projects/CelebA.html. As in the experiment of GFRSNs, we choose 65,536 training images and 16,384 testing images. For the training process, we divide the training set into two subsets of the same size, namely, A and B, to meet the unique circular training process. By training CycleGAN through 32,768 images in domain A and 32,768 images in domain B, we can calculate Train PSNR and Train SSIM. For the testing process, we also divide the testing set into two subsets of the same size, namely, A and B, to meet the unique circular training process. By training CycleGAN through 8,192 images in domain A and 8,192 images in domain B, we can calculate Test PSNR and Test SSIM. It can be known from the experimental process that in order to calculate the PSNR and SSIM values of the training data set and the testing data set, there are several characteristics when using CycleGAN:

(1) The training and testing processes are performed separately; that is, the trained generator of CycleGAN is not used in the testing process, since CycleGAN performs the task of image-to-image translation or style transfer (Gatys et al., 2016). In

<footnote n="2">²https://github.com/facebookresearch/pytorch_GAN_zoo</footnote>

<footnote n="3">³https://github.com/junyanz/pytorch-CycleGAN-and-pix2pix</footnote>
<table><caption>**TABLE 6** | Quantitative results of CycleGAN and the proposed GFRSNs with fractional orders $\alpha_1 = 0.4$, $\alpha_2 = 1$ on the CelebA dataset.</caption><thead><tr><th></th><th>Train PSNR</th><th>TestPSNR</th><th>Train SSIM</th><th>Test SSIM</th></tr></thead><tbody><tr><td>Cycle GAN</td><td>30.8059</td><td>32.6890</td><td>0.9824</td><td>0.9822</td></tr><tr><td>Ours</td><td>27.9721</td><td>21.4631</td><td>0.9629</td><td>0.9350</td></tr></tbody></table>

order to get the Test PSNR and Test SSIM of the testing images, we still need to train CycleGAN with the testing images. For example, as is shown in https://mmlab.ie.cuhk.edu.hk/projects/CelebA.html, the generator $G_{AB}$ takes an image from domain A, and then tries to do an image-to-image translation, so that the output will be a fake image with a style similar to domain B. However, there is only one style of images in the CelebA dataset; therefore, the generator will learn the same image as the input. That is, it is unfair to use the PSNR or SSIM score to measure the quality of CycleGAN to some extent, since CycleGAN trains the testing images.

(2) In CycleGAN, the role of the generator is not focused on generating images from noise. On the contrary, the generator takes their effort to the task of image-to-image translation. When the style of two subsets is the same, this kind of image-to-image method will undoubtedly lead to pixel-level alignment and, hence, failure of pixel error-based metrics, such as PSNR and SSIM. That is, the PSNR and SSIM scores can be seen as the upper bound of other methods.

The results of the comparison of PSNR and SSIM scores of the proposed GFRSNs with CycleGAN are shown in **Table 6**, from which we can see that the result of GFRSNs is worse than that of CycleGAN, especially on the testing set. This is not surprising, because CycleGAN implements style transfer between training data and testing data, while GFRSNs implements reconstruction from FrScatNet features to images. The PSNR and SSIM scores of CycleGAN can be seen as the upper bound of GFRSNs; that is, the proposed GFRSNs still have a lot of room for improvement.

## CONCLUSIONS

This study proposes generative fractional scattering networks (GFRSNs), which use fractional wavelet scattering networks (FrScatNets) as encoder to obtain features (or FrScatNet embeddings) and deconvolutional neural networks as decoder to generate an image. Additionally, this study develops a new feature-map fusion (FMF) method to reduce the dimensionality of FrScatNet embeddings. The impact of image fusion is also discussed in this study. The experimental results on the CIFAR-10 and CelebA datasets show that the proposed GFRSNs can lead to better generated images than the original GSNs in the testing dataset. Compared with GANs, the proposed GFRSNs

## REFERENCES

Abdi, H., and Williams, L. J. (2010). Principal component analysis. *Wiley Interdisc. Rev. Comput. Stat.* 2, 433–459. doi: 10.1002/wics.101

lack details of the generated image because of the essence of the autoencoder structure; however, the proposed GFRSNs have the following merits:

(1) They can obtain an image code (or a latent vector), which is very helpful for downstream tasks such as image classification.

(2) They use learning-free FrScatNets instead of CNNs in the encoder stage, which significantly reduces the parameters.

(3) They may have a potentially good performance in the differential privacy (DP) learning framework, since Tramer and Boneh (2021) show that ScatNet outperforms deep CNNs in differential private classifiers. We studied the image generation performance of GFRSNs under the framework of differential privacy learning. Appendix A gives some preliminary results.

## DATA AVAILABILITY STATEMENT

Publicly available datasets were analyzed in this study. This data can be found here: http://www.cs.toronto.edu/~kriz/cifar.html; http://mmlab.ie.cuhk.edu.hk/projects/CelebA.html.

## AUTHOR CONTRIBUTIONS

JW did conceptualization, methodology, writing—reviewing, and editing. XQ did validation and revised and edited the manuscript. JZ did writing—original draft preparation, software, and visualization. FW did software, validation, and data curation. YK and GY did validation and project administration. LS did formal analysis, writing—reviewing, and editing. HS did supervision and resources. All authors contributed to the article and approved the submitted version.

## FUNDING

This study was funded by the National Natural Science Foundation of China under Grant Nos: 61876037, 62171125, 31800825, 61871117, 61871124, 61773117, and 61872079 and in part by INSERM under the grant called IAL and IRP.

## ACKNOWLEDGMENTS

We thank the Big Data Computing Center of Southeast University for providing facility support on the numerical calculations in this study.

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnbot.2021.752752/full#supplementary-material

Angles, T., and Mallat, S. (2018). “Generative networks as inverse problems with scattering transforms,” in *2018 International Conference on Learning Representations* (ICLR) (Vancouver, BC).
Arjovsky, M., Chintala, S., and Bottou, L. (2017). "Wasserstein generative adversarial networks," in *2017 International Conference on Machine Learning (ICML)* (Sydney, NSW), 214–223.

Blei, D. M., Kucukelbir, A., and McAuliffe, J. D. (2017). Variational inference: a review for statisticians. *J. Am. Stat. Assoc.* 112, 859–877. doi: 10.1080/01621459.2017.1285773

Brock, A., Donahue, J., and Simonyan, K. (2019). "Large scale GAN training for high fidelity natural image synthesis," in *2019 International Conference on Learning Representations (ICLR)* (New Orleans, LA).

Bruna, J., and Mallat, S. (2013). Invariant scattering convolution networks. *IEEE Trans. Pattern Anal. Mach. Intell.* 35, 1872–1886. doi: 10.1109/TPAMI.2012.230

Chen, X., Duan, Y., Houthooft, R., Schulman, J., Sutskever, I., and Abbeel, P. (2016). "InfoGAN: interpretable representation learning by information maximizing generative adversarial nets," in *2016 Proceedings of the 30th International Conference on Neural Information Processing Systems, (NeurIPS)* (Barcelona), 29, 2180–2188.

Gatys, L. A., Ecker, A. S., and Bethge, M. (2016). "Image style transfer using convolutional neural networks," in *IEEE Conference on Computer Vision and Pattern Recognition (CVPR)* (Las Vegas, NV), 2414–2423.

Goodfellow, I., Pouget-Abadie, J., Mirza, M., Xu, B., Warde-Farley, D., Ozair, S., et al. (2014). Generative adversarial nets. *Adv. Neural Inf. Process. Syst.* 27, 2672–2680. doi: 10.3156/JSOFT.29.5_177_2

Gregor, K., Papamakarios, G., Besse, F., Buesing, L., and Weber, T. (2018). "Temporal difference variational auto-encoder," in *2018 International Conference on Learning Representations (ICLR)* (Vancouver, BC).

Gulrajani, I., Ahmed, F., Arjovsky, M., Dumoulin, V., and Courville, A. (2017). "Improved training of wasserstein GANs," in *2017 Proceedings of the 31st International Conference on Neural Information Processing Systems (NeurIPS)* (Long Beach, CA) 30, 5769–5779.

Higgins, I., Matthey, L., Pal, A., Burgess, C., Glorot, X., Botvinick, M., et al. (2017). "beta-VAE: learning basic visual concepts with a constrained variational framework," in *2017 International Conference on Learning Representations (ICLR)* (Toulon).

Hudson, D. A., and Zitnick, L. (2021). "Generative adversarial transformers," in *Proceedings of the 38th International Conference on Machine Learning (ICML)* (Virtual), 4487–4499.

Ioffe, S., and Szegedy, C. (2015). "Batch normalization: accelerating deep network training by reducing internal covariate shift," in *Proceedings of The 32nd International Conference on Machine Learning (ICML)* (Lille) 1, 448–456.

Karras, T., Aila, T., Laine, S., and Lehtinen, J. (2018). "Progressive growing of GANs for improved quality, stability, and variation," in *2018 International Conference on Learning Representations (ICLR)* (Vancouver, BC).

Karras, T., Laine, S., and Aila, T. (2020a). A style-based generator architecture for generative adversarial networks. *IEEE Trans. Pattern Anal. Mach. Intell.* 1:1. doi: 10.1109/TPAMI.2020.2970919

Karras, T., Laine, S., Aittala, M., Hellsten, J., Lehtinen, J., and Aila, T. (2020b). "Analyzing and improving the image quality of StyleGAN," in *2020 IEEE Conference on Computer Vision and Pattern Recognition (CVPR)* (Virtual), 8110–8119.

Kingma, D. P., and Ba, J. L. (2015). "Adam: a method for stochastic optimization," in *2015 International Conference on Learning Representations (ICML)* (Lille).

Kingma, D. P., and Welling, M. (2014). "Auto-encoding variational Bayes," in *2014 International Conference on Learning Representations (ICLR)* (Banff).

Krizhevsky, A. (2009). Learning Multiple Layers of Features from Tiny Images.

Liu, L., Wu, J., Li, D., Senhadji, L., and Shu, H. (2019). Fractional wavelet scattering network and applications. *IEEE Trans. Biomed. Eng.* 66, 553–563. doi: 10.1109/TBME.2018.2850356

Liu, Y., Chen, X., Ward, R. K., and Wang, Z. J. (2016). Image fusion with convolutional sparse representation. *IEEE Signal Process. Lett.* 23, 1882–1886. doi: 10.1109/LSP.2016.2618776

Liu, Z., Luo, P., Wang, X., and Tang, X. (2015). "Deep learning face attributes in the wild," in *2015 IEEE International Conference on Computer Vision (ICCV) (Santiago),* 3730–3738.

Long, J., Shelhamer, E., and Darrell, T. (2015). "Fully convolutional networks for semantic segmentation," in *2015 IEEE Conference on Computer Vision and Pattern Recognition (CVPR)*, 3431–3440.

Mallat, S. (2012). Group invariant scattering. *Commun. Pure Appl. Math.* 65, 1331–1398. doi: 10.1002/cpa.21413

Mirza, M., and Osindero, S. (2014). Conditional generative adversarial nets. *arXiv.*

Miyato, T., Kataoka, T., Koyama, M., and Yoshida, Y. (2018). "Spectral normalization for generative adversarial networks," in *2018 International Conference on Learning Representations (ICLR)* (Vancouver, BC).

Nair, V., and Hinton, G. E. (2010). "Rectified linear units improve restricted boltzmann machines," in *Proceedings of the 27th International Conference on Machine Learning (ICML)* (Haifa), 807–814.

Nowozin, S., Cseke, B., and Tomioka, R. (2016). "f-GAN: training generative neural samplers using variational divergence minimization," in *2016 Proceedings of the 30th International Conference on Neural Information Processing Systems (NeurIPS)* (Barcelona) 29, 271–279.

Oord, A., van den, Vinyals, O., and kavukcuoglu, koray. (2017). "Neural discrete representation learning," in *2017 Proceedings of the 31st International Conference on Neural Information Processing Systems (NeurIPS)* (Long Beach, CA) 30, 6306–6315.

Qi, G.-J. (2020). Loss-sensitive generative adversarial networks on lipschitz densities. *Int. J. Comput. Vis.* 128, 1118–1140. doi: 10.1007/s11263-019-01265-2

Radford, A., Metz, L., and Chintala, S. (2016). "Unsupervised representation learning with deep convolutional generative adversarial networks," in *2016 International Conference on Learning Representations (ICLR)* (San Juan).

Razavi, A., Oord, A., van den, and Vinyals, O. (2019). "Generating diverse high-fidelity images with VQ-VAE-2," in *2019 Proceedings of the 33rd International Conference on Neural Information Processing Systems (NeurIPS)* (Vancouver, BC) 32, 14837–14847.

Rezende, D., and Mohamed, S. (2015). "Variational inference with normalizing flows," in *2015 Proceedings of the 32nd International Conference on Machine Learning (ICML)* (Lille), 1530–1538.

Rezende, D. J., Mohamed, S., and Wierstra, D. (2014). "Stochastic backpropagation and approximate inference in deep generative models," in *2014 Proceedings of the 31st International Conference on Machine Learning (ICML)* (Beijing), 1278–1286.

Rusu, A. A., Rabinowitz, N. C., Desjardins, G., Soyer, H., Kirkpatrick, J., Kavukcuoglu, K., et al. (2017). *Progressive neural networks*. *arXiv.*

Salimans, T., Kingma, D., and Welling, M. (2015). "Markov Chain Monte Carlo and variational inference: bridging the gap." in *Proceedings of the 32nd International Conference on Machine Learning (ICML)* (Lille) 37, 1218–1226.

Shi, J., Chi, Y., and Zhang, N. (2010). Multichannel sampling and reconstruction of bandlimited signals in fractional fourier domain. *IEEE Signal Process. Lett.* 17, 909–912. doi: 10.1109/LSP.2010.2071383

Simonovsky, M., and Komodakis, N. (2018). "GraphVAE: towards generation of small graphs using variational autoencoders," in *27th International Conference on Artificial Neural Networks (ICANN) (Rhodes)*, 412–422.

Sohn, K., Yan, X., and Lee, H. (2015). "Learning structured output representation using deep conditional generative models," in *2015 Proceedings of the 29th International Conference on Neural Information Processing Systems, (NeurIPS) (Montreal, QC) 28, 3483–3491.

Sonderby, C. K., Raiko, T., Maaløe, L., Sonderby, S. K., and Winther, O. (2016). "Ladder variational autoencoders," in *2016 30th Annual Conference on Neural Information Processing Systems Conference (NeurIPS)* (Barcelona) 29, 3738–3746.

Tramer, F., and Boneh, D. (2021). "Differentially private learning needs better features (or much more data)," in *2021 International Conference on Learning Representations (ICLR)* (Virtual).

Wang, Y., Li, J., Yi, L., Yao, F., and J., Q. (2003). "Image quality evaluation based on image weighted separating block peak signal to noise ratio," in *Proceedings of the 2003 International Conference on Neural Networks and Signal Processing (Nanjing) 2*, 994–97.

Wang, Z., Bovik, A. C., Sheikh, H. R., and Simoncelli, E. P. (2004). Image quality assessment: from error visibility to structural similarity. *IEEE Trans. Image Proc.* 13, 600–612. doi: 10.1109/TIP.2003.819861

Wang, Z., She, Q., and Ward, T. E. (2021). Generative adversarial networks in computer vision: a survey and taxonomy. *ACM Comput. Surveys* 54, 1–38. doi: 10.1145/3459992

Yang, B., Zhong, J., Li, Y., and Chen, Z. (2017). Multi-focus image fusion and super-resolution with convolutional neural network. *Int. J. Wavelets Multiresol. Inform. Proc.* 15:1750037. doi: 10.1142/S0219691317500370
Zhao, H., Gallo, O., Frosio, I., and Kautz, J. (2017). Loss functions for image restoration with neural networks. *IEEE Transac. Comput. Imaging* 3, 47-57. doi: 10.1109/TCI.2016.2644865

Zhu, J. Y., Park, T., Isola, P., and Efros, A. A. (2017). "Unpaired image-to-image translation using cycle-consistent adversarial networks," in *2017 IEEE International Conference on Computer Vision (ICCV)*, 2242-2251.

Zou, D., and Lerman, G. (2019). "Encoding robust representation for graph generation," in *2019 International Joint Conference on Neural Networks (IJCNN)* (Budapest), 1-9.

**Conflict of Interest:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

**Publisher's Note:** All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Copyright © 2021 Wu, Qiu, Zhang, Wu, Kong, Yang, Senhadji and Shu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.