How to get geotiff of single wrapped (and unwrapped) Interferograms from MintPy

[Chui Samathys]

Hi all

I'm using Mintpy after generating a stack of Interferograms with ISCE. I'm still quite new to it and I'm wondering how I can get single Interferograms in a GeoTiff format? I sometimes need them for data analysis (without the unwrapping) and I can't figure out how to get that. I see that it's possible to get all wrapped Inteferograms automatically (unwrapPhase_wrap.png) but I would like to get them as single GeoTiffs. 

Any help would be greatly appreciated. Also if someone know how I could get those from the ISCE output directly. Ive been struggling since a long time to get useable geocoded (.geo) files from my stackSentinel.py outputs. 

Thanks so much already :)

Chui

[Manudeo Singh]

You can use Mintpy's geocode.py to geocode ISCE products directly. You’ll need to use lat.rdr and lon.rdr from merged/geom_reference for --lat-file and --lon-file 

Example for height data:

cd merged/geom_reference

geocode.py hgt.rdr --lat-file lat.rdr --lon-file lon.rdr -o geo_height.h5

save_gdal.py geo_height.h5 -o geo_height.tif

Best,

Manu

_________________________________
Dr Manudeo Singh
Royal Society - Newton International Fellow
Department of Geography & Earth Sciences
Aberystwyth University, Wales, UK

On Mon, 9 Sep 2024 at 14:53, Chui Samathys <chuisa...@gmail.com> wrote:

Hi all

I'm using Mintpy after generating a stack of Interferograms with ISCE. I'm still quite new to it and I'm wondering how I can get single Interferograms in a GeoTiff format? I sometimes need them data analysis (without the unwrapping) and I can't figure out how to get that. I see that it's possible to get all wrapped Inteferograms automatically (unwrapPhase_wrap.png) but I would like to get them as single GeoTiffs. 

Any help would be greatly appreciated. Also if someone know how I could get those from the ISCE output directly. Ive been struggling since a long time to get useable geocoded (.geo) files from my stackSentinel.py outputs. 

Thanks so much already :)

Chui

--
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 on the web visit https://groups.google.com/d/msgid/mintpy/b2cc4c29-6c7f-46a6-a0b2-cd06862edfa3n%40googlegroups.com.

[Chui Samathys]

OMG thanks so much!!!! Saving me so much time :) Super grateful for the fast response! Thanks a lot!

All my best

Chui

[Rob Sharifi]

Hi Chui,

Just out of curiosity, are there any reasons that you don't do it on ISCE products using gdal_translate?

Sohrab (Rob)

[Chui Samathys]

Hi Rob

I just honestly couldn't find out how to do it with gdal_translate only since my end product from ISCE is not geocoded. I'm using stackSentinel.py, so the final files that I get are the ones you can see in the screenshot attacked. But any help with this is greatly appreciated :) If you know another way, I'm keen to learn!

Cheers

Chui