在python中获取栅格化边框线内栅格文件的像元索引

from osgeo import gdal

# 打开栅格文件
dataset = gdal.Open('path/to/raster.tif')

# 获取栅格文件的边界
transform = dataset.GetGeoTransform()
xmin = transform[0]
ymax = transform[3]
xmax = xmin + transform[1] * dataset.RasterXSize
ymin = ymax + transform[5] * dataset.RasterYSize

# 获取栅格文件的像元大小
pixel_width = transform[1]
pixel_height = transform[5]

# 计算栅格文件的行列数
num_rows = dataset.RasterYSize
num_cols = dataset.RasterXSize

# 计算栅格文件的像元索引
row_indices = range(int((ymax - ymin) / pixel_height), int((ymax - ymin) / pixel_height) + num_rows)
col_indices = range(int((xmin - xmax) / pixel_width), int((xmin - xmax) / pixel_width) + num_cols)

# 输出像元索引
for row in row_indices:
    for col in col_indices:
        print(f'Row: {row}, Col: {col}')

# 关闭栅格文件
dataset = None