vtk 鼠标截取数据 使用 vtkClipPolyData 裁剪删除选中的矩形区域

vtkVolumeMapper类中提供了两种裁剪技术，分别为Cropping和Clipping

按键盘 R 可以进行矩形区域裁剪，再按R返回正常状态

效果

：

#include "InteractorStyle.hpp"

#pragma once
#include"vtkContourFilter.h"
#include"vtkAreaPicker.h"

#include "vtkInteractorStyleRubberBandPick.h"
#include "vtkClipPolyData.h"
#include "vtkSmartPointer.h"
#include "vtkDataSetMapper.h"

#include 

#include 
#include 
#include "vtkRenderWindowInteractor.h"
#include "vtkRenderer.h"
#include "vtkRenderWindow.h"
#include "vtkMarchingCubes.h"
#include "vtkStripper.h"
#include "vtkActor.h"
#include "vtkPolyDataMapper.h"
#include "vtkSmoothPolyDataFilter.h"
#include "vtkPolyDataNormals.h"
#include "vtkDecimatePro.h"
#include "vtkProperty.h"
#include 
#include "itkImageToVTKImageFilter.h"
#include "itkVTKImageToImageFilter.h"
#include "itkMeanImageFilter.h"

#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include "vtkRendererCollection.h"
 
 #include 
#include 
#include 


#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 

#include 

#define VTKISRBP_ORIENT 0
#define VTKISRBP_SELECT 1
class InteractorStyle : public vtkInteractorStyleRubberBandPick //重载vtkInteractorStyleRubberBandPick

{
public:
	 static InteractorStyle* New();
	 vtkTypeMacro(InteractorStyle, vtkInteractorStyleRubberBandPick);

	InteractorStyle()
	{
		selectedMapper = vtkSmartPointer::New();
		selectedActor = vtkSmartPointer::New();
		renderer = vtkSmartPointer::New();
	}

	virtual void OnLeftButtonUp()//重写左键按下消息
	{
		// Forward events  
		vtkInteractorStyleRubberBandPick::OnLeftButtonUp();
		if (this->CurrentMode == VTKISRBP_SELECT)
		{
			vtkPlanes* frustum = static_cast(this->GetInteractor()->GetPicker())->GetFrustum();//获得鼠标框选矩形

			vtkClipPolyData *clipper = vtkClipPolyData::New();//裁剪polydata
			clipper->SetInputData(this->Data);
			clipper->SetClipFunction(frustum);//!!!!很重要的一步，添加自定义隐函数
			clipper->GenerateClipScalarsOn();
			clipper->GenerateClippedOutputOn();
		 
			clipper->SetValue(0.5);
			this->selectedMapper->SetInputConnection(clipper->GetOutputPort());
			this->selectedMapper->ScalarVisibilityOff();
			this->selectedActor->SetMapper(selectedMapper);
			this->selectedActor->GetProperty()->SetColor(1.0, 0.0, 0.0); //(R,G,B)  
			this->selectedActor->GetProperty()->SetRepresentationToWireframe();

		 //	vtkSmartPointer vtkWriter = vtkSmartPointer::New();
			//vtkWriter->SetInputConnection(clipper->GetOutputPort());
			//vtkWriter->SetFileName("D:\\test.vtk");
			//vtkWriter->Write(); 

			renderer->SetBackground(0.6, 0.8, 0.8); // Blue
			renderer->AddActor(selectedActor);
			renderer->SetViewport(0.5, 0, 1, 1);

			this->Interactor->GetRenderWindow()->AddRenderer(renderer);
			//	this->GetInteractor()->GetRenderWindow()->Render();
			this->Interactor->GetRenderWindow()->GetRenderers()->GetFirstRenderer()->Render();
		}
	}

  
	vtkSmartPointer renderer;
	vtkSmartPointer Data;
	vtkSmartPointer selectedMapper;
	vtkSmartPointer selectedActor;
};
vtkStandardNewMacro(InteractorStyle);

 

int showVTKDataTTT(vtkSmartPointer vtkData)
{
	vtkSmartPointer m_skinExtractor2 = vtkSmartPointer::New();//生成等值面/线

	m_skinExtractor2->SetInputConnection(vtkData->GetOutputPort());
	m_skinExtractor2->SetValue(0, 400);//调整阈值，400以上是骨组织


	//进行精简reduce the number of triangles in a mesh
	vtkSmartPointer deci = vtkSmartPointer::New();
	deci->SetInputConnection(m_skinExtractor2->GetOutputPort());
	deci->SetTargetReduction(0.3);//Specify the desired reduction in the total number of polygons (e.g., if TargetReduction is set to 0.9, this filter will try to reduce the data set to 10% of its original size).
	deci->PreserveTopologyOn();//Turn on/off whether to preserve the topology of the original mesh. If on, mesh splitting and hole elimination will not occur. This may limit the maximum reduction that may be achieved.

	//设置优化 用拉普拉斯平滑来调整点的位置
	vtkSmartPointer smoother = vtkSmartPointer::New();
	smoother->SetInputConnection(deci->GetOutputPort());
	smoother->SetNumberOfIterations(50);//Specify the number of iterations for Laplacian smoothing

	//计算多边形网格的法线
	vtkSmartPointer normals = vtkSmartPointer::New();
	normals->SetInputConnection(smoother->GetOutputPort());
	normals->FlipNormalsOn();


	vtkSmartPointer mapper =
		vtkSmartPointer::New();
	mapper->SetInputConnection(normals->GetOutputPort());


	vtkSmartPointer actor =
		vtkSmartPointer::New();
	actor->SetMapper(mapper);
	actor->GetProperty()->SetColor(1.0, 1.0, 0.0);
	actor->GetProperty()->SetRepresentationToWireframe();

	vtkSmartPointer renderer =
		vtkSmartPointer::New();
	vtkSmartPointer renderer2 =
		vtkSmartPointer::New();


	vtkSmartPointer renderWindow =
		vtkSmartPointer::New();
	renderWindow->AddRenderer(renderer);

	vtkSmartPointer areaPicker = vtkSmartPointer::New();
	vtkSmartPointer renderWindowInteractor =
		vtkSmartPointer::New();
	renderWindowInteractor->SetRenderWindow(renderWindow);
	renderWindowInteractor->SetPicker(areaPicker);
	renderWindowInteractor->Initialize();

	// Set the custom stype to use for interaction.
	vtkSmartPointer style =
		vtkSmartPointer::New();

	style->Data = normals->GetOutput();

	renderWindowInteractor->SetInteractorStyle(style);

	renderWindow->SetSize(640, 480);

	renderer->AddActor(actor);
	renderer->ResetCamera();
	renderer->SetViewport(0, 0, 0.5, 1);
	renderer->SetBackground(0.5, 0.5, 0.5); // Blue
	//renderer->SetBackground(.2, .3, .4);
	renderWindow->Render();
	renderWindowInteractor->Start();

	return EXIT_SUCCESS;
}