技术标签: VTK笔记-裁剪分割
VTK加载二维数据和三维数据,我们通常观察到的是三视图图像,或者是vtkImageReslice获取到的转换矩阵对应的观察点三视图或者是斜平面图像,想要看到图像中间部分的渲染结果就要借助“手术刀”工具对图像数据进行分割裁剪,去掉不重要数据,将重要的脏器或者病灶暴露出来;
与抽取轮廓类似,裁剪也是指获取数据集中的部分数据。不过裁剪是从一个完全不同的角度来处理的。裁剪可以控制数据获取的区域,而抽取轮廓做不到这一点。
裁剪操作需要的一把“刀”——隐函数,用来完成切割的操作,并且要指定裁剪的位置——裁剪值。给定这个两个参数后,就确定了切面——隐函数对应的的函数曲面。裁剪的处理过程是首先将数据集中的点坐标带入隐函数中进行计算,然后将计算出的函数值等于裁剪值得部分(包括属性数据)取出来,即已知函数值,求自变量。可以看出裁剪是基于坐标的操作。注意:数据被“刀”分成了两部分:裁下的部分和裁掉的部分。默认只产生裁下的部分这一个输出。
vtkClipPolyData使用用户指定的隐式函数(user-specified implicit function)或输入标量数据(input scalar data)剪裁多边形数据;
vtkClipPolyData是一种Filter,接受vtkImplicitFunction派生出的任何子类或者标量数据集,vtkClipPolyData实现的剪切功能实际上“剪切”数据集的单元格,返回指定隐式函数(或大于标量值)内的所有内容,包括单元格的“片段”(将此与vtkExtractGeometry进行比较,该几何图形将拉出整个未切割单元格。)此筛选器的输出为多边形数据(vtkPolyData)。
要使用此Filter,必须决定是否使用隐式函数剪裁,还是使用输入标量数据。
如果要使用隐式函数进行剪裁,必须:
1)定义隐式函数
2)使用SetClipFunction方法设置它
3)如果未指定ClipFunction,则使用方法GenerateClipScalarsOn开启,或者是GenerateClipScalars关闭(默认值),则将使用输入的标量数据来剪裁多数据。
可以指定标量值,该值用于决定隐式函数的内部和外部是什么。
可以通过设置不同输出,选择那部分数据是被裁剪下的(切割算法通过计算隐式函数值或使用数据集中每个点的输入标量数据来进行。这与用于确定内部/外部的标量值进行比较。),对外显示出的效果就截然不同;
vtkClipPolyData可以配置为计算第二个输出。第二个输出是被剪裁的多边形数据。如果要访问此输出数据,需要将布尔值GenerateClippedData设置为on。
注意:为了在多边形数据中剪切所有类型的单元格,vtkClipPolyData将一些单元格三角化,然后剪切结果的简化(即点、线和三角形)。这意味着,结果输出可能包含不同于输入数据的单元格类型。在之前VTK版本中,需要手动添加vtkTriangleFilter来忽略多边形始末端重复的点,并对其三角化;
SetClipFunction方法接收vtkImplicitFunction类型的派生子类指针;如果未定义隐式函数,则将使用输入的标量数据进行剪裁。
vtkImplicitFunction* ClipFunction;
virtual void SetClipFunction(vtkImplicitFunction*);
vtkGetObjectMacro(ClipFunction, vtkImplicitFunction);
隐式函数的剪裁值(如果使用隐式函数剪裁时,该值为剪裁值)或标量值(如果使用标量剪裁,该值为标量值)。默认值为0.0。
double Value;
vtkSetMacro(Value, double);
vtkGetMacro(Value, double);
该值被设置为禁用时,如果隐式函数的值大于值ivar,则会在隐式函数中考虑顶点。启用“内部输出”时,如果隐式函数的隐式函数值小于或等于值ivar,则会在隐式函数中考虑顶点。默认情况下,InsideOut处于关闭状态。
vtkTypeBool InsideOut;
vtkSetMacro(InsideOut, vtkTypeBool);
vtkGetMacro(InsideOut, vtkTypeBool);
vtkBooleanMacro(InsideOut, vtkTypeBool);
如果启用此标志,则输出标量值将从隐式函数值中插值,而不是从输入标量数据中插值。
如果启用此标志但不提供隐式函数,则会报告错误。默认情况下,GenerateClipScalars处于禁用状态。
vtkTypeBool GenerateClipScalars;
vtkSetMacro(GenerateClipScalars, vtkTypeBool);
vtkGetMacro(GenerateClipScalars, vtkTypeBool);
vtkBooleanMacro(GenerateClipScalars, vtkTypeBool);
GenerateClippedOutput用于控制是否生成第二个输出。第二个输出包含被剪裁掉的多边形数据。默认情况下,GenerateClippedOutput处于禁用状态。
GetClippedOutputPort和GetClippedOutput用来输出第二个输出,是被剪裁掉的多边形数据。
GetOutput方法和GetOutputPort方法用来输出裁剪后剩余的多边形数据,是第一个输出。
vtkTypeBool GenerateClippedOutput;
vtkSetMacro(GenerateClippedOutput, vtkTypeBool);
vtkGetMacro(GenerateClippedOutput, vtkTypeBool);
vtkBooleanMacro(GenerateClippedOutput, vtkTypeBool);
vtkPolyData* GetClippedOutput();
vtkAlgorithmOutput* GetClippedOutputPort() {
return this->GetOutputPort(1); }
OutputPointsPrecision记录了对外输出数据精度类型,对应vtkAlgorithm::DesiredOutputPrecision中的SINGLE_PRECISION、
DOUBLE_PRECISION以及DEFAULT_PRECISION。
int OutputPointsPrecision;
vtkSetMacro(OutputPointsPrecision, int);
vtkGetMacro(OutputPointsPrecision, int);
场景中有三个Actor,分别是裁剪后牛剩余部分,裁剪掉的部分以及裁剪隐式函数的刨面部分;
这里使用的是一个裁剪平面,是vtkPlane,vtkPlane继承自vtkImplicitFunction,定义了一个平面;
使用vtkClipPolyData类的SetClipFunction方法设定裁剪平面。
#pragma once
#include "vtk_include.h"
#include <vtkActor.h>
#include <vtkCamera.h>
#include <vtkClipPolyData.h>
#include <vtkDataSetMapper.h>
#include <vtkFeatureEdges.h>
#include <vtkNamedColors.h>
#include <vtkNew.h>
#include <vtkPlane.h>
#include <vtkPolyData.h>
#include <vtkPolyDataMapper.h>
#include <vtkProperty.h>
#include <vtkRenderWindow.h>
#include <vtkRenderWindowInteractor.h>
#include <vtkRenderer.h>
#include <vtkSmartPointer.h>
#include <vtkSphereSource.h>
#include <vtkStripper.h>
#include <vtkXMLPolyDataReader.h>
// Readers
#include <vtkBYUReader.h>
#include <vtkOBJReader.h>
#include <vtkPLYReader.h>
#include <vtkPolyDataReader.h>
#include <vtkSTLReader.h>
#include <vtkXMLPolyDataReader.h>
namespace {
vtkSmartPointer<vtkPolyData> ReadPolyData(std::string const& fileName)
{
vtkSmartPointer<vtkPolyData> polyData;
std::string extension = "";
if (fileName.find_last_of(".") != std::string::npos){
extension = fileName.substr(fileName.find_last_of("."));
}
// Make the extension lowercase
std::transform(extension.begin(), extension.end(), extension.begin(),::tolower);
if (extension == ".ply"){
vtkNew<vtkPLYReader> reader;
reader->SetFileName(fileName.c_str());
reader->Update();
polyData = reader->GetOutput();
}
else if (extension == ".vtp"){
vtkNew<vtkXMLPolyDataReader> reader;
reader->SetFileName(fileName.c_str());
reader->Update();
polyData = reader->GetOutput();
}
else if (extension == ".obj"){
vtkNew<vtkOBJReader> reader;
reader->SetFileName(fileName.c_str());
reader->Update();
polyData = reader->GetOutput();
}
else if (extension == ".stl"){
vtkNew<vtkSTLReader> reader;
reader->SetFileName(fileName.c_str());
reader->Update();
polyData = reader->GetOutput();
}
else if (extension == ".vtk"){
vtkNew<vtkPolyDataReader> reader;
reader->SetFileName(fileName.c_str());
reader->Update();
polyData = reader->GetOutput();
}
else if (extension == ".g"){
vtkNew<vtkBYUReader> reader;
reader->SetGeometryFileName(fileName.c_str());
reader->Update();
polyData = reader->GetOutput();
}
else{
// Return a polydata sphere if the extension is unknown.
vtkNew<vtkSphereSource> source;
source->SetThetaResolution(20);
source->SetPhiResolution(11);
source->Update();
polyData = source->GetOutput();
}
return polyData;
}
} // namespace
class Test_Clip_cow {
public:
static void Test() {
vtkNew<vtkNamedColors> colors;
auto backgroundColor = colors->GetColor3d("steel_blue");
auto boundaryColor = colors->GetColor3d("banana");
auto clipColor = colors->GetColor3d("tomato");
// PolyData to process
auto polyData = ReadPolyData("G:\\Data\\cow.g");
vtkNew<vtkPlane> plane;
plane->SetOrigin(polyData->GetCenter());
plane->SetNormal(1.0, -1.0, -1.0);
vtkNew<vtkClipPolyData> clipper;
clipper->SetInputData(polyData);
clipper->SetClipFunction(plane);
clipper->GenerateClipScalarsOn();
clipper->GenerateClippedOutputOn();
clipper->SetValue(0.5);
clipper->Update();
polyData = clipper->GetOutput();
vtkNew<vtkDataSetMapper> clipMapper;
clipMapper->SetInputData(polyData);
clipMapper->ScalarVisibilityOff();
vtkNew<vtkActor> clipActor;
clipActor->SetMapper(clipMapper);
clipActor->GetProperty()->SetDiffuseColor(clipColor.GetData());
clipActor->GetProperty()->SetInterpolationToFlat();
clipActor->GetProperty()->EdgeVisibilityOn();
// Now extract feature edges
vtkNew<vtkFeatureEdges> boundaryEdges;
boundaryEdges->SetInputData(polyData);
boundaryEdges->BoundaryEdgesOn();
boundaryEdges->FeatureEdgesOff();
boundaryEdges->NonManifoldEdgesOff();
boundaryEdges->ManifoldEdgesOff();
vtkNew<vtkStripper> boundaryStrips;
boundaryStrips->SetInputConnection(boundaryEdges->GetOutputPort());
boundaryStrips->Update();
// Change the polylines into polygons
vtkNew<vtkPolyData> boundaryPoly;
boundaryPoly->SetPoints(boundaryStrips->GetOutput()->GetPoints());
boundaryPoly->SetPolys(boundaryStrips->GetOutput()->GetLines());
vtkNew<vtkPolyDataMapper> boundaryMapper;
boundaryMapper->SetInputData(boundaryPoly);
vtkNew<vtkActor> boundaryActor;
boundaryActor->SetMapper(boundaryMapper);
boundaryActor->GetProperty()->SetDiffuseColor(boundaryColor.GetData());
vtkNew<vtkPolyDataMapper> restMapper;
restMapper->SetInputConnection(clipper->GetClippedOutputPort());
restMapper->ScalarVisibilityOff();
vtkNew<vtkActor> restActor;
restActor->SetMapper(restMapper);
restActor->GetProperty()->SetRepresentationToWireframe();
// Create graphics stuff
vtkNew<vtkRenderer> renderer;
renderer->SetBackground(backgroundColor.GetData());
renderer->UseHiddenLineRemovalOn();
vtkNew<vtkRenderWindow> renderWindow;
renderWindow->AddRenderer(renderer);
renderWindow->SetSize(640, 480);
vtkNew<vtkRenderWindowInteractor> interactor;
interactor->SetRenderWindow(renderWindow);
// Add the actors to the renderer, set the background and size
renderer->AddActor(clipActor);
renderer->AddActor(boundaryActor);
renderer->AddActor(restActor);
// Generate an interesting view
renderer->ResetCamera();
renderer->GetActiveCamera()->Azimuth(30);
renderer->GetActiveCamera()->Elevation(30);
renderer->GetActiveCamera()->Dolly(1.2);
renderer->ResetCameraClippingRange();
renderWindow->Render();
renderWindow->SetWindowName("CapClip");
renderWindow->Render();
interactor->Start();
}
};
文章浏览阅读3.4k次,点赞8次,收藏42次。一、什么是内部类?or 内部类的概念内部类是定义在另一个类中的类;下面类TestB是类TestA的内部类。即内部类对象引用了实例化该内部对象的外围类对象。public class TestA{ class TestB {}}二、 为什么需要内部类?or 内部类有什么作用?1、 内部类方法可以访问该类定义所在的作用域中的数据,包括私有数据。2、内部类可以对同一个包中的其他类隐藏起来。3、 当想要定义一个回调函数且不想编写大量代码时,使用匿名内部类比较便捷。三、 内部类的分类成员内部_成员内部类和局部内部类的区别
文章浏览阅读118次。分布式系统要求拆分分布式思想的实质搭配要求分布式系统要求按照某些特定的规则将项目进行拆分。如果将一个项目的所有模板功能都写到一起,当某个模块出现问题时将直接导致整个服务器出现问题。拆分按照业务拆分为不同的服务器,有效的降低系统架构的耦合性在业务拆分的基础上可按照代码层级进行拆分(view、controller、service、pojo)分布式思想的实质分布式思想的实质是为了系统的..._分布式系统运维工具
文章浏览阅读174次。1.数据源准备2.数据处理step1:数据表处理应用函数:①VLOOKUP函数; ② CONCATENATE函数终表:step2:数据透视表统计分析(1) 透视表汇总不同渠道用户数, 金额(2)透视表汇总不同日期购买用户数,金额(3)透视表汇总不同用户购买订单数,金额step3:讲第二步结果可视化, 比如, 柱形图(1)不同渠道用户数, 金额(2)不同日期..._exce l趋势分析数据量
文章浏览阅读3.3k次。堡垒机可以为企业实现服务器、网络设备、数据库、安全设备等的集中管控和安全可靠运行,帮助IT运维人员提高工作效率。通俗来说,就是用来控制哪些人可以登录哪些资产(事先防范和事中控制),以及录像记录登录资产后做了什么事情(事后溯源)。由于堡垒机内部保存着企业所有的设备资产和权限关系,是企业内部信息安全的重要一环。但目前出现的以下问题产生了很大安全隐患:密码设置过于简单,容易被暴力破解;为方便记忆,设置统一的密码,一旦单点被破,极易引发全面危机。在单一的静态密码验证机制下,登录密码是堡垒机安全的唯一_horizon宁盾双因素配置
文章浏览阅读7.7k次,点赞4次,收藏16次。Chrome作为一款挺不错的浏览器,其有着诸多的优良特性,并且支持跨平台。其支持(Windows、Linux、Mac OS X、BSD、Android),在绝大多数情况下,其的安装都很简单,但有时会由于网络原因,无法安装,所以在这里总结下Chrome的安装。Windows下的安装:在线安装:离线安装:Linux下的安装:在线安装:离线安装:..._chrome linux debian离线安装依赖
文章浏览阅读153次。中国发达城市榜单每天都在刷新,但无非是北上广轮流坐庄。北京拥有最顶尖的文化资源,上海是“摩登”的国际化大都市,广州是活力四射的千年商都。GDP和发展潜力是衡量城市的数字指...
文章浏览阅读3.3k次。前言spark在java使用比较少,多是scala的用法,我这里介绍一下我在项目中使用的代码配置详细算法的使用请点击我主页列表查看版本jar版本说明spark3.0.1scala2.12这个版本注意和spark版本对应,只是为了引jar包springboot版本2.3.2.RELEASEmaven<!-- spark --> <dependency> <gro_使用java调用spark注册进去的程序
文章浏览阅读4.8k次。汽车零部件开发工具巨头V公司全套bootloader中UDS协议栈源代码,自己完成底层外设驱动开发后,集成即可使用,代码精简高效,大厂出品有量产保证。:139800617636213023darcy169_uds协议栈 源代码
文章浏览阅读4.6k次,点赞20次,收藏148次。AUTOSAR基础篇之OS(下)前言首先,请问大家几个小小的问题,你清楚:你知道多核OS在什么场景下使用吗?多核系统OS又是如何协同启动或者关闭的呢?AUTOSAR OS存在哪些功能安全等方面的要求呢?多核OS之间的启动关闭与单核相比又存在哪些异同呢?。。。。。。今天,我们来一起探索并回答这些问题。为了便于大家理解,以下是本文的主题大纲:[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-JCXrdI0k-1636287756923)(https://gite_autosar 定义了 5 种多核支持类型
文章浏览阅读2.2k次,点赞6次,收藏14次。原因:自己写的头文件没有被加入到方案的包含目录中去,无法被检索到,也就无法打开。将自己写的头文件都放入header files。然后在VS界面上,右键方案名,点击属性。将自己头文件夹的目录添加进去。_vs2013打不开自己定义的头文件
文章浏览阅读3.3w次,点赞80次,收藏342次。此时,可以将系统中所有用户的 Session 数据全部保存到 Redis 中,用户在提交新的请求后,系统先从Redis 中查找相应的Session 数据,如果存在,则再进行相关操作,否则跳转到登录页面。此时,可以将系统中所有用户的 Session 数据全部保存到 Redis 中,用户在提交新的请求后,系统先从Redis 中查找相应的Session 数据,如果存在,则再进行相关操作,否则跳转到登录页面。当数据量很大时,count 的数量的指定可能会不起作用,Redis 会自动调整每次的遍历数目。_redis命令
文章浏览阅读449次,点赞3次,收藏3次。URP的设计目标是在保持高性能的同时,提供更多的渲染功能和自定义选项。与普通项目相比,会多出Presets文件夹,里面包含着一些设置,包括本色,声音,法线,贴图等设置。全局只有主光源和附加光源,主光源只支持平行光,附加光源数量有限制,主光源和附加光源在一次Pass中可以一起着色。URP:全局只有主光源和附加光源,主光源只支持平行光,附加光源数量有限制,一次Pass可以计算多个光源。可编程渲染管线:渲染策略是可以供程序员定制的,可以定制的有:光照计算和光源,深度测试,摄像机光照烘焙,后期处理策略等等。_urp渲染管线