AlexeyPechnikov / pygmtsar

PyGMTSAR (Python InSAR): Powerful and Accessible Satellite Interferometry
http://insar.dev/
BSD 3-Clause "New" or "Revised" License
417 stars 91 forks source link

[Help]: Gaps in the displacement files #147

Open georgeboldeanu opened 1 month ago

georgeboldeanu commented 1 month ago

Hello! I do a sbas analysis over some area of Philippines. After I computed the displacement with the correction for trend and 2 turbulance phases I noticed some gaps in the resulted data (rectangular gaps, circled in red). It's this a known issue? I didn't get any error message or something like that. Screenshot 2024-07-12 104803

For a smaller area using SBAS and PS combination the results were clear without any gaps. It can be something resource related? The resulted displacement stack is not that big, with around 17GB. image

Thank you!

AlexeyPechnikov commented 1 month ago

It seems you have filtered out some low coherence areas like water surfaces, but not completely, as there are still invalid borders around them.

georgeboldeanu commented 1 month ago

@AlexeyPechnikov I forgot to mention that I am referring to the rectangular gaps. Sorry for the omission. Also, I mention that all the data that enters in the computation of the displacement doesn't have those gaps, so I guess is from the processing of the displacement.

AlexeyPechnikov commented 1 month ago

Trend correction is per-chunk in the spatial domain, so if you do not have enough pixels in the chunks, the correction is not possible for some chunks (prezented as rectangular gaps on the plots). It’s better to use a single chunk for this correction, and you can filter out low-correlation values for faster processing on large grids.

georgeboldeanu commented 1 month ago

@AlexeyPechnikov But those areas that have gaps are in an urban area with high coherence, so no filtering is done. How is possible that there are not sufficient pixels in an area like this? If you look at the sample image, in the upper part is understandable, but in the middle part I don't understand the behavior.

AlexeyPechnikov commented 1 month ago

You need to increase the chunk size and compare the results.