12 | | The spectrogram is computed as follows: a subwindow of N samples of the selected trace window is transformed to frequency domain using the FFT algorithm (see Numerical Recipes, Press et al., Cambridge University Press, 1988). N is determined by the configuration parameter spectrogram_width and must be a power of 2. Reasonable values are e.g. 512 or 1024. Before FFT the input trace is tapered in time domain using a cosine taper (2*cos(x)-1, -pi < x < pi). The squared frequency amplitudes are displayed in a colour-coded vertical bar/line in the output window. The next subwindow for FFT is found by shifting the last one by K samples to the right. K is determined by the configuration parameter spectrogram_step. The first subwindow starts with the first sample of the selected trace part, the last subwindow ends between 0 and K-1 samples before the end of the selected trace part. The output matrix contains frequency over time. The header info of the output file contains DELTA (sample distance in time, s), START (absolute start time of trace), LENGTH (number of frequency samples on each time step), LINES (number of time steps), DF (frequency sample distance in Hz). Following are LENGTH*LINES numbers giving frequency amplitudes over the time steps. The first LENGTH samples refer to time sample 1, the second LENGTH samples refer to time sample 2 and so on. |
| 14 | SHX: |
| 15 | Creates png figure in specifed output file. Opened if /SHOW qualifier |
| 16 | specified or no output file is given. |
| 17 | |
| 18 | SH: |
| 19 | The spectrogram is computed as follows: a subwindow of N samples of the |
| 20 | selected trace window is transformed to frequency domain using the FFT |
| 21 | algorithm (see Numerical Recipes, Press et al., Cambridge University Press, |
| 22 | 1988). N is determined by the configuration parameter spectrogram_width |
| 23 | and must be a power of 2. Reasonable values are e.g. 512 or 1024. Before |
| 24 | FFT the input trace is tapered in time domain using a cosine taper |
| 25 | (2*cos(x)-1, -pi < x < pi). The squared frequency amplitudes are displayed |
| 26 | in a colour-coded vertical bar/line in the output window. The next |
| 27 | subwindow for FFT is found by shifting the last one by K samples to the |
| 28 | right. K is determined by the configuration parameter spectrogram_step. |
| 29 | The first subwindow starts with the first sample of the selected trace |
| 30 | part, the last subwindow ends between 0 and K-1 samples before the end of |
| 31 | the selected trace part. The output matrix contains frequency over time. |
| 32 | The header info of the output file contains DELTA (sample distance in time, |
| 33 | s), START (absolute start time of trace), LENGTH (number of frequency |
| 34 | samples on each time step), LINES (number of time steps), DF (frequency |
| 35 | sample distance in Hz). Following are LENGTH*LINES numbers giving frequency |
| 36 | amplitudes over the time steps. The first LENGTH samples refer to time |
| 37 | sample 1, the second LENGTH samples refer to time sample 2 and so on. |
| 38 | |
| 39 | SHX: |
| 40 | Use obspy trace spectrogram method. |
| 41 | |
| 42 | === qualifiers === |
| 43 | |
| 44 | * /SHOW (SHX only) [[br]] |
| 45 | Open and show specified output file. |