geopsy.org

Posted: **Sat Nov 28, 2009 10:57 am**

Good morning Marc,

I haven't find an automatic way to have the value of the peak value (amplification) associated with the frequency. The file that can be saved shows only, if I am not wrong, the frequency values.

Is there a way to also have an automatic extraction ? For the moment, I am just using the tool provided to look on the display, but that is not very precise.

Thanks in advance
Corine

Posted: **Thu Dec 03, 2009 9:29 am**

Hi Corine,

Correct me if I misunderstood your question. Do you want to extract the H/V amplitude at the peak frequency? In the saved file, you can extract the amplitude at the peak frequency.

Code: Select all

# GEOPSY output version 1.1                          
# Number of windows = 8                       
# f0 from average       2.07265
# Number of windows for f0 = 8
# f0 from windows       2.05321 1.9398  2.16663
# Frequency     Average Min       Max
0.01            1.82063 0.480536  6.89794
0.0108984       1.82464 0.47767   6.96989
0.0118775       1.82822 0.474326  7.04659
0.0129446       1.83116 0.470404  7.12825 
...

Run this script (edit file name or modify script interface):

Code: Select all

#!/bin/bash

FILE=S10A.hv
FREQUENCY=$(grep "f0 from windows" $FILE | awk '{print $5}')
gpcurve -value $FREQUENCY $FILE

Which gives as output:

Code: Select all

# GEOPSY output version 1.1
# Number of windows = 8
# f0 from average       2.07265
# Number of windows for f0 = 8
# f0 from windows       2.05321 1.9398  2.16663
# Frequency     Average Min     Max
2.0532099999999999795 9.5992748074446918594 8.1223796780990280553 11.348137726793865454

Keep in mind that H/V amplitude is not a straight forward proxy to real amplification...

geopsy.org

A values with geopsy

A values with geopsy