Problem Plotting standard deviation in matlab

I am having problem with plotting standard deviation. how can i solve this problem any idea?

Calculating standard deviation of mean velocity…
4 ref PS selected
The logical indices in position 1 contain a true value outside of the array bounds.

Error in ps_mean_v>@(x)single(lscov(G(x,:),ph_bit(x,ix_non_nan),ifg_cov(ix_non_nan,ix_non_nan))) (line 280)
[mean_v_dist_temp,boot_ix] = bootstrp(n_boot, @(x)
single(lscov(G(x,:),ph_bit(x,ix_non_nan),ifg_cov(ix_non_nan,ix_non_nan))), [1:N]);

Error in bootstrp (line 178)
bootstat = feval(bootfun,bootargs{:});

Error in ps_mean_v (line 280)
[mean_v_dist_temp,boot_ix] = bootstrp(n_boot, @(x)
single(lscov(G(x,:),ph_bit(x,ix_non_nan),ifg_cov(ix_non_nan,ix_non_nan))), [1:N]);

Error in ps_plot (line 1835)

Caused by:
Unable to evaluate BOOTFUN with the supplied arguments.

i also have the question,do you have resaulted it?please tell me .thank you

not yet!!
It works with one dataset and not with another, may there is some fault in data. I re processed the multiple times with StaMPS but still the same error arises. I find a lot of NaN values in the exported time series data too. These NaN values might be the problem but i dont know how to sort it out.

I also have many Nan values. I deleted them. If you solve them, teach me. Thank you

Sure @1027042853

I have the same problem, how should I solve it.

I also encountered such a problem, with a lot of NaN values.

Just replace your stamps code, with below given StaMPS code

It’s working good. Try once.

Still it is giving error, please process once again mt_prep_snap and stamps(1,8)

Good luck.

1 Like

Dear @ABraun,
Which file (after successful completion of Stamps in Matlab) contains the Standard deviation values for the LOS velocity plot?

you have to create this file first by plotting the data


In your case it would be ps_plot('vs',-1)
Source StaMPS/ps_plot.m at master · dbekaert/StaMPS · GitHub

1 Like

Dear @ABraun ,
I just started using STaMPS and this is the rate graph I got with PSI after isce2 preprocessing, I am not able to explain this result( i.e. the settling point is scattered). PSI result figure v v-dao vs-dao…
Is this normal? :raising_hand_man:
When I processed with SBAS (multilook20X5) there were not so many discrete settling points.
Help me :pray: :pray: :pray:
Thanks in advance




could you please describe what you mean here?

The result looks okay at first sight, but of course you have to know the area to determine if these are feasible.
vs is the standard devisation of displacement which is lowest in the center and highest at the edges - this looks a bit too systematic to me but I can’t tell where it comes from. Could it have been introduced after applying GACOS?

Dear @ABraun ,
Thank you very much for your reply!
First of all, my AOI about 1500 square kilometers, located in Guiyang City, China, and I processed 72 S1 , during three years. full resolution,not multilook,about 270000 point.The current results do not introduce reference points.
From v-dao results, the subsidence area is not aggregated, but scattered, so from my current experience, I am not sure about the accuracy of this map.
VS-dao was drawn after step 7, I don’t know if this is correct.
Forgive my ignorance,Thanks in advance.

-vdao means you remove DEM errors as well as atmospheric and oribtal trends. Depending on how well this works, your data can be improved but in some cases new errors are introduced.

I cannot tell you which of the results is most feasible because I don’t know the study area. Explanation and interpretation of the result is your duty :slight_smile:

Dear @ABraun ,
Thank you!
next step, I will use GNSS observation data to check…

Maybe you are lucky here: MAGNET + Global GPS Network Map

Dear @ABraun ,
Thank you sir for your help! :pray:
I understood this information earlier, but my AOI does not cover,
I will come back to share good news when I have it.

I am still struggling with the methodology on how to apply such corrections either to the raw interferograms before exporting them into StaMPS processing or alongside with it after applying phase unwrapping. I was wondering if you have any tutorial to guide users through that correction process or if you could direct me into a step-by-step manual for how to use GACOS’s Binary-grid or Geotiff products with Interferograms. Most of what I know for now is that we need to do time differing between two delays maps of two dates and space differing for the resulted Zenith Path Delay Difference Map (ZPDDM) by subtracting the delay value of a particular reference lat/long location on (ZPDDM) from the remaining phase delay values from the same entire map. Hence, the differenced in time and space map should be either subtracted or added to the interferograms depending on whether later or earlier is the master date for each interferogram. The problem is that I do not know how to practically apply such process using matlab nor SNAP toolbox if possible, and based on what may a choose the previously mentioned reference point ?