kmaterna / Strain_2D

Strain rate modeling from GNSS velocity fields
MIT License
44 stars 19 forks source link

Comparison_strain result no contains visr figure #35

Open ditafaith opened 9 months ago

ditafaith commented 9 months ago

Hello @kmaterna Thanks for the codes, and helping me for the beginner level. If possible, it would be glad if there is sample for InSAR as well.

Furthermore, following your guidance, I got strain from all methods, including visr result, however when I ran to compare various, only 4 methods available. Also, how I can minimize the distance of colorbar with the figure?

strain_2d

kmaterna commented 9 months ago

Hi @ditafaith For the strain rate comparison, I would suggest to look into the config file and see the parameter input_dirs. This is a list of directories for comparison, separated by colons. Have you made sure that points to the visr directory? I suspect it's still showing the example case, and that the following change needs to be made:

input_dirs = Output/delaunay:Output/gpsgridder:Output/loc_avg_grad:Output/geostats to input_dirs = Output/visr:Output/gpsgridder:Output/loc_avg_grad:Output/geostats

As for the exact look of the visuals, our intention for this library was always to provide a starting point from which users can create their own publication-ready plots in other ways. I don't really want to start modifying or maintaining these plots to perfection. For templates or inspiration, the code starting point is here: https://github.com/kmaterna/Strain_2D/blob/f87a4baa4c732a61ac805ac9b63da8b934def5e5/Strain_Tools/strain/pygmt_plots.py#L259 and the line in question is here: https://github.com/kmaterna/Strain_2D/blob/f87a4baa4c732a61ac805ac9b63da8b934def5e5/Strain_Tools/strain/pygmt_plots.py#L280

(You can also make the desired source code change, e.g. +o0.7i/-0.5i to +o0.2i/-0.5i, on your own machine, and then python setup.py install to implement them on the local system. That's probably the easiest way for you to get what you want. If it works and you like it better, I'm happy to take your PR's contributed back to the repo.)

ditafaith commented 9 months ago

Thanks @kmaterna for the response and may be I am not clear in the previous question. For sure, I made the input_dirs = Output/delaunay:Output/gpsgridder:Output/loc_avg_grad:Output/geostats:Output/visr I assumed it could produce 5 figures since there are 5 input dirs. However, only 4 figures (without visr). This is my previous question, why didn't get 5 figures which is same number of inputdir.

kmaterna commented 9 months ago

I think the plot is hard-coded to plot only the first 4 on the list. My recommendation would be to discard the Delaunay one because it’s the least sophisticated method and produces the most obvious artifacts.

The velmap method may not be completely integrated into the software yet, but we hope to have it in the future.

On Wed, Dec 13, 2023 at 10:00 PM ditafaith @.***> wrote:

Thanks @kmaterna https://github.com/kmaterna for the response and may be I am not clear in the previous question. For sure, I made the input_dirs = Output/delaunay:Output/gpsgridder:Output/loc_avg_grad:Output/geostats:Output/visr I assumed it could produce 5 figures since there are 5 input dirs. However, only 4 figures (without visr). This is my previous question.

Can you explain why no results come for velmap method ?

— Reply to this email directly, view it on GitHub https://github.com/kmaterna/Strain_2D/issues/35#issuecomment-1855192694, or unsubscribe https://github.com/notifications/unsubscribe-auth/AC7QOVAGK3OVOB76X6BCOKDYJKIWZAVCNFSM6AAAAABASV2UZOVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMYTQNJVGE4TENRZGQ . You are receiving this because you were mentioned.Message ID: @.***>

ditafaith commented 3 months ago

Hello @kmaterna @jlmaurer After successfully run the example data, however, when I used my own, the visr result didn't draw entire map. Here, I'd like to hear from your side, what parameters do I need to change for the 3 methods below? I don't see any compression-extension, where is the error?

Selection_010