Python GDAL：使用一系列像素值计算光栅图像的像素数

Question

我有一些光栅文件，我想计算每个文件中像素值落在某个特定范围之间的像素数。我尝试了以下代码，虽然它运行但产生了错误的像素数。

代码：

import os, sys
import numpy as np
import gdalnumeric
path = ('dir')
os.chdir(path)
for rasters in os.listdir(path):
    raster_file = gdalnumeric.LoadFile(rasters)
    pixel_count = (raster_file =>0 & <= 25).sum()
    print(pixel_count)

Answer 1

你的语法有点不对劲——我不认为它可以编译，但是有

((raster_file >= 0) & (raster_file <= 25)).sum()

你应该得到正确的结果