Research Article  Open Access
Identification of Weakly PitchShifted Voice Based on Convolutional Neural Network
Abstract
Pitch shifting is a common voice editing technique in which the original pitch of a digital voice is raised or lowered. It is likely to be abused by the malicious attacker to conceal his/her true identity. Existing forensic detection methods are no longer effective for weakly pitchshifted voice. In this paper, we proposed a convolutional neural network (CNN) to detect not only strongly pitchshifted voice but also weakly pitchshifted voice of which the shifting factor is less than ±4 semitones. Specifically, linear frequency cepstral coefficients (LFCC) computed from power spectrums are considered and their dynamic coefficients are extracted as the discriminative features. And the CNN model is carefully designed with particular attention to the input feature map, the activation function and the network topology. We evaluated the algorithm on voices from two datasets with three pitch shifting software. Extensive results show that the algorithm achieves high detection rates for both binary and multiple classifications.
1. Introduction
Voice disguising [1] is commonly used in forensic scenario as an effective mean of concealing the identity of the speaker. And it can be divided into two categories, nonelectronic disguising and electronic disguising. Nonelectronic disguising voice is usually obtained by pinching the nose, covering the mouth, pulling the check, etc., which is easy to be noticed under human supervision. Electronic disguising is achieved by using electronic devices or software to modify the voice pitch and format.
The simplest way of electronic disguising is to change the playback speed of the target voice. Although the speaker’s identity could be concealed, the rhythm of the disguised voice generated in this way is relatively unnatural and is not often adopted by the attackers in practice. Pitch shifting is a typical electronic disguising technique in which the pitch of the voice is changed while keeping the duration unchanged. Generally, the pitchshifted voice is more natural in terms of timbre, tone, etc., and difficult to be detected. In this paper, we mainly focus on identification of pitchshifted voices.
Clark [2] studied the ability of the human to distinguish the electronic disguised voice, and quantitatively analyzed the different effects of the different pitched voice on the human hearing. Wu et al. [3–5] studied the mechanism of pitch shifting and constructed a pitch shifting dataset with various voice software/tools. The final detection accuracy of their method can reach up to 90%, while keeps the false alarm rate less than 10%. However, the performance on weakly pitchshifted voices is relatively poor. Especially for the voices shifted with semitones, the detection rates drop lower than 90%. In [6], environment noise is considered in identifying the pitch shifting. The experimental results show that the features extracted from linear frequency cepstral coefficients (LFCC) and formant can be effectively discriminant the natural and pitchshifted voice. However, the experimental results on weakly pitchshifted voices have not been given in [6].
Recently, some studies on the detection of weakly pitchshifted voices have been reported. Based on [5], Liang et al. [7] focused on voice with the shifting factors of semitones, but the promotion is limited. Singh [8] compared performance of different classifiers on the voice shifted with semitones from to . However, the result performed on a dataset with dozens of voice samples is not consistent.
Convolutional Neural Networks (CNN) [9] have achieved stateoftheart performance in computer vision, data mining, as well as automatic speaker verification. And CNN have been adopted to audio forensics as well [10, 11]. Chen et al. [12] identified various audio post processing operations by a CNN. Especially for small size voice samples, the network achieves significant improvement comparing with other works. In [13], unlike other handcrafted features, a CNN is adopted to capture the steganographic modifications adaptively and outperform the traditional methods.
Although many methods have been proposed for pitch shifting identification, there is still room to improve the performance especially when the suspected voices are weaklyshifted. In this paper, a CNN model for pitch shifting detection is proposed. By analyzing the principle of voice pitch shifting, LFCC and the first derivative coefficients are used as identification features. Comparing to other related works, the proposed CNN achieves remarkable performance in both binary and multiple classifications. The main contributions of our work are summarized as follows.(i)High accuracy is achieved on identifying weakly pitchshifted voice. Since the difference between the original voice and the weakly pitchshifted voice is little, the identification is a challenging task in previous work.(ii)Utilizing CNN architecture to identify the pitchshifting voice, which improve the performance compared to the previous work. And the proposed network architecture is carefully devised.(iii)Massive experiments are conducted on two dataset and three pitch shifting software, which indicates the proposed method achieved great robustness.
The remainder of the paper is organized as follows. In Section 2, we briefly introduce the principle of voice pitch shifting. Section 3 presents the identification features and describes the proposed CNN topology. In Section 4, a series of experiment results are given. Finally, the paper is concluded in Section 5.
2. Voice Pitch Shifting
Voice pitch shifting can be performed in either timedomain or frequency domain. Timedomain Pitch Synchronous Overlap Add (TDPSOLA) is a commonly used approach which works by windowing [14]. Upsampling achieves pitch shifting by moving the segments further apart and downsampling achieves by moving closer together. Upsampling can achieve the compression of the spectrum, which lowers the pitch. Downsampling can achieve the expansion of the spectrum, thus raise the pitch. In real scenarios, more stateofart voice synthesis algorithms are applied in audio editing software. These algorithms have better performance in timbre and rhythm. In our work, Audition [15], GoldWave [16] and Audacity [17] are considered as pitch shifting methods.
In this paper, we use semitone to measure the pitch of shifted voice. A semitone is the smallest interval between two tones. It is defined as the interval between two adjacent notes in a 12tone scale [18], which means the frequency between two adjacent semitones has an equal ratio of . In other words, if the voice frequency is raised or lowered by times, the pitch can be raised or lowered by one semitone. Let be the frequency of original voice, and the frequency of pitchshifted voice is given by the following formula
where represents the semitones of pitchshifted voice compared to original one. A positive means raising the pitch of voice and a negative one means lowering the pitch of voice. In this paper, we use as a shifting factor which denotes the pitchshifted voice.
3. Identification Algorithm Based on CNN
3.1. Feature Extraction
We randomly choose a voice sample from the TIMIT [19] dataset and shift the voice by setting in Equation (1) to and respectively. The waveform and spectrogram of original and pitchshifted voice are shown in Figure 1. As we can see, the shifting operation changes the waveform little while leaves traces on the frequency domain. Thus, acoustic features which characterize frequency domain can be applied to the proposed algorithm.
(a)
(b)
LFCC is a cepstral feature widely used in voice identification and achieves significantly performance [20]. Recent works [21] show that LFCC can more effectively captures the lower as well as higher frequency characteristics than other cepstral coefficients. Hence, in this work, LFCC is considered to extract the identification feature. The extraction procedure of LFCC is as follows.
The voice signal is first preprocessed with preemphasized and then windowed. Let be the preprocessed voice signal and , where is the duration of the signal. Suppose the frequency spectrum of the th voice frame is calculated by shorttime Fourier transform (STFT), refers to the th spectrum. Then the power spectrum filtered by a set of linearlyspaced triangular filters can be defined by
where is the number of filters and is the number of frames in a voice sample. is defined as
where , and are the lowest frequency, central frequency, and the highest frequency of th filter, respectively. The adjacent filters have .
Finally, the DCT is applied to the Logpower of the filters to calculate the LFCC of
where is LFCC of th frame, and is the index of DCT coefficients.
Since most of the pitch shifting techniques do not fully model temporal characteristics of voice [22], the dynamic coefficients, such as the first and second derivatives, could be useful in identifying pitchshifted voice. In this work, we take the first derivative into consideration and it could be given by
is the first derivative coefficient of th frame, which computed in term of the static coefficients to . A typical value for is 2.
3.2. Proposed CNN Architecture
3.2.1. Network Topology
Convolution neural networks have shown remarkable performance in various classification tasks. It generally consists of an input layer, multiple hidden layers and an output layer. The hidden layers are crucial to the network performance, which typically are combination of different kinds of layers such as convolutional layers, pooling layers and full connected layers [9].
The proposed network architecture is shown in Figure 2. The input of the network is the matrix, and the output is a predict label, which indicates the suspected voice is pitch shifted or not. The entire network consists of three convolutional groups, a fully connected layer and a softmax layer. In the training stage, after extracting features of voice segments, the feature matrix is fed into the network. The specific size of matrix depends on length of each frame and number of filters. Then it undergoes three convolutional groups which are stacked one after another. Next, the feature map of the last convolutional group is fed into the fully connected layer. All the weight values in the network will be updated via back propagation. The testing stage is mostly as same as the training stage. The feature matrix of the suspected voice is first extracted and undergoes the whole network. A softmax is used as the classifier at the end of network.
3.2.2. Convolutional Group
In our network, each convolutional group includes two convolutional layers and a pooling layer. The convolution layer consists of a set of linear convolutional filters which can generate local feature maps. Twodimensional convolutional layer preforms a convolution on the input feature map with a specific kernel size. Let be the input feature map of the th neuron at layer , output feature map is computed as
where is the output map of the th neuron at layer , and is the weight value between the th neuron at layer and the th neuron at the previous layer . All convolutional layers use the same kernel size and number of stride (5 ^{∗} 5 size, 1 ^{∗} 1 stride). Since the feature map is a twodimensional matrix, the first convolutional layer in the first group has one input channel and 64 output channels, while the other convolutional layers have both input channels and output channels with number of 64. Nonlinear activation functions can enhance the mapping capacity of the model by introducing nonlinearity into the network.
Pooling layers are adopted after convolutional layers which can obtain more global information by combining the feature information extracted from the convolution layer. Max pooling is commonly used in the pooling layer. It is a downsampling operation, which chooses the maximum value within a local window is taken as the output
where is the pooling region in feature map. The region is defined by pool size and number of strides. Pooling layers reduce the number of parameters in the network significantly and have little effect on input feature map, thus decrease the computational cost and prevent overfitting. All maxpooling layers use the same pool size and number of stride (2 ^{∗} 2 size, 2 ^{∗} 2 stride).
3.2.3. Rest Part of Network
After three convolutional groups, the fully connected layer acts as a “classification” map in the network, which can do the highlevel reasoning and learn distributed feature representation. Neurons in fully connected (FC) layer are connected to all activation functions in the previous layer. However, overly complex networks will reduce the generalization of the model. Dropout is a simple and effective regularization technique to prevent overfitting [23]. Hence, in our network, we drop out half of input neurons in the FC layer.
Softmax can be considered as an effective multipleoutput competitive whose output represents the likelihood of classification. Therefore, the dimension of its output represents the number of classes. Let be the number of classes, the probabilities of input data over different classes are predicted by the softmax function
where is the output of the FC layer on each class. Finally, the predicted label depends on the largest probability .
In summary, the architecture and parameters of the proposed network are shown in Table 1.
 
