asc_desc2horz_vert & reprojection

205 views
Skip to first unread message

David Nguyen

unread,
Nov 25, 2024, 11:15:26 PM11/25/24
to MintPy
Dear all,

I've obtained interferograms for both ascending & descending paths as follows:
- Ascending path:  Full-scene interferograms processed via ASF's HyP3 service (asf.alaska.edu).
- Descending path: Interferograms produced using hyp3-isce2 because my study area is covered by two separate full-scene SLCs, so I utilized insar_tops_burst.

However, I encountered an issue where the ascending products are in UTM Zone 51N, while the descending products are in UTM Zone 50N. I also uploaded the text files containing the output of info.py velocity.h5 for your reference.

I’m wondering if there are any available functions in the MintPy package that can help align the coordinates of these two datasets.

Alternatively, if I need to write a custom script to perform the reprojection, I would greatly appreciate any hints or suggestions to get started.

Thank you for your time and assistance!


descending_velocity.txt
ascending_velocity.txt

David Nguyen

unread,
Dec 1, 2024, 10:00:12 PM12/1/24
to MintPy
I figured out the solution myself.

When you want to perform the asc_desc2horz_vert, you need to check the information of input .tif files before loading your data into MintPy. You can apply gdalinfo to do so, and check if the coordinate system between ascending & descending files are similar or not.

If they are not the same (as in my case), you can write a program to apply gdalwarp to reproject the coordinate systems of the .tif files, then load them into MintPy.

I uploaded my naïve Python program to perform this task below for your reference.
reproject_tif.py

Gabriela Quintana Sánchez

unread,
Jan 30, 2025, 1:06:02 PM1/30/25
to MintPy
Hi David, I have a similar problem, I processed my data with HyP3, when I obtained my velocity map it is in .h5 format, thus, I cannot apply asc_desc2horz_vert.py code, because  I cannot geocoded my input!!!... my velocity.h5 in ascending and descending...

Sayyed Mohammad Javad Mirzadeh

unread,
Jan 30, 2025, 2:04:46 PM1/30/25
to MintPy
Hi Gabriela,

You can easily use the "CutGeotiffs.py" code from the Hyp3lib library to convert all HyP3 service products from the UTM to the WGS84 coordinate system and go ahead with other steps as same as what you did with the ISCE-processed input files.

Sincerely,
Sayyed



Gabriela Quintana Sánchez

unread,
Jan 31, 2025, 6:10:10 PM1/31/25
to MintPy
Thanks Mr. Sayyed

Esubalew Adem

unread,
Feb 9, 2025, 3:01:24 AM2/9/25
to min...@googlegroups.com
Dear David,
I have a similar issue: my study area is covered by two separate full scenes, and I want to use hyp3-isce2. I am new to hyp3-isce2; could you please guide me on how to download bursts, merge them, and process the time series?

--
You received this message because you are subscribed to the Google Groups "MintPy" group.
To unsubscribe from this group and stop receiving emails from it, send an email to mintpy+un...@googlegroups.com.
To view this discussion visit https://groups.google.com/d/msgid/mintpy/0f468b88-aa3f-498a-8c47-852c10efb7f0n%40googlegroups.com.


--

Esubalew Adem

Reply all
Reply to author
Forward
0 new messages