BREAKING NEWS
S transform

## Summary

S transform as a time–frequency distribution was developed in 1994 for analyzing geophysics data.[1][2] In this way, the S transform is a generalization of the short-time Fourier transform (STFT), extending the continuous wavelet transform and overcoming some of its disadvantages. For one, modulation sinusoids are fixed with respect to the time axis; this localizes the scalable Gaussian window dilations and translations in S transform. Moreover, the S transform doesn't have a cross-term problem and yields a better signal clarity than Gabor transform. However, the S transform has its own disadvantages: the clarity is worse than Wigner distribution function and Cohen's class distribution function.[citation needed]

A fast S transform algorithm was invented in 2010.[3][4] It reduces the computational complexity from O[N2·log(N)] to O[N·log(N)] and makes the transform one-to-one, where the transform has the same number of points as the source signal or image, compared to storage complexity of N2 for the original formulation.[4][5] An implementation is available to the research community under an open source license.[6]

A general formulation of the S transform[4] makes clear the relationship to other time frequency transforms such as the Fourier, short time Fourier, and wavelet transforms.[4]

## Definition

There are several ways to represent the idea of the S transform. In here, S transform is derived as the phase correction of the continuous wavelet transform with window being the Gaussian function.

• S-Transform
${\displaystyle S_{x}(t,f)=\int _{-\infty }^{\infty }x(\tau )|f|e^{-\pi (t-\tau )^{2}f^{2}}e^{-j2\pi f\tau }\,d\tau }$
• Inverse S-Transform
${\displaystyle x(\tau )=\int _{-\infty }^{\infty }\left[\int _{-\infty }^{\infty }S_{x}(t,f)\,dt\right]\,e^{j2\pi f\tau }\,df}$

## Modified form

• Spectrum Form

The above definition implies that the s-transform function can be express as the convolution of ${\displaystyle (x(\tau )e^{-j2\pi f\tau })}$  and ${\displaystyle (|f|e^{-\pi t^{2}f^{2}})}$ .
Applying the Fourier transform to both ${\displaystyle (x(\tau )e^{-j2\pi f\tau })}$  and ${\displaystyle (|f|e^{-\pi t^{2}f^{2}})}$  gives

${\displaystyle S_{x}(t,f)=\int _{-\infty }^{\infty }X(f+\alpha )\,e^{-\pi \alpha ^{2}/f^{2}}\,e^{j2\pi \alpha t}\,d\alpha }$ .
• Discrete-time S-transform

From the spectrum rorm of S-transform, we can derive the discrete-time S-transform.
Let ${\displaystyle t=n\Delta _{T}\,\,f=m\Delta _{F}\,\,\alpha =p\Delta _{F}}$ , where ${\displaystyle \Delta _{T}}$  is the sampling interval and ${\displaystyle \Delta _{F}}$  is the sampling frequency.
The Discrete time S-transform can then be expressed as:

${\displaystyle S_{x}(n\Delta _{T}\,,m\Delta _{F})=\sum _{p=0}^{N-1}X[(p+m)\,\Delta _{F}]\,e^{-\pi {\frac {p^{2}}{m^{2}}}}\,e^{\frac {j2pn}{N}}}$

## Implementation of discrete-time S-transform

Below is the Pseudo code of the implementation.

  Step1.Compute ${\displaystyle X[p\Delta _{F}]\,}$
loop over m (voices){
Step2.Compute ${\displaystyle e^{-\pi {\frac {p^{2}}{m^{2}}}}}$ for ${\displaystyle f=m\Delta _{F}}$
Step3.Move ${\displaystyle X[p\Delta _{F}]}$  to ${\displaystyle X[(p+m)\Delta _{F}]}$
Step4.Multiply Step2 and Step3 ${\displaystyle B[m,p]=X[(p+m)\Delta _{F}]\cdot e^{-\pi {\frac {p^{2}}{m^{2}}}}}$
Step5.IDFT(${\displaystyle B[m,p]}$ ).
Repeat.}


## Comparison with other time–frequency analysis tools

### Comparison with Gabor transform

The only difference between the Gabor transform (GT) and the S transform is the window size. For GT, the windows size is a Gaussian function ${\displaystyle (e^{-\pi (t-\tau )^{2}})}$ , meanwhile, the window function for S-Transform is a function of f. With a window function proportional to frequency, S Transform performs well in frequency domain analysis when the input frequency is low. When the input frequency is high, S-Transform has a better clarity in the time domain. As table below.

 Low-frequency Bad clarity in time domain Good clarity in frequency domain High-frequency Bad clarity in frequency domain Good clarity in time domain

This kind of property makes S-Transform a powerful tool to analyze sound because human is sensitive to low frequency part in a sound signal.

### Comparison with Wigner transform

The main problem with the Wigner Transform is the cross term, which stems from the auto-correlation function in the Wigner Transform function. This cross term may cause noise and distortions in signal analyses. S-transform analyses avoid this issue.

### Comparison with the short-time Fourier transform

We can compare the S transform and short-time Fourier transform (STFT).[2][7] First, a high frequency signal, a low frequency signal, and a high frequency burst signal are used in the experiment to compare the performance. The S transform characteristic of frequency dependent resolution allows the detection of the high frequency burst. On the other hand, as the STFT consists of a constant window width, it leads to the result having poorer definition. In the second experiment, two more high frequency bursts are added to crossed chirps. In the result, all four frequencies were detected by the S transform. On the other hand, the two high frequencies bursts are not detected by STFT. The high frequencies bursts cross term caused STFT to have a single frequency at lower frequency.

## Applications

• Signal filterings[8]
• Magnetic resonance imaging (MRI)[9]
• Power system disturbance recognition
• S transform has been proven to be able to identify a few types of disturbances, like voltage sag, voltage swell, momentary interruption, and oscillatory transients.[10]
• S transform also be applied for other types of disturbances such as notches, harmonics with sag and swells etc.
• S transform generates contours which are suitable for simple visual inspection. However, wavelet transform requires specific tools like standard multiresolution analysis.
• Geophysical signal analysis

## References

