图形处理之抽取轮廓

最新推荐文章于 2024-07-29 22:22:26 发布
最新推荐文章于 2024-07-29 22:22:26 发布
阅读量7.1k 收藏 13
点赞数 1
分类专栏： Study-VTK 文章标签： vtk 等值面 轮廓抽取 等值线
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/webzhuce/article/details/69092342
版权

等值面（线）

        抽取轮廓的操作对象是标量数据。其思想是：将数据集中标量值等于某一指定恒量值的部分提取出来。对于3D的数据集而言，产生的是一个等值面；对于2D的数据集而言，产生的是一个等值线。其典型的应用有气象图中的等温线、地形图中的等高线。对于医学数据而言，不同的标量值代表的是人体的不同部分，因而可以分别提取出人的皮肤或骨头。
　　抽取轮廓的功能是由一个过滤器实现的，如vtkContourFilter 。vtkContourFilter 可以接受任意数据集类型作为输入，因而具有 一般性。使用vtkContourFilter 时，除了需要设置输入数据集外，还需要指定一个或多个用于抽取的标量值。可用如下两种方法进行设置。

  1. 使用方法SetValue()逐个设置抽取值。该方法有个两个参数：第一个参数是抽取值的索引号，表示第几个 抽取值。索引号从0开始计数；第二个参数就是指定的抽取值。
  2. 使用方法GenerateValues()自动产生一系列抽取值。该方法有三个参数：第一个参数是抽取值的个数，后面两个参数是抽取值的取值范围。例如：
vtkSmartPointer<vtkContourFilter> contours = vtkSmartPointer<vtkContourFilter>::New();
contours->GenerateValues(5.0, 0.0, 1.2);

最终产生的抽取值为0.0， 0.3， 0.6， 0.9和1.0，这是一个等差数列。

vtkContourFilter生成等值线示例

CMakeLists.txt文件代码如下：

CMAKE_MINIMUM_REQUIRED(VERSION 2.8)
PROJECT(Contour)
FIND_PACKAGE(VTK REQUIRED)
INCLUDE(${VTK_USE_FILE})
ADD_EXECUTABLE(Contour	Contour.cpp)
TARGET_LINK_LIBRARIES(Contour ${VTK_LIBRARIES})

Contour.cpp文件代码如下：

#include <vtkActor.h>
#include <vtkCamera.h>
#include <vtkCellArray.h>
#include <vtkFloatArray.h>
#include <vtkPointData.h>
#include <vtkPoints.h>
#include <vtkPolyData.h>
#include <vtkPolyDataMapper.h>
#include <vtkRenderWindow.h>
#include <vtkRenderer.h>
#include <vtkLookupTable.h>
#include <vtkContourFilter.h>
#include <vtkRenderWindowInteractor.h>
#include <vtkProperty.h>
#include <vtkSmartPointer.h>

