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
