gdal
在讀取tif的時候，實際上行列讀取順序是和數組相反的，所以賦值的時候一定是行列互換，否則讀取的時候就會越界

###分塊讀取數據的生成器, gdal 行和列與數組相反！！！！！！！！
def readTifAsBlockGenerator(im,n_band,row,col,input_size,width,height):
for i in range(row):
for j in range(col):
if j != col - 1 and i != row - 1:

xoff = int(0.5 * j * input_size)yoff = int(0.5 * i * input_size)# print(str(im))im1 = mytest.readTifAsBlock(im, n_band=n_band, xoff=xoff, yoff=yoff,xsize=input_size,ysize=input_size)yield im1elif j != col - 1 and i == row - 1:##########gdal不支持負號的indexxoff = int(0.5 * j * input_size)yoff = height-input_sizeim1 = mytest.readTifAsBlock(im, n_band=n_band, xoff=xoff, yoff=yoff, xsize=input_size, ysize=input_size)yield im1# X_tmp[tmp_num] = x[-input_size:,int(0.5 * j * input_size):int(0.5 * (j + 2) * input_size), :]# Y_tmp[tmp_num] = y[-input_size:,int(0.5 * j * input_size):int(0.5 * (j + 2) * input_size), :]elif j == col - 1 and i != row - 1:##########gdal不支持負號的indexxoff = width-input_sizeyoff = int(0.5 * i * input_size)im1 = mytest.readTifAsBlock(im, n_band=n_band, xoff=xoff, yoff=yoff, xsize=input_size, ysize=input_size)yield im1# X_tmp[tmp_num] = x[int(0.5 * i * input_size):int(0.5 * (i + 2) * input_size),-input_size:, :]# Y_tmp[tmp_num] = y[int(0.5 * i * input_size):int(0.5 * (i + 2) * input_size),-input_size:, :]else:##########gdal不支持負號的indexxoff = width-input_sizeyoff = height - input_sizeim1 = mytest.readTifAsBlock(im, n_band=n_band, xoff=xoff, yoff=yoff, xsize=input_size, ysize=input_size)yield im1

例如，要讀取從圖1中所示的第6000行和第1400列開始的六列和三行，你可以執行以下操作：

data = band.ReadAsArray(1400, 6000, 6, 3)

總結

以上是生活随笔為你收集整理的gdal的知识的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。