#include <itksys/SystemTools.hxx>
-#ifndef _WIN32
+#ifdef _WIN32
#define _USE_MATH_DEFINES
#include <cmath>
+ const double M_PI = std::acos(-1.0);
#endif
// clitk include
mMesh->SetPoints(mPoints);
vtkIdType ids[2];
for (unsigned int idx=0 ; idx<mNbOfPoints ; idx++) {
- double pointIn[4], pointOut[4];
+ double pointIn[4];
for (unsigned int j=0 ; j<3; ++j)
pointIn[j] = mData->GetPoint(idx)[j];
pointIn[4] = 1.0;
- /*mTransformMatrix->MultiplyPoint(pointIn, pointOut);
+ /*double pointOut[4];
+ mTransformMatrix->MultiplyPoint(pointIn, pointOut);
std::cout << pointOut[0] << " " << pointOut[1] << " " << pointOut[2] << " " << pointOut[3] << std::endl;
mMesh->GetPoints()->InsertNextPoint(pointOut[0],
pointOut[1],
int n = input->GetLargestPossibleRegion().GetNumberOfPixels();
std::vector<int> should_keep(values.size());;
std::fill(should_keep.begin(), should_keep.end(), 0);
- while (current<max and i<n) { // loop by decreasing pixel values
+ while (current<max && i<n) { // loop by decreasing pixel values
current += values[indices[i]];
should_keep[indices[i]] = 1.0;
++i;
int n = input->GetLargestPossibleRegion().GetNumberOfPixels();
std::vector<int> should_keep(values.size());;
std::fill(should_keep.begin(), should_keep.end(), 0);
- while (current<max and i<n) { // loop by decreasing pixel values
+ while (current<max && i<n) { // loop by decreasing pixel values
current += values[indices[i]];
should_keep[indices[i]] = 1.0;
++i;