int main()
{
	int i;
	//定义网格的顶点坐标
	vtkSmartPointer<vtkPoints> points = vtkSmartPointer<vtkPoints>::New();
	//第1行5个点
	points->InsertNextPoint(-2.0, 2.0, 0.0);
	points->InsertNextPoint(-1.0, 2.0, 0.0);
	points->InsertNextPoint(0.0, 2.0, 0.0);
	points->InsertNextPoint(1.0, 2.0, 0.0);
	points->InsertNextPoint(2.0, 2.0, 0.0);
	//第2行5个点
	points->InsertNextPoint(-2.0, 1.0, 0.0);
	points->InsertNextPoint(-1.0, 1.0, 0.0);
	points->InsertNextPoint(0.0, 1.0, 0.0);
	points->InsertNextPoint(1.0, 1.0, 0.0);
	points->InsertNextPoint(2.0, 1.0, 0.0);
	//第3行5个点
	points->InsertNextPoint(-2.0, 0.0, 0.0);
	points->InsertNextPoint(-1.0, 0.0, 0.0);
	points->InsertNextPoint(0.0, 0.0, 0.0);
	points->InsertNextPoint(1.0, 0.0, 0.0);
	points->InsertNextPoint(2.0, 0.0, 0.0);
	//第4行5个点
	points->InsertNextPoint(-2.0, -1.0, 0.0);
	points->InsertNextPoint(-1.0, -1.0, 0.0);
	points->InsertNextPoint(0.0, -1.0, 0.0);
	points->InsertNextPoint(1.0, -1.0, 0.0);
	points->InsertNextPoint(2.0, -1.0, 0.0);
	//第5行5个点
	points->InsertNextPoint(-2.0, -2.0, 0.0);
	points->InsertNextPoint(-1.0, -2.0, 0.0);
	points->InsertNextPoint(0.0, -2.0, 0.0);
	points->InsertNextPoint(1.0, -2.0, 0.0);
	points->InsertNextPoint(2.0, -2.0, 0.0);

	//定义单元，每个顶点建立一个四边形单元，共计个单元
	static vtkIdType pts[16][4] = { { 0, 5, 6, 1 }, { 1, 6, 7, 2 }, { 2, 7, 8, 3 },
	{ 3, 8, 9, 4 }, { 5, 10, 11, 6 }, { 6, 11, 12, 7 }, { 7, 12, 13, 8 }, { 8, 13, 14, 9 },
	{ 10, 15, 16, 11 }, { 11, 16, 17, 12 }, { 12, 17, 18, 13 }, { 13, 18, 19, 14 },
	{ 15, 20, 21, 16 }, { 16, 21, 22, 17 }, { 17, 22, 23, 18 }, { 18, 23, 24, 19 } };

	//创建对象
	vtkSmartPointer<vtkPolyData> pGrid = vtkSmartPointer<vtkPolyData>::New();
	vtkSmartPointer<vtkCellArray> polys = vtkSmartPointer<vtkCellArray>::New();
	for (i = 0; i < 16; i++)
	{
		polys->InsertNextCell(4, pts[i]);
	}
	vtkSmartPointer<vtkFloatArray> scalars = vtkSmartPointer<vtkFloatArray>::New();
	//第1行
	scalars->InsertTuple1(0, 0);
	scalars->InsertTuple1(1, 1);
	scalars->InsertTuple1(2, 1);
	scalars->InsertTuple1(3, 3);
	scalars->InsertTuple1(4, 2);
	//第2行
	scalars->InsertTuple1(5, 1);
	scalars->InsertTuple1(6, 3);
	scalars->InsertTuple1(7, 6);
	scalars->InsertTuple1(8, 6);
	scalars->InsertTuple1(9, 3);
	//第3行
	scalars->InsertTuple1(10, 3);
	scalars->InsertTuple1(11, 7);
	scalars->InsertTuple1(12, 9);
	scalars->InsertTuple1(13, 7);
	scalars->InsertTuple1(14, 3);
	//第4行
	scalars->InsertTuple1(15, 2);
	scalars->InsertTuple1(16, 7);
	scalars->InsertTuple1(17, 8);
	scalars->InsertTuple1(18, 6);
	scalars->InsertTuple1(19, 2);
	//第5行
	scalars->InsertTuple1(20, 1);
	scalars->InsertTuple1(21, 2);
	scalars->InsertTuple1(22, 3);
	scalars->InsertTuple1(23, 4);
	scalars->InsertTuple1(24, 3);

	pGrid->SetPoints(points);
	pGrid->SetPolys(polys);
	pGrid->GetPointData()->SetScalars(scalars);

	//定义颜色映射表
	vtkSmartPointer<vtkLookupTable> pColorTable = vtkSmartPointer<vtkLookupTable>::New();
	pColorTable->SetHueRange(0.667, 0.887);
	pColorTable->SetNumberOfColors(255);
	pColorTable->Build();

	vtkSmartPointer<vtkPolyDataMapper> cubeMapper = vtkSmartPointer<vtkPolyDataMapper>::New();
	cubeMapper->SetInputData(pGrid);
	cubeMapper->SetScalarRange(0, 9);
	cubeMapper->SetLookupTable(pColorTable);
	vtkSmartPointer<vtkActor> cubeActor = vtkSmartPointer<vtkActor>::New();
	cubeActor->SetMapper(cubeMapper);

	//提取轮廓线
	vtkSmartPointer<vtkContourFilter> pContourFilter = vtkSmartPointer<vtkContourFilter>::New();
	pContourFilter->SetValue(0, 5.0);
	pContourFilter->SetInputData(pGrid);
	pContourFilter->Update();

	vtkSmartPointer<vtkPolyDataMapper> pContourMapper = vtkSmartPointer<vtkPolyDataMapper>::New();
	pContourMapper->SetInputData(pContourFilter->GetOutput());

	pContourMapper->ScalarVisibilityOff();
	vtkSmartPointer<vtkActor> pContourActor = vtkSmartPointer<vtkActor>::New();
	pContourActor->SetMapper(pContourMapper);
	pContourActor->GetProperty()->SetColor(1.0, 1.0, 0.0);
	pContourActor->GetProperty()->SetLineWidth(5.0);

	vtkSmartPointer<vtkRenderer> renderer = vtkSmartPointer<vtkRenderer>::New();
	renderer->AddActor(pContourActor);
	renderer->AddActor(cubeActor);
	renderer->SetBackground(1.0, 1.0, 1.0);

	vtkSmartPointer<vtkRenderWindow> renWin = vtkSmartPointer<vtkRenderWindow>::New();
	renWin->SetSize(300, 300);
	renWin->AddRenderer(renderer);
	renWin->SetWindowName("vtkContourFilter");

	vtkSmartPointer<vtkRenderWindowInteractor> interactor = vtkSmartPointer<vtkRenderWindowInteractor>::New();
	interactor->SetRenderWindow(renWin);

	interactor->Initialize();
	interactor->Start();
}

