mImage->GetTransform()->GetMatrix()->SetElement(i,3,10*mImage->GetTransform()->GetMatrix()->GetElement(i,3));
mImage->GetTransform()->Inverse();
- mImage->UpdateReslice();
}
}
}
f.close();
itk::Matrix<double, 4, 4> itkMat = clitk::ReadMatrix3D(filename);
+
+ vtkSmartPointer<vtkMatrix4x4> matrix = vtkSmartPointer<vtkMatrix4x4>::New();
+ matrix->Identity();
for(int j=0; j<4; j++)
for(int i=0; i<4; i++)
- mImage->GetTransform()->GetMatrix()->SetElement(j,i,itkMat[j][i]);
- mImage->UpdateReslice();
+ matrix->SetElement(j,i,itkMat[j][i]);
+
+ mImage->GetTransform()->PostMultiply();
+ mImage->GetTransform()->Concatenate(matrix);
+ mImage->GetTransform()->Update();
}
}
//------------------------------------------------------------------------------