^{1} depends on specific the number of classes. 
3.3. Proposed Identification Algorithm for PitchShifted Voice
The proposed identification algorithm is based on the first derivative of LFCC and CNN classifier. With a group of equaling distributed triangular filters, LFCC can capture more characteristics both in low frequency and high frequency comparing with other acoustics features such as MFCC. Thus, the difference between the original voice and the pitchshifting voice are easier to be distinguished. CNN is considered to have better performance in classification task for multilayers process with less time and subsampling layers give better feature extraction. The proposed algorithm consists of training and testing stages, as shown in Figure 3.
In the training stage, the voice pitchshifted different factors and the original voice are considered as separate classes. After extracting the first derivative of LFCC based on Equation (5), feature map together with labels are fed into the network for training.
In the testing stage, the first derivative of LFCC are first extracted and then fed into the trained CNN model. The probability given by softmax in Equation (8) reveals the voice is more likely to be the original one or shifted with which semitone.
4. Results and Discussion
4.1. Experiment Setup
In the experiments, the proposed algorithm is evaluated on TIMIT [19] and UME [24]. TIMIT consists of 6300 voice samples from 630 speakers with the average duration of 3 s. And it is turned into three different subdatasets using Audition, GoldWave, and Audacity respectively, each of which contains sixteen shifting factors from semitones to semitones. Hence, there are totally 100800 voice samples in each subdataset of TIMIT. Similarly, UME consists of 4040 voice samples from 202 speakers with the average duration of 5 s. TIMIT and UME are turned into three subdatasets respectively, each of which composed of 64640 voice samples. In each subdataset, 60% of voice samples are selected randomly into training dataset, 20% sample into validation dataset and the remaining 20% sample into testing dataset. Speaker identity is not considered while splitting, and two datasets are from different speakers. Thus, the datasets are supposed to be speaker independent. Those voice samples with the shifting factor less than semitones are considered as weakly pitchshifted, while others are strongly pitchshifted. All the voice samples from both datasets are WAV, 16 KHz sampling rate, 16bit quantization and mono.
For each voice sample, 20dimensional LFCC feature map is extracted by setting the length of frame to 256 and the number of filters to 20 in Equation (2). In [6], LFCC with SVM classifier achieves great robustness detecting disguised voice in noisy environment. In our work, the GMM classifier is used as a comparison, among which the number of GMM kernels is set to 256.
The detection rate is used to evaluate the performance of the proposed network. Let be the number of pitchshifted voice samples and be the number of original voice samples. Assuming that and are the voice samples from pitchshifted voices and original voices which are identified as pitchshifted. The detection rate is defined as . Meanwhile, a false alarm is the most serious of the voiceprint authentication system errors to some extent. Therefore, in addition to using the detection rate to assess the proposed algorithm, we also considered the False Alarm Rate (FAR) is the testing stage. The FAR is defined as .
4.2. CNN Training
In this paper, TanH is utilized as activation function in the proposed network. We use Adam algorithm [25] with an initial learning rate of 0.0001 to accelerate the training. The proposed network is trained for 2000 iterations with the batch size of 32. The training process is presented in Figure 4, which shows the proposed network is neither overfitting nor underfitting.
(a)
(b)
tDistributed Stochastic Neighbor Embedding (tSNE) is a dimensionality reduction method which tries to place the objects in a lowdimensional space so as to optimally preserve neighbourhood identity. And it is particularly suitable for the visualization of highdimensional data [26] such as the output feature maps of convolutional layers.
We randomly choose 100 voice samples from each subdataset of TIMIT which shifted with shifting factors from semitones to semitones by Audition. Each sample is fed into the trained the network respectively, and the output feature maps of convolution layers are recorded. Figure 5 shows the visualization results of four feature maps using tSNE. The process from Figures 5(a) to 5(d) demonstrates the proposed network can capture the difference between the original voice and voice pitchshifted with different factors. In Figure 5(a), all voice samples are mixed together at first, which indicates that the characteristic represented by first derivative of LFCC is more related to voice itself rather than pitchshifting factors. In Figure 5(d), samples from same class clustered well, which indicated that the trained network can achieve both binary and multiple classifications.
(a)
(b)
(c)
(d)
4.3. Strongly PitchShifted
In this case, as a comparison to [6] and [8], we focus on the voice strongly shifted with factors from ±5 to ±8 semitones. Firstly, we try to identify whether the suspected voice is original or pitchshifted one. All the pitchshifted voice (shifted ±5 to ±8 semitones) are taken as negative samples in binary classification. In real forensic scenarios, the pitchshifted voice can be recorded by variety of devices in different environments. Hence, crossdataset experiments are necessary and important. The detection rates and FARs of this case are presented in Table 2.
 
