I would like to use
scipy.interpolate.griddata to interpolate my Sentinel-1 products (see two sample images below) taken from different orbits resulting in having the uniform datapoint coordinates.
I would like to use a similar approach as describe in Smoothing/interpolating raster in Python using GDAL?
However, I am not able to figure it how to get points at which to interpolate data from a
geotiff. So far, I have this for values and for basic geographic description:
from osgeo import gdal # Loading terrain corrected image master = gdal.Open("S1A_TC.tif") # Data values of master image master_array = master.ReadAsArray() # Getting geographic data (e.g. pixelSizeX) gt = master.GetGeoTransform()
For getting array of data point coordinates or points at which to interpolate data, I found this topic: https://stackoverflow.com/questions/43370069/converting-xy-co-ords-of-geotiff-to-numpy-array-positions-in-python
The GDAL link is dead and I am not able to found something useful.
How can I read the data point coordinates as