[发明专利]一种用于HCC多期相数据的窗宽窗位自适应调节方法

权利要求书

1.一种用于HCC感兴趣区域多期相数据的窗宽窗位自适应调节方法，其特征在于，包括：

步骤1：获取多期相医学图像数据；

步骤2：获取图像肿瘤中心坐标和肿瘤外接框；

步骤3：根据肿瘤中心坐标及肿瘤外接框，以肿瘤中心坐标为中心，以1.3倍于肿瘤外接框的扩展外接框为裁剪边缘，裁剪所有期相中的肿瘤，作为感兴趣肿瘤体积组；

步骤4：分别获取步骤3感兴趣肿瘤体积组中肿瘤体积的z轴向的中间切片，作为该肿瘤体积的感兴趣切片组；

步骤5：遍历步骤4中得到的所有感兴趣切片组，得到该感兴趣切片组所对应的最佳窗宽和最佳窗位；

步骤6：重复步骤4～5，直至得到所有期相中所有肿瘤体积的最佳窗宽和最佳窗位；

步骤7：利用步骤6中计算出的最佳窗宽和最佳窗位来裁剪所对应的肿瘤体积，并存入感兴趣肿瘤体积组，作为所对应期相的默认窗宽和窗位；

所述的步骤5的具体过程如下：

步骤5.1：遍历感兴趣切片组内所有感兴趣切片；

步骤5.1.1：计算感兴趣切片的HU最大值Gray_max和HU最小值Gray_min；

步骤5.1.2：固定窗宽为约定窗宽Gray_fix-ww，将移动窗位Gray_move-wl从HU最小值Gray_min到HU最大值Gray_max进行遍历；

步骤5.1.2.1：每次对窗位进行移动时，计算约定窗宽Gray_fix-ww和该时刻的移动窗位Gray_move-wl下图像的HU值直方图；所述的HU值直方图去除掉最小HU值和最大HU值所对应的统计值；

步骤5.1.2.2：使用直方图双峰法计算并记录步骤5.1.2.1所述HU值直方图中的双峰差值；

步骤5.1.2.3：对窗位遍历时步长设置，具体过程如下：

当步骤5.1.2.2所述双峰差值小于约定的经验阈值T1时，移动窗位Gray_move-wl遍历HU的步长为100，当所述双峰差值大于约定的经验阈值T1时，步长减半，最小步长为10；当所述双峰差值连续增大3次则恢复步长为100，至HU最大值Gray_max停止；

步骤5.1.2.4：统计所有双峰差值，将双峰差值最大的移动窗位Gray_move-wl作为最佳窗位Gray_best-wl；

步骤5.1.3：计算步骤5.1.2.4所述最佳窗位Gray_best-wl下的最大窗宽，使用该最大窗宽作为窗宽最大值Gray_ww-max，使用10作为窗宽最小值Gray_ww-min；

步骤5.1.4：将最佳窗位Gray_best-wl作为固定窗位，将移动窗宽Gray_move-ww从窗宽最小值Gray_ww-min到窗宽最大值Gray_ww-max进行遍历；

步骤5.1.4.1：每次对窗宽进行移动时，计算最佳窗位Gray_best-wl和该时刻的移动窗宽Gray_move-ww下图像的HU值直方图；所述的HU值直方图去掉灰度为0和灰度为255所对应的统计值；

步骤5.1.4.2：使用直方图双峰法计算并记录步骤5.1.4.1所述HU值直方图中的双峰差值；

步骤5.1.4.3：对窗宽遍历时步长设置，具体过程如下：

当步骤5.1.4.2所述双峰差值小于约定的经验阈值T2时，移动窗宽Gray_move-ww遍历HU的步长为50，当所述双峰差值大于约定的经验阈值T2时，步长减半，最小步长为10；当所述双峰差值连续增大3次则恢复步长为50，至窗宽最大值停止；

步骤5.1.4.4：统计所有双峰差值，将双峰差值最大的移动窗宽Gray_move-ww作为最佳窗宽Gray_best-ww；

步骤5.2：使用最佳窗宽Gray_best-ww和最佳窗位Gray_best-wl裁剪步骤3所述肿瘤体积；得到最佳窗宽和最佳窗位。

2.根据权利要求1所述的一种用于HCC感兴趣区域多期相数据的窗宽窗位自适应调节方法，其特征在于，所述直方图双峰法计算的计算方式如下：

将HU值直方图中的最高点作为第一峰值GH_first＝hist(r)，对应的灰度值为r；

第二峰的取值公式为此时的x即为第二峰的灰度值GH_second；

双峰差值为GH_first-GH_second。