[vtkusers] vtkgdcmimagereader for volume rendering

[Taha Ceritli]

Hi,

I want to read a series of 2D dicom images. Then, I'll display it as a 3D
volume. For that, I followed this steps:

1- ipp sorter
2- setfilenames
3- vtkimagechangeinformation

After these, I tried to visualize them using vtkimageviewer2.
I'm sending almost all the code I wrote. I'm pretty sure it loads and sorts
the images. But it gives an error like

Sorting succeeded:0.5

Sorting succeeded:1.5

The program has unexpectedly finished.



vtkGDCMImageReader *dicomReader =
vtkGDCMImageReader::New();

std::vector<std::string> filenames;
const char *filename = "/home/telemed/Downloads/Dentascan";

gdcm::Directory d;
d.Load( filename, false);

filenames = d.GetFilenames();

gdcm::IPPSorter s;
s.SetComputeZSpacing( true );
s.SetZSpacingTolerance( 1e-3 );
bool b = s.Sort( filenames );
if( !b )
{
std::cerr << "Failed to sort:" << "s" << std::endl;
return 1;
}
std::cout << "Sorting succeeded:" << s.GetZSpacing() << std::endl;
//s.Print( std::cout );
const std::vector<std::string> & sorted = s.GetFilenames();

vtkStringArray *files =
vtkStringArray::New();
std::vector< std::string >::const_iterator it = sorted.begin();
for( ; it != sorted.end(); ++it)
{
const std::string &f = *it;

files->InsertNextValue( f.c_str() );
}
dicomReader->SetFileNames( files );
//dicomReader->Update();

vtkImageChangeInformation *changer =
vtkImageChangeInformation::New();
changer->SetInputConnection(dicomReader->GetOutputPort());
changer->SetOutputSpacing(VTK_DOUBLE_MAX,VTK_DOUBLE_MAX, 2.0);
changer->Update();


vtkInteractorStyleImage *imageStyle =
vtkInteractorStyleImage::New();
vtkRenderWindowInteractor *interactor =
vtkRenderWindowInteractor::New();
// Calculate the center of the volume

dicomReader->GetOutput()->UpdateInformation();
int extent[6];
double spacing[3];
double origin[3];
dicomReader->GetOutput()->GetWholeExtent(extent);
dicomReader->GetOutput()->GetSpacing(spacing);
dicomReader->GetOutput()->GetOrigin(origin);
spacing[2] = 1.5;

vtkImageData* image = vtkImageData ::New();
image->SetSpacing(spacing);

std::cout << "Sorting succeeded:" << spacing[2] << std::endl;
double center[3];
center[0] = origin[0] + spacing[0] * 0.5 * (extent[0] + extent[1]);
center[1] = origin[1] + spacing[1] * 0.5 * (extent[2] + extent[3]);
center[2] = origin[2] + spacing[2] * 0.5 * (extent[4] + extent[5]);

// Matrices for axial, coronal, sagittal, oblique view orientations
static double axialElements[16] = {
1, 0, 0, center[0],
0, 1, 0, center[1],
0, 0, 1, center[2],
0, 0, 0, 1 };

static double coronalElements[16] = {
1, 0, 0, center[0],
0, 0, 1, center[1],
0,-1, 0, center[2],
0, 0, 0, 1 };

static double sagittalElements[16] = {
0, 0,-1, center[0],
1, 0, 0, center[1],
0,-1, 0, center[2],
0, 0, 0, 1 };

// Set the slice orientation
vtkMatrix4x4 *resliceAxes = vtkMatrix4x4::New();
resliceAxes->DeepCopy(coronalElements);


vtkImageReslice *Reslice = vtkImageReslice::New();
Reslice->SetInputConnection(changer->GetOutputPort());
Reslice->SetOutputDimensionality(3);
Reslice->SetResliceAxes(resliceAxes);
Reslice->SetInterpolationModeToLinear();

vtkImageViewer2 *viewer = vtkImageViewer2::New();
viewer->SetupInteractor(interactor);
viewer->SetSize(500,500);
viewer->SetColorWindow(1024);
viewer->SetColorLevel(800);
viewer->SetInputConnection(Reslice->GetOutputPort());

// Create a greyscale lookup table
vtkLookupTable *table = vtkLookupTable::New();
table->SetRange(0, 200); // image intensity range
table->SetValueRange(0.0, 1.0); // from black to white
table->SetSaturationRange(0.0, 0.0); // no color saturation
table->SetRampToLinear();
table->Build();

// Map the image through the lookup table
vtkImageMapToColors *color = vtkImageMapToColors::New();
color->SetLookupTable(table);
color->SetInputConnection(Reslice->GetOutputPort());

// Display the image
vtkImageActor *actor = vtkImageActor::New();
actor->SetInput(color->GetOutput());
vtkRenderer *renderer = vtkRenderer::New();
renderer->AddActor(actor);
vtkRenderWindow *window = vtkRenderWindow::New();
window->AddRenderer(renderer);


interactor->SetInteractorStyle(imageStyle);
window->SetInteractor(interactor);
window->Render();
interactor->Start();

Thanks for your help.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.vtk.org/pipermail/vtkusers/attachments/20130125/5f2c9b16/attachment.htm>