-
-
Notifications
You must be signed in to change notification settings - Fork 88
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
[Help]: PSI and SBAS geocoding #101
Comments
You open internal processing files which have a special format. Geocode and export the files as explained in the examples. |
So, is my spatial resolution output correct? |
This is the rasters pixel spacing, not meters. |
Well @AlexeyPechnikov I managed to get the GoldenValley notebook working. However, I have a question. I want to separate my sbas and psi results, that is, my LOS changes, as vertical and east-west. Which function will I do this with? |
@AlexeyPechnikov Also, after a while, I get this error in the middle of the code and I cannot proceed to the result. What is the reason for this? I use Colab Pro, but I encounter this error when processing big data. Problems with the session usually start after this code: ERROR:asyncio:Task exception was never retrieved |
'detrend' is your unwrapped and detrended phase. You can miss the detrending step if you don't need it for your case. 'asyncio' is just a web sockets communication library and you can ignore the message. |
@AlexeyPechnikov ,I have one last question, in this notebook I produced the LOS displacement map using my images in the ascending direction. How do I decide if this data is an up-down or east-west change? I tried parsing ud - ew using methods in other notebooks, but the results were too big and ridiculous.
|
To obtain the actual values, you need to sum the west-east or vertical displacements from two orbits, or employ an alternative method to determine the projections. The 'detrend' variable contains the detrended, unwrapped phase. You can utilize the unwrapped phase directly or subtract the linear regression trend from the phase, among other techniques. |
`--------------------------------------------------------------------------- 5 frames ValueError: dimension 'pair' already exists as a scalar variable` @AlexeyPechnikov ,What is the reason for this error? |
I have developed a vertical and east-west decomposition based on an algorithm from published literature. It uses a search radius to combine points onto a lower resolution grid and decomposes the projections based on local incidence angle derived from a DEM which takes into account the slope. It does however depend on SNAP to create the local incidence angle, and using other tools such as GDAL did not give me the expected results. It would be interesting to see how you would approach this. |
The task should be resolved by using proper geocoding on the same grid instead of later interpolations. PyGMTSAR produces geocoded results aligned with DEM pixels, meaning that when using the same DEM, we have identical output grids for two orbits. Simply sum the grids to produce the actual up-down and east-west components. |
The 'stack' dimension is incorrect. It appears that you are applying 3D processing code to a single interferogram, or vice versa. |
Hello @AlexeyPechnikov,
I produced displacement maps in a region for PSI and SBAS in a notebook called GoldenValley. However, my output results do not have a coordinate system, they appear as row and column indexes, as in the image. How can I do that. I need georeferencing for my disp_ps.grd and disp_sbas.grd result files. Another question I have is that when I open the output results, the spatial resolution is 1 meter.
Also, the PSI result image appears differently, I cannot understand whether there is an error.
The text was updated successfully, but these errors were encountered: