QT C++ vtk 加载二维dicom 三维渲染 3DMPR

最新推荐文章于 2024-06-02 09:45:57 发布
最新推荐文章于 2024-06-02 09:45:57 发布
阅读量1.8k 收藏 9
点赞数
分类专栏： tianya.cn/blog-5287971-1.shtml blog.sina.com.cn/u/2090406854 文章标签： servlet jar java
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/ermu114/article/details/127193229
版权

 

 

 

 

 

 

 

 

 

 

#include "widget3d.h"
#include "tabwidget.h"
#include <vtkRendererCollection.h>
#include <vtkGenericOpenGLRenderWindow.h>
#include <QKeyEvent>
#include <QtConcurrent/qtconcurrentrun.h>

asclepios::gui::Widget3D::Widget3D(QWidget* parent)
    : WidgetBase(parent)
{
    initData();
    initView();
    createConnections();
    m_tabWidget = parent;
}

//-----------------------------------------------------------------------------
void asclepios::gui::Widget3D::render()
{
    if (!m_image)
    {
        throw std::runtime_error("No image for widget 3d!");
    }
    try
    {
        m_toolbar->getUI().toolButtonCrop->setVisible(false);
        m_toolbar->getUI().comboBoxFilters->setVisible(false);
        startLoadingAnimation();
        m_vtkWidget->setImage(m_image);
        m_vtkWidget->setSeries(m_series);
        m_vtkWidget->setInteractor(m_qtvtkWidget->GetInteractor());
        m_future = QtConcurrent::run(onRenderAsync, this);
        Q_UNUSED(connect(this, &Widget3D::finishedRenderAsync,
            this, &Widget3D::onFinishedRenderAsync));
    }
    catch (const std::exception& ex)
    {
        //todo log
    }
}

//-----------------------------------------------------------------------------
bool asclepios::gui::Widget3D::eventFilter(QObject* watched, QEvent* event)
{
    if (event->type() == QEvent::KeyPress)
    {
        auto* const keyEvent = dynamic_cast<QKeyEvent*>(event);
        const int key = keyEvent->key();
        auto* const combo = m_toolbar->getUI().comboBoxFilters;
        switch (key)
        {
        case Qt::Key_Left:
        {
            const int currentIndex = combo->currentIndex();
            combo->setCurrentIndex(!combo->currentIndex()
                ? combo->count() - 1
                : currentIndex - 1);
            break;
        }
        case Qt::Key_Right:
        {
            const int currentIndex = combo->currentIndex();
            combo->setCurrentIndex(currentIndex == combo->count() - 1
                ? 0
                : currentIndex + 1);
            break;
        }
        default:
            break;
        }
    }
    return QWidget::eventFilter(watched, event);
}

//-----------------------------------------------------------------------------
void asclepios::gui::Widget3D::onfilterChanged(const QString& t_filter) const
{
    if (m_qtvtkWidget && m_vtkWidget)
    {
        m_vtkWidget->setFilter(t_filter);
        m_qtvtkWidget->GetRenderWindow()->Render();
    }
}

//-----------------------------------------------------------------------------
void asclepios::gui::Widget3D::onCropPressed(const bool& t_pressed) const
{
    m_vtkWidget->activateBoxWidget(t_pressed);
    m_qtvtkWidget->GetRenderWindow()->Render();
}

//-----------------------------------------------------------------------------
void asclepios::gui::Widget3D::onActivateWidget(const bool& t_flag)
{
    if (t_flag)
    {
        auto* event = new QFocusEvent(QEvent::FocusIn,
            Qt::FocusReason::MouseFocusReason);
        focusInEvent(event);
        delete event;
    }
}

//-----------------------------------------------------------------------------
void asclepios::gui::Widget3D::onSetMaximized() const
{
    if (m_tabWidget)
    {
        dynamic_cast<TabWidget*>
            (m_tabWidget)->onMaximize();
    }
}

