comparison data/model/FFTModel.cpp @ 280:daf89d31f45c

* Restore Scale setting which had mistakenly disappeared from spectrum
* Some tweaks to spectrum peaks
* Make spectrum default to 4096/75%... or is that overdoing it?
* Show harmonic cursor even during playback if in measure mode and mouse in
  relevant pane
* Avoid crash when trying to measure a non-layer
* Various adjustments to some icons to make them look better in white-on-black
  schemes

 
     size_t sampleRate = getSampleRate();
 
     std::deque<float> window;
     std::vector<size_t> inrange;
-    size_t medianWinSize = getPeakPickWindowSize(type, sampleRate, ymin);
+    float dist = 0.5;
+    size_t medianWinSize = getPeakPickWindowSize(type, sampleRate, ymin, dist);
     size_t halfWin = medianWinSize/2;
 
     size_t binmin;
     if (ymin > halfWin) binmin = ymin - halfWin;
     else binmin = 0;

 
         float value = values[bin];
 
         window.push_back(value);
 
-        medianWinSize = getPeakPickWindowSize(type, sampleRate, bin);
+        // so-called median will actually be the dist*100'th percentile
+        medianWinSize = getPeakPickWindowSize(type, sampleRate, bin, dist);
         halfWin = medianWinSize/2;
 
         while (window.size() > medianWinSize) window.pop_front();
 
         if (type == MajorPitchAdaptivePeaks) {

             else binmax = values.size()-1;
         }
 
         std::deque<float> sorted(window);
         std::sort(sorted.begin(), sorted.end());
-        float median = sorted[sorted.size()/2];
+        float median = sorted[int(sorted.size() * dist)];
 
         if (value > median) {
             inrange.push_back(bin);
         }
 

 
     return peaks;
 }
 
 size_t
-FFTModel::getPeakPickWindowSize(PeakPickType type, size_t sampleRate, size_t bin) const
-{
+FFTModel::getPeakPickWindowSize(PeakPickType type, size_t sampleRate,
+                                size_t bin, float &percentile) const
+{
+    percentile = 0.5;
     if (type == MajorPeaks) return 10;
     if (bin == 0) return 3;
+
     size_t fftSize = m_server->getFFTSize() >> m_yshift;
     float binfreq = (sampleRate * bin) / fftSize;
     float hifreq = Pitch::getFrequencyForPitch(73, 0, binfreq);
+
     int hibin = lrintf((hifreq * fftSize) / sampleRate);
     int medianWinSize = hibin - bin;
     if (medianWinSize < 3) medianWinSize = 3;
+
+    percentile = 0.5 + (binfreq / sampleRate);
+
     return medianWinSize;
 }
 
 FFTModel::PeakSet
 FFTModel::getPeakFrequencies(PeakPickType type, size_t x,