hist.py

hist

# get counts for raster using rio

from collections import Counter
import numpy as np
import rasterio as rio
from sys import argv

_, filepath = argv

with rio.open(filepath) as src:
    bands = src.read()
    height = len(bands[0])
    width = len(bands[0][0])
    print("height:", height, width)
    counts = Counter()
    for row in bands[0]:
      counts.update(row)
    print("counts:", counts.most_common(5))