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 <vtkImageCast.h>

#include <vtkSmartPointer.h>
#include <vtkMetaImageReader.h>
#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 <vtkInteractorStyleImage.h>
#include "itkImageToVTKImageFilter.h"
#include "itkVTKImageToImageFilter.h"
#include "itkMeanImageFilter.h"

#include <vtkSmartPointer.h>
#include <vtkGPUVolumeRayCastMapper.h>
#include <vtkColorTransferFunction.h>
#include <vtkPiecewiseFunction.h>
#include <vtkRenderer.h>
#include <vtkRenderWindow.h>
#include <vtkPlanes.h>
#include <vtkObjectFactory.h>
#include "vtkRendererCollection.h"
 
 #include <vtkPolyDataToImageStencil.h>
#include <vtkImageStencil.h>
#include <vtkPointData.h>


#include <vtkVersion.h>
#include <vtkSmartPointer.h>
#include <vtkPolyData.h>
#include <vtkImageData.h>
#include <vtkSphereSource.h>
#include <vtkMetaImageWriter.h>
#include <vtkPolyDataToImageStencil.h>
#include <vtkImageStencil.h>
#include <vtkPointData.h>
#include <vtkCutter.h>
#include <vtkPlane.h>
#include <vtkStripper.h>
#include <vtkLinearExtrusionFilter.h>
#include <vtkXMLPolyDataWriter.h>
#include <vtkImageData.h>
#include <vtkMetaImageReader.h>
#include <vtkSmartPointer.h>
#include <vtkInteractorStyleImage.h>
#include <vtkRenderer.h>
#include <vtkImageActor.h>
#include <vtkImageMapper3D.h>
#include <vtkRenderWindow.h>
#include <vtkRenderWindowInteractor.h>
#include <vtkPolyDataWriter.h>

#include <vtkAlgorithmOutput.h>

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

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

  InteractorStyle()
  {
    selectedMapper = vtkSmartPointer<vtkDataSetMapper>::New();
    selectedActor = vtkSmartPointer<vtkActor>::New();
    renderer = vtkSmartPointer<vtkRenderer>::New();
  }

  virtual void OnLeftButtonUp()//重写左键按下消息
  {
    // Forward events  
    vtkInteractorStyleRubberBandPick::OnLeftButtonUp();
    if (this->CurrentMode == VTKISRBP_SELECT)
    {
      vtkPlanes* frustum = static_cast<vtkAreaPicker*>(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<vtkPolyDataWriter> vtkWriter = vtkSmartPointer<vtkPolyDataWriter>::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<vtkRenderer> renderer;
  vtkSmartPointer<vtkPolyData> Data;
  vtkSmartPointer<vtkDataSetMapper> selectedMapper;
  vtkSmartPointer<vtkActor> selectedActor;
};
vtkStandardNewMacro(InteractorStyle);

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

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


  //进行精简reduce the number of triangles in a mesh
  vtkSmartPointer<vtkDecimatePro> deci = vtkSmartPointer<vtkDecimatePro>::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<vtkSmoothPolyDataFilter> smoother = vtkSmartPointer<vtkSmoothPolyDataFilter>::New();
  smoother->SetInputConnection(deci->GetOutputPort());
  smoother->SetNumberOfIterations(50);//Specify the number of iterations for Laplacian smoothing

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


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


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

  vtkSmartPointer<vtkRenderer> renderer =
    vtkSmartPointer<vtkRenderer>::New();
  vtkSmartPointer<vtkRenderer> renderer2 =
    vtkSmartPointer<vtkRenderer>::New();


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

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

  // Set the custom stype to use for interaction.
  vtkSmartPointer<InteractorStyle> style =
    vtkSmartPointer<InteractorStyle>::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;
}

注意：

记得添加

vtkSmartPointer<vtkAreaPicker> areaPicker = vtkSmartPointer<vtkAreaPicker>::New();
 
    renderWindowInteractor->SetPicker(areaPicker);