//-----------------------------------------------------------------------------
void asclepios::gui::Widget3D::onFinishedRenderAsync()
{
    auto* const renderWindow =
        m_qtvtkWidget->GetRenderWindow();
    renderWindow->AddRenderer(m_vtkWidget->
        getRenderWindows()[0]->GetRenderers()->
        GetFirstRenderer());
    renderWindow->Render();
    onfilterChanged(m_toolbar->getUI()
        .comboBoxFilters->itemData(0).toString());
    stopLoadingAnimation();
    disconnect(this, &Widget3D::finishedRenderAsync,
        this, &Widget3D::onFinishedRenderAsync);
    m_toolbar->getUI().toolButtonCrop->setVisible(true);
    m_toolbar->getUI().comboBoxFilters->setVisible(true);
    installEventFilter(this);
}

//-----------------------------------------------------------------------------
void asclepios::gui::Widget3D::onRenderAsync(Widget3D* t_self)
{
    t_self->m_vtkWidget->render();
    emit t_self->finishedRenderAsync();
}

//-----------------------------------------------------------------------------
void asclepios::gui::Widget3D::initView()
{
    m_ui.setupUi(this);
    setAttribute(Qt::WA_DeleteOnClose);
    auto* const layout = new QVBoxLayout(this);
    layout->setSpacing(0);
    layout->setMargin(0);
    layout->addWidget(m_toolbar);
    layout->addWidget(m_qtvtkWidget);
    setLayout(layout);
}

//-----------------------------------------------------------------------------
void asclepios::gui::Widget3D::initData()
{
    m_qtvtkWidget = new QVTKOpenGLNativeWidget(this);
    m_qtvtkWidget->SetRenderWindow(vtkNew<vtkGenericOpenGLRenderWindow>());
    m_qtvtkWidget->GetRenderWindow()->SetDoubleBuffer(true);
    m_vtkWidget = std::make_unique<vtkWidget3D>();
    m_toolbar = new ToolbarWidget3D(this);
    m_vtkEvents = std::make_unique<vtkEventFilter>(this);
    setWidgetType(WidgetType::widget3d);
}

//-----------------------------------------------------------------------------
void asclepios::gui::Widget3D::createConnections()
{
    if (m_toolbar)
    {
        Q_UNUSED(connect(m_toolbar,
            &ToolbarWidget3D::filterChanged, this,
            &Widget3D::onfilterChanged));
        Q_UNUSED(connect(m_toolbar, &ToolbarWidget3D::cropPressed,
            this, &Widget3D::onCropPressed));
    }
    setFocusPolicy(Qt::FocusPolicy::WheelFocus);
    m_qtvtkWidget->installEventFilter(m_vtkEvents.get());
    Q_UNUSED(connect(m_vtkEvents.get(),
        &vtkEventFilter::activateWidget,
        this, &Widget3D::onActivateWidget));
    Q_UNUSED(connect(m_vtkEvents.get(),
        &vtkEventFilter::setMaximized,
        this, &Widget3D::onSetMaximized));
}

//-----------------------------------------------------------------------------
void asclepios::gui::Widget3D::startLoadingAnimation()
{
    m_loadingAnimation = std::make_unique<LoadingAnimation>(this);
    m_loadingAnimation->setWindowFlags(Qt::Widget);
    layout()->addWidget(m_loadingAnimation.get());
    m_loadingAnimation->show();
}

#include "vtkwidget3d.h"
#include <vtkRenderer.h>
#include <vtkDICOMMetaData.h>
#include <vtkWidgetRepresentation.h>
#include <vtkVolumeProperty.h>
#include <vtkRenderWindowInteractor.h>
#include <vtkContourValues.h>
#include <QString>