运行结果

这里写图片描述

示例原理

  1. 定义规则网格数据顶点的属性值，示例中设置标量值。
  2. 确定提取区域的属性值。
  3. 用提取区域的属性值和每个网格数据的每个单元进行比对，确定其是否在该单元内，如果在该单元内，要用插值的方法确定其和该单元边界的交点，然后继续比对下一个单元，如果不在该单元内，继续比对下一个单元。
  4. 将各个单元的交点相连成线，即为边界区域。
    ##原理示意图：
    这里写图片描述
阿兵-AI医疗
关注
专栏目录
基于道格拉斯普克算法的轮廓点简化
03-21 579
道格拉斯普克算法（DP）对元素进行简化，这里介绍对轮廓点进行简化，如下图展示了轮廓点提取、DP保留的关键点、关键点生成的矢量线： 有序轮廓点 DP简化轮廓点后的关键点 关键点相连生成的轮廓线局部放大图 1、算法原理 道格拉斯普克算法（Douglas-Peukcer）算法是一种简化线状要素的经典算法。先介绍其原始的计算原理：对每一条曲线的首末点虚连一条直线，求所有点与直线的距离,并找出最大距离值dmax，用dmax与限差D相比。若dmax<D，这条曲线上的中间点全部舍去；若dmax ≥D，..
VTK笔记-CT图像获取皮肤等值面-vtkContourFilter类的使用
黑山老妖的博客
05-09 3860
系列文章目录 提示：这里可以添加系列文章的所有文章的目录，目录需要自己手动添加 例如：第一章 Python 机器学习入门之pandas的使用 提示：写完文章后，目录可以自动生成，如何生成可参考右边的帮助文档 文章目录系列文章目录前言一、pandas是什么？二、使用步骤1.引入库2.读入数据总结 前言 提示：这里可以添加本文要记录的大概内容： 例如：随着人工智能的不断发展，机器学习这门技术也越来越重要，很多人都开启了学习机器学习，本文就介绍了机器学习的基础内容。 提示：以下是本篇文章正文内容，下面案例
[OpenGL] 基于屏幕空间的轮廓线提取
(ง •̀_•́)ง
01-13 3280
  概念引入          《Real-time Rendering》一书中介绍了多种轮廓线提取的方式，根据对比以及综合考虑，最终使用了图像空间的方法来实现实时提取轮廓。该方法需要先将场景中的物体信息先渲染到纹理中，对使用延迟渲染的框架非常友好，并且弊端也是多种轮廓线提取方法中比较少的，具体内容可以参照原书的NPR章节。 方法概述           (1) 准备2个颜色缓冲区，1...
