geopsy.org

Posted: **Fri Mar 18, 2022 10:13 am**

Hi Marc,

I'm trying to use the geopsy-hv to estimate H/V hourly,
I have some problems when I use the -load-windows option.

The terminal shows:
The ASCII format without header imply that you fill in the correct information yourself in a table. Choose menu "view/new table" and change the SamplingPeriod (or SamplingFrequency) to allow correct plotting in a graphic
geopsy-hv: Loading (My data path)
geopsy-hv: Selected groups:
Option '-load-windows' requires an argument, see -help for details.

I also tried to edit the FROM_TIME_TEXT and TO_TIME_TEXT in the param file to achieve the purpose, but it didn't work.

Could you provide an example of the command line?

Sincerely,
Uki

Posted: **Tue Mar 22, 2022 8:46 am**

Hi Uki,

Here is an example of a bash script to compute H/V on a hourly basis.
Note that geopsy-hv accepts directly signal files from version 3.4.0, not before.
In this code, the files names are supposed to end with YEAR.DOY (Day Of Year).

Code: Select all

(
cat <<END
FROM_TIME_TYPE=Absolute
TO_TIME_TYPE=Delta
TO_TIME_TEXT=1h
MINIMUM_FREQUENCY=1
MAXIMUM_FREQUENCY=20
WINDOW_LENGTH_TYPE=Exactly
WINDOW_MIN_LENGTH(s)=25
WINDOW_MAX_LENGTH(s)=25
WINDOW_MAX_COUNT=0
END
) > hv.param

FILE="STAT01.HH?.2022.034"
YEAR=${FILE:0-8:4}
DEBUT=$(echo ${FILE:0-3} | awk '{print ($1-1)*86400}' | gptime -i "seconds since ${YEAR}0101")
DEBUT_T=$(echo $DEBUT | gptime -i "yyyy-MM-dd hh:mm:ssz" -o time_t)

for START_T in $(echo $DEBUT_T | awk '{for(t=0; t<86400; t+=3600) print $1+t}'); do
  START_TIME=$(echo $START_T | gptime -i time_t -o "yyyyMMddhhmmss")
  echo $START_TIME
  geopsy-hv ${FILE} -param hv.param -set FROM_TIME_TEXT=$START_TIME -o RESULTS_${STATION}
done

Posted: **Tue Mar 29, 2022 10:46 am**

Hi Marc,

Thank you for the information.
I'm able to run it now.

Best,
Uki

geopsy.org

Questions about geopsy-hv usage

Questions about geopsy-hv usage

Re: Questions about geopsy-hv usage

Re: Questions about geopsy-hv usage