GACOS CORRECTION

[kquresh2]

Hi,
I have been running script for multiple areas and it worked fine. For one area, it has the issue of file sizes, so I was able to get help from previously answered questions. Right now, I am facing this issue that after GACOS correction (step-0-3), slc.mli.par file is not generated due to which the next step of unwrapping (step-0-4) also dont work. I checked and found that file size is not an issue. I even tried to re-download dataset but the issue is not resolved. I have attached error screenshot that appears after the GACOS correction completes. Also, I have attached screenshot of nan values for some Ifg in GACOS text files. Can you please help with this issue? Thanks

[kquresh2]

@yumorishita Looking for your help related to the above query. Thanks

[yumorishita]

I am not sure but the interferograms with nan in GACOS_info.txt might be corrupted. A similar problem has occurred before: #155 (comment). I suggest removing the suspicious interferograms and retrying step03.

[kquresh2]

@yumorishita Thanks for the quick reply. I was thinking the same to remove bad interferograms using GACOS info file. I did removed all bad ifg's from GEOC folder and even deleted GEOCml1 folder and started again from step 0-2. After doing so, I did not get any error or NaN values after GACOS correction (step 0-3). Now, the process is going smooth. Cheers!

[kquresh2]

@yumorishita I have all my results but they are raster images. I want to convert final displacement to feature layer so that I can do hotspot analysis with vector data. Is there any script or workout for that? Thanks

[yumorishita]

What do you mean by "feature layer"? Maybe LiCSBAS_flt2geotiff.py and gdal tools (e.g., gdal2xyz.py) can help.

[kquresh2]

@yumorishita by feature layer, I mean conversion from raster to point feature layer. Once I have the point data with displacement, I can run hotspot analysis tool to see the displacement pattern. I used to do SBAS in sarscape which give me feature layer with velocity point rather than raster layer. LiCSBAS_flt2geotiff.py output is single raster image which can be open in ArcGIS Pro but it is not what I am looking for.

[yumorishita]

You should be able to convert the raster data into vector data (i.e., feature layer with velocity point) in ArcGIS, though I am not familiar with ArcGIS.
https://desktop.arcgis.com/en/arcmap/latest/extensions/3d-analyst/converting-raster-surfaces-to-features.htm
https://desktop.arcgis.com/en/arcmap/latest/extensions/3d-analyst/converting-surface-data-to-vector-data.htm
https://desktop.arcgis.com/en/arcmap/latest/tools/conversion-toolbox/raster-to-point.htm

[kquresh2]

@yumorishita I did Licsbas processing for multiple scenes (around 9 ascending) when I try to mosaic them all using gdalwarp tool or in ArcGIS Pro, there is edge effect and the color of the final product is not same because every scene has different min, max, and mean values. Do you know how to fix it and bring the whole mosaic map on a same color scale? I see some people in COMET have mosaic multiple images for Anatolia fault and in Tibetan plateau. Looking forward for recommendations. Thanks

[yumorishita]

#147
#24

[kquresh2]

@yumorishita cum_filt.h5 file do not have coordinate system assigned to it. Do you know how we can assign it a coordinates to it? I am not interested in one final Geotiff image which we can do by using LiCSBAS_flt2geotiff.py tool. I want to work with cum_filt.h5 file. Please let me know if you know how to do that. Thanks