VTK讀取序列JPG圖片用MC演算法進行三維重建，JPG圖片是用底層程式碼插值生成

 1 #include <vtkRenderer.h>
 2 #include <vtkRenderWindow.h>
 3 #include <vtkRenderWindowInteractor.h>
 4 #include <vtkJPEGReader.h>
 5 #include <vtkMarchingCubes.h>
 6 #include <vtkStripper.h>
 7 #include <vtkPolyDataMapper.h>
 8
 
 #include <vtkActor.h>
 9 #include <vtkProperty.h>
10 #include <vtkOutlineFilter.h>
11 #include <vtkSmartPointer.h>
12 #include <vtkCamera.h>
13  
14 void main()
15 {
16     vtkSmartPointer<vtkRenderer>ren=vtkSmartPointer<vtkRenderer>::New();
17     vtkSmartPointer<vtkRenderWindow>renWin=vtkSmartPointer<vtkRenderWindow>::New();
 
18     renWin->AddRenderer(ren);
19     vtkSmartPointer<vtkRenderWindowInteractor>iren=vtkSmartPointer<vtkRenderWindowInteractor>::New();
20     iren->SetRenderWindow(renWin);
21     
22     vtkSmartPointer<vtkJPEGReader>Reader=vtkSmartPointer<vtkJPEGReader>::New();
23     Reader->SetDataScalarTypeToUnsignedChar();                
 
24     Reader->SetFileDimensionality(3);
25     Reader->SetFilePrefix("E:\\JPG\\rgb");    
26     Reader->SetFileNameSliceSpacing(1);
27     Reader->SetFilePattern("%s%d.jpg");                       
28     Reader->SetDataExtent( 0, 200, 0, 400, 1,9);        
29     Reader->SetDataSpacing(1,1,10);                
30     Reader->Update();
31  
32     vtkSmartPointer<vtkMarchingCubes>marchingcube=vtkSmartPointer<vtkMarchingCubes>::New();
33     marchingcube->SetInput((vtkDataSet *)Reader->GetOutput());
34     //marchingcube->SetInputConnection(Reader->GetOutputPort());      //第二種讀取資料的方法
35     marchingcube->SetValue(0,200);  
36  
37     vtkSmartPointer<vtkStripper>Stripper=vtkSmartPointer<vtkStripper>::New();
38     Stripper->SetInput( marchingcube->GetOutput()); 
39  
40     vtkSmartPointer<vtkPolyDataMapper>polyMapper=vtkSmartPointer<vtkPolyDataMapper>::New();
41     polyMapper->SetInput(Stripper->GetOutput());
42     polyMapper->ScalarVisibilityOff();
43  
44     vtkSmartPointer<vtkActor>actor=vtkSmartPointer<vtkActor>::New();
45     actor->SetMapper(polyMapper);
46     actor->GetProperty()->SetDiffuseColor(1,0.19,0.15);
47     actor->GetProperty()->SetSpecular(0.1);
48     actor->GetProperty()->SetSpecularPower(10);
49     actor->GetProperty()->SetColor(1,0,0);
50  
51     vtkSmartPointer<vtkOutlineFilter>outlinefilter=vtkSmartPointer<vtkOutlineFilter>::New();
52     outlinefilter->SetInputConnection(Reader->GetOutputPort());
53  
54     vtkSmartPointer<vtkPolyDataMapper>outlineMapper=vtkSmartPointer<vtkPolyDataMapper>::New();
55     outlineMapper->SetInputConnection(outlinefilter->GetOutputPort());
56  
57     vtkSmartPointer<vtkActor>outlineActor=vtkSmartPointer<vtkActor>::New();
58     outlineActor->SetMapper(outlineMapper);
59     outlineActor->GetProperty()->SetColor(0,0,0);
60  
61     vtkSmartPointer<vtkCamera>aCamera=vtkSmartPointer<vtkCamera>::New();
62     aCamera->SetViewUp(0,0,-1);
63     aCamera->SetPosition(0,1,0);
64     aCamera->SetFocalPoint(0,0,0);
65     aCamera->ComputeViewPlaneNormal();
66  
67     ren->AddActor(actor);
68     ren->AddActor(outlineActor);
69     ren->SetActiveCamera(aCamera);
70     ren->ResetCamera();
71     aCamera->Dolly(1.5);
72  
73     ren->SetBackground(1,1,1);
74     renWin->SetSize(640,640);
75  
76     ren->ResetCameraClippingRange();
77     iren->Initialize();
78     iren->Start();
79 }
80