search

GeoscienceAustralia / PyRate

A Python tool for estimating velocity and time-series from Interferometric Synthetic Aperture Radar (InSAR) data.
https://geoscienceaustralia.github.io/PyRate/
Apache License 2.0
200 stars 70 forks source link

conv2tif does not work with different X and Y pixel resolutions #305

Closed mcgarth closed 3 years ago

mcgarth commented 3 years ago

Hi

I am wanting to convert from the matlab-based version of pyrate to this python-based version. I am using ISCE products, and am running into issues with the conv2tif step. Is this the appropriate place to raise questions about specific errors?

Thanks for any help, and thank you for providing an awesome github repo

Originally posted by @Jesse-kearse in https://github.com/GeoscienceAustralia/PyRate/issues/40#issuecomment-716085052

mcgarth commented 3 years ago

Hi @Jesse-kearse, if you give some details here of the problem you are encountering with conv2tif we can help look in to it. Nominally, it sounds like if you've generated .rsc header files from your ISCE products it should work in PyRate

Jesse-kearse commented 3 years ago

Thanks for your response, mcgarth.

As part of the conv2tif process, the module "write_fullres_geotiff" within pyrate.core.shared.py calls "check_pixel_res_mismatch" https://github.com/GeoscienceAustralia/PyRate/blob/eb4a4440672e02fafd8820b684a7684d02f41b9f/pyrate/core/shared.py#L896

It seems to want the x-step and y-step values to be the same, yet the dem that I am using for geocoding has different x and y steps. Will I need to use a "square" dem grid instead? example of geocoded interferogram .vrt file from ISCE

<VRTDataset rasterXSize="23179" rasterYSize="10601">
    <SRS>EPSG:4326</SRS>
    <GeoTransform>171.06490487583588, 0.00010261051442909, 0.0, -42.99785174151206, 0.0, -0.0001027265486936</GeoTransform>
    <VRTRasterBand dataType="Float32" band="1" subClass="VRTRawRasterBand">
        <SourceFilename relativeToVRT="1">filt_topophase.unw.geo</SourceFilename>
        <ByteOrder>LSB</ByteOrder>
        <ImageOffset>0</ImageOffset>
        <PixelOffset>4</PixelOffset>
        <LineOffset>185432</LineOffset>
    </VRTRasterBand>
    <VRTRasterBand dataType="Float32" band="2" subClass="VRTRawRasterBand">
        <SourceFilename relativeToVRT="1">filt_topophase.unw.geo</SourceFilename>
        <ByteOrder>LSB</ByteOrder>
        <ImageOffset>92716</ImageOffset>
        <PixelOffset>4</PixelOffset>
        <LineOffset>185432</LineOffset>
    </VRTRasterBand>
</VRTDataset>
Jesse-kearse commented 3 years ago

sorry the file did not print well in the above comment:

EPSG:4326 171.06490487583588, 0.00010261051442909, 0.0, -42.99785174151206, 0.0, -0.0001027265486936 filt_topophase.unw.geo LSB 0 4 185432 filt_topophase.unw.geo LSB 92716 4 185432
mcgarth commented 3 years ago

I edited your comment, and wrapped the .vrt content in triple backquotes, which solves the display issue

Jesse-kearse commented 3 years ago

Thankyou Matt

I hope the issue I was trying to explain makes sense to you.

Cheers

On Mon, Oct 26, 2020 at 5:55 PM Matt Garthwaite notifications@github.com wrote:

I edited your comment, and wrapped the .vrt content in , which solves the display issue

— You are receiving this because you were mentioned. Reply to this email directly, view it on GitHub https://github.com/GeoscienceAustralia/PyRate/issues/305#issuecomment-716302360, or unsubscribe https://github.com/notifications/unsubscribe-auth/AQHOJLRCTOW6VRIKWMTAD2DSMT6LRANCNFSM4S62J5IQ .

--

Jesse Kearse PhD Student Victoria University of Wellington

mcgarth commented 3 years ago

Yes I understand what the issue is @Jesse-kearse and will look in to it with the team. Cheers

mcgarth commented 3 years ago

Hi @Jesse-kearse, I have raised PR https://github.com/GeoscienceAustralia/PyRate/pull/306 which I hope will sort things out for you. It would be great if you could test that for me and let me know if you are able to successfully get through the conv2tif step.

mcgarth commented 3 years ago

https://github.com/GeoscienceAustralia/PyRate/pull/306 has been merged in to master branch

Jesse-kearse commented 3 years ago

Hi Matt

I will explore these changes in the next day or so and get back to you. Thanks for your help

On Wed, 28 Oct 2020, 11:04 Matt Garthwaite, notifications@github.com wrote:

306 https://github.com/GeoscienceAustralia/PyRate/pull/306 has been

merged in to master branch

— You are receiving this because you were mentioned. Reply to this email directly, view it on GitHub https://github.com/GeoscienceAustralia/PyRate/issues/305#issuecomment-717567518, or unsubscribe https://github.com/notifications/unsubscribe-auth/AQHOJLSXSWUCCWHLHQCT6ZTSM47YHANCNFSM4S62J5IQ .

mcgarth commented 3 years ago

Hi @Jesse-kearse, did you have a chance to check whether the changes I made sorted your problem?

Jesse-kearse commented 3 years ago

Hi Matt

Sorry, I have been caught up in other research - I will run this shortly and get back to you

On Mon, Nov 16, 2020 at 6:22 PM Matt Garthwaite notifications@github.com wrote:

Hi @Jesse-kearse https://github.com/Jesse-kearse, did you have a chance to check whether the changes I made sorted your problem?

— You are receiving this because you were mentioned. Reply to this email directly, view it on GitHub https://github.com/GeoscienceAustralia/PyRate/issues/305#issuecomment-727741199, or unsubscribe https://github.com/notifications/unsubscribe-auth/AQHOJLTN3FNMVHSTOLTUVOLSQCZH7ANCNFSM4S62J5IQ .

--

Jesse Kearse PhD Student Victoria University of Wellington