yumorishita / LiCSBAS

LiCSBAS: InSAR time series analysis package using LiCSAR products
https://doi.org/10.3390/RS12030424
GNU General Public License v3.0
229 stars 110 forks source link

"All pixels are masked" error #415

Open geomaticss opened 1 day ago

geomaticss commented 1 day ago

first of all hello there. i was going to get my geoclip file as a result of a very long time process but i realized that all pixels were masked. you can see the relevant part of the code below. what could be the reason? how can i fix it? also do i need to redo the whole process? it really took a lot of time...

".... Elapsed time for 1th patch: 106 sec

Output png images...

Elapsed time: 00h 01m 46s

LiCSBAS14_vel_std.py Successfully finished!!

Output directory: TS_GEOCml1clip

LiCSBAS15_mask_ts.py ver1.8.1 20200911 Y. Morishita LiCSBAS15_mask_ts.py -t TS_GEOCml1clip

All pixels would be masked with n_ifg_noloop thre of 50 Automatically change the thre to 52.0 (ceil of min value)

Noise index : Threshold (rate to be masked)

Masked pixels : 63509/63510 (100.0%) Kept pixels : 1/63510 (0.0%)

All pixels are masked!! Try again with different threshold."

I worked with a dataset that is almost 10 years old so there are a lot of image pairs etc. I had previously experimented with 6-7 year old datasets in different regions and did not encounter such a problem. thank you in advance for your reply.

yumorishita commented 18 hours ago

Try increasing p15_n_ifg_noloop_thre to >52 and rerun from step15. https://github.com/yumorishita/LiCSBAS/blob/master/batch_LiCSBAS.sh#L57