在Python中创建堆叠栅格图像的numpy数组

Question

I'm trying to read a stacked raster image as a numpy array. Actually I'd like to read all bands. But when I check array dimension, there is anything just 0. Please tell me why? and what should I do for this?

>>> import numpy, gdal, gdalconst, os, pandas
>>> os.chdir(r'N:\PhenologyMOD09\ImagesUsed\Condition\ndvi\Filter_tiff')
>>> gdal.AllRegister()
>>> inRas = gdal.Open("2000_allbands.tif", gdalconst.GA_ReadOnly)
>>> myarray = numpy.array([gdal.Open("2000_allbands.tif", gdalconst.GA_ReadOnly).ReadAsArray()
                for name, descr in inRas.GetSubDatasets()])

>>> myarray.shape
(0,)

Thanks in advance

Answer 1

Here is a simplified version of what you want to do:

from osgeo import gdal
inRas = gdal.Open("2000_allbands.tif")
myarray = inRas.ReadAsArray()
print(myarray.shape)  # (num_bands, y_size, x_size)