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
203 stars 71 forks source link

Metadata addition: LOS direction in GeoTIFF metadata #359

Closed adeane-ga closed 2 years ago

adeane-ga commented 2 years ago

Problem When users are interpreting PyRate velocity outputs, it would be useful to quickly be able to query in which direction surface motion can/cannot be observed based on the direction of the LOS observations. Currently a user has to query the azimuth.tif file to obtain this information. It would be convenient if one could access the median azimuth angle in the metadata of the GeoTIFFs within the velocity_dir.

Request

velocity_dir/linear_rate.tif velocity_dir/linear_error.tif velocity_dir/linear_rsquared.tif velocity_dir/linear_samples.tif velocity_dir/linear_error.tif

But potentially in the time series files within timeseries_dir directory too.

mcgarth commented 2 years ago

I believe the average incidence angle is already saved to the interferogram geotiff metadata?

adeane-ga commented 2 years ago

I believe the average incidence angle is already saved to the interferogram geotiff metadata?

Oh, here I meant the average azimuth angle (azimuth direction of target to sensor).

My reasoning for the additional metadata here is just so non-experts can have something to refer to independently when investigating the velocity file and have specific infrastructure in question where they may already know the expected direction of motion and whether it may or may not be observed in their velocity file due to the sensing geometry.

adeane-ga commented 2 years ago

Closing this issue because it was fixed with this merge into develop: https://github.com/GeoscienceAustralia/PyRate/pull/369