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

如果你想统计矢量数据中栅格像元的个数，可以使用以下步骤：

1. 首先，将矢量数据转化为栅格数据。这可以通过使用 GIS 软件（如 ArcGIS、QGIS 等）中的 "矢量转栅格" 工具来实现。在转换过程中，你需要指定栅格分辨率和像元大小。

2. 一旦你获得了栅格数据，你可以使用 Python 中的 GDAL 库来读取栅格数据。你可以使用以下代码来加载栅格数据：

   from osgeo import gdal

   # 打开栅格数据
   raster = gdal.Open('raster.tif')

3. 接下来，你可以使用 GDAL 库中的 `ReadAsArray()` 方法来读取栅格数据中每个像元的值。例如，以下代码将读取栅格数据中所有像元的值，并将它们存储在一个 NumPy 数组中：

   import numpy as np

   # 读取栅格数据中所有像元的值
   band = raster.GetRasterBand(1)
   arr = band.ReadAsArray()

   # 统计像元个数
   num_pixels = np.count_nonzero(arr)

在上面的代码中，`GetRasterBand()` 方法用于获取栅格数据的一个波段，`ReadAsArray()` 方法用于读取该波段的像元值。最后，`count_nonzero()` 方法用于计算非零像元的数量。

4. 最后，你可以将像元个数输出到控制台或将其保存到一个文件中，以便后续使用。例如，以下代码将像元个数输出到控制台：

   print('Number of pixels:', num_pixels)

希望这可以帮助你解决问题！