3D Slicer auto W/L实现

3D Slicer auto W/L实现

参考博客：
3Dslicer1:入门及基本控制
自动窗宽窗位的一些思路
How auto W/L is implemented in 3DSlicer?

python代码：
根据How auto W/L is implemented in 3DSlicer?实现

// auto W/L
def auto_wl(img, low=0.1, high=0.99):
    imhist, bins = np.histogram(img.flatten(), int(math.ceil(img.max()) - np.min(img[img > -math.inf]) + 1),
                                range=(np.min(img[img > -math.inf]), math.ceil(img.max())), density=True)
    pre_sum = np.zeros(shape=imhist.shape)
    temp_sum = 0
    for i in range(pre_sum.shape[0]):
        temp_sum += imhist[i]
        pre_sum[i] = temp_sum

    f = interp1d(pre_sum, bins[:-1], bounds_error=False, fill_value=0)
    low_value, high_value = f(low), f(high)

    assert high_value > 0

    l = math.ceil(low_value + high_value / 2)
    w = math.floor(high_value - low_value)

    return w, l

效果：
Fixed W/L

Auto W/L

可以看到对比度观感有明显提升