search

ZeyuJin / geodetic_inversion

homogeneous/layered inversion using InSAR/GPS
30 stars 15 forks source link

The weight of GPS #12

Open lightwjq opened 1 month ago

lightwjq commented 1 month ago

Dear Dr. Zou,

Could you please let me know if the program includes a module for searching the optimal weight values between GNSS and InSAR using the search method? If so, which function file is it? I would appreciate your guidance.

Best, Jingqi

x3zou commented 1 month ago

Dear Jingqi,

We don't currently have a module for searching for the optimal weight but it should be easy to do. To search for the optimal relative weight, you may want to keep the weight of one kind of data set to be 1 and vary the weight of another dataset. Then compute the RMS misfit of the model for the two datasets respectively and find a trade-off value. You can refer to the figure 7 of this paper to get an idea Wang2014.pdf

Thanks, Xiaoyu

lightwjq commented 1 month ago

Dear Dr. Zou,

Thank you very much for your reply! I have another question to ask you. Is it possible to restrict the rake angle during the inversion? If so, could you please let me know where I can modify the relevant settings?

Thanks, Jingqi

x3zou commented 1 month ago

Hi Jingqi,

In our codes, we don't specify a rank angle, but instead, we constrain the dip-slip and strike-slip, as rake is the arctan of dip-slip/strike-slip. You can go to bounds_new.m to see the constraints. The 'Con' parameter controls the sign of strike-slip and dip-slip. Also, you can adjust the upper/lower bounds of strike-slip and dip-slip inside bounds_new.m

Xiaoyu

lightwjq commented 1 month ago

Dear Dr. Zou,

Thank you very much for your guidance! I will modify the corresponding file according to your prompts. In addition, I want to use the hierarchical earth model to calculate the Green's function. We find that there are related calculation functions in the Layered_ EDCMP folder. In the founction of ''build_green_from_edcmp'', we have not found the case of the input file, can you provide relevant cases for our reference?

Thanks, Jingqi

x3zou commented 1 month ago

Hello Jingqi,

We are using EDGRN to calculate the layered Green's function. To do that, 1) please click here to download and install EDCMP-EDGRN 2) After that please run the EDGRN and make sure you have created the correct directories for EDGRN outputs (usually name it as 'edgrnfcts') 3) Then run the getedgrn() MATLAB function to generate EDGRN input for MATLAB. 4) Load your fault geometry or okada model and the input data created in your previous steps and run the calc_green_LOS_edcmp(). Please see the calc_green_layered.m in the edcmp_edgrn_util directory. I have also uploaded a sample input file ('edgrnPamir') for EDGRN within that directory but you should change the layered velocity structure based on your region of interest.

I know it might look a bit complicated so please let me know if you encounter any problems.

Xiaoyu

lightwjq commented 1 month ago

Hi Xiaoyu,

Thank you for your patient guidance! I have completed the first three steps according to your instructions. However, I find that the 'izmhs.cl' file is different from the other files obtained by EDGRN outputs. And I modify the getedgrn() function to ensure the code running smoothly. But I do not know how to make fault geometry by using the files obtained from getedgrn() function. I have found the functions in the edcmp_edgrn_util, Layered_EDCMP and other folders but no solution has been found. I hope to get your help. Thanks for your time!

Jingqi

x3zou commented 1 month ago

Hi Jingqi,

The fault geometry is made in your previous steps. Have you tried running the Okada model already? If so you can just save your Okada slip model and load it into the calc_green_layered.m, and it will generate layered_green.mat in your data directries. Then you can switch the inversion type to 'layered' and it can make layered inversion. Let me know if that works. Thanks

Xiaoyu