Ray casting算法--CT三维重建简例

#include<vtkRenderer.h>
#include<vtkRenderWindow.h>
#include<vtkRenderWindowInteractor.h>
#include<vtkVolume16Reader.h>
#include<vtkPolyDataMapper.h>
#include<vtkActor.h>
#include<vtkCamera.h>
#include<vtkStripper.h>
#include<vtkContourFilter.h>
#include<vtkPiecewiseFunction.h>
#include<vtkColorTransferFunction.h>
#include<vtkVolumeRayCastCompositeFunction.h>
#include<vtkVolumeRayCastMIPFunction.h>
#include<vtkVolumeRayCastMapper.h>
#include<vtkVolumeProperty.h>
#include<vtkVolume.h>
#include<vtkProperty.h>
#include<vtkImporter.h> 
#include<vtkTransform.h>


int main(int argc, char* argv[])
{

        
        vtkRenderWindow *Renwin=vtkRenderWindow::New(); 
        vtkRenderer *aRenderer=vtkRenderer::New();
        Renwin->AddRenderer(aRenderer);

vtkRenderWindowInteractor *iren=vtkRenderWindowInteractor::New();

iren->SetRenderWindow(Renwin);




        vtkVolume16Reader *reader=vtkVolume16Reader::New();
        reader->SetDataDimensions(64,64);
        reader->SetDataByteOrderToLittleEndian();
        reader->SetFilePrefix("E:\\SRTP\\VTK_ITKandCMake\\vtkdata-5.0.3\\Data\\headsq\\quarter");
        reader->SetImageRange(1,93);
        reader->SetDataSpacing(3.2,3.2,1.5);
reader->SetDataOrigin(-32*3.2,-32*3.2,-46*1.5);



vtkTransform *transform = vtkTransform::New();
// transform->RotateX(90);
reader->SetTransform(transform);
reader->Update();

        

        vtkContourFilter *skinextractor=vtkContourFilter::New();
        skinextractor->SetInput((vtkDataObject*)reader->GetOutput());
        skinextractor->SetValue(0,500.0);


        vtkStripper *skinstripper=vtkStripper::New();
        skinstripper->SetInput(skinextractor->GetOutput());