#include "itkPoint.h"
#include "clitkImageCommon.h"
#include "clitkCommon.h"
-
+#define VTK_EXCLUDE_STRSTREAM_HEADERS
+#include <vtkMatrix4x4.h>
+#include <vtkSmartPointer.h>
namespace clitk
{
return matrix;
}
- inline itk::Matrix<double, 4, 4> ReadMatrix3D(std::string fileName)
- {
+ itk::Matrix<double, 4, 4> ReadMatrix3D(std::string fileName);
+
+ inline vtkMatrix4x4* ReadVTKMatrix3D(std::string fileName) {
// read input matrix
std::ifstream is;
openFileForReading(is, fileName);
while (!is.eof()) {
nb.push_back(x);
skipComment(is);
- is >> x;
+ is >> x;
}
-
- //copy it to the matrix
- itk::Matrix<double, 4, 4> matrix;
+
+ vtkSmartPointer<vtkMatrix4x4> matrix = vtkSmartPointer<vtkMatrix4x4>::New();
unsigned int index=0;
for (unsigned int i=0;i<4;i++)
for (unsigned int j=0;j<4;j++)
- matrix[i][j]=nb[index++];
- return matrix;
+ matrix->SetElement(i,j, nb[index++]);
+
+ return matrix;
}
-
+
inline itk::Matrix<double, 3, 3> ReadMatrix2D(std::string fileName)
{
// read input matrix