X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=common%2FvvImage.cxx;h=f75d3212267441b0ee75e2a0ffe221d00ad9066a;hb=49169efa175f39cb8052cea08186c032835bb2f2;hp=c7e13cc60b12fe9b98193efb232c8b1792802371;hpb=8b0f32e46083e5fbc2efe128c9eae002f129cbd8;p=clitk.git diff --git a/common/vvImage.cxx b/common/vvImage.cxx index c7e13cc..f75d321 100644 --- a/common/vvImage.cxx +++ b/common/vvImage.cxx @@ -261,6 +261,47 @@ const std::vector< vtkSmartPointer >& vvImage::GetTransform() //-------------------------------------------------------------------- +//-------------------------------------------------------------------- +void vvImage::InitializeTransform() +{ + for (int i = 0; i < mTransform.size(); i++) + mTransform[i]->Identity(); +} +//-------------------------------------------------------------------- + + +//-------------------------------------------------------------------- +std::vector< std::vector > vvImage::GetDirection() +{ + int dim = this->GetNumberOfDimensions(); + vtkSmartPointer transform = vtkSmartPointer::New(); + transform->Identity(); + for (int i = 0; i < mTransform.size(); i++) + transform->Concatenate(this->mTransform[i]); + + vtkSmartPointer matrix = transform->GetMatrix(); + matrix->Invert(); + std::vector > direction0; + for (int i = 0; i < dim; i++) { + if (i != 3) { + std::vector direction1; + for (int j = 0; j < dim; j++) { + if (j != 3) +#if VTK_MAJOR_VERSION <= 6 + direction1.push_back((*matrix)[i][j]); +#else + direction1.push_back((*matrix).GetElement(i,j)); +#endif + } + direction0.push_back(direction1); + } + } + return direction0; +} +//-------------------------------------------------------------------- + + + //-------------------------------------------------------------------- bool vvImage::HaveSameSizeAndSpacingThan(vvImage * other) {