计算图像的自信息和信息熵，需包含下述步骤： (1) 读入一幅图像 (2) 计算图中每个灰度级的频数、频率 (3) 从频率出发计算该图像的自信息、信息熵 (4) 可视化显示原图、灰度级频数、自信息和信息熵 计算两只骰子点数组合(随机事件)的熵，需包含下述步骤： (1) 用列表类型保存一只骰子(色子)的可能结果； (2) 表示并保存两只骰子的可能点数组合；提示(不强求)：可以使用itertools.product()函数产生两个列表的笛卡尔积。 (3) 将两只骰子的点数之和作为随机事件，例如两只骰子分别取点1和点2，与分别取点2和点1视为同一种组合，计算出每种组合的数量(频数)，存入一个列表结构。 (4) 由每种组合的频数，计算该组合的频率(概率)。 (5) 计算两只骰子点数和的概率分布的熵。

计算图像的自信息和信息熵：

import numpy as np
import matplotlib.pyplot as plt

# 读入一幅图像
img = plt.imread('image.png')

# 计算图中每个灰度级的频数、频率
hist, bins = np.histogram(img, bins=256, range=(0, 1))
freq = hist.astype(float) / np.sum(hist)

# 计算该图像的自信息、信息熵
log_freq = np.log2(freq)
entropy = -np.sum(freq * log_freq)
self_info = -log_freq

# 可视化显示原图、灰度级频数、自信息和信息熵
fig, axs = plt.subplots(2, 2, figsize=(10, 10), sharex=True)

axs[0, 0].imshow(img, cmap='gray')
axs[0, 0].set_title('Original Image')

axs[0, 1].bar(bins[:-1], hist, width=1)
axs[0, 1].set_title('Histogram')

axs[1, 0].bar(bins[:-1], self_info, width=1)
axs[1, 0].set_title('Self Information')

axs[1, 1].text(0.5, 0.5, 'Entropy: {:.3f}'.format(entropy),
               horizontalalignment='center', verticalalignment='center', fontsize=20)
axs[1, 1].axis('off')
axs[1, 1].set_title('Information Entropy')

plt.show()

计算两只骰子点数组合(随机事件)的熵：

import itertools

# 用列表类型保存一只骰子(色子)的可能结果
dice = [1, 2, 3, 4, 5, 6]

# 表示并保存两只骰子的可能点数组合
combinations = list(itertools.product(dice, dice))

# 将两只骰子的点数之和作为随机事件
sums = [sum(c) for c in combinations]

# 计算出每种组合的数量(频数)，存入一个列表结构
hist, bins = np.histogram(sums, bins=11, range=(2, 12))
freq = hist.astype(float) / np.sum(hist)

# 由每种组合的频数，计算该组合的频率(概率)
log_freq = np.log2(freq)
entropy = -np.sum(freq * log_freq)

# 计算两只骰子点数和的概率分布的熵
plt.bar(bins[:-1], freq, width=1)
plt.title('Probability Distribution')
plt.show()

print('Entropy: {:.3f}'.format(entropy))