Question

我想更好地了解dicom体积渲染。

我有一组dicom图像，我可以从中提取轴向，冠状和矢状切口，如下所示：

我首先想从头开始生成3D模型，但这似乎太难了。

所以我听说过VTK / ITK，我一直在使用这段代码从我的图像集生成一个.vtk文件：

http://www.itk.org/Doxygen46/html/IO_2DicomSeriesReadImageWrite2_8cxx-example.html

它有效，但我需要一些解释：

当我使用ParaView打开此文件时，我得到以下结果：

首先，这可能是一个愚蠢的问题，但为什么它是蓝色的？

有没有办法切割和看到模型的内部？

我的目标不是使用ParaView，而且我想制作我自己的.vtk读卡器，我发现这段代码我不记得哪里，我认为它可以起作用，但我得到的只是绿色背景，仅此而已：

#include <vtkPolyDataReader.h>
#include <vtkSmartPointer.h>
#include <vtkPolyDataMapper.h>
#include <vtkActor.h>
#include <vtkRenderWindow.h>
#include <vtkRenderer.h>
#include <vtkRenderWindowInteractor.h>

int main ( int argc, char *argv[] ) {

  // Parse command line arguments                                                                     
  if (argc != 2) {
    std::cerr << "Usage: " << argv[0] << " Filename(.vtk)" << std::endl;
    return EXIT_FAILURE;
  }

  std::string filename = argv[1];

  // Read all the data from the file                                                                  
  vtkSmartPointer<vtkPolyDataReader> reader = vtkSmartPointer<vtkPolyDataReader>::New();
  reader->SetFileName(filename.c_str());
  reader->Update();

  // Visualize                                                                                        
  vtkSmartPointer<vtkPolyDataMapper> mapper = vtkSmartPointer<vtkPolyDataMapper>::New();
  mapper->SetInputConnection(reader->GetOutputPort());

  vtkSmartPointer<vtkActor> actor = vtkSmartPointer<vtkActor>::New();
  actor->SetMapper(mapper);

  vtkSmartPointer<vtkRenderer> renderer = vtkSmartPointer<vtkRenderer>::New();
  vtkSmartPointer<vtkRenderWindow> renderWindow = vtkSmartPointer<vtkRenderWindow>::New();
  renderWindow->AddRenderer(renderer);

  vtkSmartPointer<vtkRenderWindowInteractor> renderWindowInteractor =
    vtkSmartPointer<vtkRenderWindowInteractor>::New();
  renderWindowInteractor->SetRenderWindow(renderWindow);

  renderer->AddActor(actor);
  renderer->SetBackground(.3, .6, .3); // Background color green                                      

  renderWindow->Render();
  renderWindowInteractor->Start();

  return EXIT_SUCCESS;
}

任何想法为什么？我在ParaView中看到我必须激活“音量”模式才能看到我的模型，这里有类似处理吗？

最后一件事，非常重要：是否可以修改.vtk文件中的3D卷？例如，如果我想改变模型特定部分的颜色，VTK是否提供了允许它的工具？

Answer 1

这里有很多问题！以下是一些答案。

渲染为蓝色，因为查询表是paraview中的默认值（从蓝色到红色）。您可以使用color map editor
确实有某种方法来削减＆＃34;在数据内部，请查看之后发布的示例。
该示例不起作用，因为您尝试使用polydata（网格）阅读器加载vtkImageData
这正是你必须选择＆＃34; Volume＆＃34;在Paraview中，这是因为你的数据是一个体积（体素的三维数组）。您可以使用VTK中提供的任何卷映射器执行相同的操作。
是的，你可以编辑卷的价值，但这会让我们现在有点走远，如果需要，我们会在稍后讨论;）

这是一个完整的示例，它读取目录下的所有DICOM文件，构建卷，使用卷渲染渲染它，并使框小部件能够以交互方式剪切卷。

#include "vtkBoxRepresentation.h"
#include "vtkBoxWidget2.h"
#include "vtkCamera.h"
#include "vtkColorTransferFunction.h"
#include "vtkCommand.h"
#include "vtkDICOMImageReader.h"
#include "vtkGPUVolumeRaycastMapper.h"
#include "vtkImageData.h"
#include "vtkInteractorStyle.h"
#include "vtkInteractorStyleTrackballCamera.h"
#include "vtkMath.h"
#include "vtkPiecewiseFunction.h"
#include "vtkPlanes.h"
#include "vtkRenderWindow.h"
#include "vtkRenderWindowInteractor.h"
#include "vtkRenderer.h"
#include "vtkVolume.h"
#include "vtkVolumeProperty.h"

// Box interaction callback
class vtkBoxCallback : public vtkCommand
{
public:
    static vtkBoxCallback *New(){ return new vtkBoxCallback; }
    vtkGPUVolumeRayCastMapper* m_mapper;
    vtkPlanes* m_planes;

    virtual void Execute( vtkObject* a_caller, unsigned long, void* ){
        vtkBoxWidget2* l_box_wdget = vtkBoxWidget2::SafeDownCast( a_caller );
        ( (vtkBoxRepresentation*)l_box_wdget->GetRepresentation() )->GetPlanes( m_planes );
        this->m_mapper->SetClippingPlanes( m_planes );
    }

