Useful for unit testing
Write a Numpy array to GeoTIFFf file.
from osgeo import gdal, osr
from osgeo import gdalnumeric
import numpy as np
# get the osr spatial reference
srs = osr.SpatialReference()
# we want meters
srs.ImportFromEPSG(3857)
# set pixel dimensions
nx = 32
ny = 32
# we create a geotiff with 1 band
dst_ds = gdal.GetDriverByName('GTiff').Create('/path/to/dummy.tif', nx, ny, 1, gdal.GDT_Float32)
# create a random grid of 32 by 32 pixels
arr = np.random.rand(32, 32)
# write the array to the geotiff
dst_ds.GetRasterBand(1).WriteArray(arr)
# optionally set some nodata
dst_ds.GetRasterBand(1).SetNoDataValue(-9999)
# set the projection
dst_ds.SetProjection(srs.ExportToWkt())
#/* top left x */
#/* w-e pixel resolution */
#/* 0 */
# /* top left y */
# /* 0 */
# /* n-s pixel resolution (negative value) */
dst_ds.SetGeoTransform([32, 1, 0, 32, 0, -1])
# FlushChache to write
dst_ds.FlushCache()
Thank you so much for your huge amounts of work. These days, I really study many things form your materials (especially Google Earth Engine). I’m expecting for your update! 🙂 Again, thank you for your great works !!
LikeLike