void asclepios::gui::vtkWidget3D::initWidget()
{
    m_renderWindows[0] = vtkSmartPointer<vtkRenderWindow>::New();
    m_mapper = vtkSmartPointer<vtkSmartVolumeMapper>::New();
    m_renderer = vtkSmartPointer<vtkRenderer>::New();
    m_volume = vtkSmartPointer<vtkVolume>::New();
    m_transferFunction = std::make_unique<TransferFunction>();
}

//-----------------------------------------------------------------------------
void asclepios::gui::vtkWidget3D::initBoxWidget()
{
    if (!m_renderer || !m_volume)
    {
        return;
    }
    m_boxWidget = vtkSmartPointer<vtkBoxWidget2>::New();
    m_boxWidget->SetInteractor(m_renderer->GetRenderWindow()->GetInteractor());
    m_boxWidget->CreateDefaultRepresentation();
    m_boxWidget->GetRepresentation()->SetPlaceFactor(1);
    m_boxWidget->GetRepresentation()->PlaceWidget(m_volume->GetBounds());
    initBoxWidgetCallback();
}

//-----------------------------------------------------------------------------
void asclepios::gui::vtkWidget3D::initBoxWidgetCallback()
{
    m_boxWidgetCallback = vtkSmartPointer<vtkBoxWidget3DCallback>::New();
    m_boxWidgetCallback->setVolume(m_volume);
    m_boxWidget->AddObserver(vtkCommand::InteractionEvent, m_boxWidgetCallback);
}

//-----------------------------------------------------------------------------
void asclepios::gui::vtkWidget3D::initInteractorStyle()
{
    m_interactorStyle =
        vtkSmartPointer<vtkWidget3DInteractorStyle>::New();
    m_interactorStyle->setWidget(this);
    m_interactorStyle->setTransferFunction(m_transferFunction.get());
    m_interactor->SetInteractorStyle(m_interactorStyle);
}

//-----------------------------------------------------------------------------
void asclepios::gui::vtkWidget3D::setVolumeMapperBlend() const
{
    m_mapper->SetBlendMode(vtkVolumeMapper::COMPOSITE_BLEND);
}

//-----------------------------------------------------------------------------
std::tuple<int, int> asclepios::gui::vtkWidget3D::getWindowLevel() const
{
    const auto imageReader =
        m_image->getImageReader();
    return std::make_tuple<int, int>(imageReader->
                                     GetMetaData()->Get(DC::WindowCenter).AsInt(),
                                     imageReader->GetMetaData()->Get(DC::WindowWidth).AsInt());
}

//-----------------------------------------------------------------------------
void asclepios::gui::vtkWidget3D::setFilter(const QString& t_filePath) 
{
    try
    {
        if (t_filePath == "MIP")
        {
            m_transferFunction.reset();
            m_transferFunction = std::make_unique<TransferFunction>();
            m_mapper->SetBlendMode(vtkVolumeMapper::MAXIMUM_INTENSITY_BLEND);
            m_transferFunction->setMaximumIntensityProjectionFunction(0, 0);
            const auto [window, level] = getWindowLevel();
            m_transferFunction->updateWindowLevel(window, level);
        }
        else
        {
            m_mapper->SetBlendMode(vtkVolumeMapper::COMPOSITE_BLEND);
            m_transferFunction->loadFilterFromFile(t_filePath);
        }
        updateFilter();
        m_interactorStyle->setTransferFunction(m_transferFunction.get());
        
    }
    catch (const std::exception& ex)
    {
        //todo log
    }
}

