search

paudetseis / SplitPy

Teleseismic shear-wave splitting analysis
https://paudetseis.github.io/SplitPy/
MIT License
20 stars 9 forks source link

pickle object to text file #10

Open bosby74 opened 3 years ago

bosby74 commented 3 years ago

Hi, I'm trying the use of your SplitPy software and I'm finding some very interesting characteristics as the use of open platform (python) or the the automatization of the process that take into account the downloading of the data. I tested your example of TGTN station without error so everything seem to be installed properly. Anyway, since I'm not a very proficient user of python, I'm finding some difficulties to "read humanly" the results saved because the pickle object resulted at the end of the analysis is not easily manageable (at least for me). Are there any scripts to convert the pickle object of the results into a csv ot txt file? If not, do you plan to write down some in the next future? Or to provide example in your website? Thanks in advance, Simone

paudetseis commented 3 years ago

Hi and thanks for your feedback. It's nice to hear that you find this code useful. Our original intent was to save disk space and essentially save a single file per object - admittedly this part is a bit clunky and could be improved. It is certainly possible to add flexibility to save the individual results into a .csv file and we will look into this in a future release. I will note that it is possible to get the estimates in a .csv table once you run the script split_average.py.

In the mean time, you can read the individual pickle files by navigating to the folder of interest (e.g., in the DATA/NY.TGTN/ folder), open a python window and type:

>>> import pickle
>>> file = open('Split_results_manual.pkl', 'rb')
>>> split_SC = pickle.load(file)
>>> split_RC = pickle.load(file)
>>> split_null = pickle.load(file)
>>> split_quality = pickle.load(file)

and then you check the content of each split variable and extract the value that you want, e.g.:

>>> split_SC.__dict__  # This shows the object split_SC as a dictionary with keys and values
>>> dt_SC = split_SC.dtt
>>> dt_SC_err = split_SC.edtt
>>> phi_SC = split_SC.phi
>>> phi_SC_err = split_SC.ephi
...

...and so on. You can also get the value of the Null (True or False) as well as the quality of the estimate:

>>> print(split_null)
>>> print(split_quality)
bosby74 commented 3 years ago

Thank you, these commands will help me on read the pickles file and to write out the result on a text file. Ciao, Simone