vtkResliceCursorWidget实现多平面重建效果

admin

vtkResliceCursorWidget和vtkResliceCursor类实现交互式图像切分，vtkResliceCursorWidget对象中需要定义相对应的vtkResliceCursor对象。vtkResliceCursorWidget定义十字坐标轴为用户提供便捷的切分和交互方式，支持坐标轴的旋转和平移；当坐标轴发生改变时，调用vtkResliceCursor来进行图像切分并更新到vtkRenderer对象中。

调整层厚的效果

1. // SPDX-FileCopyrightText: Copyright (c) Ken Martin, Will Schroeder, Bill
   
2. // Lorensen SPDX-License-Identifier: BSD-3-Clause
   
3. #include "vtkRegressionTestImage.h"
   
4. #include "vtkSmartPointer.h"
   
5. 
   
6. #include "vtkActor.h"
   
7. #include "vtkBiDimensionalWidget.h"
   
8. #include "vtkCamera.h"
   
9. #include "vtkCellPicker.h"
   
10. #include "vtkCommand.h"
    
11. #include "vtkDICOMImageReader.h"
    
12. #include "vtkImageActor.h"
    
13. #include "vtkImageData.h"
    
14. #include "vtkImageMapToColors.h"
    
15. #include "vtkImagePlaneWidget.h"
    
16. #include "vtkImageReader.h"
    
17. #include "vtkImageReslice.h"
    
18. #include "vtkInteractorEventRecorder.h"
    
19. #include "vtkInteractorStyleImage.h"
    
20. #include "vtkLookupTable.h"
    
21. #include "vtkOutlineFilter.h"
    
22. #include "vtkPlane.h"
    
23. #include "vtkPlaneSource.h"
    
24. #include "vtkPointData.h"
    
25. #include "vtkPolyDataMapper.h"
    
26. #include "vtkProperty.h"
    
27. #include "vtkRenderWindow.h"
    
28. #include "vtkRenderWindowInteractor.h"
    
29. #include "vtkRenderer.h"
    
30. #include "vtkResliceCursor.h"
    
31. #include "vtkResliceCursorActor.h"
    
32. #include "vtkResliceCursorLineRepresentation.h"
    
33. #include "vtkResliceCursorPolyDataAlgorithm.h"
    
34. #include "vtkResliceCursorThickLineRepresentation.h"
    
35. #include "vtkResliceCursorWidget.h"
    
36. #include "vtkVolume16Reader.h"
    
37. 
    
38. #include "vtkTestUtilities.h"
    
39. 
    
40. //------------------------------------------------------------------------------
    
41. class vtkResliceCursorCallback3 : public vtkCommand
    
