Question

我想将我的dicom图像的窗口级别从肺窗改为胸窗。我知道窗口调平需要的值。但是如何在python中实现它？或者任何人都可以向我提供这个过程的详细描述将受到高度赞赏。

Answer 1

我已经在Python中实现了这一点。查看dicomparser库中dicompyler-core模块中的GetImage函数。

基本上它遵循kritzel_sw的建议。

Answer 2

以下开源代码可以实现骨窗口。

def get_pixels_hu(slices):
    image = numpy.stack([s.pixel_array for s in slices])
    image = image.astype(numpy.int16)
    image[image == -2000] = 0
    for slice_number in range(len(slices)):
        intercept = slices[slice_number].RescaleIntercept
        slope = slices[slice_number].RescaleSlope
        if slope != 1:
            image[slice_number] = slope * image[slice_number].astype(numpy.float64)
            image[slice_number] = image[slice_number].astype(numpy.int16)
        image[slice_number] += numpy.int16(intercept)

    return numpy.array(image, dtype=numpy.int16)

我添加了以下代码

image[slice_number] = image[slice_number]*3.5 + mean2*0.1

后 image [slice_number] + = numpy.int16（拦截）可以将骨窗改为脑组织窗。关键是设置参数3.5和0.1。我只是尝试了这两个适合脑组织窗口的参数。也许你可以调整胸部窗口。

python中的窗口调平

2 个答案: