** This is part of the CTSim program
** Copyright (C) 1983-2000 Kevin Rosenberg
**
-** $Id: imagefile.cpp,v 1.9 2000/07/15 08:36:13 kevin Exp $
+** $Id: imagefile.cpp,v 1.10 2000/08/22 16:49:56 kevin Exp $
**
** This program is free software; you can redistribute it and/or modify
** it under the terms of the GNU General Public License (version 2) as
}
mean /= (nx * ny);
- static const int nbin = 256;
+ static const int nbin = 1024;
int hist[ nbin ] = {0};
double spread = max - min;
mode = 0;
}
mode = (max_binindex * spread / (nbin - 1)) + min;
-
- median = 0.;
+
+ int nPixels = nx * ny;
+ slist<double> vecImage;
+ for (int ix = 0; ix < nx; ix++)
+ for (int iy = 0; iy < ny; iy++)
+ vecImage.push_front (v[ix][iy]);
+ vecImage.sort();
+ slist<double>::const_iterator iter = vecImage.begin();
+ for (int i = 0; i < (nPixels / 2) - 1; i++)
+ iter++; // Advance iterator to (nPixels / 2) - 1;
+
+ if (nPixels % 2) { // Odd
+ iter++;
+ median = *iter;
+ } else // Even
+ median = (*iter++ + *iter) / 2;
}