[vtkusers] Re: how to access 3D volume

[Divya Rathore]

Search in the examples directory if you are using the source code of VTK.
If you are using the binary (installable) version then also it is a matter of a simple file searching.

Windows/Linux?
  ----- Original Message ----- 
  From: Ravi Samala 
  To: vtkusers at vtk.org 
  Sent: Tuesday, October 04, 2005 10:00 PM
  Subject: Re: [vtkusers] Re: how to access 3D volume


  hi Rathore,

  I am kind of new to VTK, can you tell me where can I find ImagePlaneWidget.cpp file.

  thank you,

  Ravi

  Divya Rathore <divyaksr at iitk.ac.in> wrote:
    Check out ImagePlaneWidget.cpp (I would like to thank Dr. Nagarajan, raman_nagarajan at vsnl.net, a member of this list).

    regards,
    D. Rathore
      ----- Original Message ----- 
      From: Ravi Samala 
      To: vtkusers at vtk.org 
      Sent: Monday, October 03, 2005 7:02 PM
      Subject: [vtkusers] how to access 3D volume


      hi vtkusers,

       I am reading 2D slices and constructing a 3D volume, using the following code....

      vtkVolume16Reader *v16Img = vtkVolume16Reader::New();

      v16Img->SetDataDimensions (sizex,sizey);

      v16Img->SetImageRange (atoi(argv[3]),atoi(argv[4])); 

      v16Img->SetHeaderSize(HeadS);

      v16Img->SetFilePrefix (argv[2]); 

      v16Img->SetDataSpacing (0.1, 0.1, 0.5);

      v16Img->SetDataOrigin(0, 0, 0);

      cout<<"Image Volume Read...."<<endl;




      //CompositeFunction with option

      vtkVolumeRayCastCompositeFunction *compositeFunction = vtkVolumeRayCastCompositeFunction::New();

      compositeFunction->SetCompositeMethodToInterpolateFirst(); 



      //Mapper with selected RayCastFunction

      vtkVolumeRayCastMapper *volumeMapper = vtkVolumeRayCastMapper::New();

      volumeMapper->SetInput(v16Img -> GetOutput());

      volumeMapper->SetVolumeRayCastFunction(compositeFunction);



      // Create transfer mapping scalar value to opacity

      vtkPiecewiseFunction *pieceFunc = vtkPiecewiseFunction::New();

      pieceFunc->AddPoint(0, 0.0); // alpha value for the lower intensity limit

      pieceFunc->AddPoint(15000, 0.5); // alpha value for the upper intensity limit



      //Set Volume Property

      vtkVolumeProperty *volumeProperty = vtkVolumeProperty::New();

      volumeProperty->SetInterpolationTypeToLinear();

      volumeProperty->SetScalarOpacity(pieceFunc);

      volumeProperty->ShadeOn();



      //The special actor with properties

      vtkVolume *vol = vtkVolume::New();

      vol->SetMapper(volumeMapper);

      vol->SetProperty(volumeProperty);



      // An outline provides context around the data.

      vtkOutlineFilter *outlineData = vtkOutlineFilter::New();

      outlineData->SetInput((vtkDataSet *) v16Img->GetOutput());

      vtkPolyDataMapper *mapOutline = vtkPolyDataMapper::New();

      mapOutline->SetInput(outlineData->GetOutput());

      vtkActor *outline = vtkActor::New();

      outline->SetMapper(mapOutline);

      outline->GetProperty()->SetColor(0,0,0); 




      //Starting to Render

      vtkRenderer *renderer=vtkRenderer::New();

      renderer->AddVolume(vol); 

      renderer->SetBackground(1,1,1);



      vtkRenderWindow *renWin=vtkRenderWindow::New();

      renWin->SetSize(500,500);

      renWin->DoubleBufferOn(); 

      renWin->AddRenderer(renderer);


      vtkRenderWindowInteractor *iren=vtkRenderWindowInteractor::New();

      iren->SetRenderWindow(renWin);

      iren->Initialize(); 

      iren->Start();



      my question is : is there any way I can access individual slices and a point in the 3 dimensional space (like p(x,y,z).



      thanks in advance,



      Ravi.



--------------------------------------------------------------------------
      Yahoo! for Good
      Click here to donate to the Hurricane Katrina relief effort. 


--------------------------------------------------------------------------


      _______________________________________________
      This is the private VTK discussion list. 
      Please keep messages on-topic. Check the FAQ at: http://www.vtk.org/Wiki/VTK_FAQ
      Follow this link to subscribe/unsubscribe:
      http://www.vtk.org/mailman/listinfo/vtkusers

    _______________________________________________
    This is the private VTK discussion list. 
    Please keep messages on-topic. Check the FAQ at: http://www.vtk.org/Wiki/VTK_FAQ
    Follow this link to subscribe/unsubscribe:
    http://www.vtk.org/mailman/listinfo/vtkusers



------------------------------------------------------------------------------
  Yahoo! for Good
  Click here to donate to the Hurricane Katrina relief effort. 


------------------------------------------------------------------------------


  _______________________________________________
  This is the private VTK discussion list. 
  Please keep messages on-topic. Check the FAQ at: http://www.vtk.org/Wiki/VTK_FAQ
  Follow this link to subscribe/unsubscribe:
  http://www.vtk.org/mailman/listinfo/vtkusers
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.vtk.org/pipermail/vtkusers/attachments/20051004/0754bdaa/attachment.htm>