Error al utilizar GDAL para exportar un raster a GeoTIFF

Question

Estoy corrigiendo y recortando datos GeoTIFF DEM utilizando GDAL 2.4.0 para obtener el raster de la imagen. Después de hacer mis manipulaciones, tengo el raster que quiero (comprobado al trazar con plt.imshow()), pero al exportar/guardar a un GeoTIFF las últimas cien filas o más se cortan. Aquí está la sección de exportación de mi script, hágame saber si usted tiene alguna idea de por qué esto está sucediendo.

# Export and Save
print("Exporting file...")
crop_tiff = driver.Create(filename, ds_array.shape[1], ds_array.shape[0], 1, gdal.GDT_Float32) #Create new GeoTiff
crop_tiff.GetRasterBand(1).WriteArray(ds_array) #Write data
crop_tiff.GetRasterBand(1).SetNoDataValue((noData))
crop_tiff.SetGeoTransform(new_gt) 
crop_tiff.SetProjection(proj) 
crop_tiff.FlushCache() #Flush and save.
print("Saving complete.")

Mi trama (correcta):

El geoTIFF exportado con error:

Answer 1

¡Acabé descubriendo esto! Al parecer, el método FlushCache() de GDAL no siempre es suficiente por sí solo para guardar, así que todo lo que tenía que hacer era añadir

crop_tiff = None

al final para borrar correctamente la memoria, y todo se arregló. ¡Frustrante, pero sencillo!