1. ^ Stockwell, RG; Mansinha, L; Lowe, RP (1996). "Localization of the complex spectrum: the S transform". IEEE Transactions on Signal Processing. 44 (4): 998–1001. Bibcode:1996ITSP...44..998S. CiteSeerX 10.1.1.462.1500. doi:10.1109/78.492555.
2. ^ a b Stockwell, RG (1999). S-transform analysis of gravity wave activity from a small scale network of airglow imagers. PhD thesis, University of Western Ontario, London, Ontario, Canada.
3. ^ Brown, RA; Frayne, R (2008). A fast discrete S-transform for biomedical signal processing. Conf Proc IEEE Eng Med Biol Soc. Vol. 2008. pp. 2586–9. doi:10.1109/IEMBS.2008.4649729. ISBN 978-1-4244-1814-5. PMID 19163232. S2CID 29974786.
4. ^ a b c d Brown, Robert A.; Lauzon, M. Louis; Frayne, Richard (January 2010). "A General Description of Linear Time-Frequency Transforms and Formulation of a Fast, Invertible Transform That Samples the Continuous S-Transform Spectrum Nonredundantly". IEEE Transactions on Signal Processing. 58 (1): 281–290. Bibcode:2010ITSP...58..281B. doi:10.1109/tsp.2009.2028972. ISSN 1053-587X. S2CID 16074001.
5. ^ Kelly Sansom, "Fast S Transform", University of Calgary, https://www.ucalgary.ca/news/utoday/may31-2011/computing
6. ^ "Fast S-Transform download | SourceForge.net".
7. ^ E. Sejdić, I. Djurović, J. Jiang, "Time-frequency feature representation using energy concentration: An overview of recent advances," Digital Signal Processing, vol. 19, no. 1, pp. 153-183, January 2009.
8. ^ Ditommaso, R, Mucciarelli M, Ponzo FC (2012). ANALYSIS OF NONSTATIONARY STRUCTURAL SYSTEMS BY USING A BAND-VARIABLE FILTER. Bulletin of Earthquake Engineering. doi:10.1007/s10518-012-9338-y. See also MATLAB file
9. ^ Hongmei Zhu, and J. Ross Mitchell, "The S Transform in Medical Imaging," University of Calgary Seaman Family MR Research Centre Foothills Medical Centre, Canada.
10. ^ Prakash K. Ray, et al. "Coherency determination in grid-connected distributed generation based hybrid system under islanding scenarios." Power and Energy (PECon), 2010 IEEE International Conference on. IEEE, 2010,doi:10.1109/PECON.2010.5697562
• Rocco Ditommaso, Felice Carlo Ponzo, Gianluca Auletta (2015). Damage detection on framed structures: modal curvature evaluation using Stockwell Transform under seismic excitation. Earthquake Engineering and Engineering Vibration. June 2015, Volume 14, Issue 2, pp 265–274.
• Rocco Ditommaso, Marco Mucciarelli, Felice C. Ponzo (2010). S-Transform based filter applied to the analysis of non-linear dynamic behaviour of soil and buildings. 14th European Conference on Earthquake Engineering. Proceedings Volume. Ohrid, Republic of Macedonia. August 30 – September 3, 2010. (downloadable from http://roccoditommaso.xoom.it)
• M. Mucciarelli, M. Bianca, R. Ditommaso, M.R. Gallipoli, A. Masi, C Milkereit, S. Parolai, M. Picozzi, M. Vona (2011). FAR FIELD DAMAGE ON RC BUILDINGS: THE CASE STUDY OF NAVELLI DURING THE L’AQUILA (ITALY) SEISMIC SEQUENCE, 2009. Bulletin of Earthquake Engineering. doi:10.1007/s10518-010-9201-y.
• J. J. Ding, "Time-frequency analysis and wavelet transform course note," the Department of Electrical Engineering, National Taiwan University (NTU), Taipei, Taiwan, 2007.
• Jaya Bharata Reddy, Dusmanta Kumar Mohanta, and B. M. Karan, "Power system disturbance recognition using wavelet and s-transform techniques," Birla institute of Technology, Mesra, Ranchi-835215, 2004.
• B. Boashash, "Notes on the use of the wigner distribution for time frequency signal analysis", IEEE Trans. on Acoust. Speech. and Signal Processing, vol. 26, no. 9, 1987
• R. N. Bracewell, The Fourier Transform and Its Applications, McGrawHill Book Company, New York, 1978
• E. O. Brigham, The Fast Fourier Transform, Prentice-Hall Inc., Englewood Cliffs, New Jersey, 1974
• Cohen, L. (1989). "Time-frequency distributions—A review". Proc. IEEE. 77 (7): 941–981. CiteSeerX 10.1.1.1026.2853. doi:10.1109/5.30749.
• I. Daubechies, "The wavelet transform, time-frequency localization and signal analysis", IEEE Trans. on Information Theory, vol. 36, no. 5, Sept. 1990
• Farge, M. (1992). "Wavelet transforms and their application to turbulence". Annual Review of Fluid Mechanics. 24: 395–457. doi:10.1146/annurev.fluid.24.1.395.
• D. Gabor, "Theory of communication", J. Inst. Elect. Eng., vol. 93, no. 3, pp. 429–457, 1946
• Goupillaud, P.; Grossmann, A.; Morlet, J. (1984). "Cycle-octave and related transforms in seismic analysis". Geoexploration. 23: 85–102. doi:10.1016/0016-7142(84)90025-5.
• F. Hlawatsch and G. F. Boudreuax-Bartels, 1992 "Linear and quadratic timefrequency signal representations", IEEE SP Magazine, pp. 21–67
• Rioul, O.; Vetterli, M. (1991). "Wavelets and signal processing" (PDF). IEEE SP Magazine. 8 (4): 14–38. Bibcode:1991ISPM....8...14R. doi:10.1109/79.91217. S2CID 13266737.
• R. K. Young, Wavelet Theory and its Applications, Kluwer Academic Publishers, Dordrecht,1993