- for (int j=0; j<3; ++j) {
- if (0 < extUpdateTemp[0][j]) {
- extOverlayUpdate[2*j] = itk::Math::Round<double>(extUpdateTemp[0][j]);
+ for (int j=0; j<3; ++j) { //Rotation
+ if (extOverlayUpdateTemp[1][j] - extOverlayUpdateTemp[0][j] > w_ext[j]-1) {
+ extOverlayUpdate[2*j] = 0;
+ extOverlayUpdate[2*j+1] = w_ext[j]-1;
+ } else {
+ extOverlayUpdate[2*j] = itk::Math::Round<double>(extOverlayUpdateTemp[0][j]);
+ extOverlayUpdate[2*j+1] = itk::Math::Round<double>(extOverlayUpdateTemp[1][j]);
+ }
+ }
+
+ //Compute translation
+ double pointOrigin[4], pointOriginUpdate[4];
+ for (int j=0; j<3; ++j)
+ {
+ pointOrigin[j] = 0 - centre[j];
+ }
+ pointOrigin[3] = 0.0;
+ matrix->MultiplyPoint(pointOrigin, pointOriginUpdate);
+ for (int j=0; j<3; ++j)
+ {
+ pointOriginUpdate[j] = (pointOriginUpdate[j] + centre[j])/mCurrentSlicerManager->GetImage()->GetSpacing()[j];
+ translation[j] = matrix->GetElement(j,3) - pointOriginUpdate[j]*mCurrentSlicerManager->GetImage()->GetSpacing()[j];
+ pointOrigin[j] = translation[j];
+ }
+ pointOrigin[3] = 0.0;
+ matrixTranspose->MultiplyPoint(pointOrigin, pointOriginUpdate);
+ for (int j=0; j<3; ++j)
+ {
+ translation[j] = pointOriginUpdate[j]/mCurrentSlicerManager->GetImage()->GetSpacing()[j];
+ }
+
+ for (int j=0; j<3; ++j) { //Translation
+ if (0 < extOverlayUpdateTemp[0][j] - translation[j]) {
+ extOverlayUpdate[2*j] = itk::Math::Round<double>(extOverlayUpdateTemp[0][j] - translation[j]);