dicom 图像感兴趣区域内对比度计算 代码

以下是Python代码，用于计算DICOM图像感兴趣区域内的对比度：

import pydicom
import numpy as np

# 加载DICOM图像
ds = pydicom.dcmread('example.dcm')

# 获取感兴趣区域（ROI）的坐标
x1, y1 = 100, 100 # 左上角坐标
x2, y2 = 200, 200 # 右下角坐标

# 获取ROI内的像素值
roi = ds.pixel_array[x1:x2, y1:y2]

# 计算ROI内的平均像素值
mean_roi = np.mean(roi)

# 计算ROI内像素值的标准差
std_roi = np.std(roi)

# 计算对比度
contrast = std_roi/mean_roi

print('对比度为：', contrast)

请注意，此代码仅用于演示目的，实际使用时需要根据DICOM图像的特定格式进行修改。此外，还需要注意DICOM图像中的灰度值范围和颜色映射，以正确计算像素值和对比度。