esrivas17 / gacos_correction_for_GMTSAR

Atmospheric GACOS correction adapted for GMTSAR
17 stars 6 forks source link

How do I use GACOS_correction.csh in sbas processing? #15

Closed MinervaLee0429 closed 4 months ago

MinervaLee0429 commented 1 year ago

Dear professor, I have a question to ask, because I am doing SBAS-INSAR processing, I want to use the script GACOS_correction.csh for correction, but because I am doing it under the processing flow of sbas, and I also see some opinions in the forum that this script can only be used after unwrap. At the same time, SBAS is unwrapped after merging multiple subwaths. I would like to know whether GACOS_correction.csh can correct the interference pattern in the merge folder after merging and unwrapping. If yes, do I need to change the merge folder? What changes, if any, should be made? I would really appreciate it if you could answer me.

esrivas17 commented 1 year ago

Hi, yes you must apply the correction to unwrapped interferograms. if you have merged multiple subswaths, you should point to the folder where your unwrap.grd files are. I think that originally the code was not meant to be applied to merged subswaths (I cant remember what was the problem) but I think it should work.

MinervaLee0429 commented 1 year ago

Thank you very much for your reply. I will try it.

MinervaLee0429 commented 1 year ago

Dear Professor, I tried to run the GACOS_correction.csh script in the merge folder but the following error occurred. 截图 2023-11-24 18-01-28

this is my merge folder and one of the merged interferogram folders. 截图 2023-11-24 18-21-00 截图 2023-11-24 18-40-57

There are also the topo folders created by following the instructions, 截图 2023-11-24 18-40-57

and the folders where the GACOS data resides. 截图 2023-11-24 18-26-05 Please guide me on how to solve this error. I have not been able to find the example of applying GACOS_correction.csh script in sbas processing, so please help me solve this problem and make the script work.

esrivas17 commented 1 year ago

Hi, could you show the content of your intflist file?

esrivas17 commented 1 year ago

ALso, you are pointing wrongly to the GACOS directory

MinervaLee0429 commented 1 year ago

@esrivas17 Hello professor, thank you for your reply, this is the content of my intflist, 截图 2023-11-24 21-19-39 the path of GACOS is changed after I failed to run the script, the path is actually no problem when the screenshot. Therefore, I really do not know where the problem is, and I would like to ask you to answer this question for me. 截图 2023-11-24 21-22-42

esrivas17 commented 1 year ago

Hi, checking your log files seems like you need to add some files to your "topo" directory... make sure that the file master.PRM is there. if not, try to get it from your non-merged directory. Also, copy the "date.SLC" files from your original directory to your merge directory. These files are used to get the dates of each SLC

MinervaLee0429 commented 1 year ago

@esrivas17 Thank you for your reply, Professor. Regarding the topo folder, I saw the script prompt that requires the master.PRM file and the corresponding. LED file and dem.grd file, I merge the supermaster.PRM generated during the process as master.PRM, and find the corresponding. LED file, placed in the topo folder.

In addition, you said "opy the "date.SLC" files from your original directory to your merge directory" is the corresponding date in each subtape. Are SLCS copied to the corresponding date interferogram folder in the merge directory? SLC, S1_20180110_100451_F1.SLC, S1_20180110_100451_F2.SLC, S1_20180110_100451_F3.SLC, S1_20180215_100450_F1.SLC, S1_20180215_1 Are 00450_F2.SLC and S1_20180215_100450_F3.SLC moved to 2018009_2018045 in the merge directory?

MinervaLee0429 commented 1 year ago

Professor, this is my topo folder. 截图 2023-11-24 21-51-29

esrivas17 commented 1 year ago

yes, please rename your supermaster.PRM to master.PRM or you could edit the code to call the same file. And yes, but please you would need to paste the SLC files in its corresponding directory in the merge folder

MinervaLee0429 commented 1 year ago

@esrivas17 Hello Professor, thank you for your tireless guidance, I followed your guidance, the script did start to run, but did not produce the final result. At runtime I get this error, 截图 2023-11-25 00-11-07 截图 2023-11-25 00-30-32 I get this result in one of my folders, 截图 2023-11-25 00-32-11 This is my log file, GACOSlog.txt Then I tried to use only one subband. the SLC file, while resolving “the number of SLC files is inconsistent ”at the beginning, 截图 2023-11-25 00-43-11 截图 2023-11-25 00-43-27 but can not produce the final result. This is my second log file GACOSlog2.txt

MinervaLee0429 commented 1 year ago

@esrivas17 Hello Professor, you are experiencing a situation where phasefilt_gacOS_Correct.grd and phasefilt_GACOS_corrected_detrended.grd cannot be generated These two results have problems after the generation of szpddm_phase.grd and szpddm_phase_LOS.grd. I don't know how to solve them. Please kindly give me your guidance.

esrivas17 commented 1 year ago

Looking at the folder you shared. You have 6 .SLC files. The scripts are designed to work with one subswath therefore it's expecting only two of these files. (1 per date). That's why you are getting the message that the SLC files are not consistent. I should work on including the option to process more than one subswath but for the moment you could leave only two SLC files (either F1, F2 or F3) and see if you can get the final grids. Also make sure your dem is larger than the coverage of GACOS. Maybe you can comment out in operacion.csh lines 121 to 125, run the whole script again and send an image of the one of the folders to see what files are generated

MinervaLee0429 commented 1 year ago

@esrivas17 Thank you for your response, Professor. I'll try to deal with it and see what happens.