python 读取geotiff_GeoTiff图像在numpy数组中被读取为全零

你看图片的方式似乎还可以。但是，您不能仅通过打印来验证图像是否为全零。它只打印角点，因此您只能看到零，因为您的图像的所有角点都是零。在

可以使用numpy验证是否存在大于0的值。将您的代码更改为以下内容。请注意，现在它将打印39935，这是光栅的最大值(不是零)from osgeo import gdal

import numpy as np

gtif = gdal.Open( "test.tif" )

band = gtif.GetRasterBand(1)

bandArray = band.ReadAsArray()

print(bandArray)

print(np.max(bandArray))

为了获得特定像素的坐标，可以使用与获取光栅右下角坐标相同的方法。在

^{pr2}$

在上面的代码中，ulx和{}是左上角的坐标。这可以直接从光栅中提取。xres和{}是每个像素的大小。鉴于此，右下角的坐标lrx和{}可以根据前面提到的值来计算。在

现在为了找到图像中任何像素的坐标，您只需将gtif.RasterXSize和{}替换为需要坐标的像素项。在ulx, xres, xskew, uly, yskew, yres = gtif.GetGeoTransform()

xCoord = ulx + (pixelX * xres)

yCoord = uly + (pixelY * yres)

所以为了找到每个像素的坐标，可以执行以下操作ulx, xres, xskew, uly, yskew, yres = gtif.GetGeoTransform()

for x in range(0, gtif.RasterXSize):

for y in range(0, gtif.RasterYSize):

xCoord = ulx + (x * xres)

yCoord = uly + (y * yres)

# now do something with the coordinate