Geographic Information Systems Stack Exchange is a question and answer site for cartographers, geographers and GIS professionals. It only takes a minute to sign up.
Sign up to join this community
When you download a large image from Google Earth Engine, it is in multiple small images. Each sub-image name contains the x_min and y_min coordinates of that image.
It is described here - https://developers.google.com/earth-engine/image_upload#tfrecord
How to join all these .tif files into one .tif file in right order by using a Python script.
Merging tiles can be accomplished using:
It is not necessary to order the tiles in a specific way... they are non-overlapping.
If you want to order the tif files by coordinate pair this will do it.
from osgeo import gdal
tiffiles=["basefilename-99-100.tif","basefilename-98.7-200.tif","basefilename-101-300.tif","basefilename-100.25-400.tif","basefilename-99.925-500.tif","basefilename-97.552-100.tif"]
coordarray=[]
tifdict={}
for tf in tiffiles:
ds = gdal.Open(tf)
w = ds.RasterXSize
h = ds.RasterYSize
gtrans = ds.GetGeoTransform()
minx = gtrans[0]
miny = gtrans[3] + w*gtrans[4] + h*gtrans[5]
# splittifname=tf.split("-")
coords=(minx,miny)
coordarray.append(coords)
tifdict[coords]=tf
coordarray=sorted(coordarray,key=lambda k:(k[0],k[1]))
for coord in coordarray:
### DO STUFF TO COMBINE THE RASTER TO A STACK
print(tifdict[coord])
I don't think I explain myself clearly, this will put your tif in order by coordinate from min x min y of the set, then you can merge or place the file names in order on .txt file and batch process. I change the code to use the coord value in the tif file instead of value in file name but I have not had a chance to test you may need to declare the spatial reference somewhere.
