python统计矢量内栅格像元个数

在Python中实现矢量内栅格像元个数的计算可以通过以下步骤：

1.导入需要的库

import geopandas as gpd
import rasterio

2.读取矢量和栅格数据

vector_data = gpd.read_file('vector.shp')
raster_data = rasterio.open('raster.tif')

3.获取栅格数据的分辨率和范围

width = raster_data.meta['width']
height = raster_data.meta['height']
xmin, ymin, xmax, ymax = raster_data.bounds
resolution = raster_data.meta['transform'][0]

4.创建一个像元计数器

pixel_count = 0

5.循环遍历矢量数据中的每个要素，将其转换为栅格坐标系下的范围，然后统计在该范围内的像素个数

for index, row in vector_data.iterrows():
    xmin_r, ymin_r, xmax_r, ymax_r = raster_data.index(row['geometry'].bounds)
    col_start = int((xmin_r - xmin) / resolution)
    col_end = int((xmax_r - xmin) / resolution) + 1
    row_start = int((ymax - ymax_r) / resolution)
    row_end = int((ymax - ymin_r) / resolution) + 1
    window = ((row_start, row_end), (col_start, col_end))
    data = raster_data.read(1, window=window)
    pixel_count += (data > 0).sum()

6.输出计算的总像元数

print(pixel_count)

这样就完成了Python统计矢量内栅格像元个数的过程。