3D空间坐标转换
在冠状图、矢状面、横截面等创建的角度组件的三个端点坐标,不能直接用在3D视图中。这是因为2D切片的坐标是基于像素的,而3D空间的坐标可能是基于实际物理尺寸的。
解决方案是使用2D点的坐标、切片的物理位置和方向以及可能的其他参数(例如切片厚度或间距),来计算空间坐标。将计算好的空间坐标p1、p2、p3传入下面的函数中,就可以在3D视图中画出角度组件。文章来源:https://www.toymoban.com/news/detail-681799.html
例程
get3DViewAngle(double* p1, double* p2, double* mid, int* indexVal) {
double modifiedPoints[3][3];
int currentAngleIndex = *indexVal;
auto iter = threeDActors.find(currentAngleIndex);
if (iter != threeDActors.end()) {
for (auto act : iter->second) {
threeDRenderer->RemoveActor(act);
}
threeDActors.erase(iter);
}
double* pts[] = {p1, p2, mid};
for(int i = 0; i < 3; i++) {
modifiedPoints[i][0] = pts[i][0] * xRange[1] / imageDims[0];
modifiedPoints[i][1] = pts[i][1] * yRange[1] / imageDims[1];
modifiedPoints[i][2] = pts[i][2] * zRange[1] / imageDims[2];
}
vtkSmartPointer<vtkSphereSource> src1 = vtkSmartPointer<vtkSphereSource>::New();
src1->SetCenter(modifiedPoints[0]);
src1->SetRadius(2);
vtkSmartPointer<vtkPolyDataMapper> map1 = vtkSmartPointer<vtkPolyDataMapper>::New();
map1->SetInputConnection(src1->GetOutputPort());
vtkSmartPointer<vtkActor> act1 = vtkSmartPointer<vtkActor>::New();
act1->SetMapper(map1);
vtkSmartPointer<vtkSphereSource> src2 = vtkSmartPointer<vtkSphereSource>::New();
src2->SetCenter(modifiedPoints[1]);
src2->SetRadius(2);
vtkSmartPointer<vtkPolyDataMapper> map2 = vtkSmartPointer<vtkPolyDataMapper>::New();
map2->SetInputConnection(src2->GetOutputPort());
vtkSmartPointer<vtkActor> act2 = vtkSmartPointer<vtkActor>::New();
act2->SetMapper(map2);
vtkSmartPointer<vtkSphereSource> src3 = vtkSmartPointer<vtkSphereSource>::New();
src3->SetCenter(modifiedPoints[2]);
src3->SetRadius(2);
vtkSmartPointer<vtkPolyDataMapper> map3 = vtkSmartPointer<vtkPolyDataMapper>::New();
map3->SetInputConnection(src3->GetOutputPort());
vtkSmartPointer<vtkActor> act3 = vtkSmartPointer<vtkActor>::New();
act3->SetMapper(map3);
vtkSmartPointer<vtkLineSource> lnSrc1 = vtkSmartPointer<vtkLineSource>::New();
lnSrc1->SetPoint1(modifiedPoints[0]);
lnSrc1->SetPoint2(modifiedPoints[2]);
vtkSmartPointer<vtkPolyDataMapper> lnMap1 = vtkSmartPointer<vtkPolyDataMapper>::New();
lnMap1->SetInputConnection(lnSrc1->GetOutputPort());
vtkSmartPointer<vtkActor> lnAct1 = vtkSmartPointer<vtkActor>::New();
lnAct1->SetMapper(lnMap1);
vtkSmartPointer<vtkLineSource> lnSrc2 = vtkSmartPointer<vtkLineSource>::New();
lnSrc2->SetPoint1(modifiedPoints[1]);
lnSrc2->SetPoint2(modifiedPoints[2]);
vtkSmartPointer<vtkPolyDataMapper> lnMap2 = vtkSmartPointer<vtkPolyDataMapper>::New();
lnMap2->SetInputConnection(lnSrc2->GetOutputPort());
vtkSmartPointer<vtkActor> lnAct2 = vtkSmartPointer<vtkActor>::New();
lnAct2->SetMapper(lnMap2);
threeDRenderer->AddActor(act1);
threeDRenderer->AddActor(act2);
threeDRenderer->AddActor(act3);
threeDRenderer->AddActor(lnAct1);
threeDRenderer->AddActor(lnAct2);
std::vector<vtkSmartPointer<vtkActor>> acts = {act1, act2, act3, lnAct1, lnAct2};
threeDActors[currentAngleIndex] = acts;
for (auto act : acts) {
threeDRenderer->AddActor(act);
}
renderWidget[3]->interactor()->Render();
}
}
心得
- 在创建角度组件后,不能直接读取组件的端点坐标来用于计算空间坐标,因为此时的端点坐标还是空。需要等到鼠标左键点击三次后,才会得到完整的端点坐标。这个时候可以写在角度组件的回调函数中,然后使用成员变量来获取创建好的端点坐标。
- 在角度组件的回调函数中,可以使用信号槽机制来保证拖动组件后能够实时更新角度组件位置。
- 如果是创建线条组件,需要手动触发一次回调函数,因为线条组件的端点位置是可以直接读取的。
lineCallback->Execute(lineWidget, vtkCommand::EndInteractionEvent, nullptr); //手动触发回调
4.如果想要保证实时删除2D和3D视图中的角度组件,可以使用map,只需要记住创建的索引就能直接删除组件。这样可以与页面上的组件数量和组件索引独立开来.文章来源地址https://www.toymoban.com/news/detail-681799.html
std::map<int,vtkSmartPointer<vtkAngleWidget>> ngleWidgets;
到了这里,关于VTK——angleWidget的3D转换的文章就介绍完了。如果您还想了解更多内容,请在右上角搜索TOY模板网以前的文章或继续浏览下面的相关文章,希望大家以后多多支持TOY模板网!