【VTK 图形基本操作进阶_表面重建技术（等值面提取）】
最新发布
weixin_42483745的博客
07-29 936
由于这一方法原理简单，易于实现，目前已经得到了较为广泛的应用，称为三维数据等值面生成的经典算法。等值面提取类根据数据类型的不同而有所侧重。首先通过一个reader对象来读取一幅图像，并将输入到vtkMarchingCubes中，提取等值面时，最重要的是要设置等值面的数值，SetValue()函数用于设置等值面的值，其第一个参数表示等值面的序号，因此可以通过这个函数设置多个等值面值来提取多个等值面。此时可通过等值面提取技术，仅提取感兴趣的一个或者几个组织轮廓，并生成网格模型以供后续的处理和研究。
VTK轮廓提取
OneSea的专栏
06-06 4193
参考资料 [1]VTK教程之九 可视化基础算法-二维轮廓线（等值线）提取
vtk 曲线 样式_VTK教程之九 可视化基础算法-二维轮廓线（等值线）提取
weixin_39591720的博客
12-21 1142
二维轮廓线(等值线)提取属于标量可视化技术范围，常用于提取一个研究区域的轮廓边界，如在气象中应用中，常根据气象卫星或数值模拟的结果，按照每个地理位置点确定的属性值，分析区域内的降雨情况，提取不同降雨量的区域边界(如大雨、中雨、小雨的区域边界)，在有限元分析中，常用于提取位移等值线等，在VTK中轮廓线提取功能主要面对的是规则网格数据集，对于非规则网格数据，必须要进行插值计算，生成规则网格数据集，由于...
bs.rar_图形图像处理_C/C++_
08-12
总结来说，C/C++在图形图像处理方面具有强大的功能，通过使用OpenCV等库，开发者可以实现条纹图分析、二值化、细化、判读和三维轮廓图生成等一系列复杂任务。对于这些技术的深入理解和熟练运用，对于提升图像处理...
实时图像轮廓抽取算法研究.docx
06-21
3. **图像演示软件**：作者们开发的软件可能包含了一个图形用户界面(GUI)，用于显示和交互式地处理图像，它能实时显示轮廓抽取的结果，帮助用户理解和评估算法性能。 4. **代数特征抽取**：在相关文献中提到，可能...
基于边分类的三维模型外轮廓线抽取算法 (2005年)
04-27
在三维图形处理领域，三维模型的轮廓线抽取是一项重要的任务，尤其是在渲染、可视化以及几何处理等应用中。传统的轮廓线抽取方法往往计算量大且容易受到模型精度的影响，特别是在处理粗糙或低精度模型时。因此，开发...
Opencv python 图像轮廓处理
weixin_52662649的博客
04-29 1514
Opencv Python 轮廓处理 Opencv python轮廓处理 本人大一数学系学渣菜鸟一枚，课程之余学习了OpenCV，把写的关于轮廓处理的知识分享一下，欢迎在评论区里给出宝贵意见！ 我把代码放到一起了，要是运行的话最好复制到新的python文件里运行哦 以下为BB子整理的OpenCV中基本的轮廓处理部分 #program 1.绘制轮廓 import numpy as np import cv2 as cv #读取图像 img = cv.imread("G:\sundries\CVpicture
VTK：对输入的三维模型在某个方向等间距提取模型的切面轮廓线
HW140701的博客
12-04 5481
测试的三维模型为bunny，即斯坦福兔子 代码示例： #include <vtkOBJReader.h> #include <vtkSmartPointer.h> #include <vtkCutter.h> #include <vtkPolyDataMapper.h> #include <vtkPlane.h> #inclu...
VTK vtkContourWidget 提取感兴趣区域
q610098308的专栏
06-19 951
主要用来表示一个点集，点集之间有连线,相邻点，而非任意点。最后一个点的选取决定了，轮廓是闭合的还是是开放的。 通常可用于兴趣区域提取； 见图 官方 demo
从分析blender轮廓提取技术中学习opengl如何绘制三维模型轮廓
微型蚂蚁的专栏
03-17 7051
项目需要对blender轮廓提取功能进行分析和提取：从blender的基本操作入手学习如何建模和修改，发现在模型（object）状态下选中一个三维模型系统会自动产生该三维模型的投影轮廓，并且移动旋转视图等操作，轮廓图像都会实时修改，于是选择最有可能发现轮廓提取算法所在的选择视图操作相关代码进行跟踪学习。发现：绘制轮廓的入口点在：BL_src : drawobject.cstatic void dr
VTK:轮廓最小凸包提取与显示
cfqcfqcfqcfqcfq的博客
11-20 4461
vtk提取凸包使用类 PointsPrjectHull。该类可以获取任意点集的最小凸包。该类输入为点集。输出为包围该点集的最小凸包轮廓点集。前面一章讲了如何寻找图像轮廓并对其进行连通性分析(连通区域分析)。现在以某一连通域的轮廓点集作为输入提取其最小凸闭包。将提取到的凸闭包用vtkPlolyLine来保存显示 关键代码： vtkSmartPointer point0=
VTK学习笔记（三十一）vtkCutter类对多边形取切面轮廓
juluwangriyue的博客
04-13 1890
VTK学习笔记（三十一）vtkCutter类对多边形取切面轮廓1、tkCutter类对多边形取切面轮廓1.1、vtkCutter类1.2、接口2、python 调用 1、tkCutter类对多边形取切面轮廓 1.1、vtkCutter类 vtkCutter可以使用使用用户指定的隐式函数切割vtkDataSet。 vtkCutter是一个过滤器，可以使用vtkImplicitFunction 的任何子类来切割数据。也就是说，对应于隐式函数 F(x,y,z) = value(s) 创建了一个多边形表面，可以在其
VTK 二维 外部轮廓线绘制
q610098308的专栏
12-06 803
有时我们需求，三维的轮廓线 可以用来做什么，自己想一下；
VTK：轮廓 Glow Pass用法实战
希望我的博客，能帮上你解决学习中工作中所遇到的问题
03-03 383
VTK：轮廓 Glow Pass用法实战程序输出程序完整源代码 程序输出 程序完整源代码 #include <vtkActor.h> #include <vtkArrowSource.h> #include <vtkCamera.h> #include <vtkNamedColors.h> #include <vtkNew.h> #include <vtkOpenGLRenderer.h> #include <vtkOutlin
VTK教程之十 可视化基础算法-三维轮廓面（等值面）提取
云海唯C的专栏
10-31 7389
VTK教程之十 可视化基础算法-三维轮廓面（等值面）提取 在VTK中，轮廓面提取功能主要面对的是三维规则网格数据集，三维规则网格数据集的单元一般为立方体单元，轮廓面的提取主要采用移动立方体算法(Marching Cubes)，VTK提供的类对该算法进行了封装，下面给出三维轮廓面提取的示例程序。 #include "stdafx.h" #include #inc
[VTK]基于VTK的任意平面切割
aees4684的博客
02-17 623
// 先贴码 以后再。。。// 切割介绍// 对于一个模型的切割需要怎么办呢，想想切西瓜就知道，首先需要有一个模型、然后有一个切割平面// 接着对于每个切割操作来更新模型，这样就可以得到切割的效果了#include "vtkPlanes.h" #include "vtkProperty.h" #include "vtkRenderer.h" #include "vtkRende...
评论 4
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值