//-----------------------------------------------------------------------------
void asclepios::gui::vtkWidget3D::render()
{
    m_renderWindows[0]->OffScreenRenderingOn();
    setVolumeMapperBlend();
    const auto [window, level] = getWindowLevel();
    const auto reader = m_image && m_image->getIsMultiFrame()
                            ? m_image->getImageReader()
                            : m_series->getReaderForAllSingleFrameImages();
    m_mapper->SetInputConnection(reader->GetOutputPort());
    m_transferFunction->updateWindowLevel(window, level);
    m_volume->SetMapper(m_mapper);
    m_renderer->AddActor(m_volume);
    m_renderWindows[0]->AddRenderer(m_renderer);
    m_renderWindows[0]->Render();
    m_renderWindows[0]->OffScreenRenderingOff();
    auto* const extend = m_volume->GetBounds();
    m_volume->SetOrigin(extend[0] + (extend[1] - extend[0]) / 2,
                        extend[2] + (extend[3] - extend[2]) / 2, 0);
    initInteractorStyle();
}

//-----------------------------------------------------------------------------
void asclepios::gui::vtkWidget3D::activateBoxWidget(const bool& t_flag)
{
    if (!m_boxWidget)
    {
        initBoxWidget();
    }
    m_boxWidget->SetEnabled(t_flag);
}

//-----------------------------------------------------------------------------
void asclepios::gui::vtkWidget3D::updateFilter() const
{
    m_volume->GetProperty()->SetInterpolationTypeToLinear();
    m_volume->GetProperty()->SetScalarOpacity(m_transferFunction->getOpacityFunction());
    m_volume->GetProperty()->SetAmbient(m_transferFunction->getAmbient());
    m_volume->GetProperty()->SetColor(m_transferFunction->getColorFunction());
    m_volume->GetProperty()->SetDiffuse(m_transferFunction->getDiffuse());
    m_volume->GetProperty()->SetSpecular(m_transferFunction->getSpecular());
    m_volume->GetProperty()->SetSpecularPower(m_transferFunction->getSpecularPower());
    m_volume->GetProperty()->GetIsoSurfaceValues()->SetValue(0, 0);
    (m_transferFunction->getHasShade()) ? m_volume->GetProperty()->ShadeOn() : m_volume->GetProperty()->ShadeOff();
}
{
  "name":"Dark Bone",
  "color": [
    {
      "value": 46.733612060546875,
      "red": 0.0,
      "green": 0.0,
      "blue": 0.0
    },
    {
      "value": 134.97621154785156,
      "red": 1,
      "green": 1,
      "blue": 1
    },
    {
      "value": 244.72689819335938,
      "red": 0.20000000298023224,
      "green": 0.20000000298023224,
      "blue": 0.20000000298023224
    },
    {
      "value": -812.04962158203125,
      "red": 0.0,
      "green": 1,
      "blue": 1
    },
    {
      "value": -622.0498046875,
      "red": 0.0,
      "green": 1,
      "blue": 1
    },
    {
      "value": -420.04998779296875,
      "red": 0.0,
      "green": 1,
      "blue":  1
    },
    {
      "value": -262.84738159179688,
      "red": 0,
      "green": 1,
      "blue": 1
    }
  ],
  "opacity": [
    {
      "value": 46.733612060546875,
      "alpha": 0.0
    },
    {
      "value": 134.97621154785156,
      "alpha": 0.25999999046325684
    },
    {
      "value": 244.72689819335938,
      "alpha": 0.5300024151802063
    },
    {
      "value": -812.04962158203125,
      "alpha": 0.0
    },
    {
      "value": -622.0498046875,
      "alpha": 0.1643165796995163
    },
    {
      "value": -420.04998779296875,
      "alpha": 0.36469146609306335
    },
    {
      "value": -262.84738159179688,
      "alpha": 0.0
    }
  ],
  "ambient":
    {
      "value": 0.4
    },
  "diffuse":
    {
      "value": 0.7
    },
  "specular":
    {
      "value": 1.0
    },
  "specularpower":
    {
      "value": 64
    },
  "shade":
    {
      "value": 1
    }
}

灿烂李
关注
  • 0
    点赞
  • 踩
  • 9
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 5
    评论
