Я хотел бы лучше понять объемный рендеринг 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, что мне нужно было активировать режим «Объем», чтобы увидеть мою модель, здесь есть что-то похожее?
Последнее, что очень важно: можно ли изменить 3D-объем внутри файла .vtk? Например, если я хочу изменить цвет определенной части модели, предоставляет ли VTK инструменты, позволяющие это сделать?