使用itkwidgets的view查看医学CT 3D对象时透明立方体问题排查

• 传输函数与CT灰度值范围不匹配：医学CT的灰度值为HU单位，常规范围在-1000（空气）到3000（骨组织）区间，itkwidgets默认的体渲染透明度传输函数是适配0-255或0-1的归一化图像的，大部分CT的HU值落在默认传输函数的全透明区间，仅保留了图像外框的立方体渲染结果。
• 空间间距参数未正确配置：CT数据通常为各向异性分辨率，z轴层厚远大于x/y轴的平面分辨率，如果没有显式传入spacing参数，工具默认按各向同性采样渲染，会导致体数据采样异常，无法识别有效组织灰度。
• 维度顺序不匹配：itkwidgets要求3D输入数组的维度顺序为(z, y, x)，如果传入的numpy数组维度顺序不对，也会触发渲染异常。

• 手动指定灰度渲染范围：先统计CT数据的有效灰度区间，可通过numpy的分位计算排除极值后，将vmin、vmax参数传入渲染函数，示例代码如下：

# 取1%和99%分位值作为渲染上下限，避免极值干扰
vmin = np.percentile(ct_array, 1)
vmax = np.percentile(ct_array, 99)
view(image=ct_array, vmin=vmin, vmax=vmax)

如果仍显示不清晰，可额外调整透明度映射参数，自定义不同灰度对应的透明度。

• 显式传入空间间距参数：根据CT元数据获取各轴的实际分辨率，按(z轴间距, y轴间距, x轴间距)的顺序传入spacing参数，示例如下：
  view(image=ct_array, vmin=-500, vmax=1500, spacing=(1.5, 0.5, 0.5))
• 调整数组维度顺序：如果确认维度顺序不匹配，用numpy的transpose方法转置数组到要求的(z, y, x)格式即可。

内容的提问来源于stack exchange，提问作者Vivek Sthanam