Problem Plotting standard deviation in matlab

s1tbx StaMPS
#1

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

image
image919×391 4.89 KB

ps_plot(‘vs’)
Calculating standard deviation of mean velocity…
4 ref PS selected
test
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)
fig_name_tca=ps_mean_v(ifg_list,1500,value_type(4:end),use_small_baselines,aps_flag);

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

#3

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

#4

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.

#5

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

#6

Sure @1027042853

#7

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

#8

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

#9

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
#10

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

#11

you have to create this file first by plotting the data

grafik

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

1 Like
#12

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

v

grafik
grafik1920×970 203 KB

v-dao(gacos)

grafik
grafik1920×970 207 KB

vs-dao(gacos)

grafik
grafik1920×970 203 KB

#13

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?

#14

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.

#15

-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:

#16

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

#17

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

#18

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.