Bold values represent the best performance under same circumstances (in same row) of the three methods. For criteria detection rate (Rate), higher is better. For criteria false alarm rate (FAR), lower is better. 
It can be seen that, all the detection methods achieve a detection rate higher than 95% and FAR lower than 2%. The method in [6] performs best in binary classification, for it achieves the highest detection rate and lowest FAR in most cases. Although the proposed method does not perform as well as [6] and [8], the gaps in both detection rates and FARs are less than 1%. These minor differences may have little effect on the detection performance.
Compared with binary classification, multiple classification is more practical for real forensic application. In this case, we not only recognize whether the suspected voice is pitchshifted, but also determine the specific shifting factor. The results are presented in Figure 6. First, as we can see from Figure 6, the detection rates of voice shifted with negative factors are higher than those with positive factors. The main reason for this phenomenon is that, downsampling (raising the pitch) will amplify the spectrum which brings more noise, while upsampling compress the spectrum. Second, different pitch shifting software have an impact on detection performance. The proposed method remains generally steady while others fluctuate greatly. Finally, the detection rates drop obviously in the crossdataset evaluation, especially for a few specific semitones are lower than 50% in [6] and [8]. And it can be seen that; the detection rates of proposed method remain higher than 60% in every case when crossing training set and testing set. Hence, for those strongly pitchshifted voice, compared with exist methods, the proposed method achieves generally the same the performance in binary classification and show more generalization ability in multiple classification.
4.4. Weakly PitchShifted
In this case, we focus on weakly pitchshifted samples shifted from to semitones which are more challenging to detect. Like Section 4.2, the binary classification is evaluated first as using all the pitchshifted voice as negative samples. The detection rates and FARs are shown in Table 3. Compared with those strongly pitchshifted voice, performance of all detection methods dropped. However, unlike Table 2, the proposed method performs best in Table 3. It achieves the highest detection rate and lowest FAR in most cases. Though the performance drops a little in intradataset, the proposed method achieves a significant improvement in crossdataset evaluation. The detection rates remain higher than 93% in every case while others drop lower than 88%. This phenomenon can be attributed to the factor that, both LFCC and MFCC mainly focus on the static features which are more related to the voice characteristic, while captures dynamic features which are more related to the shifting trace.

