from PIL import Image
import numpy as np


def compute_histogram(
    image_path, bins: int, value_range=(0, 256), normalize: bool = True
):
    image = Image.open(image_path)
    image = np.array(image)

    histogram, _ = np.histogramdd(
        image.reshape(-1, 3), bins=bins, range=[value_range, value_range, value_range]
    ).astype(np.float64)

    if normalize:
        histogram = histogram / np.sum(histogram)

    return histogram