Comparison_strain result no contains visr figure

[ditafaith]

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?

[kmaterna]

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:

Strain_2D/Strain_Tools/strain/pygmt_plots.py

Line 259 in f87a4ba

def plot_method_differences(strain_values_ds, average_strains, region, outdir, outfile):

and the line in question is here:

Strain_2D/Strain_Tools/strain/pygmt_plots.py

Line 280 in f87a4ba

fig.colorbar(position="JCR+w4.0i+v+o0.7i/-0.5i", cmap=outdir+"/mycpt.cpt", truncate="-300/300",

(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]

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]

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 <#35 (comment)>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/AC7QOVAGK3OVOB76X6BCOKDYJKIWZAVCNFSM6AAAAABASV2UZOVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMYTQNJVGE4TENRZGQ> . You are receiving this because you were mentioned.Message ID: ***@***.***>