problems with gphistogram and RTBF geopsy 3.4.1.

[luigiV]

Hi Marc,
today I have processed a minishark dataset (7 stations 3C) with FK tool (RTBF). The .max file was obtained correctly, but using gphistogram I cannot use any print options (-p). If I don't type anything after gphistogram I can open the file.max and Rayleigh's slowness-frequency plot appears. If I type -p R (or L, or -ell etc) this message appears:
"bin $ gphistogram -p R
Reading file 'C: /Lavori/2021/IND00/fk.max' ...
Slowness max file with pattern 'R'
Loading: ----
[1ALoading: 0
No samples available
bin $ "
If I run gpviewmax I can see all plots perfectly (R, L, ell, azimuth, noise, power).
Also if I run RTBF radial (which should be the rerun for Love, is that correct? Because other options are not there to get L dispersion, as you have recently suggested) the program crashes.
I have also been able to verify that gphistogram of 3.3.6 works perfectly with RTBF max files of 3.3.6., but it does not read files.max of 3.4.1.
RTBF Radial of 3.3.6 runs correctly.

Regards
Luigi

ps: I would like to know why there is such an obvious difference, not only cosmetic, for the "passive RTBF" outputs generated by 3.3.6 and 3.4.1 on the same data set and with exactly the same parameters

[admin]

Hi Luigi,

In the original implementation of RTBF, Love wave were computed at the same time. The .max files were containing both types of waves with a column dedicated to a flag: "Vertical", "Rayleigh" or "Love". This discussion may bring interesting explanations.

"RTBF Radial" is still a RTBF process that produces Rayleigh dispersion curves. The "Radial" option acts on the type of power. "RTBF" uses power P^{prime}_{R_s} as defined in Wathelet et al. (2018) which is theoretically insensitive to incoherent noise. "RTBF Radial" uses power P^{prime}_{R_h}. The latter is biased by the incoherent noise but it is more stable for ellipticities approaching zero. It is good to run both options.

I have also been able to verify that gphistogram of 3.3.6 works perfectly with RTBF max files of 3.3.6., but it does not read files.max of 3.4.1.

On branch 3.3, there was the additional column for the type of wave. The contrary should work. If not there are scripts to translate them.

ps: I would like to know why there is such an obvious difference, not only cosmetic, for the "passive RTBF" outputs generated by 3.3.6 and 3.4.1 on the same data set and with exactly the same parameters

The output of 3.3.6 looks effectively strange. I worked on RTBF implementation between 3.3 and 3.4 to fix various issues (details can be found be checking the differences in git repository, but I can't remember). I would better trust the results of 3.4.1. Currently working on 3.5, I found again minor errors in 3.4.1. They are fixed in 3.4.2 to be publish soon.

[luigiV]

In your answer to Giuseppe you confirmed the need to re-run FK to get L dispersion with the script:
"CODE: SELECT ALL
geopsy-fk [database / group / signal specification] -set PROCESS_TYPE = CaponTransverse"

Using the Windows version I cannot use this process_type and therefore I have to conclude that this is currently not available from console (passive RTBF, radial RTBF, passive conventional are the only options).
Will we have to wait for the next releases?

regards
Luigi

[admin]

Release 3.4.2 will be ready for Window within a few minutes. It contains CaponTransverse inside the graphical interface. The command line should work also for Windows (e.g. under GeopsyLand). Why do you wrtie that you cannot under Windows? Is there any error message?

[luigiV]

It is surely due to an error in the my script. The data to be processed is in a gpy file that I copied, to simplify, to the geopsy bin directory.
geopsy-fk file.gpy -set PROCESS_TYPE = CaponTransverse
geopsy-fk: no database specified, see-help
What's wrong?
Luigi

[admin]

Code: Select all

geopsy-fk -db file.gpy -set PROCESS_TYPE=CaponTransverse

• Remove the spaces around the "="
• Insert option "-db" before the database file name

I checked under 3.4.2 that it works even without "-db". Are you already under 3.4.2?

[luigiV]

Sorry, one time more because I was working inadvertently on GeopsyLand 3.3.6. and so, I redid the script on GeopsyLand 3.4.1.
bin$ geopsy-fk file.gpy -set PROCESS_TYPE=CaponTransverse

The process started and runned till the end but did not save the .max file:
geopsy-fk: Process ended at 2021-07-30 18:55:34
geopsy-fk: Process run in 00:00:32
geopsy-fk: error writing to file 'capontransverse-importedsignals.max'
geopsy-fk: Error saving results to 'capontransverse-importedsignals.max'.

i have tried on 3.4.2, but the result is the same:
geopsy-fk: Process ended at 2021-08-01 10:15:45
geopsy-fk: Process run in 00:00:37
geopsy-fk: error writing to file 'capontransverse-importedsignals.max
geopsy-fk: Error saving results to 'capontransverse-importedsignals.max'.


gphistogram of 3.4.2 has the same problems as gphistogram 3.4.1.
You can only see the slowness-frequency graph. The others with gpviewmax.
bin $ gphistogram -p R
Reading file 'C: /Works/2021/MYDIR/fk_RTBF.max' ...
Slowness max file with pattern 'R'
Loading: -
[1ALoading: 0
No samples available

bin $ gphistogram -p R -ell
Reading file 'C: /Works/2021/MYDIR/fk_RTBF.max' ...
Slowness max file with pattern 'R
Loading: -
[1ALoading: 0
Loaded 0 samples
Filtering Y <-10
remaining 0 samples (nan%)
Filtering Y> 10
remaining 0 samples (nan%)
No samples available

Last question about SPAC:
if I set in Spac toolbox the option "auto" for the rings, the target output file is not readible by Spac2disp because:
bin$ spac2disp
Unknown tag at line 23 for context TargetList: EllipticityPeakTarget
bin$ spac2disp
Unknown tag at line 2351 for context TargetList: EllipticityPeakTarget

If I run spac2disp from the bin directory an error message in the white box is showed:
No frequency samples found in target

How can I use the option "auto"?

Regards
Luigi

[admin]

i have tried on 3.4.2, but the result is the same:
geopsy-fk: Process ended at 2021-08-01 10:15:45
geopsy-fk: Process run in 00:00:37
geopsy-fk: error writing to file 'capontransverse-importedsignals.max
geopsy-fk: Error saving results to 'capontransverse-importedsignals.max'.

You have to check if you have enough permissions to write in the current directory. I know that permissions under Windows are a bit obscure. Try to add option "-o /path/where/you/can/write". "-o" specifies the base name (and optionally the path) of the output file. Did you try to run "CaponTransverse" in the graphical interface?

A simple way to check if it is a permission related issue, run the following commands:

Code: Select all

$ echo "Ciao mondo!" > toto
$ cat toto
Ciao mondo!

gphistogram of 3.4.2 has the same problems as gphistogram 3.4.1.
You can only see the slowness-frequency graph. The others with gpviewmax.

Do not use the "-p" option with the new .max files. The pattern "R" is not found in the new .max files and as expected there are no samples to load.

Last question about SPAC:
if I set in Spac toolbox the option "auto" for the rings, the target output file is not readible by Spac2disp because:
bin$ spac2disp
Unknown tag at line 23 for context TargetList: EllipticityPeakTarget
bin$ spac2disp
Unknown tag at line 2351 for context TargetList: EllipticityPeakTarget

The "Auto" button generates automatically the ring list but if you create it yourself manually, it does not change the output file. I'm puzzled about the error sent by spac2disp. In any case there is always a section "EllipticityPeakTarget" in all target files, which is empty for SPAC target files. I checked the attached file under Windows (geopsypack-3.4.2) with no error. Can you try to load it in your spa2disp? Is it spac2disp from geopsypack-3.4.2 (check it with option "-version"), I guess so, but let's check everything? Can you post your unreadable .target file?

[luigiV]

The "Auto" button generates automatically the ring list but if you create it yourself manually, it does not change the output file. I'm puzzled about the error sent by spac2disp. In any case there is always a section "EllipticityPeakTarget" in all target files, which is empty for SPAC target files. I checked the attached file under Windows (geopsypack-3.4.2) with no error. Can you try to load it in your spa2disp? Is it spac2disp from geopsypack-3.4.2 (check it with option "-version"), I guess so, but let's check everything? Can you post your unreadable .target file?

Spac2disp reads your file correctly (I had already successfully replicated the Spac on the passive linear array file).
I am attaching the file generated with the "auto" option on a triangular array with 7 stations. Please note that in spac toolbox / rings no ring appears when I activate the "auto" option.

You have to check if you have enough permissions to write in the current directory. I know that permissions under Windows are a bit obscure. Try to add option "-o /path/where/you/can/write". "-o" specifies the base name (and optionally the path) of the output file. Did you try to run "CaponTransverse" in the graphical interface?

A simple way to check if it is a permission related issue, run the following commands:
CODE: SELECT ALL

$ echo "Ciao mondo!" > toto
$ cat toto
Ciao mondo!

With minishell the script doesn't work. This is the answer:
bin $ echo "Hello world"> toto
C: / Program Files (x86) /geopsy3.4.2/bin/toto: Cannot write to file
bin $ cat toto
cat: toto: No such a file
bin $
The -o command is not recognized by minishell and I don't know how to use a Linux terminal.
Running the script: geopsy-fk -set PROCESS_TYPE = CaponTranserve no CaponTransverse.max file is saved
How I have to do to run the CaponTransverse graphical interface?

Do not use the "-p" option with the new .max files. The pattern "R" is not found in the new .max files and as expected there are no samples to load.

So I can only view with gpviewmax the other RTBF outputs if no print options are accepted...

Regards
Luigi

[admin]

I am attaching the file generated with the "auto" option on a triangular array with 7 stations. Please note that in spac toolbox / rings no ring appears when I activate the "auto" option.

Effectively if the ring list is empty, the .target file will be also empty. The automatic definition of rings is based on some reasonable but arbitrary conditions that must be all satisfied:

• There must be at least 2 pairs of stations in the whole dataset.
• The angle step between pairs in a ring must be less than 40 deg. (only for non linear arrays).
• There must be at least 6 pairs in a ring.
• The relative thickness of a ring must be less than 20% of its average radius (average of its minimum and maximum radii, not the average of the pair radii).

The limits of rings are rounded to a cm precision.

Thus for a triangular array, there is a maximum of three pairs. The automatic process fails to define a ring. The manual definition is recommended.

bin $ echo "Hello world"> toto
C: / Program Files (x86) /geopsy3.4.2/bin/toto: Cannot write to file

You do not have the permission to write into such location. That's the reason why geopsy-fk did not save the output file. Use

Code: Select all

cd /your/working/directory

to change the current directory that can be checked with "pwd" (print working directory). Alternatively add "-o" option to geopsy-fk.

Code: Select all

geopsy-fk -db file.gpy -set PROCESS_TYPE=CaponTransverse -o /your/working/directory/toto

How I have to do to run the CaponTransverse graphical interface?

Go into the FK toolbox, with 3.4.2, the list of available process types has been updated. If you provide horizontal or three-component signals, "CaponTransverse" is available. Not with only vertical components.

So I can only view with gpviewmax the other RTBF outputs if no print options are accepted...

Option "-p" stands for "pattern", it rejects all lines that do not contain the pattern provided as argument after the option "-p". Just run gphistogram without the option "-p" and it will show the results as expected.

Best regards,

Marc

[luigiV]

Thus for a triangular array, there is a maximum of three pairs. The automatic process fails to define a ring. The manual definition is recommended.

Ok, but for a triangular array with 6 stations+1 central I think there are more than 3 pairs. Right?

How I have to do to run the CaponTransverse graphical interface?
Go into the FK toolbox, with 3.4.2, the list of available process types has been updated. If you provide horizontal or three-component signals, "CaponTransverse" is available. Not with only vertical components

Sorry, but as you can see, in the release 3.4.2 downloaded, CaponTransverse is not in the process list.


Thanks a lot for the other answers

Best regards
Luigi