GDAL create Grid trae una Grid inválida en la primera ejecución

Question

Estoy jugando con gdal y numpy para crear una cuadrícula aleatoria en QGIS. El siguiente fragmento de código crea una cuadrícula e intenta cargarla en QGIS. Pero la capa resultante no es válida. Si ejecuto el script por segunda vez, justo después de la primera ejecución, la capa se crea y se carga sin ningún problema... ¿me falta algo en mi sniplet de código?

from osgeo import gdal
import numpy as np

fn = "d:/newgrid.tif"
rasterband = np.random.random((10,10))
driver = gdal.GetDriverByName("GTiff")
ds = driver.Create(fn, xsize=10, ysize=10, bands=1, eType=gdal.GDT_Float32)
ds.GetRasterBand(1).WriteArray(rasterband)
rlayer = iface.addRasterLayer(fn)

El mensaje de error en la primera ejecución es: Capa inválida : proveedor de GDAL No se puede abrir el conjunto de datos GDAL d:/newgrid.tif: `d:/newgrid.tif' no se reconoce como un formato de archivo compatible. El proveedor de la capa raster no es válido (proveedor: gdal, URI: d:/newgrid.tif

Answer 1

En comparación con https://pcjericks.github.io/py-gdalogr-cookbook/raster_layers.html#create-raster-from-array Falta FlushCache().