    vtkBoxCallback(){}
};



int main( int argc, char *argv[] ){

    // Read volume
    vtkDICOMImageReader* l_reader = vtkDICOMImageReader::New();
    l_reader->SetDirectoryName( "C:/PathToDicomFiles/" );
    l_reader->Update();

    // Setup rendering stuff
    vtkRenderer* l_renderer = vtkRenderer::New();
    l_renderer->SetBackground( 0.3, 0.3, 0.3 );

    vtkRenderWindow* l_render_windows = vtkRenderWindow::New();
    l_render_windows->AddRenderer( l_renderer );
    l_render_windows->SetSize( 900, 900 );

    vtkInteractorStyleTrackballCamera* l_trackball = vtkInteractorStyleTrackballCamera::New();

    vtkRenderWindowInteractor* l_iren = vtkRenderWindowInteractor::New();
    l_iren->SetInteractorStyle( l_trackball );
    l_iren->SetRenderWindow( l_render_windows );
    l_iren->GetInteractorStyle()->SetDefaultRenderer( l_renderer );
    l_iren->SetDesiredUpdateRate( 15 );

    // Make sure we have an opengl context
    l_render_windows->Render();


    // Setup GPU volume raycast mapper
    vtkGPUVolumeRayCastMapper* l_gpu_mapper = vtkGPUVolumeRayCastMapper::New();
    l_gpu_mapper->SetInputConnection( l_reader->GetOutputPort() );

    // Setup Volume property
    // Window/Level
    double wl = 260;
    double ww = 270;

    // Color function
    vtkColorTransferFunction* l_color = vtkColorTransferFunction::New();
    l_color->SetColorSpaceToRGB();
    l_color->AddRGBPoint( wl - ww / 2, 0, 0, 0 );
    l_color->AddRGBPoint( wl - ww / 2 + 94 * ( ww / 255.0 ), 1., 21. / 255.0, 27. / 255.0 );
    l_color->AddRGBPoint( wl - ww / 2 + 147 * ( ww / 255.0 ), 1., 176. / 255.0, 9. / 255.0 );
    l_color->AddRGBPoint( wl - ww / 2 + 201 * ( ww / 255.0 ), 1., 241. / 255.0, 39. / 255.0 );
    l_color->AddRGBPoint( wl - ww / 2 + 255 * ( ww / 255.0 ), 1, 1, 1. );
    l_color->Build();

    // Opacity function
    vtkPiecewiseFunction* l_opacity = vtkPiecewiseFunction::New();
    l_opacity->AddPoint( wl - ww / 2, 0 );
    l_opacity->AddPoint( wl + ww / 2, 1 );

    // Volume property, light, shading
    vtkVolumeProperty* l_volume_property = vtkVolumeProperty::New();
    l_volume_property->SetColor( l_color );
    l_volume_property->SetScalarOpacity( l_opacity );
    l_volume_property->SetInterpolationTypeToLinear();
    l_volume_property->ShadeOn();
    l_volume_property->SetAmbient( 0.15 );
    l_volume_property->SetDiffuse( 0.8 );
    l_volume_property->SetSpecular( 0.25 );
    l_volume_property->SetSpecularPower( 40 );

    // Put everything together
    vtkVolume* l_volume = vtkVolume::New();
    l_volume->SetProperty( l_volume_property );
    l_volume->SetMapper( l_gpu_mapper );
    l_renderer->AddVolume( l_volume );
    l_renderer->ResetCamera();

    // setup Box interactive widget
    vtkBoxRepresentation* l_box_rep = vtkBoxRepresentation::New();
    l_box_rep->SetInsideOut( true );

    vtkBoxWidget2* l_voi_widget = vtkBoxWidget2::New();
    l_voi_widget->SetRepresentation( l_box_rep );
    l_voi_widget->SetInteractor( l_iren );
    l_voi_widget->GetRepresentation()->SetPlaceFactor( 1. );
    l_voi_widget->GetRepresentation()->PlaceWidget( l_reader->GetOutput()->GetBounds() );
    l_voi_widget->SetEnabled( true );

    vtkPlanes* l_planes = vtkPlanes::New();

    vtkBoxCallback* l_callback = vtkBoxCallback::New();
    l_callback->m_mapper = l_gpu_mapper;
    l_callback->m_planes = l_planes;
    l_voi_widget->AddObserver( vtkCommand::InteractionEvent, l_callback );

    // Go rendering !
    l_iren->Start();

    // Memory cleanup
    l_reader->Delete();
    l_renderer->Delete();
    l_render_windows->Delete();
    l_trackball->Delete();
    l_iren->Delete();
    l_gpu_mapper->Delete();
    l_color->Delete();
    l_opacity->Delete();
    l_volume_property->Delete();
    l_volume->Delete();
    l_voi_widget->Delete();
    l_planes->Delete();
    l_callback->Delete();
}

作为一般性建议，我建议您阅读VTK examples，它可以帮助您了解所有VTK功能。

希望有所帮助：）

从.dcm文件生成和读取.vtk文件，以及修改模型

1 个答案: