rotate spectrum

[gdg]

Hi Marc
my aim is to get (and save as ascci file) rotated spectra for regular bins for the components of motion. I mean not H/V but the spectra of the single component.

I have a vague remind that it was possible with the very old release of geopsy.

Now I use the tool Spectrum Rotate, but it seems to return a single horizontal spectrum. What is exactly this result? I suppose a kind of mean between rotated NS and EW components.
How can I save separately the rotated spectra of the single component (I mean N+10, N+20 etc) ?

Thanks so much
Giuseppe

[admin]

Hi Giuseppe,

I've just tried... Simply run

Code: Select all

geopsy-hv -db ~/sandbox/bestrtbf/Mirandola.gpy -group CN01-3C -rotate

It generates file 'CN01.hv' which is a grid file that can be viewed with

Code: Select all

cat CN01.hv | figue -g

The text file contains 3 columns: frequency, azimuth, h/v value

For spectra:

Code: Select all

geopsy-hv -db ~/sandbox/bestrtbf/Mirandola.gpy -group CN01-3C -rotate -spectrum

And to plot:

Code: Select all

cat CN01_Z.spec | figue -g

The naming could be improved. I do not know why there is such 'Z'

[gdg]

Thanks Marc, the script works; tested with the 3.4.1 release.
For the Spectra, in the output *Z.spec file for each angle I have a spectrum. But I would like to individually save the rotated EW and NS components. Or in other words, the script on rotate spectra works using 3C data as input, but I cannot find the way to recover separately the rotated NS and EW component.
My best. Giuseppe

[admin]

Hi Giuseppe,

I'm not sure I completely got your question.
In the result file, there is one spectrum per direction from 0 to 180 deg. (excluded), with a default step of 10 deg. (adjustable in the parameters). If you want the rotated NS and EW, let's say rotated by 20 deg., you can run

Code: Select all

awk '{if($2==20) print $1, $3}' CN01_Z.spec > NS_20.spec
awk '{if($2==100) print $1, $3}' CN01_Z.spec > EW_20.spec

If you want the rotated signals, you can use rotate from Waveform menu. It is also accessible from the command line via the waveform console.

Read you

Marc

[gdg]

Perfect Marc. I had misinterpreted the output *spec file, now clear thanks you answer.
solved !!
Giuseppe

[luigiV]

Hi Marc, I submitted you (07/16) an approved question that seems to me similar to those posed by Giuseppe and to which you have not answered. The problem is the following:
"Hi Marc,
I want to ask you if it is possible to save the result of an H / V curve obtained with the directional energy that can be used for the calculation of the power spectrum and which represents the curves of the single components that lead to that spectral ratio . In my attempts, for example saving the H / V result of a 30 ° rotation of the horizontal components (with directional energy option) and loading the relative log file for the calculation of the power spectrum, I realized that it is identical to the one of the spectral ratio obtained with squared energy. I also tried to rotate the components with the rotate components option in the Waveform section and then calculate H / V with the squared energy, but the result is different from the one using the directional energy option.
Is there something wrong with my request?
I am attaching a plot with 4 well-differentiated spectral ratio curves, which correspond to four log files that provide the same power spectrum.
Thanks
Luigi "
With Windows can I simply think about rotating the components of the desired angle (in waveform/rotate components) and then directly calculate the spectrum and save the two components?
Why using Waveform / rotate components and calculating H / V I get a different result from the one with the directional energy option of the same angle?
Thank you
Luigi