如果电脑没有安装OpenCV,那么请问该如何仅用VTK实现图形任意平面的切割?求大佬代码指点!

网上找了一些关于基于VTK的任意平面切割,但是代码都是那一个一样的基于OpenCV
的代码,而假如我不想使用OpenCV,(而且我还找不到头文件"VTKReBuild.h")那么这个程序改怎么写呢?而且我想切割出来的是一个有一定厚度的立体图,相当于一个圆环类似(以人体模型为例,若选择从腹部切割开来的话,保留的便是以切割平面为基准大概1cm上下的切割图)恳求大佬指点!

3个回答

duanbanyu
duanbanyu 然后错误是在vtkTextRender.h里报错引发了未经处理的异常:读取访问权限冲突。 **this** 是 nullptr。
一年多之前 回复
duanbanyu
duanbanyu 不胜感激!
一年多之前 回复
duanbanyu
duanbanyu 麻烦大佬再为我解答一下,下面的问题,我就改了输入我的dcm文档和 删去了vs2017无法识别的int retVal = vtkRegressionTestImage(renWin); if ( retVal == vtkRegressionTester::DO_INTERACTOR),程序为什么就运行不了呢?
一年多之前 回复

十分感谢你的回答,但是我想请问下头文件vtkTestUtilities.h和vtkRegressionTestImage.h是哪来的呢?请问是在安装VTK的时候我漏装了环境而使我没有这两个头文件吗?

#include "vtkAutoInit.h"

VTK_MODULE_INIT(vtkRenderingOpenGL2);
VTK_MODULE_INIT(vtkInteractionStyle);
#include "vtkSmartPointer.h"
#include "vtkRegressionTestImage.h"
#include "vtkActor.h"
#include "vtkCamera.h"
#include "vtkCellPicker.h"
#include "vtkCommand.h"
#include "vtkImageActor.h"
#include "vtkImageReslice.h"
#include "vtkInteractorStyleImage.h"
#include "vtkImageMapToColors.h"
#include "vtkImagePlaneWidget.h"
#include "vtkImageReader.h"
#include "vtkInteractorEventRecorder.h"
#include "vtkLookupTable.h"
#include "vtkOutlineFilter.h"
#include "vtkDICOMImageReader.h"
#include "vtkPolyDataMapper.h"
#include "vtkProperty.h"
#include "vtkRenderWindow.h"
#include "vtkRenderWindowInteractor.h"
#include "vtkRenderer.h"
#include "vtkVolume16Reader.h"
#include "vtkImageData.h"
#include "vtkPointData.h"
#include "vtkPlaneSource.h"
#include "vtkPlane.h"
#include "vtkResliceCursorActor.h"
#include "vtkResliceCursorPolyDataAlgorithm.h"
#include "vtkResliceCursor.h"
#include "vtkResliceCursorWidget.h"
#include "vtkResliceCursorLineRepresentation.h"
#include "vtkBiDimensionalWidget.h"

#include "vtkTestUtilities.h"

