[vtkusers] vtkvolume to be used in dealunary3D

[InfoSeekerr]

Hi all,

I am trying to use the following steps to create tetrahdral mesh as appear
in the photo ins this link :
http://www.iue.tuwien.ac.at/phd/fleischmann/img346.gif 

1. I read the 2D slices using vtkDICOMImageReader 
2. apply delaunry3D or vtkDataSetTriangleFilte for the vtkdicomimagereader
output.
3. defien iso value to render the skull 
4. create a volume and rendering it. 


there're two error appears , which are :
1.vtkstreamingdemandpibline error. 
2. vtkfixedpointvolumeraycastmapper error : No input.




I think the problem in use the output of vtkdicomimagereader in delaunary3d
any suggestion please 






the code is :



//Variables
			array<double> ^metaScalarRange = gcnew array<double>(2) { 0, 0 };
			vtk:: vtkRenderWindow renwM = gcnew vtk::vtkRenderWindow();

			vtk:: vtkRenderer ren1 = gcnew vtkRenderer();
			vtk:: vtkRenderWindowInteractor iren = gcnew
vtk::vtkRenderWindowInteractor();

			vtk:: vtkVolume volume = gcnew vtk::vtkVolume();
			renwM::AddRenderer(ren1);


			iren::SetRenderWindow(renwM);

			// Read a DICOM series path:

			vtk:: vtkDICOMImageReader VDR = gcnew vtk::vtkDICOMImageReader();
			VDR::SetDirectoryName("H:\\work\\Master
Degree\\DataSet\\case1-\\DICOM\\PA1\\ST1\\SE1");

			VDR::SetDataOrigin(0, 0, 0);

			VDR->Update();


			/////
			vtk:: vtkImageChangeInformation VIC = gcnew
vtk::vtkImageChangeInformation();

			VIC::SetInput(VDR::GetOutput());

			VIC::CenterImageOn();


			VIC->Update();

			///////// Applying the mesh

			vtkDataSetTriangleFilter ^VDST = gcnew vtkDataSetTriangleFilter();
			VDST->SetInputConnection(VDR::GetOutputPort());


//            
//            vtkDelaunay3D VDST = new vtkDelaunay3D();
//
//            VDST.SetInputConnection(VDR.GetOutputPort()); 

			// Load the data on vtkimagedata.

			vtk:: vtkImageData VoxelData = VIC::GetOutput();
			// Start adjust the rendering

			vtkVolumeProperty ^volumeProperty = gcnew vtkVolumeProperty();
			metaScalarRange = VoxelData::GetScalarRange();


			int isovalue = Convert::ToInt32(0.8 * (metaScalarRange[1] +
metaScalarRange[0]));
			/////////////////////////////////////////////

			array<double> ^inputRange = VoxelData::GetScalarRange();
			vtk:: vtkPiecewiseFunction opacityIsoTransferFunction = gcnew
vtk::vtkPiecewiseFunction();

			double per = ((inputRange[1] - inputRange[0]) / 100);
			if (isovalue > (inputRange[0] + per))

			{

				opacityIsoTransferFunction::AddPoint(inputRange[0], 0);


				opacityIsoTransferFunction::AddPoint(isovalue - (per), 0);

			}


			opacityIsoTransferFunction::AddPoint(isovalue, 1);

			opacityIsoTransferFunction::AddPoint(inputRange[1], 1);


			vtk:: vtkColorTransferFunction colorTransferFunction = gcnew
vtk::vtkColorTransferFunction();
			colorTransferFunction::AddRGBPoint(inputRange[0], 1.0, 1.0, 1.0);


			colorTransferFunction::AddRGBPoint(inputRange[1], 1.0, 1.0, 1.0);

			volumeProperty->SetColor(colorTransferFunction);


			volumeProperty->SetScalarOpacity(opacityIsoTransferFunction);

			volumeProperty->ShadeOn();


			volumeProperty->SetInterpolationTypeToLinear();

			volume::SetProperty(volumeProperty);


			///////////////////////////
			vtk:: vtkFixedPointVolumeRayCastMapper volumeMapper = gcnew
vtk::vtkFixedPointVolumeRayCastMapper();

			volumeMapper::SetInputConnection(VDST->GetOutputPort());

			volumeMapper::AutoAdjustSampleDistancesOn();


			volumeMapper->Update();

			volume::SetMapper(volumeMapper);


			volume->Update();

			//////////////////////////////////

			vtkActor ^actor = gcnew vtkActor();
			actor->GetProperty()->SetRepresentationToWireframe();


			ren1::AddVolume(volume);

			ren1::AddActor(actor);


			ren1::ResetCamera();

//            ren1.GetActiveCamera().Azimuth(0);
//            ren1.GetActiveCamera().Roll(0);


			ren1::GetActiveCamera()->Elevation(90);

			renwM::Render();


			iren->Initialize();
			iren::Enable();



best regards 

InfoSeekerr wrote:
> 
> I have a vtkvolume into  vtkdelaunary3d?
> 
> 
> by another words:\how to convert this volume to polydata, unstructured
> grid,... to be suitable to enter to the filter 
> 

-- 
View this message in context: http://www.nabble.com/vtkvolume-to-be-used-in-dealunary3D-tp25966433p26034468.html
Sent from the VTK - Users mailing list archive at Nabble.com.