- vecImage.push_front (v[ix5][iy]);
- vecImage.sort();
- std::list<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;
+ vecImage.push_back (v[ix5][iy]);
+ std::sort(vecImage.begin(), vecImage.end());\r
+\r
+ if (nPixels % 2) // Odd
+ median = vecImage[((nPixels - 1) / 2)];\r
+ else // Even
+ median = (vecImage[(nPixels / 2) - 1] + vecImage[nPixels / 2]) / 2;