[vtkusers] Visualize 3D image (.mhd)
Urvoas, Paul
paul.urvoas at philips.com
Mon Jul 12 08:20:40 EDT 2010
It change anything.
I have just found a program to read Meta Image. But I get the following error (during the execution):
NDims required and not defined.
MetaObject: Read: MET_Read Failed
MetaImage: M_Read: Error parsing file
MetaImage: Read: Cannot parse file
NDims required and not defined.
MetaObject: Read: MET_Read Failed
MetaImage: M_Read: Error parsing file
MetaImage: Read: Cannot parse file
My input is of course a 3d mhd file...
Thanks in advance,
Paul
This is the program I am trying to use:
int main( int argc, char * argv[] )
{
char * in_file_name = argv[1];
char * out_file_name = argv[2];
// View up
double vux = atof( argv[3] );
double vuy = atof( argv[4] );
double vuz = atof( argv[5] );
// Camera position direction from fp
double px = atof( argv[6] );
double py = atof( argv[7] );
double pz = atof( argv[8] );
// MIP or not
int mip = atoi( argv[9] );
// Focal point set to center.
double fpx;
double fpy;
double fpz;
// Graphics Factory
vtkGraphicsFactory * graphics_factory
= vtkGraphicsFactory::New();
graphics_factory->SetOffScreenOnlyMode( 1);
graphics_factory->SetUseMesaClasses( 1 );
// Imaging Factory
vtkImagingFactory * imaging_factory
= vtkImagingFactory::New();
imaging_factory->SetUseMesaClasses( 1 );
// Image Reader
//vtkXMLImageDataReader * reader
// = vtkXMLImageDataReader::New();
vtkMetaImageReader * reader = vtkMetaImageReader::New();
reader->SetFileName( in_file_name );
reader->ReleaseDataFlagOn();
// Get center
reader->Update();
vtkImageData * image = reader->GetOutput();
double * image_spacing = image->GetSpacing();
double * image_origin = image->GetOrigin();
int * whole_extent = image->GetWholeExtent();
int image_dim[3];
image_dim[0] = whole_extent[1] + 1;
image_dim[1] = whole_extent[3] + 1;
image_dim[2] = whole_extent[5] + 1;
fpx = image_origin[0]
+ image_dim[0] / 2 * image_spacing[0];
fpy = image_origin[1]
+ image_dim[1] / 2 * image_spacing[1];
fpz = image_origin[2]
+ image_dim[2] / 2 * image_spacing[2];
// Opacity
vtkPiecewiseFunction * opacity
= vtkPiecewiseFunction::New();
opacity->AddPoint( 0, 0.0 );
opacity->AddPoint( 255, 0.8 );
opacity->AddPoint( 256, 0.0 );
opacity->AddPoint( 511, 0.8 );
opacity->AddPoint( 512, 0.0 );
opacity->AddPoint( 767, 1.0 );
// Color
vtkColorTransferFunction * color
= vtkColorTransferFunction::New();
color->AddRGBPoint( 0, 1.0, 0.0, 0.0 );
color->AddRGBPoint( 255, 1.0, 0.0, 0.0 );
color->AddRGBPoint( 256, 0.0, 1.0, 0.0 );
color->AddRGBPoint( 511, 0.0, 1.0, 0.0 );
color->AddRGBPoint( 512, 1.0, 1.0, 0.0 );
color->AddRGBPoint( 767, 1.0, 1.0, 0.0 );
// Volume Property
vtkVolumeProperty * vol_prop = vtkVolumeProperty::New();
vol_prop->SetColor( color );
if ( mip == 0 )
{
vol_prop->SetScalarOpacity( opacity );
}
//vol_prop->SetInterpolationTypeToLinear();
vol_prop->SetInterpolationTypeToNearest();
// Mapper
vtkFixedPointVolumeRayCastMapper * mapper
= vtkFixedPointVolumeRayCastMapper::New();
if ( mip == 0 )
{
mapper->SetBlendModeToComposite();
}
else
{
mapper->SetBlendModeToMaximumIntensity();
}
mapper->SetInputConnection( reader->GetOutputPort() );
// Volume
vtkVolume * volume = vtkVolume::New();
volume->SetMapper( mapper );
volume->SetProperty( vol_prop );
// Camera
vtkCamera * camera = vtkCamera::New();
camera->ParallelProjectionOn();
camera->SetViewUp (vux, vuy, vuz);
camera->SetPosition (
fpx + px,
fpy + py,
fpz + pz );
camera->SetFocalPoint (fpx, fpy, fpz);
camera->ComputeViewPlaneNormal();
// Scale PolyData
vtkPoints * scale_points = vtkPoints::New();
scale_points->InsertNextPoint(-10.0, 0.0, 0.0);
scale_points->InsertNextPoint(-10.0, 100.0, 0.0);
vtkCellArray * scale_cells = vtkCellArray::New();
scale_cells->InsertNextCell( 2 );
scale_cells->InsertCellPoint( 0 );
scale_cells->InsertCellPoint( 1 );
vtkPolyData * scale_poly = vtkPolyData::New();
scale_poly->SetPoints( scale_points );
scale_poly->SetLines( scale_cells );
vtkPolyDataMapper * scale_mapper = vtkPolyDataMapper::New();
scale_mapper->SetInput( scale_poly );
vtkActor * scale_actor = vtkActor::New();
scale_actor->GetProperty()->SetColor(1.0, 1.0, 1.0);
scale_actor->SetMapper( scale_mapper );
// Renderer
vtkRenderer * renderer = vtkRenderer::New();
renderer->SetBackground( 0, 0, 0 );
renderer->AddVolume( volume );
//renderer->AddActor( scale_actor );
renderer->SetActiveCamera(camera);
renderer->ResetCamera();
camera->Zoom(1.3);
// Render Window
vtkRenderWindow * render_window = vtkRenderWindow::New();
render_window->SetSize( WSIZE, WSIZE );
render_window->SetOffScreenRendering( 1 );
render_window->AddRenderer( renderer );
// ***** STEREO *****
//render_window->StereoRenderOn();
// Window to Image
vtkWindowToImageFilter * win_2_image
= vtkWindowToImageFilter::New();
win_2_image->SetInput( render_window );
// PNG Writer
vtkPNGWriter * writer = vtkPNGWriter::New();
writer->SetFileName( out_file_name );
writer->SetInputConnection( win_2_image->GetOutputPort() );
render_window->Render();
writer->Write();
return 0;
}
-----Original Message-----
From: Bill Lorensen [mailto:bill.lorensen at gmail.com]
Sent: 2010 Jul 12 1:59 PM
To: Urvoas, Paul
Cc: vtkusers at vtk.org
Subject: Re: [vtkusers] Visualize 3D image (.mhd)
Updates are usually done on the filters, not the outpus of the filter.Try:
reader->Update();
and
caster->Update();
On Mon, Jul 12, 2010 at 7:08 AM, Urvoas, Paul <paul.urvoas at philips.com> wrote:
> Hi,
>
> I would like to visualize 3d mhd image with VTK, thanks to another program,
> I wrote the following code.
>
> It crashes just after cout<<"caster ok" and I don't find my mistakes.
>
>
>
> Thanks in advance,
>
> Paul
>
>
>
>
>
> int main(int argc, char* argv)
>
> {
>
>
>
> vtkImageData** ImageOriginale = new vtkImageData*[10];
>
>
>
> vtkViewImage3D* view;
>
>
>
> cout << "start" <<endl;
>
>
>
> vtkMetaImageReader* reader = vtkMetaImageReader::New();
>
> reader->SetFileName ("Synthetic-04.mhd");
>
> reader->GetOutput()->Update();
>
> cout << "reader ok" <<endl;
>
>
>
> vtkImageCast* caster = vtkImageCast::New();
>
> caster->SetInput (reader->GetOutput());
>
> caster->SetOutputScalarTypeToUnsignedShort ();
>
> caster->GetOutput()->Update();
>
>
>
> cout << "caster ok" <<endl;
>
>
>
> ImageOriginale[0]->DeepCopy(caster->GetOutput()) ;
>
>
>
> cout << "Copy volume to ImageOriaginale" <<endl;
>
>
>
> caster->Delete();
>
> reader->Delete();
>
>
>
> cout << "reader caster deleted" <<endl;
>
>
>
> int dims[3];
>
> ImageOriginale[0]->GetDimensions(dims);
>
> double spacing[3];
>
> ImageOriginale[0]->GetSpacing(spacing);
>
> cout << "spacing & dimension" <<endl;
>
>
>
>
>
> view->SetImage(ImageOriginale[0]);
>
>
>
> view->SetRenderingModeToVR();
>
> view->SetVolumeRayCastFunctionToMIP();
>
>
>
> vtkRenderer *ren = vtkRenderer::New();
>
> vtkRenderWindow *renWin = vtkRenderWindow::New();
>
> renWin->AddRenderer(ren);
>
> vtkRenderWindowInteractor *iren = vtkRenderWindowInteractor::New();
>
> iren->SetRenderWindow(renWin);
>
>
>
>
>
> view->SetImage (ImageOriginale[0]);
>
> view->SetRenderingModeToVR();
>
> view->SetVolumeRayCastFunctionToMIP();
>
>
>
> cout << "befor render" <<endl;
>
> renWin->Render();
>
>
>
> }
>
> ________________________________
> The information contained in this message may be confidential and legally
> protected under applicable law. The message is intended solely for the
> addressee(s). If you are not the intended recipient, you are hereby notified
> that any use, forwarding, dissemination, or reproduction of this message is
> strictly prohibited and may be unlawful. If you are not the intended
> recipient, please contact the sender by return e-mail and destroy all copies
> of the original message.
>
> _______________________________________________
> Powered by www.kitware.com
>
> Visit other Kitware open-source projects at
> http://www.kitware.com/opensource/opensource.html
>
> Please keep messages on-topic and check the VTK FAQ at:
> http://www.vtk.org/Wiki/VTK_FAQ
>
> Follow this link to subscribe/unsubscribe:
> http://www.vtk.org/mailman/listinfo/vtkusers
>
>
The information contained in this message may be confidential and legally protected under applicable law. The message is intended solely for the addressee(s). If you are not the intended recipient, you are hereby notified that any use, forwarding, dissemination, or reproduction of this message is strictly prohibited and may be unlawful. If you are not the intended recipient, please contact the sender by return e-mail and destroy all copies of the original message.
More information about the vtkusers
mailing list