//----------------------------------------------------------------------------
class vtkResliceCursorCallback : public vtkCommand
{
public:
static vtkResliceCursorCallback *New()
{
return new vtkResliceCursorCallback;
}

void Execute(vtkObject *caller, unsigned long /*ev*/,
    void *callData) override
{
    vtkImagePlaneWidget* ipw =
        dynamic_cast< vtkImagePlaneWidget* >(caller);
    if (ipw)
    {
        double* wl = static_cast<double*>(callData);

        if (ipw == this->IPW[0])
        {
            this->IPW[1]->SetWindowLevel(wl[0], wl[1], 1);
            this->IPW[2]->SetWindowLevel(wl[0], wl[1], 1);
        }
        else if (ipw == this->IPW[1])
        {
            this->IPW[0]->SetWindowLevel(wl[0], wl[1], 1);
            this->IPW[2]->SetWindowLevel(wl[0], wl[1], 1);
        }
        else if (ipw == this->IPW[2])
        {
            this->IPW[0]->SetWindowLevel(wl[0], wl[1], 1);
            this->IPW[1]->SetWindowLevel(wl[0], wl[1], 1);
        }
    }

    vtkResliceCursorWidget *rcw = dynamic_cast<
        vtkResliceCursorWidget * >(caller);
    if (rcw)
    {
        vtkResliceCursorLineRepresentation *rep = dynamic_cast<
            vtkResliceCursorLineRepresentation * >(rcw->GetRepresentation());
        vtkResliceCursor *rc = rep->GetResliceCursorActor()->
            GetCursorAlgorithm()->GetResliceCursor();
        for (int i = 0; i < 3; i++)
        {
            vtkPlaneSource *ps = static_cast< vtkPlaneSource * >(
                this->IPW[i]->GetPolyDataAlgorithm());
            ps->SetNormal(rc->GetPlane(i)->GetNormal());
            ps->SetCenter(rc->GetPlane(i)->GetOrigin());

            // If the reslice plane has modified, update it on the 3D widget
            this->IPW[i]->UpdatePlacement();

            //std::cout << "Updating placement of plane: " << i << " " <<
            //  rc->GetPlane(i)->GetNormal()[0] << " " <<
            //  rc->GetPlane(i)->GetNormal()[1] << " " <<
            //  rc->GetPlane(i)->GetNormal()[2] << std::endl;
            //this->IPW[i]->GetReslice()->Print(cout);
            //rep->GetReslice()->Print(cout);
            //std::cout << "---------------------" << std::endl;
        }
    }

    // Render everything
    this->RCW[0]->Render();
}

vtkResliceCursorCallback() = default;
vtkImagePlaneWidget* IPW[3];
vtkResliceCursorWidget *RCW[3];

};

