Python计算CT数据灰度直方图

计算 CT 数据灰度直方图的方法与计算普通图像的灰度直方图类似，只是需要先将 CT 数据转换为灰度图像。以下是 Python 代码实现：

import pydicom
import numpy as np
import matplotlib.pyplot as plt

# 读取 CT 数据
dcm = pydicom.dcmread('ct.dcm')
img = dcm.pixel_array.astype(float)

# 根据窗宽窗位进行 CT 值的映射
ww = dcm.WindowWidth
wl = dcm.WindowCenter
img = np.clip((img - (wl - 0.5)) / (ww - 1.0) + 0.5, 0, 1)
img = (img * 255).astype(np.uint8)

# 统计灰度值的频数
hist, bins = np.histogram(img.ravel(), 256, [0, 256])

# 绘制直方图
plt.bar(bins[:-1], hist, width=1)
plt.xlim([0, 256])
plt.show()

其中，`pydicom.dcmread()` 函数用于读取 CT 数据，`np.clip()` 函数用于根据窗宽窗位进行 CT 值的映射，`np.histogram()` 函数用于统计灰度值的频数，`plt.bar()` 函数用于绘制直方图。