VTK中Dicom影像显示的当前缩放倍率获取

Posted 2022-05-18 farmer-d

找了好久发现找不到类似的接口，只能自己研究了，主要还是要对VTK中的坐标系系统要了解。这里涉及到两个坐标系一个是世界坐标系（World），还有一个就是显示坐标系（Display）。世界坐标系为VTK中唯一的坐标系不会改变，显示坐标系则为我们眼睛看到的坐标系了。所以我们要将世界坐标系中的图像尺寸映射到显示坐标系，然后和原始尺寸进行对比，就可以得到当前的缩放倍率了。

 1 // vtkSmartPointer<vtkImageActor> dcmActor = vtkSmartPointer<vtkImageActor>::New();
 2 
 3 // 获取显示的边界
 4 double* pBounds = dcmActor->GetBounds();
 5 
 6 // 获取右下角坐标
 7 vtkSmartPointer<vtkCoordinate> coordMax = vtkSmartPointer<vtkCoordinate>::New();
 8 coordMax->SetCoordinateSystemToWorld();
 9 coordMax->SetValue(pBounds[1], pBounds[3]);
10 int* pDisp = coordMax->GetComputedDisplayValue(render);
11 
12 // 获取左上角坐标
13 vtkSmartPointer<vtkCoordinate> coordMin = vtkSmartPointer<vtkCoordinate>::New();
14 coordMin->SetCoordinateSystemToWorld();
15 coordMin->SetValue(pBounds[0], pBounds[2]);
16 int* pDisp1 = coordMin->GetComputedDisplayValue(render);
17 
18 // 计算倍率
19 int nWidth = pDisp[0] - pDisp1[0];
20 double scale = nWidth * 1.0 / 512;