//----------------------------------------------------------------------------
int main(int argc, char *argv[])
{
vtkSmartPointer reader =
vtkSmartPointer::New();
reader->SetDataDimensions(64, 64);
reader->SetDataByteOrderToLittleEndian();
reader->SetImageRange(1, 93);
reader->SetDataSpacing(3.2, 3.2, 1.5);
reader->SetFilePrefix("D:/programe/VTK/mhd/肝脏/1/062d1");
reader->SetFilePattern("%s%30a.dcm");
reader->ReleaseDataFlagOn();
reader->SetDataMask(0x7fff);
reader->Update();

vtkSmartPointer<vtkOutlineFilter> outline =
    vtkSmartPointer<vtkOutlineFilter>::New();
outline->SetInputConnection(reader->GetOutputPort());

vtkSmartPointer<vtkPolyDataMapper> outlineMapper =
    vtkSmartPointer<vtkPolyDataMapper>::New();
outlineMapper->SetInputConnection(outline->GetOutputPort());

vtkSmartPointer<vtkActor> outlineActor =
    vtkSmartPointer<vtkActor>::New();
outlineActor->SetMapper(outlineMapper);

vtkSmartPointer<vtkRenderer> ren[4];

vtkSmartPointer<vtkRenderWindow> renWin =
    vtkSmartPointer<vtkRenderWindow>::New();
renWin->SetMultiSamples(0);

for (int i = 0; i < 4; i++)
{
    ren[i] = vtkSmartPointer<vtkRenderer>::New();
    renWin->AddRenderer(ren[i]);
}

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

vtkSmartPointer<vtkCellPicker> picker =
    vtkSmartPointer<vtkCellPicker>::New();
picker->SetTolerance(0.005);

vtkSmartPointer<vtkProperty> ipwProp =
    vtkSmartPointer<vtkProperty>::New();


//assign default props to the ipw's texture plane actor
vtkSmartPointer<vtkImagePlaneWidget> planeWidget[3];
int imageDims[3];
reader->GetOutput()->GetDimensions(imageDims);

for (int i = 0; i < 3; i++)
{
    planeWidget[i] = vtkSmartPointer<vtkImagePlaneWidget>::New();
    planeWidget[i]->SetInteractor(iren);
    planeWidget[i]->SetPicker(picker);
    planeWidget[i]->RestrictPlaneToVolumeOn();
    double color[3] = { 0, 0, 0 };
    color[i] = 1;
    planeWidget[i]->GetPlaneProperty()->SetColor(color);
    planeWidget[i]->SetTexturePlaneProperty(ipwProp);
    planeWidget[i]->TextureInterpolateOff();
    planeWidget[i]->SetResliceInterpolateToLinear();
    planeWidget[i]->SetInputConnection(reader->GetOutputPort());
    planeWidget[i]->SetPlaneOrientation(i);
    planeWidget[i]->SetSliceIndex(imageDims[i] / 2);
    planeWidget[i]->DisplayTextOn();
    planeWidget[i]->SetDefaultRenderer(ren[3]);
    planeWidget[i]->SetWindowLevel(1358, -27);
    planeWidget[i]->On();
    planeWidget[i]->InteractionOn();
}

planeWidget[1]->SetLookupTable(planeWidget[0]->GetLookupTable());
planeWidget[2]->SetLookupTable(planeWidget[0]->GetLookupTable());


vtkSmartPointer<vtkResliceCursorCallback> cbk =
    vtkSmartPointer<vtkResliceCursorCallback>::New();

// Create the reslice cursor, widget and rep

vtkSmartPointer< vtkResliceCursor > resliceCursor =
    vtkSmartPointer< vtkResliceCursor >::New();
resliceCursor->SetCenter(reader->GetOutput()->GetCenter());
resliceCursor->SetThickMode(0);
resliceCursor->SetThickness(10, 10, 10);
resliceCursor->SetImage(reader->GetOutput());

vtkSmartPointer< vtkResliceCursorWidget > resliceCursorWidget[3];
vtkSmartPointer< vtkResliceCursorLineRepresentation > resliceCursorRep[3];

double viewUp[3][3] = { { 0, 0, -1 },{ 0, 0, 1 },{ 0, 1, 0 } };
for (int i = 0; i < 3; i++)
{
    resliceCursorWidget[i] = vtkSmartPointer< vtkResliceCursorWidget >::New();
    resliceCursorWidget[i]->SetInteractor(iren);

    resliceCursorRep[i] =
        vtkSmartPointer< vtkResliceCursorLineRepresentation >::New();
    resliceCursorWidget[i]->SetRepresentation(resliceCursorRep[i]);
    resliceCursorRep[i]->GetResliceCursorActor()->
        GetCursorAlgorithm()->SetResliceCursor(resliceCursor);
    resliceCursorRep[i]->GetResliceCursorActor()->
        GetCursorAlgorithm()->SetReslicePlaneNormal(i);

    const double minVal = reader->GetOutput()->GetScalarRange()[0];
    if (vtkImageReslice *reslice =
        vtkImageReslice::SafeDownCast(resliceCursorRep[i]->GetReslice()))
    {
        reslice->SetBackgroundColor(minVal, minVal, minVal, minVal);
    }

    resliceCursorWidget[i]->SetDefaultRenderer(ren[i]);
    resliceCursorWidget[i]->SetEnabled(1);

    ren[i]->GetActiveCamera()->SetFocalPoint(0, 0, 0);
    double camPos[3] = { 0,0,0 };
    camPos[i] = 1;
    ren[i]->GetActiveCamera()->SetPosition(camPos);

    ren[i]->GetActiveCamera()->ParallelProjectionOn();
    ren[i]->GetActiveCamera()->SetViewUp(
        viewUp[i][0], viewUp[i][1], viewUp[i][2]);
    ren[i]->ResetCamera();
    //ren[i]->ResetCameraClippingRange();

    // Tie the Image plane widget and the reslice cursor widget together
    cbk->IPW[i] = planeWidget[i];
    cbk->RCW[i] = resliceCursorWidget[i];
    resliceCursorWidget[i]->AddObserver(
        vtkResliceCursorWidget::ResliceAxesChangedEvent, cbk);

    // Initialize the window level to a sensible value
    double range[2];
    reader->GetOutput()->GetScalarRange(range);
    resliceCursorRep[i]->SetWindowLevel(range[1] - range[0], (range[0] + range[1]) / 2.0);
    planeWidget[i]->SetWindowLevel(range[1] - range[0], (range[0] + range[1]) / 2.0);

    // Make them all share the same color map.
    resliceCursorRep[i]->SetLookupTable(resliceCursorRep[0]->GetLookupTable());
    planeWidget[i]->GetColorMap()->SetLookupTable(resliceCursorRep[0]->GetLookupTable());
}


// Add the actors
//
ren[0]->SetBackground(0.3, 0.1, 0.1);
ren[1]->SetBackground(0.1, 0.3, 0.1);
ren[2]->SetBackground(0.1, 0.1, 0.3);
ren[3]->AddActor(outlineActor);
ren[3]->SetBackground(0.1, 0.1, 0.1);
renWin->SetSize(600, 600);
//renWin->SetFullScreen(1);

ren[0]->SetViewport(0, 0, 0.5, 0.5);
ren[1]->SetViewport(0.5, 0, 1, 0.5);
ren[2]->SetViewport(0, 0.5, 0.5, 1);
ren[3]->SetViewport(0.5, 0.5, 1, 1);

// Set the actors' positions
//
renWin->Render();

ren[3]->GetActiveCamera()->Elevation(110);
ren[3]->GetActiveCamera()->SetViewUp(0, 0, -1);
ren[3]->GetActiveCamera()->Azimuth(45);
ren[3]->GetActiveCamera()->Dolly(1.15);
ren[3]->ResetCameraClippingRange();

vtkSmartPointer< vtkInteractorStyleImage > style =
    vtkSmartPointer< vtkInteractorStyleImage >::New();
iren->SetInteractorStyle(style);

iren->Initialize();


    iren->Start();


return 0;

}
按照上面的指示,我只是轻微的修改了细微的几个地方而已,
http://www.vtk.org/gitweb?p=VTK.git;a=blob_plain;f=Interaction/Widgets/Testing/Cxx/TestResliceCursorWidget2.cxx
这是官网的代码,可是我的代码却运行不了呢?麻烦大家解答。

