利用vtkwidget 实现剪裁
想利用vtkBoxWidget和vtkExtractVOI实现剪裁。之前介绍过vtkBoxWidget，主要是说基本功能和接口，没有贴出完整实例。
Study-VTK：vtkWidget 分割/配准类之 正交六面体3D小部件 vtkBoxWidget
干脆直接介绍如何利用vtkwidget 实现剪裁。本文代码基本就是 cp 一遍vtk相关接口。
利用vtkwidget 实现剪裁
- 1 利用vtkwidget 剪裁步骤
- 2 利用vtkwidget 剪裁模型展示
- 3 利用vtkwidget 剪裁模型代码
- 4 利用vtkwidget 剪裁影像展示
- 5 利用vtkwidget 剪裁影像代码
1 利用vtkwidget 剪裁步骤
实现剪裁需要三步：
-
输入被剪裁模型；
被剪裁数据一般有：
影像（vti/vtr格式 dcm/图片/矩阵 等数据）、
模型（vtp/vtu格式 stl/obj/体网格 等数据）、
场景 -
剪切区域选择（交互）；
vtk提供大量现成的交互vtkWidget。可以自己开发交互，我只搞过一个，很麻烦。现成的虽然丑点但是是很多前辈积累的，比较实用。一般不建议自己开发。跟剪切相关的，遵循kiss原则，大家直接想到的无非就是：
平面上任意点两个点，模型分成两个部分
平面上任意点多个点，练成闭合样条曲线，分成曲线内、曲线外
空间放置一个球，分成球内求外
空间放置一个立方体，截取立方体内外 -
剪切算法实现。
模型的基本上是自己搞一个无限长的平面，利用在利用现成的与或非
影像的话vtkExtractVOI可以实现任意剪裁
2 利用vtkwidget 剪裁模型展示
临时搞一个demo，可以先看下效果。
被剪裁模型：
- 本地stl文件
剪切区域选择：
- 球、
- box、
- 平面任意两点、
- 平面任意多点
剪切算法：
- vtkClipPolyData实现两个模型取并、
- vtkClipPolyData实现两个模型取异、
- vtkCutter实现两个模型去交线、
- vtkBoxClipDataSet生成无限平面后去并
stl、box、取交集
stl、box、取外部
stl、box、取交线
stl、ball、取交集
stl、ball、取外部
stl、平面直线、取最大连通域（直线设置比较细，录屏压缩后看不到了）
stl、平面多点、取最大连通域（顺时针内、逆时针外）
3 利用vtkwidget 剪裁模型代码
vtkWidget 大部分使用方法基本一样，vtkBoxWidget举例
1 初始化
if (widget_type_ == BOX) {
if (clip_box_widget_ == nullptr) {
clip_box_widget_ = vtkSmartPointer<vtkBoxWidget>::New();
clip_box_widget_->GetFaceProperty()->SetColor(0.6, 0.6, 0.2);
clip_box_widget_->GetFaceProperty()->SetOpacity(0.25);
clip_box_widget_->SetInteractor(
vmtk_renderer_->GetRenderWindowInteractor());
}
}
1 设置rans
clip_box_widget_->GetTransform(transform_);
2 开启交互
if (widget_type_ == BOX) {
clip_box_widget_->SetInputData(surface_);
clip_box_widget_->PlaceWidget();
} else if (widget_type_ == SPHERE) {
} else if (widget_type_ == LINE) {
}
if (transform_ && widget_type_ == BOX) {
clip_box_widget_->SetTransform(transform_);
clip_box_widget_->On();
}
3 交互完成后获取数据
vtkSmartPointer<vtkClipPolyData> clipper = vtkSmartPointer<vtkClipPolyData>::New();
vtkSmartPointer<vtkPlanes> clip_planes_function;
vtkSmartPointer<vtkSphere> clip_sphere_function;
clipper->SetInputData(surface_);
clipper->GenerateClippedOutputOn();
clipper->SetInsideOut(inside_out_);
if (widget_type_ == BOX) {
clip_planes_function = vtkSmartPointer<vtkPlanes>::New();
clipper->SetClipFunction(clip_planes_function);
} else if (widget_type_ == SPHERE) {
clip_sphere_function = vtkSmartPointer<vtkSphere>::New();
clipper->SetClipFunction(clip_sphere_function);
}
vtkSmartPointer<vtkCutter> cutter = vtkSmartPointer<vtkCutter>::New();
cutter->SetInputData(surface_);
if (widget_type_ == BOX) {
cutter->SetCutFunction(clip_planes_function);
} else if (widget_type_ == SPHERE) {
cutter->SetCutFunction(clip_sphere_function);
}
clipped_surface_ = vtkSmartPointer<vtkPolyData>::New();
cut_lines_ = vtkSmartPointer<vtkPolyData>::New();
if (widget_type_ == BOX) {
clip_box_widget_->GetPlanes(clip_planes_function);
} else if (widget_type_ == SPHERE) {
clip_sphere_widget_->GetSphere(clip_sphere_function);
}
clipper->Update();
if (widget_type_ == BOX) {
clip_box_widget_->Off();
} else if (widget_type_ == SPHERE) {
clip_sphere_widget_->Off();
}
cutter->Update();
4 关闭交互
if (widget_type_ == BOX) {
clip_box_widget_->Off();
} else if (widget_type_ == SPHERE) {
clip_sphere_widget_->Off();
}
完整代码
我这里把vtkSeedWidget、vtkBoxWidget、vtkSphereWidget封装在一起。
#ifndef VMTKSURFACECLIPPER_H
#define VMTKSURFACECLIPPER_H
#include #include "vmtksurfaceclipper.h"
#include 4 利用vtkwidget 剪裁影像展示
影像剪裁一般在平面选择切割方向让后利用vtkExtractVOI重新切割影像，网友问vtkExtractVOI + vtkBoxWidget应该是在空间剪裁。方法跟上边模型切割一样，周末不上边了补上。
5 利用vtkwidget 剪裁影像代码
Study-VTK：三维影像实现任意方向、大小的切割