42. {
    
43. public:
    
44.   static vtkResliceCursorCallback3* New()
    
45.   {
    
46.     return new vtkResliceCursorCallback3;
    
47.   }
    
48. 
    
49.   void Execute(vtkObject* caller, unsigned long /*ev*/, void* callData) override
    
50.   {
    
51.     vtkImagePlaneWidget* ipw = dynamic_cast<vtkImagePlaneWidget*>(caller);
    
52.     if (ipw)
    
53.     {
    
54.       double* wl = static_cast<double*>(callData);
    
55. 
    
56.       if (ipw == this->IPW[0])
    
57.       {
    
58.         this->IPW[1]->SetWindowLevel(wl[0], wl[1], 1);
    
59.         this->IPW[2]->SetWindowLevel(wl[0], wl[1], 1);
    
60.       }
    
61.       else if (ipw == this->IPW[1])
    
62.       {
    
63.         this->IPW[0]->SetWindowLevel(wl[0], wl[1], 1);
    
64.         this->IPW[2]->SetWindowLevel(wl[0], wl[1], 1);
    
65.       }
    
66.       else if (ipw == this->IPW[2])
    
67.       {
    
68.         this->IPW[0]->SetWindowLevel(wl[0], wl[1], 1);
    
69.         this->IPW[1]->SetWindowLevel(wl[0], wl[1], 1);
    
70.       }
    
71.     }
    
72. 
    
73.     vtkResliceCursorWidget* rcw = dynamic_cast<vtkResliceCursorWidget*>(caller);
    
74.     if (rcw)
    
75.     {
    
76.       vtkResliceCursorLineRepresentation* rep =
    
77.           dynamic_cast<vtkResliceCursorLineRepresentation*>(
    
78.               rcw->GetRepresentation());
    
79.       vtkResliceCursor* rc = rep->GetResliceCursorActor()
    
80.                                  ->GetCursorAlgorithm()
    
81.                                  ->GetResliceCursor();
    
82.       for (int i = 0; i < 3; i++)
    
83.       {
    
84.         vtkPlaneSource* ps =
    
85.             static_cast<vtkPlaneSource*>(this->IPW[i]->GetPolyDataAlgorithm());
    
86.         ps->SetNormal(rc->GetPlane(i)->GetNormal());
    
87.         ps->SetCenter(rc->GetPlane(i)->GetOrigin());
    
88. 
    
89.         // If the reslice plane has modified, update it on the 3D widget
    
90.         this->IPW[i]->UpdatePlacement();
    
91. 
    
92.         // std::cout << "Updating placement of plane: " << i << " " <<
    
93.         //  rc->GetPlane(i)->GetNormal()[0] << " " <<
    
94.         //  rc->GetPlane(i)->GetNormal()[1] << " " <<
    
95.         //  rc->GetPlane(i)->GetNormal()[2] << std::endl;
    
96.         // this->IPW[i]->GetReslice()->Print(cout);
    
97.         // rep->GetReslice()->Print(cout);
    
98.         // std::cout << "---------------------" << std::endl;
    
99.       }
    
100.     }
     
101. 
     
102.     // Render everything
     
103.     this->RCW[0]->Render();
     
104.   }
     
105. 
     
106.   vtkResliceCursorCallback3() = default;
     
107.   vtkImagePlaneWidget* IPW[3];
     
108.   vtkResliceCursorWidget* RCW[3];
     
109. };
     
110. 
     
111. //------------------------------------------------------------------------------
     
112. int main(int argc, char* argv[])
     