Csdn user default icon
上传中...
上传图片
插入图片
抄袭、复制答案,以达到刷声望分或其他目的的行为,在CSDN问答是严格禁止的,一经发现立刻封号。是时候展现真正的技术了!
其他相关推荐
opencv,cmake,vtk相关使用问题
-
opencv求助,安装问题?!
-
opencv怎样实现图像投影
-
opencv 使用问题
-
怎样实现c++利用opencv实现人脸检测与识别
-
opencv安装不成功,求大神们指点!!
-
opencv运动物体追踪的实现
-
如何使用conda更换opencv的版本
-
opencv如何不配置直接使用
-
opencv houghlines出错
-
linux安装opencv是make出错
-
用anaconda安装的opencv能用于c语言编程吗?
-
opencv3.1.0中没有opencv2\legacy\compat.hpp
-
ubuntu18.04安装Opencv3.10的问题?
-
如何在Android中用OPENCV实现人脸检测,人脸识别等功能,如何实现,用什么方法?
-
使用opencv4Android如何实现自动聚焦
-
ARM-Linux环境下安装OpenCV
-
前台c# 算法基于c++的opencv ide:vs2015 如何打包成一个交付产品,在其他电脑上安装部署?
-
如何快速查阅openCV的函数作用和使用方法?
-
记一道字节跳动的算法面试题
点击蓝色“五分钟学算法”关注我哟加个“星标”,天天中午 12:15,一起学算法作者 | 帅地来源公众号 | 苦逼的码农前几天有个朋友去面试字节跳动,面试官问了他一道链表相...
程序员真是太太太太太有趣了!!!
网络上虽然已经有了很多关于程序员的话题,但大部分人对这个群体还是很陌生。我们在谈论程序员的时候,究竟该聊些什么呢?各位程序员大佬们,请让我听到你们的声音!不管你是前端开发...
史上最详细的IDEA优雅整合Maven+SSM框架(详细思路+附带源码)
网上很多整合SSM博客文章并不能让初探ssm的同学思路完全的清晰,可以试着关掉整合教程,摇两下头骨,哈一大口气,就在万事具备的时候,开整,这个时候你可能思路全无 ~中招了咩~ ,还有一些同学依旧在使用eclipse或者Myeclipse开发,我想对这些朋友说IDEA 的编译速度很快,人生苦短,来不及解释了,直接上手idea吧。这篇文章每一步搭建过程都测试过了,应该不会有什么差错。本文章还有个比较优秀的特点,就是idea的使用,基本上关于idea的操作都算是比较详细的,所以不用太担心不会撸idea!最后,本文
吃人的那些 Java 名词:对象、引用、堆、栈
作为一个有着 8 年 Java 编程经验的 IT 老兵,说起来很惭愧,我被 Java 当中的四五个名词一直困扰着:**对象、引用、堆、栈、堆栈**(栈可同堆栈,因此是四个名词,也是五个名词)。每次我看到这几个名词,都隐隐约约觉得自己在被一只无形的大口慢慢地吞噬,只剩下满地的衣服碎屑(为什么不是骨头,因为骨头也好吃)。
LeetCode解题汇总目录
此篇为LeetCode刷题的汇总目录,方便大家查找,一起刷题,一起PK交流! 已解题目 考点 LeetCode 1. 两数之和(哈希) LeetCode 2. 两数相加(单链表反转) LeetCode 9. 回文数 LeetCode 11. 盛最多水的容器(双指针) LeetCode 15. 三数之和 LeetCode 17. 电话号码的字母组合(回溯...
我花了一夜用数据结构给女朋友写个H5走迷宫游戏
起因 又到深夜了,我按照以往在csdn和公众号写着数据结构!这占用了我大量的时间!我的超越妹妹严重缺乏陪伴而 怨气满满! 而女朋友时常埋怨,认为数据结构这么抽象难懂的东西没啥作用,常会问道:天天写这玩意,有啥作用。而我答道:能干事情多了,比如写个迷宫小游戏啥的! 当我码完字准备睡觉时:写不好别睡觉! 分析 如果用数据结构与算法造出东西来呢? ...
不识 Pandas,纵是老手也枉然?
作者 |周志鹏 责编 | 郭 芮 这段时间和一些做数据分析的同学闲聊,我发现数据分析技能入门阶段存在一个普遍性的问题,很多凭着兴趣入坑的同学,都能够很快熟悉Python基础语法,然后不约而同的一头扎进《利用Python进行数据分析》这本经典之中,硬着头皮啃完之后,好像自己什么都会了一点,然而实际操作起来既不知从何操起,又漏洞百出。 至于原因嘛,理解不够,实践不够是两条老牌的拦路...
接班马云的为何是张勇?
上海人、职业经理人、CFO 背景,集齐马云三大不喜欢的张勇怎么就成了阿里接班人? 作者|王琳 本文经授权转载自燃财经(ID:rancaijing) 9月10日,张勇转正了,他由阿里巴巴董事局候任主席正式成为阿里巴巴董事局主席,这也意味着阿里巴巴将正式开启“逍遥子时代”。 从2015年接任CEO开始,张勇已经将阿里巴巴股价拉升了超过200%。但和马云强大的个人光环比,张勇显得尤其...
14 个实用的数据库设计技巧
点击上方“后端技术精选”,选择“置顶公众号”技术文章第一时间送达!作者:echozhjuejin.im/post/5d5b4c6951882569eb570958原始单据...
我在快手认识了 4 位工程师,看到了快速发展的公司和员工如何彼此成就!
作者 | 胡巍巍 出品 | CSDN(ID:CSDNnews) 从西二旗地铁站B口出来,步行700多米可以看到一个工业建筑风格的院子。这个独立的院子和后厂村各大互联网公司的高楼林立有些不同。 院子里有7栋6层高的楼,几栋楼之间打通,可以从A栋自由穿行到F栋。这里就是快手总部。这个园区可以容纳6000多名员工,目前40%以上是研发人员。 这些研发人员维护着快手这款日活超过2亿的ap...
让程序员崩溃的瞬间(非程序员勿入)
今天给大家带来点快乐,程序员才能看懂。 来源:https://zhuanlan.zhihu.com/p/47066521 1. 公司实习生找 Bug 2.在调试时,将断点设置在错误的位置 3.当我有一个很棒的调试想法时 4.偶然间看到自己多年前写的代码 5.当我第一次启动我的单元测试时 ...
用Python分析2000款避孕套,得出这些有趣的结论
到现在为止,我们的淘宝教程已经写到了第四篇,前三篇分别是: 第一篇:Python模拟登录淘宝,详细讲解如何使用requests库登录淘宝pc端。 第二篇:淘宝自动登录2.0,新增Cookies序列化,教大家如何将cookies保存起来。 第三篇:Python爬取淘宝商品避孕套,教大家如何爬取淘宝pc端商品信息。 今天,我们来看看淘宝系列的第四篇 我们在上一篇的时候已经将淘宝数据爬取下来了,...
Spring高级技术梳理
Spring高级技术梳理 序言正文SpringDate部分Spring全家桶之SpringData——预科阶段Spring全家桶之SpringData——Spring 整合Hibernate与Hibernate JpaSpring全家桶之SpringData——Spring Data JPASpring全家桶之SpringData——SpringData RedisSpringBoot部分Sp...
如何在Windows中开启"上帝模式"
原文链接 : https://mp.weixin.qq.com/s?__biz=MzIwMjE1MjMyMw==&amp;mid=2650202982&amp;idx=1&amp;sn=2c6c609ce06db1cee81abf2ba797be1b&amp;chksm=8ee1438ab996ca9c2d0cd0f76426e92faa835beef20ae21b537c0867ec2773be...
Docker 零基础从入门到使用
诺!这只可爱的小鲸鱼就是docker了! Docker 是什么? Docker 是一个开源的应用容器引擎,让开发者可以打包他们的应用以及依赖包到一个可移植的镜像中,然后发布到任何流行的 Linux 或 Windows 机器上( 摘自百度 )。 Docker 能干什么? 在讲 Docker 能干什么之前,我们不妨先看看没有 Docker 和有Docker分别是个什么样子的? 场景一 某公司需要开发...
再见 Docker,是时候拥抱下一代容器工具了
什么是 Linux 容器?Linux 容器是由 Linux 内核所提供的具有特定隔离功能的进程,Linux 容器技术能够让你对应用及其整个运行时环境(包括全部所需文件)一...
不足20行 python 代码,高效实现 k-means 均值聚类算法
关于 k-means 均值聚类算法的原理介绍、实现代码,网上有很多,但运行效率似乎都有点问题。今天稍微有点空闲,写了一个不足20行的 k-means 均值聚类算法,1万个样本平均耗时20毫秒(10次均值)。同样的数据样本,网上流行的算法平均耗时3000毫秒(10次均值)。差距竟然达百倍以上,令我深感意外,不由得再次向 numpy 献上膝盖!
分享靠写代码赚钱的一些门路
作者 mezod,译者 josephchang10如今,通过自己的代码去赚钱变得越来越简单,不过对很多人来说依然还是很难,因为他们不知道有哪些门路。今天给大家分享一个精彩...
北漂程序员,扬帆起航的地方
随着耳畔传来“你看这碗又大又圆、你看这面又长又宽......碗大宽无影、像儿时的回忆......”听着挺带劲,于是看了一下手机,原来是吴亦凡的作品《大碗宽面》,随着入耳的旋律,脑子也不由自主的想起 10 年前,在平西府吃 5 块钱一大碗牛肉板面的情景。 平西府最有名的就是这个牌坊啦。记得每当有同事问起住哪里?都会自豪的说住在王府里;隔三差五也会邀请朋友去府上坐坐。其实打内心里讲,平西府是一个...
技术人员要拿百万年薪,必须要经历这9个段位
很多人都问,技术人员如何成长,每个阶段又是怎样的,如何才能走出当前的迷茫,实现自我的突破。所以我结合我自己10多年的从业经验,总结了技术人员成长的9个段位,希望对大家的职...
多线程编程是后台开发人员的基本功
这里先给大家分享一个小故事:在我刚开始参加工作的那年,公司安排我开发一款即时通讯软件(IM,类似于 QQ 聊天软件),在这之前我心里也知道如果多线程操作一个整型值是要加锁...
相关热词 c#panel增加滚动条 c#中生成的dll文件 c# 模板类 c# 截取txt文本内容 c# 内存 占用 c#时间格式化 不带- c#替换字符串中指定位置 c# rdlc 动态报表 c# 获取txt编码格式 c#事件主动调用