/*=========================================================================
Program: vv http://www.creatis.insa-lyon.fr/rio/vv
- Authors belong to:
+ Authors belong to:
- University of LYON http://www.universite-lyon.fr/
- Léon Bérard cancer center http://oncora1.lyon.fnclcc.fr
- CREATIS CNRS laboratory http://www.creatis.insa-lyon.fr
#include "clitkImageExtractLine_ggo.h"
#include "clitkIO.h"
#include "clitkImageCommon.h"
+#include "clitkCommon.h"
#include <itkLineConstIterator.h>
//--------------------------------------------------------------------
-int main(int argc, char * argv[]) {
+int main(int argc, char * argv[])
+{
// Init command line
GGO(clitkImageExtractLine, args_info);
typedef float PixelType;
const unsigned int Dimension=3;
typedef itk::Image<PixelType, Dimension> ImageType;
-
+
// Check options
if (args_info.firstIndex_given != Dimension) {
std::cerr << "Please give " << Dimension << "values to --firstIndex option" << std::endl;
double length = 0.0;
for(unsigned int i=0; i<Dimension; i++) {
firstIndex[i] = args_info.firstIndex_arg[i];
- lastIndex[i] = args_info.lastIndex_arg[i];
+ lastIndex[i] = args_info.lastIndex_arg[i];
if (args_info.mm_flag) {
firstIndex[i] /= spacing[i];
lastIndex[i] /= spacing[i];
++iter;
}
double step = length/values.size();
-
+
// If isocenter is used
double isoDistance = 0.0;
if (args_info.isocenter_given) { // isoCenter is in mm
IndexType isoCenter;
- for(unsigned int i=0; i<Dimension; i++) {
- isoCenter[i] = args_info.isocenter_arg[i];
+ for(unsigned int i=0; i<Dimension; i++) {
+ isoCenter[i] = args_info.isocenter_arg[i];
isoDistance += pow(isoCenter[i] - firstIndex[i]*spacing[i],2);
}
DD(isoCenter);
lg += step;
}
os.close();
- }
- else {
+ } else {
for(unsigned int i=0; i<values.size(); i++) {
os << values[i] << std::endl;
}
os.close();
}
- // this is the end my friend
+ // this is the end my friend
return 0;
} // end main