113. {
     
114.     vtkSmartPointer<vtkDICOMImageReader> reader =
     
115.       vtkSmartPointer<vtkDICOMImageReader>::New();
     
116.   reader->SetDirectoryName("D:/Data/two/SE1");
     
117.   reader->Update();
     
118. 
     
119.   vtkSmartPointer<vtkOutlineFilter> outline =
     
120.       vtkSmartPointer<vtkOutlineFilter>::New();
     
121.   outline->SetInputConnection(reader->GetOutputPort());
     
122. 
     
123.   vtkSmartPointer<vtkPolyDataMapper> outlineMapper =
     
124.       vtkSmartPointer<vtkPolyDataMapper>::New();
     
125.   outlineMapper->SetInputConnection(outline->GetOutputPort());
     
126. 
     
127.   vtkSmartPointer<vtkActor> outlineActor = vtkSmartPointer<vtkActor>::New();
     
128.   outlineActor->SetMapper(outlineMapper);
     
129. 
     
130.   vtkSmartPointer<vtkRenderer> ren[4];
     
131. 
     
132.   vtkSmartPointer<vtkRenderWindow> renWin =
     
133.       vtkSmartPointer<vtkRenderWindow>::New();
     
134.   renWin->SetMultiSamples(0);
     
135. 
     
136.   for (int i = 0; i < 4; i++)
     
137.   {
     
138.     ren[i] = vtkSmartPointer<vtkRenderer>::New();
     
139.     renWin->AddRenderer(ren[i]);
     
140.   }
     
141. 
     
142.   vtkSmartPointer<vtkRenderWindowInteractor> iren =
     
143.       vtkSmartPointer<vtkRenderWindowInteractor>::New();
     
144.   iren->SetRenderWindow(renWin);
     
145. 
     
146.   vtkSmartPointer<vtkCellPicker> picker = vtkSmartPointer<vtkCellPicker>::New();
     
147.   picker->SetTolerance(0.005);
     
148. 
     
149.   vtkSmartPointer<vtkProperty> ipwProp = vtkSmartPointer<vtkProperty>::New();
     
150. 
     
151.   // assign default props to the ipw's texture plane actor
     
152.   vtkSmartPointer<vtkImagePlaneWidget> planeWidget[3];
     
153.   int imageDims[3];
     
154.   reader->GetOutput()->GetDimensions(imageDims);
     
155. 
     
156.   for (int i = 0; i < 3; i++)
     
157.   {
     
158.     planeWidget[i] = vtkSmartPointer<vtkImagePlaneWidget>::New();
     
159.     planeWidget[i]->SetInteractor(iren);
     
160.     planeWidget[i]->SetPicker(picker);
     
161.     planeWidget[i]->RestrictPlaneToVolumeOn();
     
162.     double color[3] = {0, 0, 0};
     
163.     color[i] = 1;
     
164.     planeWidget[i]->GetPlaneProperty()->SetColor(color);
     
165.     planeWidget[i]->SetTexturePlaneProperty(ipwProp);
     
166.     planeWidget[i]->TextureInterpolateOff();
     
167.     planeWidget[i]->SetResliceInterpolateToLinear();
     
168.     planeWidget[i]->SetInputConnection(reader->GetOutputPort());
     
169.     planeWidget[i]->SetPlaneOrientation(i);
     
170.     planeWidget[i]->SetSliceIndex(imageDims[i] / 2);
     
171.     planeWidget[i]->DisplayTextOn();
     
172.     planeWidget[i]->SetDefaultRenderer(ren[3]);
     
173.     planeWidget[i]->SetWindowLevel(1358, -27);
     
174.     planeWidget[i]->On();
     
175.     planeWidget[i]->InteractionOn();
     
176.   }
     
177. 
     
178.   planeWidget[1]->SetLookupTable(planeWidget[0]->GetLookupTable());
     
179.   planeWidget[2]->SetLookupTable(planeWidget[0]->GetLookupTable());
     
180. 
     
181.   vtkSmartPointer<vtkResliceCursorCallback3> cbk =
     
182.       vtkSmartPointer<vtkResliceCursorCallback3>::New();
     
183. 
     
184.   // Create the reslice cursor, widget and rep
     
185. 
     
186.   vtkSmartPointer<vtkResliceCursor> resliceCursor =
     
187.       vtkSmartPointer<vtkResliceCursor>::New();
     
188.   resliceCursor->SetCenter(reader->GetOutput()->GetCenter());
     
189.   resliceCursor->SetThickMode(1);
     
190.   resliceCursor->SetThickness(10, 10, 10);
     
191.   resliceCursor->SetImage(reader->GetOutput());
     
192. 
     
193.   vtkSmartPointer<vtkResliceCursorWidget> resliceCursorWidget[3];
     
194.   vtkSmartPointer<vtkResliceCursorThickLineRepresentation> resliceCursorRep[3];
     
195. 
     
196.   double viewUp[3][3] = {{0, 0, -1}, {0, 0, 1}, {0, 1, 0}};
     
197.   for (int i = 0; i < 3; i++)
     
198.   {
     
199.     resliceCursorWidget[i] = vtkSmartPointer<vtkResliceCursorWidget>::New();
     
200.     resliceCursorWidget[i]->SetInteractor(iren);
     
201. 
     
202.     resliceCursorRep[i] =
     
203.         vtkSmartPointer<vtkResliceCursorThickLineRepresentation>::New();
     
204.     resliceCursorWidget[i]->SetRepresentation(resliceCursorRep[i]);
     
205.     resliceCursorRep[i]
     
206.         ->GetResliceCursorActor()
     
207.         ->GetCursorAlgorithm()
     
208.         ->SetResliceCursor(resliceCursor);
     
209.     resliceCursorRep[i]
     
210.         ->GetResliceCursorActor()
     
211.         ->GetCursorAlgorithm()
     
212.         ->SetReslicePlaneNormal(i);
     
213. 
     
214.     const double minVal = reader->GetOutput()->GetScalarRange()[0];
     
215.     if (vtkImageReslice* reslice =
     
216.             vtkImageReslice::SafeDownCast(resliceCursorRep[i]->GetReslice()))
     
217.     {
     
218.       reslice->SetBackgroundColor(minVal, minVal, minVal, minVal);
     
219.     }
     
220. 
     
221.     resliceCursorWidget[i]->SetDefaultRenderer(ren[i]);
     
222.     resliceCursorWidget[i]->SetEnabled(1);
     
223. 
     
224.     ren[i]->GetActiveCamera()->SetFocalPoint(0, 0, 0);
     
225.     double camPos[3] = {0, 0, 0};
     
226.     camPos[i] = 1;
     
227.     ren[i]->GetActiveCamera()->SetPosition(camPos);
     
228. 
     
229.     ren[i]->GetActiveCamera()->ParallelProjectionOn();
     
230.     ren[i]->GetActiveCamera()->SetViewUp(viewUp[i]);
     
231.     ren[i]->ResetCamera();
     
232.     // ren[i]->ResetCameraClippingRange();
     
233. 
     
234.     // Tie the Image plane widget and the reslice cursor widget together
     
235.     cbk->IPW[i] = planeWidget[i];
     
236.     cbk->RCW[i] = resliceCursorWidget[i];
     
237.     resliceCursorWidget[i]->AddObserver(
     
238.         vtkResliceCursorWidget::ResliceAxesChangedEvent, cbk);
     
239. 
     
240.     // Initialize the window level to a sensible value
     
241.     double range[2];
     
242.     reader->GetOutput()->GetScalarRange(range);
     
243.     resliceCursorRep[i]->SetWindowLevel(range[1] - range[0],
     
244.                                         (range[0] + range[1]) / 2.0);
     
245.     planeWidget[i]->SetWindowLevel(range[1] - range[0],
     
246.                                    (range[0] + range[1]) / 2.0);
     
247. 
     
248.     // Make them all share the same color map.
     
249.     resliceCursorRep[i]->SetLookupTable(resliceCursorRep[0]->GetLookupTable());
     
250.     planeWidget[i]->GetColorMap()->SetLookupTable(
     
251.         resliceCursorRep[0]->GetLookupTable());
     
252. 
     
253.     // clang-format off
     
254.     // Workaround VTK issue #18441
     
255.     // Make sure vtkResliceCursorActor is visible by forcing its representation to wireframe.
     
256.     // vtkResliceCursorActor is a quad with a normal parallel to the camera view up vector.
     
257.     // When represented as a surface, it has a thickness of 0 pixels. The class internally turns
     
258.     // edge visibility on to workaround the problem, which does not seem to be enough.
     
259.     resliceCursorRep[i]->GetResliceCursorActor()->GetCenterlineProperty(0)->SetRepresentationToWireframe();
     
260.     resliceCursorRep[i]->GetResliceCursorActor()->GetCenterlineProperty(1)->SetRepresentationToWireframe();
     
261.     resliceCursorRep[i]->GetResliceCursorActor()->GetCenterlineProperty(2)->SetRepresentationToWireframe();
     
262.     resliceCursorRep[i]->GetResliceCursorActor()->GetThickSlabProperty(0)->SetRepresentationToWireframe();
     
263.     resliceCursorRep[i]->GetResliceCursorActor()->GetThickSlabProperty(1)->SetRepresentationToWireframe();
     
264.     resliceCursorRep[i]->GetResliceCursorActor()->GetThickSlabProperty(2)->SetRepresentationToWireframe();
     
265.     // Workaround rendering artefacts with Intel chipsets and osmesa, where lines are rendered
     
266.     // black if perfectly aligned with the camera viewup (see #18453)
     
267.     resliceCursorRep[i]->GetResliceCursorActor()->GetCenterlineProperty(0)->RenderLinesAsTubesOn();
     
268.     resliceCursorRep[i]->GetResliceCursorActor()->GetCenterlineProperty(1)->RenderLinesAsTubesOn();
     
269.     resliceCursorRep[i]->GetResliceCursorActor()->GetCenterlineProperty(2)->RenderLinesAsTubesOn();
     
270.     resliceCursorRep[i]->GetResliceCursorActor()->GetCenterlineProperty(0)->SetLineWidth(2);
     
271.     resliceCursorRep[i]->GetResliceCursorActor()->GetCenterlineProperty(1)->SetLineWidth(2);
     
272.     resliceCursorRep[i]->GetResliceCursorActor()->GetCenterlineProperty(2)->SetLineWidth(2);
     
273.     resliceCursorRep[i]->GetResliceCursorActor()->GetThickSlabProperty(0)->RenderLinesAsTubesOn();
     
274.     resliceCursorRep[i]->GetResliceCursorActor()->GetThickSlabProperty(1)->RenderLinesAsTubesOn();
     
275.     resliceCursorRep[i]->GetResliceCursorActor()->GetThickSlabProperty(2)->RenderLinesAsTubesOn();
     
276.     resliceCursorRep[i]->GetResliceCursorActor()->GetThickSlabProperty(0)->SetLineWidth(2);
     
277.     resliceCursorRep[i]->GetResliceCursorActor()->GetThickSlabProperty(1)->SetLineWidth(2);
     
278.     resliceCursorRep[i]->GetResliceCursorActor()->GetThickSlabProperty(2)->SetLineWidth(2);
     
279.     // clang-format on
     
280.   }
     
281. 
     
282.   // Add the actors
     
283.   //
     
284.   ren[0]->SetBackground(0.3, 0.1, 0.1);
     
285.   ren[1]->SetBackground(0.1, 0.3, 0.1);
     
286.   ren[2]->SetBackground(0.1, 0.1, 0.3);
     
287.   ren[3]->AddActor(outlineActor);
     
288.   ren[3]->SetBackground(0.1, 0.1, 0.1);
     
289.   renWin->SetSize(600, 600);
     
290.   // renWin->SetFullScreen(1);
     
291. 
     
292.   ren[0]->SetViewport(0, 0, 0.5, 0.5);
     
293.   ren[1]->SetViewport(0.5, 0, 1, 0.5);
     
294.   ren[2]->SetViewport(0, 0.5, 0.5, 1);
     
295.   ren[3]->SetViewport(0.5, 0.5, 1, 1);
     
296. 
     
297.   // Set the actors' positions
     
298.   //
     
299.   renWin->Render();
     
300. 
     
301.   ren[3]->GetActiveCamera()->Elevation(110);
     
302.   ren[3]->GetActiveCamera()->SetViewUp(0, 0, -1);
     
303.   ren[3]->GetActiveCamera()->Azimuth(45);
     
304.   ren[3]->GetActiveCamera()->Dolly(1.15);
     
305.   ren[3]->ResetCameraClippingRange();
     
306. 
     
307.   vtkSmartPointer<vtkInteractorStyleImage> style =
     
308.       vtkSmartPointer<vtkInteractorStyleImage>::New();
     
309.   iren->SetInteractorStyle(style);
     
310.   iren->Initialize();
     
311.   iren->Start();
     
312. 
     
313.   return 1;
     
314. }
     

复制代码

hosstd

先看看先。。。

meddev

谢谢你的辛苦劳动了！！！

13636600141

先学习一下，以留备用

dentaldev

大哥谢了

mri

厉害

hosstd

次贴无用,检定完毕!

李孟

有吸引力，回复看看

medinfo

人生真是寂寞如雪啊。

proton

先看看在说