56 double camWorldPos[4];
60 double camWorldFocalPoint[4];
62 camWorldFocalPoint[3] = 1.0;
66 InverseVolumeMatrix->DeepCopy(volMatrix);
67 InverseVolumeMatrix->Invert();
68 InverseVolumeMatrix->MultiplyPoint(camWorldPos, camWorldPos);
69 InverseVolumeMatrix->MultiplyPoint(camWorldFocalPoint, camWorldFocalPoint);
71 this->CameraPosition =
vtkVector3d(camWorldPos[0], camWorldPos[1], camWorldPos[2]);
72 this->CameraPosition = this->CameraPosition /
vtkVector3d(camWorldPos[3]);
74 vtkVector3d camFP(camWorldFocalPoint[0], camWorldFocalPoint[1], camWorldFocalPoint[2]);
110 for (
int i = 0; i < 6; ++i)
112 int low = 2 * (i / 2);
113 bboundsP[i] = bbounds[i];
114 if (bboundsP[i] < abounds[low])
116 bboundsP[i] = abounds[low];
118 if (bboundsP[i] > abounds[low + 1])
120 bboundsP[i] = abounds[low + 1];
122 aboundsP[i] = abounds[i];
123 if (aboundsP[i] < bbounds[low])
125 aboundsP[i] = bbounds[low];
127 if (aboundsP[i] > bbounds[low + 1])
129 aboundsP[i] = bbounds[low + 1];
137 for (
int i = 0; i < 6; i += 2)
139 if (aboundsP[i] != aboundsP[i + 1])
141 dimSize[dims] = aboundsP[i + 1] - aboundsP[i];
146 degenAxes[degenDims] = i / 2;
156 if (dimSize[0] < dimSize[1])
158 if (dimSize[0] < dimSize[2])
169 if (dimSize[1] < dimSize[2])
183 double atobdir[3] = { bbounds[0] + bbounds[1] - abounds[0] - abounds[1],
184 bbounds[2] + bbounds[3] - abounds[2] - abounds[3],
185 bbounds[4] + bbounds[5] - abounds[4] - abounds[5] };
189 if (fabs(aboundsP[degenAxes[0] * 2] - bboundsP[degenAxes[0] * 2]) > 0.01 * atoblength)
196 if (this->CameraIsParallel)
203 vtkVector3d planePoint(0.25 * (aboundsP[0] + aboundsP[1] + bboundsP[0] + bboundsP[1]),
204 0.25 * (aboundsP[2] + aboundsP[3] + bboundsP[2] + bboundsP[3]),
205 0.25 * (aboundsP[4] + aboundsP[5] + bboundsP[4] + bboundsP[5]));
213 double plane[3] = { 0, 0, 0 };
214 plane[degenAxes[0]] = 1.0;
216 double dot = dir[0] * plane[0] + dir[1] * plane[1] + dir[2] * plane[2];
222 double side = plane[0] * atobdir[0] + plane[1] * atobdir[1] + plane[2] * atobdir[2];
223 return (dot * side < 0 ? 1 : -1);