Python-VTK批量读取二维切片并显示三维模型

 更新时间:2022年04月18日 12:46:52   作者:派大大大星   
这篇文章主要介绍了Python-VTK批量读取二维切片并显示三维模型,文章基于python的相关资料展开对主题的详细介绍,具有一定的参考价值,需要的小伙伴可以参考一下

前言:

使用Python-VTK批量读取二维切片,并显示三维模型

VTK,(visualizationtoolkit)是一个开放资源的免费软件系统,主要用于三维计算机图形学、图像处理和可视化。Vtk是在面向对象原理的基础上设计和实现的,它的内核是用C++构建的,包含有大约250,000行代码,2000多个类,还包含有几个转换界面,因此也可以自由的通过Java,Tcl/Tk和Python各种语言使用vtk。

主要函数介绍

vtk.vtkJPEGReader(): 读取图片的接口,可以通过该接口设置读取的图片位置、大小、维度。需要注意的是,如果需要将一列二维图片读取成一个三维图像需要将一些图片命名成自然顺序0、1、2、3···,这样VTK才能正确的读取图片。如果需要读取其他格式的图片只需要将函数中的JPEG换成对应的格式即可,如PNG。

vtk.vtkMarchingCubes(): 三维重建函数,MarchingCubes算法实际上是一个分而治之的方法,因为其将等值面的抽取分布于每一个体素(voxel)中进行。对于每个被处理的体素,以三角面片来逼近其内部的等值面。每个体素是一个小立方体(cube),在构造三角面片的处理过程中对每个体素都“扫描”一遍,就好像是一个处理器在这些体素上移动一样,也因此而得名。··在等值面抽取的过程中将一系列二维切片数据看做是一个三维的数据场,从中将具有某种阈值的物质抽取出来,以某种拓扑形式连接成三角面片,所以MC算法也被称为“等值面提取”(IsosurfaceExtraction)算法。··在医学应用上,采用MC算法可以重建人体外部轮廓、内部组织器官,使医生能够直接在三维图像上观察感兴趣的器官与周围组织的空间关系。

import vtk
# 定义渲染窗口、交互模式
aRender = vtk.vtkRenderer()
Renwin = vtk.vtkRenderWindow()
Renwin.AddRenderer(aRender)
iren = vtk.vtkRenderWindowInteractor()
iren.SetRenderWindow(Renwin)

# 定义个图片读取接口
#读取PNG图片就换成PNG_Reader = vtk.vtkPNGReader()
Jpg_Reader = vtk.vtkJPEGReader()
Jpg_Reader.SetNumberOfScalarComponents(1)
Jpg_Reader.SetFileDimensionality(3)  # 说明图像是三维的
 # 定义图像大小,本行表示图像大小为(512*512*240)
Jpg_Reader.SetDataExtent(0, 512, 0, 512, 0, 240) 
 # 设置图像的存放位置
Jpg_Reader.SetFilePrefix("E:/outer/label/5/")
 # 设置图像前缀名字
 #表示图像前缀为数字(如:0.jpg)
Jpg_Reader.SetFilePattern("%s%d.jpg")
Jpg_Reader.Update()
Jpg_Reader.SetDataByteOrderToLittleEndian()

# 计算轮廓的方法
contour = vtk.vtkMarchingCubes()  
contour.SetInputConnection(Jpg_Reader.GetOutputPort())
contour.ComputeNormalsOn()
contour.SetValue(0, 255)


mapper = vtk.vtkPolyDataMapper()
mapper.SetInputConnection(contour.GetOutputPort())
mapper.ScalarVisibilityOff()

actor = vtk.vtkActor()
actor.SetMapper(mapper)

renderer = vtk.vtkRenderer()
renderer.SetBackground([0.1, 0.1, 0.5])
renderer.AddActor(actor)

window = vtk.vtkRenderWindow()
window.SetSize(512, 512)
window.AddRenderer(renderer)


interactor = vtk.vtkRenderWindowInteractor()
interactor.SetRenderWindow(window)

# 开始显示
window.Render()
interactor.Initialize()
interactor.Start()

效果展示

1636290573(1).jpg

到此这篇关于Python-VTK批量读取二维切片并显示三维模型的文章就介绍到这了,更多相关Python-VTK读取二维切片内容请搜索脚本之家以前的文章或继续浏览下面的相关文章希望大家以后多多支持脚本之家!

相关文章

  • python实现将一维列表转换为多维列表(numpy+reshape)

    python实现将一维列表转换为多维列表(numpy+reshape)

    今天小编就为大家分享一篇python实现将一维列表转换为多维列表(numpy+reshape),具有很好的参考价值,希望对大家有所帮助。一起跟随小编过来看看吧
    2019-11-11
  • 联邦学习神经网络FedAvg算法实现

    联邦学习神经网络FedAvg算法实现

    这篇文章主要为大家介绍了联邦学习神经网络FedAvg算法实现,有需要的朋友可以借鉴参考下,希望能够有所帮助,祝大家多多进步,早日升职加薪
    2022-05-05
  • Python使用pydub模块转换音频格式以及对音频进行剪辑

    Python使用pydub模块转换音频格式以及对音频进行剪辑

    这篇文章主要给大家介绍了关于Python使用pydub模块转换音频格式以及对音频进行剪辑的相关资料pydub是python的高级一个音频处理库,可以让你以一种不那么蠢的方法处理音频。需要的朋友可以参考下
    2021-06-06
  • python监控日志中的报错并进行邮件报警

    python监控日志中的报错并进行邮件报警

    这篇文章主要为大家介绍了python监控日志中的报错并进行邮件报警的实现实例,有需要的朋友可以借鉴参考下,希望能够有所帮助,祝大家多多进步,早日升职加薪
    2022-06-06
  • 解决编码问题:UnicodeDecodeError: 'utf-8' codec can't decod

    解决编码问题:UnicodeDecodeError: 'utf-8' codec

    这篇文章主要介绍了快速解决编码问题:UnicodeDecodeError: 'utf-8' codec can't decod,具有很好的参考价值,希望对大家有所帮助。如有错误或未考虑完全的地方,望不吝赐教
    2022-05-05
  • Django ORM查询操作方式

    Django ORM查询操作方式

    Django提供了一套非常方便的类似SqlAlchemy ORM的通过对象调用的方式操作数据库表的ORM框架,,本文给大家详细介绍Django ORM查询操作方式,感兴趣的朋友一起看看吧
    2023-10-10
  • python的继承详解

    python的继承详解

    这篇文章主要介绍了Python中继承的的相关资料,文中讲解非常细致,代码帮助大家更好的理解和学习,感兴趣的朋友可以了解下
    2021-10-10
  • pycharm找不到conda可执行文件怎么解决

    pycharm找不到conda可执行文件怎么解决

    这篇文章主要给大家介绍了关于pycharm找不到conda可执行文件怎么解决的相关资料,文中通过图文介绍非常的详细,对大家的学习或者工作具有一的参考学习价值,需要的朋友可以参考下
    2023-07-07
  • pycharm实现猜数游戏

    pycharm实现猜数游戏

    这篇文章主要为大家详细介绍了pycharm实现猜数游戏,文中示例代码介绍的非常详细,具有一定的参考价值,感兴趣的小伙伴们可以参考一下
    2020-12-12
  • 解决python3中自定义wsgi函数,make_server函数报错的问题

    解决python3中自定义wsgi函数,make_server函数报错的问题

    下面小编就为大家分享一篇解决python3中自定义wsgi函数,make_server函数报错的问题,具有很好的参考价值,希望对大家有所帮助。一起跟随小编过来看看吧
    2017-11-11

最新评论