dbekaert / StaMPS

Stanford Method for Persistent Scatterers
http://homepages.see.leeds.ac.uk/~earahoo/stamps/
GNU General Public License v3.0
233 stars 120 forks source link

Step 2 in StaMPS #15

Open redagalal opened 5 years ago

redagalal commented 5 years ago

Hi, I'm working on StaMPS in step-2 I got the following error: Error in ps_est_gamma_quick (line 221) ph_grid(grid_ij(i,1),grid_ij(i,2),:)=ph_grid(grid_ij(i,1),grid_ij(i,2),:)+shiftdim(ph_weight(i,:),-1);

Screen Shot 2019-06-29 at 2 21 37 PM

I need your help to solve it.

hruol commented 2 years ago

Hi, i meet the same error in step2, had u resolve it? x_x

constantinevi6 commented 2 years ago

Hello, To deal with this problem, you have to check if the patch really have PSCs selected by mt_prep, for none PSC case, stamps will have error like this. An easy way to check PSCs is looking at the file "pscands.1.ij" in PATCH_X folder, if the file size is 0, or it has only few rows of data (less than 100 in experience), please skip this patch by edit "patch.list" in main PS folder, removing PATCH_X line depend on you case.

hruol commented 2 years ago

thanks for your reply, but my pscands.1.ij file looks normal... Screenshot from 2022-03-31 11-39-02

constantinevi6 commented 2 years ago

Hello, After that, please check geocode of your image, if you use SNAP to processing interferogram, please check the orthorectified Lon and Lat band if there are normal. Any NaN data in Lon or Lat will caused problem in my experience. If it still looks fine, I can figure out what happen in your data, It's hard to debug without having that dataset.

hruol commented 2 years ago

Hello, After that, please check geocode of your image, if you use SNAP to processing interferogram, please check the orthorectified Lon and Lat band if there are normal. Any NaN data in Lon or Lat will caused problem in my experience. If it still looks fine, I can figure out what happen in your data, It's hard to debug without having that dataset.

well I found the NaN data in my orthorectified Lon and Lat, but the distribution of these empty areas looks regular, im so curious about that.

clinomaniacop commented 2 years ago

Hey! @constantinevi6 Everything seems to be right in my dataset. However , I am still encountering the same issue. The PSCs are 3852. And The ortho rectified lat and long are normal as well. Any idea why I am still having the problem?

clinomaniacop commented 2 years ago

@constantinevi6 I have again checked the whole stack, neither the interferograms nor the orthorectified lat long has any Nan values. I am still getting this error.

constantinevi6 commented 2 years ago

Hi, if you still have problem, it should be something else I can't help, sorry. In this case, may be you should try to check every variable you use in that line to find the bug directly.

darenhusam commented 1 year ago

Hi, I'm working on StaMPS in step-2 I got the following error: Error using signal.internal.sigcheckfloattype Can't reload '/usr/local/MATLAB/R2018aaa/bin/glnxa64/libmwxcp_dwarf.so'

Error in signal.internal.sigcasttofloat (line 36) signal.internal.sigcheckfloattype(x,'', fcnName, varName, datacheckflag);

Error in gausswin (line 40) a = signal.internal.sigcasttofloat(a,'double','gausswin','ALPHA',...

Error in clap_filt (line 47) B=gausswin(7)*gausswin(7)';

Error in ps_est_gamma_quick (line 225) ph_filt(:,:,i)=clap_filt(ph_grid(:,:,i),clap_alpha,clap_beta,n_win0.75,n_win0.25,low_pass);

Error in stamps (line 327) ps_est_gamma_quick(est_gamma_parm); image

hruol commented 1 year ago

Hi, I'm working on StaMPS in step-2 I got the following error: Error using signal.internal.sigcheckfloattype Can't reload '/usr/local/MATLAB/R2018aaa/bin/glnxa64/libmwxcp_dwarf.so'

Error in signal.internal.sigcasttofloat (line 36) signal.internal.sigcheckfloattype(x,'', fcnName, varName, datacheckflag);

Error in gausswin (line 40) a = signal.internal.sigcasttofloat(a,'double','gausswin','ALPHA',...

Error in clap_filt (line 47) B=gausswin(7)*gausswin(7)';

Error in ps_est_gamma_quick (line 225) ph_filt(:,:,i)=clap_filt(ph_grid(:,:,i),clap_alpha,clap_beta,n_win_0.75,n_win_0.25,low_pass);

Error in stamps (line 327) ps_est_gamma_quick(est_gamma_parm); image

Your MATLAB version is 2018a, I think you can try to uninstall this version and install the new version.

DD-ZHAO commented 9 months ago

well I found the NaN data in my orthorectified Lon and Lat, but the distribution of these empty areas looks regular, im so curious about that.

@constantinevi6 Hello, thank you very much for your advice, I checked my pscands.1.ij file and it looks normal. I found the NaN data in my orthorectified Lon and Lat in SNAP, It has a regular distribution as shown below.

555

I would like to ask you how to solve the stamps(2,2) bug, thank you again!

constantinevi6 commented 8 months ago

well I found the NaN data in my orthorectified Lon and Lat, but the distribution of these empty areas looks regular, im so curious about that.

@constantinevi6 Hello, thank you very much for your advice, I checked my pscands.1.ij file and it looks normal. I found the NaN data in my orthorectified Lon and Lat in SNAP, It has a regular distribution as shown below. I would like to ask you how to solve the stamps(2,2) bug, thank you again!

Hello, you can try to kick out the PSCs that have NaN in Lon and Lat, simply check the "lonlat" in ps1.mat and remove them, other files and variables also need to apply the same changes.

Or, you have to prevent SNAP output elevation with such a problem, the problem is from the DEM you use. I notice that the interferogram across the coast sometimes has this problem, while the DEM has set ocean as no data, you can try to use SRTM as input and see if the problem is solved.