View difference between Paste ID: G2Hirxbb and PArvffSq

 #include "itkImage.h"
 #include "itkImageFileWriter.h"
 #include "itkImageFileReader.h"
 #include "itkMinimumMaximumImageCalculator.h"
 #include "itkImageRegionIterator.h"
 #include <iostream>
 #include <string>
 #include <stdlib.h>
 #include <fstream>
 
 using namespace std;
 
-const unsigned int RDImageDimension = ImageDimension-1;
+
 typedef double  PixelType;
 typedef itk::Image< PixelType, ImageDimension > ImageType;
-typedef itk::Image< PixelType, RDImageDimension> RDImageType;
+
 typedef itk::ImageFileWriter< ImageType > WriterType;
 
-ofstream file;
+
 {
-void writetofile(PixelType min, PixelType max)
+
 	typedef itk::MinimumMaximumImageCalculator <ImageType> 	ImageCalculatorFilterType;
-	file << min << " " << max << std::endl;
+
 
 	imageCalculatorFilter->SetImage( inputimage );
 	imageCalculatorFilter->Compute();
 
 	PixelType Imin = imageCalculatorFilter->GetMinimum();
 	PixelType Imax = imageCalculatorFilter->GetMaximum();
 	PixelType Omin = 0;
 	PixelType Omax = 1;
-	double pi = 4.0*atan(1.0);
+
-	double U = rand()/float(RAND_MAX); double V = rand()/float(RAND_MAX);
+	typedef itk::ImageRegionIterator< ImageType > IteratorType;
-	double X = sqrt(-2*log(U))*cos(2*pi*V);
+	IteratorType imageIterator(inputimage,iterregion);
-	double Y = sqrt(-2*log(U))*sin(2*pi*V);
+
-	double randomnumber = (Y*Y + 0.1)*(X/abs(X));
+
 
 	std::cout << "image iterator at end?" << imageIterator.IsAtEnd() << std::endl; //gives 1 always
 
 	while(!imageIterator.IsAtEnd())
-	std::cout << 25*randomnumber << std::endl;
+
 		PixelType inputPixel = imageIterator.Get();
-	PixelType Omin = 25*randomnumber;
+		PixelType outputPixel = (inputPixel-Imin)*((Omax-Omin)/(Imax-Imin))+Omin;
-	PixelType Omax = 255-Omin;
+
 		++imageIterator;
 	}
-	writetofile(Imin, Imax);
+
 	ImageType::Pointer outputimage = inputimage;
-	itk::ImageRegionIterator<ImageType> imageIterator(inputimage,iterregion);
+
-	PixelType outputPixel;
+
-	PixelType inputPixel;
+
 
 int main(int argc , char *argv[])
 {
 	if( argc < 3 )
 	{
 		std::cerr << "Usage: " << std::endl;
 		std::cerr << argv[1] << " inputImageFile" << std::endl;
-		inputPixel = imageIterator.Get();
+
-		outputPixel = (inputPixel-Imin)*((Omax-Omin)/(Imax-Imin))+Omin;
+
 	}
 
 	ReaderType::Pointer	reader = ReaderType::New();
 
 	reader->SetFileName(argv[1]);
 	reader->Update();
 
 	ImageType::Pointer image = reader->GetOutput();
 	ImageType::RegionType region = image->GetLargestPossibleRegion();
 	ImageType::SizeType size = region.GetSize();
 
-	//srand(0);
+
-	file.open("minmax.txt");
+
 
 	int nrofloops = size[ImageDimension-1];
 
 	size[ImageDimension-1] = 0;
 
 	region.SetSize(size);
 
 	ImageType::Pointer scaledImage;
 	scaledImage = image;
 
 	for(int i=0; i < nrofloops; i++)
 	{
 		volumestart[ImageDimension-1] = i;
 		region.SetIndex(volumestart);
 		//PixelType max = Imaxvector.get(i);
 		scaledImage = randomscaler(scaledImage, max, region, scalingback);
 	}
 
 	return 0;
 }