X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=vv%2FvvToolRigidReg.cxx;h=c02d7e6017fefca88875c19cde908fcc17e47389;hb=42df5bd85f2ad8eb5f9323fdf2441c7e36e2129f;hp=28591ba05ed1b6c6cc8af79af440ef5eaea6eb37;hpb=b449b21225742bcca0953edce7987c9fae44020b;p=clitk.git diff --git a/vv/vvToolRigidReg.cxx b/vv/vvToolRigidReg.cxx index 28591ba..c02d7e6 100644 --- a/vv/vvToolRigidReg.cxx +++ b/vv/vvToolRigidReg.cxx @@ -19,17 +19,23 @@ // vv #include "vvToolRigidReg.h" #include "vvSlicer.h" +#include // vtk #include #include +#include #include +#include +#include +#include // itk #include // clitk #include "clitkTransformUtilities.h" +#include "clitkMatrix.h" // qt #include @@ -118,8 +124,9 @@ void vvToolRigidReg::InputIsSelected(vvSlicerManager *input) //backup original matrix for(int j=0; j<4; j++) for(int i=0; i<4; i++) - mInitialMatrix->SetElement(i,j, mCurrentSlicerManager->GetImage()->GetTransform()->GetMatrix()->GetElement(i,j)); - QString origTransformString = dynamic_cast(mMainWindow)->Get4x4MatrixDoubleAsString(mInitialMatrix); + // TODO SR and BP: check on the list of transforms and not the first only + mInitialMatrix->SetElement(i,j, mCurrentSlicerManager->GetImage()->GetTransform()[0]->GetMatrix()->GetElement(i,j)); + QString origTransformString(clitk::Get4x4MatrixDoubleAsString(mInitialMatrix).c_str()); transformationLabel->setText(origTransformString); SetTransform(mInitialMatrix); @@ -242,7 +249,8 @@ void vvToolRigidReg::SpinBoxChange(double newVal) } // Compute transform and set - vtkSmartPointer transform_final=mInput->GetImage()->GetTransform(); + // TODO SR and BP: check on the list of transforms and not the first only + vtkSmartPointer transform_final=mInput->GetImage()->GetTransform()[0]; transform_final->Identity(); transform_final->PostMultiply(); @@ -294,8 +302,9 @@ void vvToolRigidReg::SaveFile() QFile file(filename); if (file.open(QFile::WriteOnly | QFile::Truncate)) { - vtkMatrix4x4* matrix = mCurrentSlicerManager->GetImage()->GetTransform()->GetMatrix(); - QString matrixStr = dynamic_cast(mMainWindow)->Get4x4MatrixDoubleAsString(matrix,16); + // TODO SR and BP: check on the list of transforms and not the first only + vtkMatrix4x4* matrix = mCurrentSlicerManager->GetImage()->GetTransform()[0]->GetMatrix(); + QString matrixStr = clitk::Get4x4MatrixDoubleAsString(matrix,16).c_str(); QTextStream out(&file); out << matrixStr; } @@ -332,7 +341,8 @@ void vvToolRigidReg::LoadFile() //------------------------------------------------------------------------------ void vvToolRigidReg::ChangeOfRotationCenter() { - SetTransform(mCurrentSlicerManager->GetImage()->GetTransform()->GetMatrix()); + // TODO SR and BP: check on the list of transforms and not the first only + SetTransform(mCurrentSlicerManager->GetImage()->GetTransform()[0]->GetMatrix()); } //------------------------------------------------------------------------------ @@ -347,7 +357,8 @@ void vvToolRigidReg::ResetTransform() void vvToolRigidReg::SetTransform(vtkMatrix4x4 *matrix) { vtkSmartPointer transform=vtkSmartPointer::New(); - mCurrentSlicerManager->GetImage()->GetTransform()->SetMatrix(matrix); + // TODO SR and BP: check on the list of transforms and not the first only + mCurrentSlicerManager->GetImage()->GetTransform()[0]->SetMatrix(matrix); transform->Update(); Render(); dynamic_cast(mMainWindow)->ImageInfoChanged(); @@ -426,13 +437,157 @@ void vvToolRigidReg::GetSlidersAndSpinBoxes(std::vector&transSliders, //------------------------------------------------------------------------------ //------------------------------------------------------------------------------ -void vvToolRigidReg::Render() +void vvToolRigidReg::ExtentMax(const double pointExtent[8][4], double maxExtent[2][3]) { - for (int i=0; iGetNumberOfSlicers(); i++) + double max, min; + for (int i=0; i<3; ++i) { + max = pointExtent[0][i]; + min = pointExtent[0][i]; + for (int j=1; j<8; ++j) { + if (pointExtent[j][i] > max) { + max = pointExtent[j][i]; + } + if (pointExtent[j][i] < min) { + min = pointExtent[j][i]; + } + } + maxExtent[0][i] = min; + maxExtent[1][i] = max; + } +} +//------------------------------------------------------------------------------ + +//------------------------------------------------------------------------------ +void vvToolRigidReg::Render() +{ //out << __func__ << endl; +#if VTK_MAJOR_VERSION > 5 +double translationValues[4], translationValuesUpdate[4]; +mCurrentSlicerManager->GetImage()->GetTransform()[0]->Print(cout); +vtkMatrix4x4* matrix = mCurrentSlicerManager->GetImage()->GetTransform()[0]->GetMatrix(); +vtkMatrix4x4* matrixTranspose = matrix->NewInstance(); +for (int i=0; i<3; ++i) { + for (int j=0; j<3; ++j) { + matrixTranspose->SetElement(i,j,matrix->GetElement(j,i)); + } +} +for (int j=0; j<3; ++j) +{ + translationValues[j] = matrix->GetElement(j,3); +} +translationValues[3] = 0.0; +matrix->MultiplyPoint(translationValues, translationValuesUpdate); +for (int i=0; i<3; ++i) { + matrixTranspose->SetElement(i,3,translationValuesUpdate[i]); +} +for (int i=0; i<4; ++i) { + matrixTranspose->SetElement(3,i,matrix->GetElement(3,i)); +} + +#endif +for (int i=0; iGetNumberOfSlicers(); i++) { +#if VTK_MAJOR_VERSION > 5 + double pointExtent[8][4], pointExtentUpdate[8][4]; + std::vector w_ext; + w_ext=mCurrentSlicerManager->GetImage()->GetSize(); + pointExtent[0][0] = 0.0; + pointExtent[0][1] = 0.0; + pointExtent[0][2] = 0.0; + pointExtent[0][3] = 1.0; + pointExtent[1][0] = w_ext[0]-1; + pointExtent[1][1] = w_ext[1]-1; + pointExtent[1][2] = w_ext[2]-1; + pointExtent[1][3] = 1.0; + pointExtent[2][0] = 0.0; + pointExtent[2][1] = w_ext[1]-1; + pointExtent[2][2] = w_ext[2]-1; + pointExtent[2][3] = 1.0; + pointExtent[3][0] = w_ext[0]-1; + pointExtent[3][1] = 0.0; + pointExtent[3][2] = w_ext[2]-1; + pointExtent[3][3] = 1.0; + pointExtent[4][0] = w_ext[0]-1; + pointExtent[4][1] = w_ext[1]-1; + pointExtent[4][2] = 0.0; + pointExtent[4][3] = 1.0; + pointExtent[5][0] = 0.0; + pointExtent[5][1] = 0.0; + pointExtent[5][2] = w_ext[2]-1; + pointExtent[5][3] = 1.0; + pointExtent[6][0] = 0.0; + pointExtent[6][1] = w_ext[1]-1; + pointExtent[6][2] = 0.0; + pointExtent[6][3] = 1.0; + pointExtent[7][0] = w_ext[0]-1; + pointExtent[7][1] = 0.0; + pointExtent[7][2] = 0.0; + pointExtent[7][3] = 1.0; + + for (int k=0; k<8; ++k) { + for (int j=0; j<3; ++j) + { + pointExtent[k][j] = mCurrentSlicerManager->GetImage()->GetSpacing()[j] * pointExtent[k][j]; + } + matrixTranspose->MultiplyPoint(pointExtent[k], pointExtentUpdate[k]); + for (int j=0; j<3; ++j) + { + pointExtentUpdate[k][j] = (pointExtentUpdate[k][j])/mCurrentSlicerManager->GetImage()->GetSpacing()[j]; + cout << pointExtentUpdate[k][j] << " "; + } + cout << endl; + } + double extUpdateTemp[2][3]; + int extUpdate[6]; + ExtentMax(pointExtentUpdate, extUpdateTemp); + for (int j=0; j<3; ++j) { + extUpdate[2*j] = 0; + extUpdate[2*j+1] = itk::Math::Round(extUpdateTemp[1][j] - extUpdateTemp[0][j]); + } + mCurrentSlicerManager->GetSlicer(i)->SetRegisterExtent(extUpdate); + extUpdate[2*mCurrentSlicerManager->GetSlicer(i)->GetOrientation()] = mCurrentSlicerManager->GetSlicer(i)->GetSlice(); + extUpdate[2*mCurrentSlicerManager->GetSlicer(i)->GetOrientation()+1] = mCurrentSlicerManager->GetSlicer(i)->GetSlice(); + + vtkSmartPointer mapperOpenGL= vtkSmartPointer::New(); + try { + mapperOpenGL = dynamic_cast(mCurrentSlicerManager->GetSlicer(i)->GetImageActor()->GetMapper()); + } catch (const std::bad_cast& e) { + std::cerr << e.what() << std::endl; + std::cerr << "Conversion error" << std::endl; + return; + } + cout << extUpdate[0] << " " << extUpdate[1] << " " << extUpdate[2] << " " << extUpdate[3] << " " << extUpdate[4] << " " << extUpdate[5] << endl; + mapperOpenGL->SetCroppingRegion(extUpdate); + + if (mCurrentSlicerManager->GetSlicer(i)->GetOverlay() && mCurrentSlicerManager->GetSlicer(i)->GetOverlayActor()->GetVisibility()) { + int extOverlayUpdate[6]; + for (int j=0; j<3; ++j) { + if (0 < extUpdateTemp[0][j]) { + extOverlayUpdate[2*j] = itk::Math::Round(extUpdateTemp[0][j]); + } else { + extOverlayUpdate[2*j] = 0; + } + if (extUpdateTemp[1][j] < w_ext[j]-1) { + extOverlayUpdate[2*j+1] = itk::Math::Round(extUpdateTemp[1][j]); + } else { + extOverlayUpdate[2*j+1] = w_ext[j]-1; + } + } + extOverlayUpdate[2*mCurrentSlicerManager->GetSlicer(i)->GetOrientation()] += mCurrentSlicerManager->GetSlicer(i)->GetSlice(); + extOverlayUpdate[2*mCurrentSlicerManager->GetSlicer(i)->GetOrientation()+1] = extOverlayUpdate[2*mCurrentSlicerManager->GetSlicer(i)->GetOrientation()]; + vtkSmartPointer mapperOpenGL= vtkSmartPointer::New(); + try { + mapperOpenGL = dynamic_cast(mCurrentSlicerManager->GetSlicer(i)->GetOverlayActor()->GetMapper()); + } catch (const std::bad_cast& e) { + std::cerr << e.what() << std::endl; + std::cerr << "Conversion error" << std::endl; + return; + } + mapperOpenGL->SetCroppingRegion(extOverlayUpdate); + } +#endif mCurrentSlicerManager->GetSlicer(i)->ForceUpdateDisplayExtent(); mCurrentSlicerManager->GetSlicer(i)->Render(); - } +} } //------------------------------------------------------------------------------