Python如何读取影像高程并划分800*800格网

Python可以使用GDAL库来读取影像高程数据，并使用numpy库来进行数据处理和划分格网。以下是一个简单的代码示例：

import gdal
import numpy as np

# 打开影像高程数据
ds = gdal.Open('elevation.tif')

# 读取数据
data = ds.ReadAsArray()

# 将数据划分为800*800的格网
grid_size = 800
rows, cols = data.shape
grid_rows = rows // grid_size
grid_cols = cols // grid_size
grid_data = np.zeros((grid_rows, grid_cols))

for i in range(grid_rows):
    for j in range(grid_cols):
        grid_data[i, j] = np.mean(data[i*grid_size:(i+1)*grid_size, j*grid_size:(j+1)*grid_size])

# 输出结果
print(grid_data)

这段代码会将影像高程数据划分为800*800的格网，并计算每个格网的平均高程值。最终输出的结果是一个800*800的矩阵，每个元素代表一个格网的平均高程值。