专栏目录
读取并可视化DICOM图像（VTK & OpenCV） - louishao的博客 - CSDN博客1
08-04
本人处理的CT图的格式是DICOM格式，这个用OpenCV不能直接读取，我采用了VTK进行读取，然后用OpenCV处理的思路（当然网上也有人说可以直接用ITK读
基于VTK和QT实现DICOM图像的像素擦除
l1783111653的博客
10-13 3361
基于VTK和QT实现DICOM图像的像素擦除
推荐：Asclepios DICOM Viewer——医疗影像的卓越视图解析工具
最新发布
gitblog_00065的博客
06-02 314
推荐：Asclepios DICOM Viewer——医疗影像的卓越视图解析工具 项目地址:https://gitcode.com/GavriloviciEduard/Asclepios-DICOM-Viewer 在医学领域，高质量的图像处理和解析是至关重要的。今天，我们要向您推荐一款名为Asclepios DICOM Viewer的开源软件，它是一个强大的C++开发的DICOM（Digital ...
VTK：利用qt实现体绘制剪裁
Beyond欣
07-29 4958
VTK：利用qt实现体绘制剪裁 做体绘制渲染一定少不了交互剪裁的功能。一般又两种方式： 搞一个包围盒，移动包围盒实现体渲染的剪裁 像小蚂蚁一样可以在屏幕画任意多边形实现剪裁 结果跟狗啃的一样主要是因为数据尺寸比较小（像素40* 40 * 40），直接抹黑像素肯定很烂 原文 VTK：利用qt实现体绘制剪裁 1 项目地址 在官方案例基础上改的 包围盒案例 https://kitware.github.io/vtk-examples/site/Cxx/Widgets/BoxWidget2/ 体渲染
VTK可视化管线之实现观察者/命令（Observer/Command）设计模式
kateyabc的博客
05-12 1394
VTK可视化管线之实现观察者/命令（Observer/Command）设计模式#include"vtkSmartPointer.h"#include"vtkConeSource.h"#include"vtkPolyDataMapper.h"#include"vtkRenderWindow.h"#include"vtkRenderWindowInteractor.h"#include"vtkCame.
6.VTK二维三维图像显示以及二维三维图像切面交互提取（回调函数、观察者-命令模式）
weixin_52194015的博客
10-09 773
遇到的一个大问题就是由于版本更新，这些资料中很多代码无法正常运行，需要进行一定的修改，所以这个专栏会记录下来我修改后的程序代码，以便于我之后温习。（3）横断面（Transverse/Axial Plane）：平行于XY平面的切面。除了以上三种特殊的切面，切面也可以是通过三维图像内部一点且平行于任意方向的平面。（1）矢状面（Sagital Plane）：平行于YZ平面的切面。（2）冠状面（Coronal Plane）：平行于XZ平面的切面。VTK中vtkImageReslice类可实现图像切面的提取。
【vtk】鼠标拖动移动Actor
weixin_41838721的博客
05-08 2169
https://blog.csdn.net/weixin_42834523/article/details/104844660 #include <vtkVersion.h> #include <vtkSmartPointer.h> #include <vtkActor.h> #include <vtkSphereSource.h> #include <vtkRendererCollection.h> #include <vtkCellArr
3D_VTK.rar_DICOM 3D_VTK 三维可视化_VTK三维视图_qt dicom_qt vtk
07-15
这是一个基于QT和VTK开发的一个三维可视化界面，可以用面绘制实现序列DICOM的三维重建。
基于Qt和VTK实现CT三维重建项目c++源码+项目说明.zip
07-24
基于Qt和VTK实现CT三维重建项目c++源码+项目说明.zip 该项目是个人毕设项目，答辩评审分达到95分，代码都经过调试测试，确保可以运行！欢迎下载使用，可用于小白学习、进阶。 该资源主要针对计算机、通信、人工智能...
VTK读取30张Dicom图片三维重建
04-20
利用VTK三维重建软件与VS2010结合，读取格式为Dicom的系列图像，并能显示三维重建结果，语言为C++语言(Using VTK 3D reconstruction software and VS2010 combination, read the format of Dicom series of images, ...
VTK.js dicom图像三维重建
04-19
VTK.js 实现网页版dicom图像三维重建
编译好的VTK+Qt5可视化库（VTK9.0）
05-08
编译好的Vtk9.0三维可视化库，c++版本，是与QT5 绑定的，可以直接使用，包括debug版本的库和release版本的库。
利用VTK的DICOM图像的显示与信息读取
11-16
利用VTK的DICOM图像的显示与信息读取
医疗三维重建技术c++,QT
02-26
总的来说，"医疗三维重建技术c++,QT"项目结合了C++的性能、QT的界面设计和VTK的可视化能力，为医疗专业人员提供了一个强大的工具，帮助他们更好地理解和分析DICOM影像中的复杂结构，从而提高诊断和治疗的精度。...
QT5.10.1+VTK8.2.0+PCL1.8.1库显示点云的编译错误
u014048569的博客
04-16 1950
编译错误1：QT5.10.1+VTK8.2.0+PCL1.8.1库显示点云，报错：“QVTKOpenGLNativeWidget” requires a 'vtkGenericOpenGLRenderWindow',‘vtkWin32OpenGLRenderWindow’is not supported. 解决方案1：1.首先在用Make欲编译PCL时，要勾选vtkGUISupportQtOpe...
【C++】将Qt界面显示到VTK三维环境中
愿岁并谢，与长友兮
11-03 2479
开发环境：win10，vs2017，VTK 8.2.0，Qt 5.11.0 项目的需求是在三维环境下实现UI，目前实现的效果如下： 左边是我的Qt窗口在界面设计器中的预览，右边是VTK三维环境下界面的显示。 （为了体现三维环境，我在右下角加了一个圆柱） 读到这里，首先你应该具备以下能力才能进行之后的步骤： 1.能将VTK和Qt编译到同一个目录下； 2.熟悉Qt和vt...
从DICOM开始入门医学图像处理-2-VTK+Qt+VS开发环境搭建
安娜的修行
11-21 3321
软件版本 安装Qt,注意选择与你安装的VS版本对应的组件 目录结构 打开CMake,如下设置 点击configure: 选中VTK节点中的VTK_Group_Qt, 点击Configure Ungrouped Entries节点,确认Qt_QMAKE_EXECUTABLE 的路径,将VTK_QT_VERSION改成5. 点击Configure Ungrouped Entries节点,确认Qt的
【vtk】renderer如何删除已添加的vtkactor
weixin_43464938的博客
06-17 2612
我在某个函数里给renderer->addactor(一个actor)。 因为renderer可以加多个actor，这里我就没有将actor写成类成员变量，只是函数里临时new出来的。 后来我想要remove掉这个actor时傻眼了。因为获取不到actor的指针了。 解决思路： actor有个类似list的容器叫actorcollection。我每次将actor加入renderer的时候，actorcollection也就additem这个对象。然后记下actor的编号到QStri..
基于VTK与Qt的体绘制程序
热门推荐
海的那边
06-10 2万+
1基于VTK与Qt的体绘制程序 1.1体绘制介绍 体绘制（Volume Rendering）是一种从离散三维数据场获得二维投影图像的技术。体绘制在医学成像和科学可视化当中有着广泛应用。用于体绘制的数据集也称为体数据（Volume Data）。典型的三维体数据是通过 CT、MRI 等成像技术采集的一组二维切片图像。通常这些图像都是按照某种规则（如每毫米一个切面）采集并重
如何用C++，QT结合vtk绘制二维的矢量箭头
04-01
以下是一个使用C++和QT结合vtk绘制二维矢量箭头的例子： ```cpp #include #include #include #include <vtkActor.h> #include <vtkRenderer.h> #include <vtkRenderWindow.h> #include #include #include #...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 5
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

灿烂李

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

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

余额充值