Like the previous section, multiple classification is adopted after the binary evaluation. The result show in Figure 7 reveals the proposed method performance on weakly pitchshifted voice form to semitones.
Generally, in Figure 7, as the same trend shown in Figure 6, raising the pitch is still difficult to detect compared with lowering the pitch. And it is noted that the fluctuation on detection rates when using different pitch shifting software is still unavoidable. The first row and the last row in Figure 7 indicate the intradataset results, the detection rates of proposed method are higher than 90% in most cases, while others are greatly affected by different pitch shifting software and even drop lower than 60%. The 2^{nd} and 3^{rd} rows show the crossdataset results, especially for a few specific semitones, both [6] and [8] lower than 20%. Proposed methods remain a steady performance with the worst case of ~60% and ~80% for most cases.
Hence, both binary and multiple classifications show that the proposed algorithm achieves good performance and has strong robustness in detecting weakly pitchshifted voice.
5. Conclusions
In this paper, an algorithm for pitchshifted voice identification is proposed. A convolutional neural network architecture is designed and adopted as the classifier to detect the pitchshifted voice while linear frequency cepstral coefficients are extracted as acoustic features. The algorithm is evaluated on two datasets and three audio editing software. Extensive results indicate that the proposed algorithm achieves much better detection rates and FARs in most cases, and the proposed network shows better generalization ability comparing to traditional classifier such as GMM. Next, network architecture which can replace handcrafted acoustic features is also one of the directions worth studying.
Data Availability
The open source databases used in this work have been listed in the reference.
Conflicts of Interest
The authors declare that they have no conflicts of interest regarding the publication of this paper.
Acknowledgments
This research was funded by the National Natural Science Foundation of China, grant numbers [61300055, 61672302]; Natural Science Foundation of Zhejiang, grant number [LY17F020010, LY20F020010]; Natural Science Foundation of Ningbo, grant number [2017A610123] and Zhejiang College Students Science and Technology Innovation Training Program, grant number [2018R405033].
References
 P. Perrot, G. Aversano, and G. Chollet, “Voice disguise and automatic detection: review and perspectives,” in Progress in Nonlinear Speech Processing (Lecture Notes in Computer Science), pp. 101–117, SpringerVerlag, New York, NY, USA, 2007. View at: Publisher Site  Google Scholar
 J. Clark and P. Foulkes, “Identification of voices in electronically disguised speech,” International Journal of Speech, Language and the Law, vol. 14, no. 2, pp. 195–221, 2007. View at: Publisher Site  Google Scholar
 Y. Wang, Y. Deng, H. Wu, and J. Huang, Blind Detection of Electronic Voice Transformation with Natural Disguise, Springer, Berlin, Heidelberg, 2013.
 H. Wu, Y. Wang, and J. Huang, “Blind detection of electronic disguised voice,” in 2013 IEEE International Conference on Acoustics, Speech and Signal Processing, pp. 3013–3017, IEEE, Vancouver, BC, Canada, 2013. View at: Publisher Site  Google Scholar
 H. Wu, Y. Wang, and J. Huang, “Identification of electronic disguised voices,” IEEE Transactions on Information Forensics and Security, vol. 9, no. 3, pp. 489–500, 2014. View at: Publisher Site  Google Scholar
 W. Cao, H. Wang, H. Zhao, Q. Qian, and S. M. Abdullahi, “Identification of electronic disguised voices in the noisy environment,” in Digital Forensics and Watermarking. IWDW 2016. Lecture Notes in Computer Science, Y. Shi, H. Kim, F. PerezGonzalez, and F. Liu, Eds., vol. 10082, pp. 75–87, Springer, Cham, 2017. View at: Publisher Site  Google Scholar
 H. Liang, X. Lin, Q. Zhang, and X. Kang, “Recognition of spoofed voice using convolutional neural networks,” in 2017 IEEE Global Conference on Signal and Information Processing (GlobalSIP), pp. 293–297, IEEE, 2017. View at: Google Scholar
 M. K. Singh, A. K. Singh, and N. Singh, “Multimedia analysis for disguised voice and classification efficiency,” Multimedia Tools and Applications, pp. 1–17, 2018. View at: Publisher Site  Google Scholar
 Y. LeCun, L. Bottou, Y. Bengio, and P. Haffner, “Gradientbased learning applied to document recognition,” Proceedings of the IEEE, vol. 86, no. 11, pp. 2278–2324, 1998. View at: Publisher Site  Google Scholar
 H. Ali, S. N. Tran, E. Benetos, and A. S. d’Avila Garcez, “Speaker recognition with hybrid features from a deep belief network,” Neural Computing Applications, vol. 29, no. 6, pp. 13–19, 2018. View at: Publisher Site  Google Scholar
 K. Sundararajan and D. L. Woodard, “Deep learning for biometrics: a survey,” ACM Computing Surveys, vol. 51, no. 3, p. 65, 2018. View at: Publisher Site  Google Scholar
 B. Chen, W. Luo, and D. Luo, “Identification of audio processing operations based on convolutional neural network,” in Proceedings of the 6th ACM Workshop on Information Hiding and Multimedia Security — IH&MMSec’18, pp. 73–77, ACM, 2018. View at: Google Scholar
 B. Chen, W. Luo, and H. Li, “Audio steganalysis with convolutional neural network,” in Proceedings of the 5th ACM Workshop on Information Hiding and Multimedia Security — IHMMSec’17, pp. 85–90, ACM, 2017. View at: Google Scholar
 F. Charpentier and M. Stella, “Diphone synthesis using an overlapadd technique for speech waveforms concatenation,” in ICASSP’86. IEEE International Conference on Acoustics, Speech, and Signal Processing, vol. 11, pp. 2015–2018, 1986. View at: Publisher Site  Google Scholar
 “Adobe Audition CC. A professional audio workstation,” Febraury 2019, https://www.adobe.com/products/audition.html. View at: Google Scholar
 “GoldWave — audio editor, recorder, converter, restoration, & analysis software,” Febraury 2019, http://www.goldwave.ca/. View at: Google Scholar
 “Audacity: free audio editor and recorder,” Febraury 2019, https://www.audacityteam.org/. View at: Google Scholar
 S. Trehub, A. Cohen, L. Thorpe, and B. Morrongiello, “Development of the perception of musical relations: semitone and diatonic structure,” Journal of Experimental Psychology: Human Perception Performance, vol. 12, no. 3, pp. 295–301, 1986. View at: Publisher Site  Google Scholar
 “Timit AcousticPhonetic Continuous Speech Corpus,” Febraury 2019, https://catalog.ldc.upenn.edu/LDC93S1. View at: Google Scholar
 F. Alegre, R. Vipperla, A. Amehraye, and N. Evans, “A new speaker verification spoofing countermeasure based on local binary patterns,” in INTERSPEECH 2013, 14th Annual Conference of the International Speech Communication Association, p. 5, 2013. View at: Google Scholar
 M. Sahidullah, T. Kinnunen, and C. Hanilçi, “A comparison of features for synthetic speech detection,” in INTERSPEECH 2015, 16th Annual Conference of the International Speech Communication Association, pp. 2087–2091, 2015. View at: Google Scholar
 F. Soong and A. Rosenberg, “On the use of instantaneous and transitional spectral information in speaker recognition,” IEEE Transactions on Acoustics, Speech, Signal Processing, vol. 36, no. 6, pp. 871–879, 1988. View at: Publisher Site  Google Scholar
 N. Srivastava, G. Hinton, A. Krizhevsky, I. Sutskever, and R. Salakhutdinov, “Dropout: a simple way to prevent neural networks from overfitting,” The Journal of Machine Learning Research, vol. 15, no. 1, pp. 1929–1958, 2014. View at: Google Scholar
 “Advanced utilization of multimedia to promote higher education reform speech database,” February 2019, http://research.nii.ac.jp/src/en/UMEERJ.html. View at: Google Scholar
 D. P. Kingma and J. Ba, “Adam a method for stochastic optimization,” 2014, https://arxiv.org/abs/1412.6980. View at: Google Scholar
 L. V. D. Maaten and G. Hinton, “Visualizing data using tSNE,” Journal of Machine Learning Research, vol. 9, pp. 2579–2605, 2008. View at: Google Scholar
Copyright
Copyright © 2020 Yongchao Ye et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.