geopsy.org

Hi Marc,

I'm creating an ipython tutorial how to use Geopsy from command and how to loop over data to demonstrate the HV variability with time. My process strategy is:
1) first to create my own PARAMETER string and fill in whatever value I want in that string.
2) then to activate the command line using this param string.

However, whatever value I input in the Frequency Rejection params, the output .HV is always the same and the maximum of windows is selected. Using the freq. rej. params in the command line does not result in a windows rejection. Because of this problem, the error in the .hv file is too high compared to manual Geopsy treatment.

Have you tested this or have you seen the same problem? I tested it by doing the exact same process in Geopsy manually and from the command line, but that gives different results. I'm using Geopsy pack win64-3.4.2.

I attach an example to clarify this problem.
Thanks for your help!
Koen Van Noten
Royal Observatory of Belgium

Hi Koen,

At first glance (without digging into the code and testing it), I would say that frequency rejection is not performed at all when running from command line. In the GUI, usual H/V must be computed first. The rejection takes place after a second start of the process. I guess that this second start is not triggered when running from command line. There is probably something to improve in this side. The simplest way to implement it is to add a mode like "-hv-rejection" (on the same level as "-hv" and "-spectrum").

Thanks for your reply Marc!
I'm looking forward to this kind of solution. In the meantime I'll play a bit with the STA/LTA params, but I found that the frequency rejection was performing better in windows rejection.
Cheers
Koen