2 #include "vtkMPR3DDataViewer.h"
5 vtkMPR3DDataViewer::vtkMPR3DDataViewer()
7 _visiblePosition[0]=false;
8 _visiblePosition[1]=false;
9 _visiblePosition[2]=false;
24 //-------------------------------------------------------------------
25 vtkMPR3DDataViewer::~vtkMPR3DDataViewer()
27 _outlineData -> Delete();
28 _mapOutline -> Delete();
31 // if (_bwLut) _bwLut -> Delete();
32 // if (_hueLut) _hueLut -> Delete();
33 // if (_satLut) _satLut -> Delete();
34 if (_ctfun) _ctfun -> Delete();
35 if (_saggitalColors) _saggitalColors -> Delete();
36 if (_saggital) _saggital -> Delete();
37 if (_axialColors) _axialColors -> Delete();
38 if (_axial) _axial -> Delete();
39 if (_coronalColors) _coronalColors -> Delete();
40 if (_coronal) _coronal -> Delete();
43 //-------------------------------------------------------------------
44 vtkActor* vtkMPR3DDataViewer::GetOutlineActor()
48 //-------------------------------------------------------------------
49 vtkImageActor* vtkMPR3DDataViewer::GetImageActor(int id)
51 vtkImageActor *tmpVtkActor=NULL;
53 tmpVtkActor = GetvtkActor_saggital();
56 tmpVtkActor = GetvtkActor_coronal();
59 tmpVtkActor = GetvtkActor_axial();
63 //-------------------------------------------------------------------
64 void vtkMPR3DDataViewer::Refresh()
66 int x = (int)(_vtkmprbasedata->GetX());
67 int y = (int)(_vtkmprbasedata->GetY());
68 int z = (int)(_vtkmprbasedata->GetZ());
73 //-------------------------------------------------------------------
74 void vtkMPR3DDataViewer::SetVtkMPRBaseData(vtkMPRBaseData *vtkmprbasedata)
76 _vtkmprbasedata=vtkmprbasedata;
78 //-------------------------------------------------------------------
79 vtkMPRBaseData* vtkMPR3DDataViewer::GetVtkMPRBaseData()
81 return _vtkmprbasedata;
83 //-------------------------------------------------------------------
85 vtkColorTransferFunction *vtkMPR3DDataViewer::GetvtkColorTransferFunction()
89 //-------------------------------------------------------------------
91 std::vector<double>* vtkMPR3DDataViewer::GetctFunVectorPoint()
93 return &_ctFunVectorPoint;
95 //-------------------------------------------------------------------
96 std::vector<double>* vtkMPR3DDataViewer::GetctFunVectorRed()
98 return &_ctFunVectorRed;
101 //-------------------------------------------------------------------
102 std::vector<double>* vtkMPR3DDataViewer::GetctFunVectorGreen()
104 return &_ctFunVectorGreen;
107 //-------------------------------------------------------------------
108 std::vector<double>* vtkMPR3DDataViewer::GetctFunVectorBlue()
110 return &_ctFunVectorBlue;
113 //-------------------------------------------------------------------
114 void vtkMPR3DDataViewer::Configure()
117 // Start by creatin a black/white lookup table.
118 _bwLut = vtkLookupTable::New();
119 _bwLut->SetTableRange (0, 2000);
120 _bwLut->SetSaturationRange (0, 0);
121 _bwLut->SetHueRange (0, 0);
122 _bwLut->SetValueRange (0, 1);
124 // Now create a lookup table that consists of the full hue circle
126 _hueLut = vtkLookupTable::New();
127 _hueLut->SetTableRange (0, 2000);
128 _hueLut->SetHueRange (0, 1);
129 _hueLut->SetSaturationRange (1, 1);
130 _hueLut->SetValueRange (1, 1);
132 // Finally, create a lookup table with a single hue but having a range
133 // in the saturation of the hue.
134 _satLut = vtkLookupTable::New();
135 _satLut->SetTableRange (0, 2000);
136 _satLut->SetHueRange (.6, .6);
137 _satLut->SetSaturationRange (0, 1);
138 _satLut->SetValueRange (1, 1);
143 _vtkmprbasedata->GetImageData()->GetScalarRange(range);
144 double max = range[1];
146 _ctFunVectorPoint.clear();
147 _ctFunVectorPoint.push_back(max*0/4);
148 _ctFunVectorPoint.push_back(max*1/4);
149 _ctFunVectorPoint.push_back(max*2/4);
150 _ctFunVectorPoint.push_back(max*3/4);
151 _ctFunVectorPoint.push_back(max*4/4);
153 _ctFunVectorRed.clear();
154 _ctFunVectorRed.push_back(0.0);
155 _ctFunVectorRed.push_back(1.0);
156 _ctFunVectorRed.push_back(0.0);
157 _ctFunVectorRed.push_back(0.0);
158 _ctFunVectorRed.push_back(0.0);
160 _ctFunVectorGreen.clear();
161 _ctFunVectorGreen.push_back(0.0);
162 _ctFunVectorGreen.push_back(0.0);
163 _ctFunVectorGreen.push_back(0.0);
164 _ctFunVectorGreen.push_back(1.0);
165 _ctFunVectorGreen.push_back(0.2);
167 _ctFunVectorBlue.clear();
168 _ctFunVectorBlue.push_back(0.0);
169 _ctFunVectorBlue.push_back(0.0);
170 _ctFunVectorBlue.push_back(1.0);
171 _ctFunVectorBlue.push_back(0.0);
172 _ctFunVectorBlue.push_back(0.0);
176 _ctfun = vtkColorTransferFunction::New();
178 _ctfun->RemoveAllPoints();
180 int i,size=_ctFunVectorPoint.size();
183 _ctfun->AddRGBPoint( _ctFunVectorPoint[i] , _ctFunVectorRed[i],_ctFunVectorGreen[i],_ctFunVectorBlue[i]);
187 _ctfun->AddRGBPoint( 100 ,0 ,0,0);
188 _ctfun->AddRGBPoint( 200 ,1 ,0,0);
189 _ctfun->AddRGBPoint( 300 ,0 ,1,0);
190 _ctfun->AddRGBPoint( 400 ,0 ,0,1);
192 // void *p=this; // JPRx
194 // Create the first of the three planes. The filter vtkImageMapToColors
195 // maps the data through the corresponding lookup table created above. The
196 // vtkImageActor is a type of vtkProp and conveniently displays an image on
197 // a single quadrilateral plane. It does this using texture mapping and as
198 // a result is quite fast. (Note: the input image has to be unsigned char
199 // values, which the vtkImageMapToColors produces.) Note also that by
200 // specifying the DisplayExtent, the pipeline requests data of this extent
201 // and the vtkImageMapToColors only processes a slice of data.
203 if(_saggitalColors==NULL){
204 _saggitalColors = vtkImageMapToColors::New();
206 _saggitalColors->RemoveAllInputs();
208 _saggitalColors->SetInput( _vtkmprbasedata->GetImageData() );
209 // _saggitalColors->SetLookupTable(_bwLut);
211 _saggitalColors->SetLookupTable(_ctfun);
216 _saggital = vtkImageActor::New();
217 _saggital->SetInput(_saggitalColors->GetOutput());
219 //_saggitalColors->Update();
223 // Create the second (axial) plane of the three planes. We use the
224 // same approach as before except that the extent differs.
227 if(_axialColors==NULL){
228 _axialColors = vtkImageMapToColors::New();
231 _axialColors->RemoveAllInputs();
232 _axialColors->SetInput( _vtkmprbasedata->GetImageData() );
233 // _axialColors->SetLookupTable(_hueLut);
234 _axialColors->SetLookupTable(_ctfun);
237 _axial = vtkImageActor::New();
238 _axial->SetInput(_axialColors->GetOutput());
241 //_axialColors->Update();
244 // Create the third (coronal) plane of the three planes. We use
245 // the same approach as before except that the extent differs.
246 if(_coronalColors==NULL){
247 _coronalColors = vtkImageMapToColors::New();
249 _coronalColors->RemoveAllInputs();
250 _coronalColors->SetInput( _vtkmprbasedata->GetImageData() );
251 // _coronalColors->SetLookupTable(_satLut);
252 _coronalColors->SetLookupTable(_ctfun);
255 _coronal = vtkImageActor::New();
256 _coronal->SetInput(_coronalColors->GetOutput());
261 // An outline provides context around the data.
263 if(_outlineData==NULL){
264 _outlineData = vtkOutlineFilter::New();
266 _outlineData->RemoveAllInputs();
267 _outlineData->SetInput((vtkDataSet *) _vtkmprbasedata->GetImageData() );
269 if(_mapOutline==NULL){
270 _mapOutline = vtkPolyDataMapper::New();
271 _mapOutline->SetInput(_outlineData->GetOutput());
273 _mapOutline->Update();
274 //_mapOutline->RemoveAllInputs();
277 _outline = vtkActor::New();
278 _outline->SetMapper(_mapOutline);
279 _outline->GetProperty()->SetColor(0,0,0);
285 _vtkmprbasedata->GetImageData()->GetExtent(ext);
292 //-------------------------------------------------------------------
293 vtkImageActor * vtkMPR3DDataViewer::GetvtkActor_saggital()
297 //-------------------------------------------------------------------
298 vtkImageActor * vtkMPR3DDataViewer::GetvtkActor_coronal()
302 //-------------------------------------------------------------------
303 vtkImageActor * vtkMPR3DDataViewer::GetvtkActor_axial()
308 //------------------------------------------------------------------------
309 void vtkMPR3DDataViewer::SetPositionX(int pos){
310 int x1,x2,y1,y2,z1,z2;
311 _vtkmprbasedata->GetDimensionExtention(&x1,&x2,&y1,&y2,&z1,&z2);
312 _saggital->SetDisplayExtent( pos , pos , y1 ,y2 , z1 , z2 );
314 //------------------------------------------------------------------------
315 void vtkMPR3DDataViewer::SetPositionY(int pos){
316 int x1,x2,y1,y2,z1,z2;
317 _vtkmprbasedata->GetDimensionExtention(&x1,&x2,&y1,&y2,&z1,&z2);
318 _coronal->SetDisplayExtent(x1,x2, pos,pos, z1,z2);
320 //------------------------------------------------------------------------
321 void vtkMPR3DDataViewer::SetPositionZ(int pos){
322 int x1,x2,y1,y2,z1,z2;
323 _vtkmprbasedata->GetDimensionExtention(&x1,&x2,&y1,&y2,&z1,&z2);
324 _axial->SetDisplayExtent(x1,x2, y1,y2, pos,pos);
326 //-------------------------------------------------------------------
327 void vtkMPR3DDataViewer::SetVisiblePosition(int idPosition, bool visible)
329 _visiblePosition[idPosition]=visible;
331 //-------------------------------------------------------------------
332 bool vtkMPR3DDataViewer::GetVisiblePosition(int idPosition)
334 return _visiblePosition[idPosition];