StaMPS-Visualizer, SNAP-StaMPS Workflow

[ABraun]

you can remove a couple of the lines from the end or upload it somewhere outside the forum.

[yzxtll]

thank you ,there are more than one hundred thousand PS points,now 30 thousand remain(remove from end)stamps_tsexportD.csv (2.7 MB)

[yzxtll]

remove from start,include the example above (eg:water):stamps_tsexportE.csv (3.2 MB)

[yzxtll]

configure: error: gdal-config not found or not executable.
ERROR: configuration failed for package ‘rgdal’
I run as @thho:

[quote=“thho, post:300, topic:9613”]
sometimes I am able to use Rstudio`s gui when installtion via command does not work…lower right window, Packages tab–>klick install–>search for rgdal–> klick install
[/quote]

[thho]

please have a look here:

StaMPS-Visualizer, SNAP-StaMPS Workflow

@CFEgildan71, I assume you are using a Unix OS, then open a terminal and do: sudo apt install libgeos-dev after that open RStudio again install.packages('rgeos') should do.

and here:

StaMPS-Visualizer, SNAP-StaMPS Workflow

For all having the same problem, use this as a starting point:

[ABraun]

you have to replace the days since year 0 by the formatted date

Then you import the points from the csv

Then you save the imported CSV as a shapefile.

Then the tool works

But I recommend to increase the number of input images to get more reliable results.
There is no clear trend in the single points and the overall image shows a strong ramp

stamps_tsexportE.zip (2.9 MB)

[yzxtll]

thnk you for your help.1/yes ,I just have a try about the scripts within 6 images to know the whole workflow .2/ compare the map between auto navi map and open street map, it seems a little shift about the each coordinate,this result in my PS points in water as showed above.3/ I download the file stamps_tsexportE,it runs successful ,very thanks for you .I’ll do further exercise to learn.

[ABraun]

a possible explanation for smaller shifts is given here:

StaMPS-Visualizer, SNAP-StaMPS Workflow

One explanation could be that your data is still in slant range geometry. What seems like a square in your SLC product actually covers a long stripe due to different pixel resolutions in azimuth and range direction. Have a look at these sources: https://earth.esa.int/web/guest/missions/esa-operational-eo-missions/ers/instruments/sar/applications/radar-courses/content-2/-/asset_publisher/qIBc6NYRXfnG/content/radar-course-2-slant-range-ground-range [image] …

[yzxtll]

hello,viStaMPS 2.1 subset progrom error :as below,I compared the two csv files structure ,one is direct ly exported from matlab2018 (stamps_tsexport.csv)the other is the example of viStaMPS 2.1(Maoxian_PS.csv).the question is what’s 2row 1:2col mean (before NaN in Maoxian_PS.csv)?

[ABraun]

I think you are mixing up the packages:

The R scripts developed by @thho are called the StaMPS Visualizer

The Matlab GUI developed by UTAD is called ViStaMPS

About the NaN column - this is explained here:

StaMPS-Visualizer, SNAP-StaMPS Workflow

The csv table is not that intuitive, I have tried to make it as compact as possible… [1,1:2] = centre of selection lon lat [1,3] = Na [1,4:ncol] = the date in days since 0 of the image, the intervals do not have to be equal, it depends on the image date. That [1,3] shows a NA value is because, the values beneath it are not part of the time series, these are the mean velocity values, used for colour coding! [2:nrow, 1:2] Measurement Points lon lat [2:nrow, 3] mean velocity [2:nrow, 4:ncol]…

[yzxtll]

thank you for your mention,I indeed mix up two plugins.Another question ,how can I sort the result of csv file or shp file ,which related to [2:nrow, 3] mean velocity ?

[ABraun]

in QGIS you can sort columsn of an attribute table by clicking on their name. https://docs.qgis.org/2.18/en/docs/user_manual/working_with_vector/attribute_table.html#sorting-columns

[yzxtll]

HI,I’ve read through https://docs.qgis.org/2.18/en/docs/user_manual/working_with_vector/attribute_table.html#sorting-columns .if there any other materials related , recommend to me please!

[ABraun]

can you describe a bit more what you mean by “sort the result of csv or shp”?

[yzxtll]

sorry for my poor english,what I want is different velocity rank in different level ,just as -100 --15,-15 --2.5,-2.5 -2.5,etc.

[yzxtll]

@thho

use a ref point, this might highly improve your results…since it is set to 0,0 as reference the overall velocity of all points is used in StaMPS as baseline…

hello!@thho:I dont know how to choose a better ref point in ps_plot frame ,that is just I wondering about ,i’ll read line by line blow,

[yzxtll]

@thho

you have to define the ref point before plotting

I’ve read throug P32 of StaMPS_Manual_v4.1b1.pdf,It seems said how to defined a ref area. I always define a ref point in ps_plot frame as random and don’t know which point is the stable one among dense candidate points.i know I totally wrong about a ref point and dont know how to correct?

[badrulhudahusain]

Hi…have a good day @thho @ABraun… I want to try to process the ALOS 2 Data for PSInsar using SNAP-StaMPS Workflow, do you know ho to prepare the ALOS 2 data…??.. is it just the same with Sentinel 1 Data Preparation…??.. Thank You in Advance…

[ABraun]

It is actually much easier, because ALOS-2 is Stripmap and does not need Split, Deburst, ect…

You simple create one stack with the coregistered products (first input for the StaMPS export) and apply interferogram generation (including topographic phase removal) on this stack. This is the second input for the StaMPS export.

[badrulhudahusain]

Thank you @ABraun …I’ll share the result later​:ok_hand: