Chris@58: /* -*- c-basic-offset: 4 indent-tabs-mode: nil -*-  vi:set ts=8 sts=4 sw=4: */
Chris@0: 
Chris@0: /*
Chris@59:     Sonic Visualiser
Chris@59:     An audio file viewer and annotation editor.
Chris@59:     Centre for Digital Music, Queen Mary, University of London.
Chris@484:     This file copyright 2006-2009 Chris Cannam and QMUL.
Chris@0:     
Chris@59:     This program is free software; you can redistribute it and/or
Chris@59:     modify it under the terms of the GNU General Public License as
Chris@59:     published by the Free Software Foundation; either version 2 of the
Chris@59:     License, or (at your option) any later version.  See the file
Chris@59:     COPYING included with this distribution for more information.
Chris@0: */
Chris@0: 
Chris@0: #include "SpectrogramLayer.h"
Chris@0: 
Chris@128: #include "view/View.h"
Chris@0: #include "base/Profiler.h"
Chris@0: #include "base/AudioLevel.h"
Chris@0: #include "base/Window.h"
Chris@24: #include "base/Pitch.h"
Chris@118: #include "base/Preferences.h"
Chris@167: #include "base/RangeMapper.h"
Chris@253: #include "base/LogRange.h"
Chris@1063: #include "base/ColumnOp.h"
Chris@376: #include "widgets/CommandHistory.h"
Chris@1078: #include "data/model/Dense3DModelPeakCache.h"
Chris@1078: 
Chris@376: #include "ColourMapper.h"
Chris@283: #include "ImageRegionFinder.h"
Chris@690: #include "PianoScale.h"
Chris@1078: #include "PaintAssistant.h"
Chris@1089: #include "Colour3DPlotRenderer.h"
Chris@0: 
Chris@0: #include <QPainter>
Chris@0: #include <QImage>
Chris@0: #include <QPixmap>
Chris@0: #include <QRect>
Chris@92: #include <QApplication>
Chris@178: #include <QMessageBox>
Chris@283: #include <QMouseEvent>
Chris@316: #include <QTextStream>
Chris@1017: #include <QSettings>
Chris@0: 
Chris@0: #include <iostream>
Chris@0: 
Chris@0: #include <cassert>
Chris@0: #include <cmath>
Chris@0: 
Chris@545: #ifndef __GNUC__
Chris@545: #include <alloca.h>
Chris@545: #endif
Chris@545: 
Chris@1090: #define DEBUG_SPECTROGRAM 1
Chris@1044: #define DEBUG_SPECTROGRAM_REPAINT 1
Chris@1025: 
Chris@1025: using namespace std;
Chris@907: 
Chris@44: SpectrogramLayer::SpectrogramLayer(Configuration config) :
Chris@0:     m_model(0),
Chris@0:     m_channel(0),
Chris@0:     m_windowSize(1024),
Chris@0:     m_windowType(HanningWindow),
Chris@97:     m_windowHopLevel(2),
Chris@0:     m_gain(1.0),
Chris@215:     m_initialGain(1.0),
Chris@37:     m_threshold(0.0),
Chris@215:     m_initialThreshold(0.0),
Chris@9:     m_colourRotation(0),
Chris@215:     m_initialRotation(0),
Chris@119:     m_minFrequency(10),
Chris@0:     m_maxFrequency(8000),
Chris@135:     m_initialMaxFrequency(8000),
Chris@1105:     m_colourScale(ColourScaleType::Log),
Chris@197:     m_colourMap(0),
Chris@1103:     m_binScale(BinScale::Linear),
Chris@1103:     m_binDisplay(BinDisplay::AllBins),
Chris@1104:     m_normalization(ColumnNormalization::None),
Chris@1104:     m_normalizeVisibleArea(false),
Chris@133:     m_lastEmittedZoomStep(-1),
Chris@390:     m_synchronous(false),
Chris@608:     m_haveDetailedScale(false),
Chris@193:     m_exiting(false),
Chris@1088:     m_fftModel(0),
Chris@1088:     m_peakCache(0),
Chris@1088:     m_peakCacheDivisor(8)
Chris@0: {
Chris@1017:     QString colourConfigName = "spectrogram-colour";
Chris@1017:     int colourConfigDefault = int(ColourMapper::Green);
Chris@1017:     
Chris@215:     if (config == FullRangeDb) {
Chris@215:         m_initialMaxFrequency = 0;
Chris@215:         setMaxFrequency(0);
Chris@215:     } else if (config == MelodicRange) {
Chris@0: 	setWindowSize(8192);
Chris@97: 	setWindowHopLevel(4);
Chris@215:         m_initialMaxFrequency = 1500;
Chris@215: 	setMaxFrequency(1500);
Chris@215:         setMinFrequency(40);
Chris@1105: 	setColourScale(ColourScaleType::Linear);
Chris@215:         setColourMap(ColourMapper::Sunset);
Chris@1103:         setBinScale(BinScale::Log);
Chris@1017:         colourConfigName = "spectrogram-melodic-colour";
Chris@1017:         colourConfigDefault = int(ColourMapper::Sunset);
Chris@224: //        setGain(20);
Chris@37:     } else if (config == MelodicPeaks) {
Chris@37: 	setWindowSize(4096);
Chris@97: 	setWindowHopLevel(5);
Chris@135:         m_initialMaxFrequency = 2000;
Chris@40: 	setMaxFrequency(2000);
Chris@37: 	setMinFrequency(40);
Chris@1103: 	setBinScale(BinScale::Log);
Chris@1105: 	setColourScale(ColourScaleType::Linear);
Chris@1103: 	setBinDisplay(BinDisplay::PeakFrequencies);
Chris@1104:         setNormalization(ColumnNormalization::Max1);
Chris@1017:         colourConfigName = "spectrogram-melodic-colour";
Chris@1017:         colourConfigDefault = int(ColourMapper::Sunset);
Chris@0:     }
Chris@110: 
Chris@1017:     QSettings settings;
Chris@1017:     settings.beginGroup("Preferences");
Chris@1017:     setColourMap(settings.value(colourConfigName, colourConfigDefault).toInt());
Chris@1017:     settings.endGroup();
Chris@1017:     
Chris@122:     Preferences *prefs = Preferences::getInstance();
Chris@122:     connect(prefs, SIGNAL(propertyChanged(PropertyContainer::PropertyName)),
Chris@122:             this, SLOT(preferenceChanged(PropertyContainer::PropertyName)));
Chris@122:     setWindowType(prefs->getWindowType());
Chris@0: }
Chris@0: 
Chris@0: SpectrogramLayer::~SpectrogramLayer()
Chris@0: {
Chris@1106:     invalidateRenderers();
Chris@1088:     invalidateFFTModel();
Chris@0: }
Chris@0: 
Chris@1105: ColourScaleType
Chris@1104: SpectrogramLayer::convertToColourScale(int value)
Chris@1104: {
Chris@1104:     switch (value) {
Chris@1105:     case 0: return ColourScaleType::Linear;
Chris@1105:     case 1: return ColourScaleType::Meter;
Chris@1105:     case 2: return ColourScaleType::Log; //!!! db^2
Chris@1105:     case 3: return ColourScaleType::Log;
Chris@1105:     case 4: return ColourScaleType::Phase;
Chris@1105:     default: return ColourScaleType::Linear;
Chris@1104:     }
Chris@1104: }
Chris@1104: 
Chris@1104: int
Chris@1105: SpectrogramLayer::convertFromColourScale(ColourScaleType scale)
Chris@1104: {
Chris@1104:     switch (scale) {
Chris@1105:     case ColourScaleType::Linear: return 0;
Chris@1105:     case ColourScaleType::Meter: return 1;
Chris@1105:     case ColourScaleType::Log: return 3; //!!! + db^2
Chris@1105:     case ColourScaleType::Phase: return 4;
Chris@1105: 
Chris@1105:     case ColourScaleType::PlusMinusOne:
Chris@1105:     case ColourScaleType::Absolute:
Chris@1104:     default: return 0;
Chris@1104:     }
Chris@1104: }
Chris@1104: 
Chris@1104: std::pair<ColumnNormalization, bool>
Chris@1104: SpectrogramLayer::convertToColumnNorm(int value)
Chris@1104: {
Chris@1104:     switch (value) {
Chris@1104:     default:
Chris@1104:     case 0: return { ColumnNormalization::None, false };
Chris@1104:     case 1: return { ColumnNormalization::Max1, false };
Chris@1104:     case 2: return { ColumnNormalization::None, true }; // visible area
Chris@1104:     case 3: return { ColumnNormalization::Hybrid, false };
Chris@1104:     }
Chris@1104: }
Chris@1104: 
Chris@1104: int
Chris@1104: SpectrogramLayer::convertFromColumnNorm(ColumnNormalization norm, bool visible)
Chris@1104: {
Chris@1104:     if (visible) return 2;
Chris@1104:     switch (norm) {
Chris@1104:     case ColumnNormalization::None: return 0;
Chris@1104:     case ColumnNormalization::Max1: return 1;
Chris@1104:     case ColumnNormalization::Hybrid: return 3;
Chris@1104: 
Chris@1104:     case ColumnNormalization::Sum1:
Chris@1104:     default: return 0;
Chris@1104:     }
Chris@1104: }
Chris@1104: 
Chris@0: void
Chris@0: SpectrogramLayer::setModel(const DenseTimeValueModel *model)
Chris@0: {
Chris@682: //    cerr << "SpectrogramLayer(" << this << "): setModel(" << model << ")" << endl;
Chris@34: 
Chris@110:     if (model == m_model) return;
Chris@110: 
Chris@0:     m_model = model;
Chris@1088:     invalidateFFTModel();
Chris@0: 
Chris@0:     if (!m_model || !m_model->isOK()) return;
Chris@0: 
Chris@320:     connectSignals(m_model);
Chris@0: 
Chris@0:     connect(m_model, SIGNAL(modelChanged()), this, SLOT(cacheInvalid()));
Chris@906:     connect(m_model, SIGNAL(modelChangedWithin(sv_frame_t, sv_frame_t)),
Chris@906: 	    this, SLOT(cacheInvalid(sv_frame_t, sv_frame_t)));
Chris@0: 
Chris@0:     emit modelReplaced();
Chris@110: }
Chris@115: 
Chris@0: Layer::PropertyList
Chris@0: SpectrogramLayer::getProperties() const
Chris@0: {
Chris@0:     PropertyList list;
Chris@87:     list.push_back("Colour");
Chris@87:     list.push_back("Colour Scale");
Chris@87:     list.push_back("Window Size");
Chris@97:     list.push_back("Window Increment");
Chris@862:     list.push_back("Normalization");
Chris@87:     list.push_back("Bin Display");
Chris@87:     list.push_back("Threshold");
Chris@87:     list.push_back("Gain");
Chris@87:     list.push_back("Colour Rotation");
Chris@153: //    list.push_back("Min Frequency");
Chris@153: //    list.push_back("Max Frequency");
Chris@87:     list.push_back("Frequency Scale");
Chris@0:     return list;
Chris@0: }
Chris@0: 
Chris@87: QString
Chris@87: SpectrogramLayer::getPropertyLabel(const PropertyName &name) const
Chris@87: {
Chris@87:     if (name == "Colour") return tr("Colour");
Chris@87:     if (name == "Colour Scale") return tr("Colour Scale");
Chris@87:     if (name == "Window Size") return tr("Window Size");
Chris@112:     if (name == "Window Increment") return tr("Window Overlap");
Chris@862:     if (name == "Normalization") return tr("Normalization");
Chris@87:     if (name == "Bin Display") return tr("Bin Display");
Chris@87:     if (name == "Threshold") return tr("Threshold");
Chris@87:     if (name == "Gain") return tr("Gain");
Chris@87:     if (name == "Colour Rotation") return tr("Colour Rotation");
Chris@87:     if (name == "Min Frequency") return tr("Min Frequency");
Chris@87:     if (name == "Max Frequency") return tr("Max Frequency");
Chris@87:     if (name == "Frequency Scale") return tr("Frequency Scale");
Chris@87:     return "";
Chris@87: }
Chris@87: 
Chris@335: QString
Chris@862: SpectrogramLayer::getPropertyIconName(const PropertyName &) const
Chris@335: {
Chris@335:     return "";
Chris@335: }
Chris@335: 
Chris@0: Layer::PropertyType
Chris@0: SpectrogramLayer::getPropertyType(const PropertyName &name) const
Chris@0: {
Chris@87:     if (name == "Gain") return RangeProperty;
Chris@87:     if (name == "Colour Rotation") return RangeProperty;
Chris@87:     if (name == "Threshold") return RangeProperty;
Chris@0:     return ValueProperty;
Chris@0: }
Chris@0: 
Chris@0: QString
Chris@0: SpectrogramLayer::getPropertyGroupName(const PropertyName &name) const
Chris@0: {
Chris@153:     if (name == "Bin Display" ||
Chris@153:         name == "Frequency Scale") return tr("Bins");
Chris@87:     if (name == "Window Size" ||
Chris@1086: 	name == "Window Increment") return tr("Window");
Chris@87:     if (name == "Colour" ||
Chris@87: 	name == "Threshold" ||
Chris@87: 	name == "Colour Rotation") return tr("Colour");
Chris@862:     if (name == "Normalization" ||
Chris@153:         name == "Gain" ||
Chris@87: 	name == "Colour Scale") return tr("Scale");
Chris@0:     return QString();
Chris@0: }
Chris@0: 
Chris@0: int
Chris@0: SpectrogramLayer::getPropertyRangeAndValue(const PropertyName &name,
Chris@216: 					   int *min, int *max, int *deflt) const
Chris@0: {
Chris@216:     int val = 0;
Chris@216: 
Chris@216:     int garbage0, garbage1, garbage2;
Chris@55:     if (!min) min = &garbage0;
Chris@55:     if (!max) max = &garbage1;
Chris@216:     if (!deflt) deflt = &garbage2;
Chris@10: 
Chris@87:     if (name == "Gain") {
Chris@0: 
Chris@0: 	*min = -50;
Chris@0: 	*max = 50;
Chris@0: 
Chris@906:         *deflt = int(lrint(log10(m_initialGain) * 20.0));
Chris@216: 	if (*deflt < *min) *deflt = *min;
Chris@216: 	if (*deflt > *max) *deflt = *max;
Chris@216: 
Chris@906: 	val = int(lrint(log10(m_gain) * 20.0));
Chris@216: 	if (val < *min) val = *min;
Chris@216: 	if (val > *max) val = *max;
Chris@0: 
Chris@87:     } else if (name == "Threshold") {
Chris@37: 
Chris@37: 	*min = -50;
Chris@37: 	*max = 0;
Chris@37: 
Chris@906:         *deflt = int(lrint(AudioLevel::multiplier_to_dB(m_initialThreshold)));
Chris@216: 	if (*deflt < *min) *deflt = *min;
Chris@216: 	if (*deflt > *max) *deflt = *max;
Chris@216: 
Chris@906: 	val = int(lrint(AudioLevel::multiplier_to_dB(m_threshold)));
Chris@216: 	if (val < *min) val = *min;
Chris@216: 	if (val > *max) val = *max;
Chris@37: 
Chris@87:     } else if (name == "Colour Rotation") {
Chris@9: 
Chris@9: 	*min = 0;
Chris@9: 	*max = 256;
Chris@216:         *deflt = m_initialRotation;
Chris@216: 
Chris@216: 	val = m_colourRotation;
Chris@9: 
Chris@87:     } else if (name == "Colour Scale") {
Chris@0: 
Chris@1099:         // linear, meter, db^2, db, phase
Chris@0: 	*min = 0;
Chris@176: 	*max = 4;
Chris@1092:         *deflt = 2;
Chris@216: 
Chris@1104: 	val = convertFromColourScale(m_colourScale);
Chris@0: 
Chris@87:     } else if (name == "Colour") {
Chris@0: 
Chris@0: 	*min = 0;
Chris@196: 	*max = ColourMapper::getColourMapCount() - 1;
Chris@216:         *deflt = 0;
Chris@216: 
Chris@216: 	val = m_colourMap;
Chris@0: 
Chris@87:     } else if (name == "Window Size") {
Chris@0: 
Chris@0: 	*min = 0;
Chris@0: 	*max = 10;
Chris@216:         *deflt = 5;
Chris@0: 	
Chris@216: 	val = 0;
Chris@0: 	int ws = m_windowSize;
Chris@216: 	while (ws > 32) { ws >>= 1; val ++; }
Chris@0: 
Chris@97:     } else if (name == "Window Increment") {
Chris@0: 	
Chris@0: 	*min = 0;
Chris@97: 	*max = 5;
Chris@216:         *deflt = 2;
Chris@216: 
Chris@216:         val = m_windowHopLevel;
Chris@0:     
Chris@87:     } else if (name == "Min Frequency") {
Chris@37: 
Chris@37: 	*min = 0;
Chris@37: 	*max = 9;
Chris@216:         *deflt = 1;
Chris@37: 
Chris@37: 	switch (m_minFrequency) {
Chris@216: 	case 0: default: val = 0; break;
Chris@216: 	case 10: val = 1; break;
Chris@216: 	case 20: val = 2; break;
Chris@216: 	case 40: val = 3; break;
Chris@216: 	case 100: val = 4; break;
Chris@216: 	case 250: val = 5; break;
Chris@216: 	case 500: val = 6; break;
Chris@216: 	case 1000: val = 7; break;
Chris@216: 	case 4000: val = 8; break;
Chris@216: 	case 10000: val = 9; break;
Chris@37: 	}
Chris@37:     
Chris@87:     } else if (name == "Max Frequency") {
Chris@0: 
Chris@0: 	*min = 0;
Chris@0: 	*max = 9;
Chris@216:         *deflt = 6;
Chris@0: 
Chris@0: 	switch (m_maxFrequency) {
Chris@216: 	case 500: val = 0; break;
Chris@216: 	case 1000: val = 1; break;
Chris@216: 	case 1500: val = 2; break;
Chris@216: 	case 2000: val = 3; break;
Chris@216: 	case 4000: val = 4; break;
Chris@216: 	case 6000: val = 5; break;
Chris@216: 	case 8000: val = 6; break;
Chris@216: 	case 12000: val = 7; break;
Chris@216: 	case 16000: val = 8; break;
Chris@216: 	default: val = 9; break;
Chris@0: 	}
Chris@0: 
Chris@87:     } else if (name == "Frequency Scale") {
Chris@0: 
Chris@0: 	*min = 0;
Chris@0: 	*max = 1;
Chris@1103:         *deflt = int(BinScale::Linear);
Chris@1093: 	val = (int)m_binScale;
Chris@0: 
Chris@87:     } else if (name == "Bin Display") {
Chris@35: 
Chris@35: 	*min = 0;
Chris@35: 	*max = 2;
Chris@1103:         *deflt = int(BinDisplay::AllBins);
Chris@216: 	val = (int)m_binDisplay;
Chris@35: 
Chris@862:     } else if (name == "Normalization") {
Chris@36: 	
Chris@862:         *min = 0;
Chris@862:         *max = 3;
Chris@1104:         *deflt = 0;
Chris@1104:         
Chris@1104:         val = convertFromColumnNorm(m_normalization, m_normalizeVisibleArea);
Chris@120: 
Chris@0:     } else {
Chris@216: 	val = Layer::getPropertyRangeAndValue(name, min, max, deflt);
Chris@0:     }
Chris@0: 
Chris@216:     return val;
Chris@0: }
Chris@0: 
Chris@0: QString
Chris@0: SpectrogramLayer::getPropertyValueLabel(const PropertyName &name,
Chris@9: 					int value) const
Chris@0: {
Chris@87:     if (name == "Colour") {
Chris@196:         return ColourMapper::getColourMapName(value);
Chris@0:     }
Chris@87:     if (name == "Colour Scale") {
Chris@0: 	switch (value) {
Chris@0: 	default:
Chris@37: 	case 0: return tr("Linear");
Chris@37: 	case 1: return tr("Meter");
Chris@215: 	case 2: return tr("dBV^2");
Chris@215: 	case 3: return tr("dBV");
Chris@119: 	case 4: return tr("Phase");
Chris@0: 	}
Chris@0:     }
Chris@862:     if (name == "Normalization") {
Chris@862:         return ""; // icon only
Chris@862:     }
Chris@87:     if (name == "Window Size") {
Chris@0: 	return QString("%1").arg(32 << value);
Chris@0:     }
Chris@97:     if (name == "Window Increment") {
Chris@0: 	switch (value) {
Chris@0: 	default:
Chris@112: 	case 0: return tr("None");
Chris@112: 	case 1: return tr("25 %");
Chris@112: 	case 2: return tr("50 %");
Chris@112: 	case 3: return tr("75 %");
Chris@112: 	case 4: return tr("87.5 %");
Chris@112: 	case 5: return tr("93.75 %");
Chris@0: 	}
Chris@0:     }
Chris@87:     if (name == "Min Frequency") {
Chris@37: 	switch (value) {
Chris@37: 	default:
Chris@38: 	case 0: return tr("No min");
Chris@37: 	case 1: return tr("10 Hz");
Chris@37: 	case 2: return tr("20 Hz");
Chris@37: 	case 3: return tr("40 Hz");
Chris@37: 	case 4: return tr("100 Hz");
Chris@37: 	case 5: return tr("250 Hz");
Chris@37: 	case 6: return tr("500 Hz");
Chris@37: 	case 7: return tr("1 KHz");
Chris@37: 	case 8: return tr("4 KHz");
Chris@37: 	case 9: return tr("10 KHz");
Chris@37: 	}
Chris@37:     }
Chris@87:     if (name == "Max Frequency") {
Chris@0: 	switch (value) {
Chris@0: 	default:
Chris@0: 	case 0: return tr("500 Hz");
Chris@0: 	case 1: return tr("1 KHz");
Chris@0: 	case 2: return tr("1.5 KHz");
Chris@0: 	case 3: return tr("2 KHz");
Chris@0: 	case 4: return tr("4 KHz");
Chris@0: 	case 5: return tr("6 KHz");
Chris@0: 	case 6: return tr("8 KHz");
Chris@0: 	case 7: return tr("12 KHz");
Chris@0: 	case 8: return tr("16 KHz");
Chris@38: 	case 9: return tr("No max");
Chris@0: 	}
Chris@0:     }
Chris@87:     if (name == "Frequency Scale") {
Chris@0: 	switch (value) {
Chris@0: 	default:
Chris@0: 	case 0: return tr("Linear");
Chris@0: 	case 1: return tr("Log");
Chris@0: 	}
Chris@0:     }
Chris@87:     if (name == "Bin Display") {
Chris@35: 	switch (value) {
Chris@35: 	default:
Chris@37: 	case 0: return tr("All Bins");
Chris@37: 	case 1: return tr("Peak Bins");
Chris@37: 	case 2: return tr("Frequencies");
Chris@35: 	}
Chris@35:     }
Chris@0:     return tr("<unknown>");
Chris@0: }
Chris@0: 
Chris@862: QString
Chris@862: SpectrogramLayer::getPropertyValueIconName(const PropertyName &name,
Chris@862:                                            int value) const
Chris@862: {
Chris@862:     if (name == "Normalization") {
Chris@862:         switch(value) {
Chris@862:         default:
Chris@862:         case 0: return "normalise-none";
Chris@862:         case 1: return "normalise-columns";
Chris@862:         case 2: return "normalise";
Chris@862:         case 3: return "normalise-hybrid";
Chris@862:         }
Chris@862:     }
Chris@862:     return "";
Chris@862: }
Chris@862: 
Chris@167: RangeMapper *
Chris@167: SpectrogramLayer::getNewPropertyRangeMapper(const PropertyName &name) const
Chris@167: {
Chris@167:     if (name == "Gain") {
Chris@167:         return new LinearRangeMapper(-50, 50, -25, 25, tr("dB"));
Chris@167:     }
Chris@167:     if (name == "Threshold") {
Chris@167:         return new LinearRangeMapper(-50, 0, -50, 0, tr("dB"));
Chris@167:     }
Chris@167:     return 0;
Chris@167: }
Chris@167: 
Chris@0: void
Chris@0: SpectrogramLayer::setProperty(const PropertyName &name, int value)
Chris@0: {
Chris@87:     if (name == "Gain") {
Chris@906: 	setGain(float(pow(10, float(value)/20.0)));
Chris@87:     } else if (name == "Threshold") {
Chris@37: 	if (value == -50) setThreshold(0.0);
Chris@906: 	else setThreshold(float(AudioLevel::dB_to_multiplier(value)));
Chris@87:     } else if (name == "Colour Rotation") {
Chris@9: 	setColourRotation(value);
Chris@87:     } else if (name == "Colour") {
Chris@197:         setColourMap(value);
Chris@87:     } else if (name == "Window Size") {
Chris@0: 	setWindowSize(32 << value);
Chris@97:     } else if (name == "Window Increment") {
Chris@97:         setWindowHopLevel(value);
Chris@87:     } else if (name == "Min Frequency") {
Chris@37: 	switch (value) {
Chris@37: 	default:
Chris@37: 	case 0: setMinFrequency(0); break;
Chris@37: 	case 1: setMinFrequency(10); break;
Chris@37: 	case 2: setMinFrequency(20); break;
Chris@37: 	case 3: setMinFrequency(40); break;
Chris@37: 	case 4: setMinFrequency(100); break;
Chris@37: 	case 5: setMinFrequency(250); break;
Chris@37: 	case 6: setMinFrequency(500); break;
Chris@37: 	case 7: setMinFrequency(1000); break;
Chris@37: 	case 8: setMinFrequency(4000); break;
Chris@37: 	case 9: setMinFrequency(10000); break;
Chris@37: 	}
Chris@133:         int vs = getCurrentVerticalZoomStep();
Chris@133:         if (vs != m_lastEmittedZoomStep) {
Chris@133:             emit verticalZoomChanged();
Chris@133:             m_lastEmittedZoomStep = vs;
Chris@133:         }
Chris@87:     } else if (name == "Max Frequency") {
Chris@0: 	switch (value) {
Chris@0: 	case 0: setMaxFrequency(500); break;
Chris@0: 	case 1: setMaxFrequency(1000); break;
Chris@0: 	case 2: setMaxFrequency(1500); break;
Chris@0: 	case 3: setMaxFrequency(2000); break;
Chris@0: 	case 4: setMaxFrequency(4000); break;
Chris@0: 	case 5: setMaxFrequency(6000); break;
Chris@0: 	case 6: setMaxFrequency(8000); break;
Chris@0: 	case 7: setMaxFrequency(12000); break;
Chris@0: 	case 8: setMaxFrequency(16000); break;
Chris@0: 	default:
Chris@0: 	case 9: setMaxFrequency(0); break;
Chris@0: 	}
Chris@133:         int vs = getCurrentVerticalZoomStep();
Chris@133:         if (vs != m_lastEmittedZoomStep) {
Chris@133:             emit verticalZoomChanged();
Chris@133:             m_lastEmittedZoomStep = vs;
Chris@133:         }
Chris@87:     } else if (name == "Colour Scale") {
Chris@0: 	switch (value) {
Chris@0: 	default:
Chris@1105: 	case 0: setColourScale(ColourScaleType::Linear); break;
Chris@1105: 	case 1: setColourScale(ColourScaleType::Meter); break;
Chris@1105: 	case 2: setColourScale(ColourScaleType::Log); break; //!!! dB^2
Chris@1105: 	case 3: setColourScale(ColourScaleType::Log); break;
Chris@1105: 	case 4: setColourScale(ColourScaleType::Phase); break;
Chris@0: 	}
Chris@87:     } else if (name == "Frequency Scale") {
Chris@0: 	switch (value) {
Chris@0: 	default:
Chris@1103: 	case 0: setBinScale(BinScale::Linear); break;
Chris@1103: 	case 1: setBinScale(BinScale::Log); break;
Chris@0: 	}
Chris@87:     } else if (name == "Bin Display") {
Chris@35: 	switch (value) {
Chris@35: 	default:
Chris@1103: 	case 0: setBinDisplay(BinDisplay::AllBins); break;
Chris@1103: 	case 1: setBinDisplay(BinDisplay::PeakBins); break;
Chris@1103: 	case 2: setBinDisplay(BinDisplay::PeakFrequencies); break;
Chris@35: 	}
Chris@862:     } else if (name == "Normalization") {
Chris@1104:         auto n = convertToColumnNorm(value);
Chris@1104:         setNormalization(n.first);
Chris@1104:         setNormalizeVisibleArea(n.second);
Chris@0:     }
Chris@0: }
Chris@0: 
Chris@0: void
Chris@1106: SpectrogramLayer::invalidateRenderers()
Chris@95: {
Chris@1044: #ifdef DEBUG_SPECTROGRAM
Chris@1106:     cerr << "SpectrogramLayer::invalidateRenderers called" << endl;
Chris@1044: #endif
Chris@1106: 
Chris@1089:     for (ViewRendererMap::iterator i = m_renderers.begin();
Chris@1089:          i != m_renderers.end(); ++i) {
Chris@1089:         delete i->second;
Chris@1089:     }
Chris@1089:     m_renderers.clear();
Chris@95: }
Chris@95: 
Chris@95: void
Chris@122: SpectrogramLayer::preferenceChanged(PropertyContainer::PropertyName name)
Chris@122: {
Chris@587:     SVDEBUG << "SpectrogramLayer::preferenceChanged(" << name << ")" << endl;
Chris@122: 
Chris@122:     if (name == "Window Type") {
Chris@122:         setWindowType(Preferences::getInstance()->getWindowType());
Chris@122:         return;
Chris@122:     }
Chris@490:     if (name == "Spectrogram Y Smoothing") {
Chris@1086:         setWindowSize(m_windowSize);
Chris@1106:         invalidateRenderers();
Chris@490:         invalidateMagnitudes();
Chris@490:         emit layerParametersChanged();
Chris@490:     }
Chris@490:     if (name == "Spectrogram X Smoothing") {
Chris@1106:         invalidateRenderers();
Chris@122:         invalidateMagnitudes();
Chris@122:         emit layerParametersChanged();
Chris@122:     }
Chris@122:     if (name == "Tuning Frequency") {
Chris@122:         emit layerParametersChanged();
Chris@122:     }
Chris@122: }
Chris@122: 
Chris@122: void
Chris@0: SpectrogramLayer::setChannel(int ch)
Chris@0: {
Chris@0:     if (m_channel == ch) return;
Chris@0: 
Chris@1106:     invalidateRenderers();
Chris@0:     m_channel = ch;
Chris@1088:     invalidateFFTModel();
Chris@9: 
Chris@0:     emit layerParametersChanged();
Chris@0: }
Chris@0: 
Chris@0: int
Chris@0: SpectrogramLayer::getChannel() const
Chris@0: {
Chris@0:     return m_channel;
Chris@0: }
Chris@0: 
Chris@1086: int
Chris@1086: SpectrogramLayer::getFFTOversampling() const
Chris@1086: {
Chris@1103:     if (m_binDisplay != BinDisplay::AllBins) {
Chris@1086:         return 1;
Chris@1086:     }
Chris@1086: 
Chris@1086:     Preferences::SpectrogramSmoothing smoothing = 
Chris@1086:         Preferences::getInstance()->getSpectrogramSmoothing();
Chris@1086:     
Chris@1086:     if (smoothing == Preferences::NoSpectrogramSmoothing ||
Chris@1086:         smoothing == Preferences::SpectrogramInterpolated) {
Chris@1086:         return 1;
Chris@1086:     }
Chris@1086: 
Chris@1086:     return 4;
Chris@1086: }
Chris@1086: 
Chris@1087: int
Chris@1087: SpectrogramLayer::getFFTSize() const
Chris@1087: {
Chris@1087:     return m_windowSize * getFFTOversampling();
Chris@1087: }
Chris@1087: 
Chris@0: void
Chris@805: SpectrogramLayer::setWindowSize(int ws)
Chris@0: {
Chris@1087:     if (m_windowSize == ws) return;
Chris@0: 
Chris@1106:     invalidateRenderers();
Chris@0:     
Chris@0:     m_windowSize = ws;
Chris@0:     
Chris@1088:     invalidateFFTModel();
Chris@9: 
Chris@9:     emit layerParametersChanged();
Chris@0: }
Chris@0: 
Chris@805: int
Chris@0: SpectrogramLayer::getWindowSize() const
Chris@0: {
Chris@0:     return m_windowSize;
Chris@0: }
Chris@0: 
Chris@0: void
Chris@805: SpectrogramLayer::setWindowHopLevel(int v)
Chris@0: {
Chris@97:     if (m_windowHopLevel == v) return;
Chris@0: 
Chris@1106:     invalidateRenderers();
Chris@0:     
Chris@97:     m_windowHopLevel = v;
Chris@0:     
Chris@1088:     invalidateFFTModel();
Chris@9: 
Chris@9:     emit layerParametersChanged();
Chris@9: 
Chris@110: //    fillCache();
Chris@0: }
Chris@0: 
Chris@805: int
Chris@97: SpectrogramLayer::getWindowHopLevel() const
Chris@0: {
Chris@97:     return m_windowHopLevel;
Chris@0: }
Chris@0: 
Chris@0: void
Chris@0: SpectrogramLayer::setWindowType(WindowType w)
Chris@0: {
Chris@0:     if (m_windowType == w) return;
Chris@0: 
Chris@1106:     invalidateRenderers();
Chris@0:     
Chris@0:     m_windowType = w;
Chris@110: 
Chris@1088:     invalidateFFTModel();
Chris@9: 
Chris@9:     emit layerParametersChanged();
Chris@0: }
Chris@0: 
Chris@0: WindowType
Chris@0: SpectrogramLayer::getWindowType() const
Chris@0: {
Chris@0:     return m_windowType;
Chris@0: }
Chris@0: 
Chris@0: void
Chris@0: SpectrogramLayer::setGain(float gain)
Chris@0: {
Chris@587: //    SVDEBUG << "SpectrogramLayer::setGain(" << gain << ") (my gain is now "
Chris@585: //	      << m_gain << ")" << endl;
Chris@55: 
Chris@40:     if (m_gain == gain) return;
Chris@0: 
Chris@1106:     invalidateRenderers();
Chris@0:     
Chris@0:     m_gain = gain;
Chris@0:     
Chris@9:     emit layerParametersChanged();
Chris@0: }
Chris@0: 
Chris@0: float
Chris@0: SpectrogramLayer::getGain() const
Chris@0: {
Chris@0:     return m_gain;
Chris@0: }
Chris@0: 
Chris@0: void
Chris@37: SpectrogramLayer::setThreshold(float threshold)
Chris@37: {
Chris@40:     if (m_threshold == threshold) return;
Chris@37: 
Chris@1106:     invalidateRenderers();
Chris@37:     
Chris@37:     m_threshold = threshold;
Chris@37: 
Chris@37:     emit layerParametersChanged();
Chris@37: }
Chris@37: 
Chris@37: float
Chris@37: SpectrogramLayer::getThreshold() const
Chris@37: {
Chris@37:     return m_threshold;
Chris@37: }
Chris@37: 
Chris@37: void
Chris@805: SpectrogramLayer::setMinFrequency(int mf)
Chris@37: {
Chris@37:     if (m_minFrequency == mf) return;
Chris@37: 
Chris@587: //    SVDEBUG << "SpectrogramLayer::setMinFrequency: " << mf << endl;
Chris@187: 
Chris@1106:     invalidateRenderers();
Chris@119:     invalidateMagnitudes();
Chris@37:     
Chris@37:     m_minFrequency = mf;
Chris@37: 
Chris@37:     emit layerParametersChanged();
Chris@37: }
Chris@37: 
Chris@805: int
Chris@37: SpectrogramLayer::getMinFrequency() const
Chris@37: {
Chris@37:     return m_minFrequency;
Chris@37: }
Chris@37: 
Chris@37: void
Chris@805: SpectrogramLayer::setMaxFrequency(int mf)
Chris@0: {
Chris@0:     if (m_maxFrequency == mf) return;
Chris@0: 
Chris@587: //    SVDEBUG << "SpectrogramLayer::setMaxFrequency: " << mf << endl;
Chris@187: 
Chris@1106:     invalidateRenderers();
Chris@119:     invalidateMagnitudes();
Chris@0:     
Chris@0:     m_maxFrequency = mf;
Chris@0:     
Chris@9:     emit layerParametersChanged();
Chris@0: }
Chris@0: 
Chris@805: int
Chris@0: SpectrogramLayer::getMaxFrequency() const
Chris@0: {
Chris@0:     return m_maxFrequency;
Chris@0: }
Chris@0: 
Chris@0: void
Chris@9: SpectrogramLayer::setColourRotation(int r)
Chris@9: {
Chris@9:     if (r < 0) r = 0;
Chris@9:     if (r > 256) r = 256;
Chris@9:     int distance = r - m_colourRotation;
Chris@9: 
Chris@9:     if (distance != 0) {
Chris@1106: //!!!	rotatePalette(-distance);
Chris@9: 	m_colourRotation = r;
Chris@9:     }
Chris@9:     
Chris@1112:     invalidateRenderers(); //!!! in theory we shouldn't have to do this... but...
Chris@1112:     
Chris@9:     emit layerParametersChanged();
Chris@9: }
Chris@9: 
Chris@9: void
Chris@1105: SpectrogramLayer::setColourScale(ColourScaleType colourScale)
Chris@0: {
Chris@0:     if (m_colourScale == colourScale) return;
Chris@0: 
Chris@1106:     invalidateRenderers();
Chris@0:     
Chris@0:     m_colourScale = colourScale;
Chris@0:     
Chris@9:     emit layerParametersChanged();
Chris@0: }
Chris@0: 
Chris@1105: ColourScaleType
Chris@0: SpectrogramLayer::getColourScale() const
Chris@0: {
Chris@0:     return m_colourScale;
Chris@0: }
Chris@0: 
Chris@0: void
Chris@197: SpectrogramLayer::setColourMap(int map)
Chris@0: {
Chris@197:     if (m_colourMap == map) return;
Chris@0: 
Chris@1106:     invalidateRenderers();
Chris@0:     
Chris@197:     m_colourMap = map;
Chris@9: 
Chris@0:     emit layerParametersChanged();
Chris@0: }
Chris@0: 
Chris@196: int
Chris@197: SpectrogramLayer::getColourMap() const
Chris@0: {
Chris@197:     return m_colourMap;
Chris@0: }
Chris@0: 
Chris@0: void
Chris@1103: SpectrogramLayer::setBinScale(BinScale binScale)
Chris@0: {
Chris@1093:     if (m_binScale == binScale) return;
Chris@0: 
Chris@1106:     invalidateRenderers();
Chris@1093:     m_binScale = binScale;
Chris@9: 
Chris@9:     emit layerParametersChanged();
Chris@0: }
Chris@0: 
Chris@1103: BinScale
Chris@1093: SpectrogramLayer::getBinScale() const
Chris@0: {
Chris@1093:     return m_binScale;
Chris@0: }
Chris@0: 
Chris@0: void
Chris@1103: SpectrogramLayer::setBinDisplay(BinDisplay binDisplay)
Chris@35: {
Chris@37:     if (m_binDisplay == binDisplay) return;
Chris@35: 
Chris@1106:     invalidateRenderers();
Chris@37:     m_binDisplay = binDisplay;
Chris@35: 
Chris@35:     emit layerParametersChanged();
Chris@35: }
Chris@35: 
Chris@1103: BinDisplay
Chris@37: SpectrogramLayer::getBinDisplay() const
Chris@35: {
Chris@37:     return m_binDisplay;
Chris@35: }
Chris@35: 
Chris@35: void
Chris@1104: SpectrogramLayer::setNormalization(ColumnNormalization n)
Chris@36: {
Chris@862:     if (m_normalization == n) return;
Chris@36: 
Chris@1106:     invalidateRenderers();
Chris@119:     invalidateMagnitudes();
Chris@862:     m_normalization = n;
Chris@36: 
Chris@36:     emit layerParametersChanged();
Chris@36: }
Chris@36: 
Chris@1104: ColumnNormalization
Chris@862: SpectrogramLayer::getNormalization() const
Chris@36: {
Chris@862:     return m_normalization;
Chris@36: }
Chris@36: 
Chris@36: void
Chris@1104: SpectrogramLayer::setNormalizeVisibleArea(bool n)
Chris@1104: {
Chris@1104:     if (m_normalizeVisibleArea == n) return;
Chris@1104: 
Chris@1106:     invalidateRenderers();
Chris@1104:     invalidateMagnitudes();
Chris@1104:     m_normalizeVisibleArea = n;
Chris@1104:     
Chris@1104:     emit layerParametersChanged();
Chris@1104: }
Chris@1104: 
Chris@1104: bool
Chris@1104: SpectrogramLayer::getNormalizeVisibleArea() const
Chris@1104: {
Chris@1104:     return m_normalizeVisibleArea;
Chris@1104: }
Chris@1104: 
Chris@1104: void
Chris@918: SpectrogramLayer::setLayerDormant(const LayerGeometryProvider *v, bool dormant)
Chris@29: {
Chris@33:     if (dormant) {
Chris@33: 
Chris@331: #ifdef DEBUG_SPECTROGRAM_REPAINT
Chris@985:         cerr << "SpectrogramLayer::setLayerDormant(" << dormant << ")"
Chris@585:                   << endl;
Chris@331: #endif
Chris@331: 
Chris@131:         if (isLayerDormant(v)) {
Chris@131:             return;
Chris@131:         }
Chris@131: 
Chris@131:         Layer::setLayerDormant(v, true);
Chris@33: 
Chris@1106: 	invalidateRenderers();
Chris@1030: 
Chris@1088:         //!!! in theory we should call invalidateFFTModel() if and
Chris@1088:         //!!! only if there are no remaining views in which we are not
Chris@1088:         //!!! dormant
Chris@33: 	
Chris@33:     } else {
Chris@33: 
Chris@131:         Layer::setLayerDormant(v, false);
Chris@33:     }
Chris@29: }
Chris@29: 
Chris@29: void
Chris@0: SpectrogramLayer::cacheInvalid()
Chris@0: {
Chris@391: #ifdef DEBUG_SPECTROGRAM_REPAINT
Chris@985:     cerr << "SpectrogramLayer::cacheInvalid()" << endl;
Chris@391: #endif
Chris@391: 
Chris@1106:     invalidateRenderers();
Chris@119:     invalidateMagnitudes();
Chris@0: }
Chris@0: 
Chris@0: void
Chris@1037: SpectrogramLayer::cacheInvalid(
Chris@1037: #ifdef DEBUG_SPECTROGRAM_REPAINT
Chris@1037:     sv_frame_t from, sv_frame_t to
Chris@1037: #else 
Chris@1037:     sv_frame_t     , sv_frame_t
Chris@1037: #endif
Chris@1037:     )
Chris@0: {
Chris@391: #ifdef DEBUG_SPECTROGRAM_REPAINT
Chris@985:     cerr << "SpectrogramLayer::cacheInvalid(" << from << ", " << to << ")" << endl;
Chris@391: #endif
Chris@391: 
Chris@1030:     // We used to call invalidateMagnitudes(from, to) to invalidate
Chris@1030:     // only those caches whose views contained some of the (from, to)
Chris@1030:     // range. That's the right thing to do; it has been lost in
Chris@1030:     // pulling out the image cache code, but it might not matter very
Chris@1030:     // much, since the underlying models for spectrogram layers don't
Chris@1030:     // change very often. Let's see.
Chris@1106:     invalidateRenderers();
Chris@391:     invalidateMagnitudes();
Chris@0: }
Chris@0: 
Chris@224: bool
Chris@224: SpectrogramLayer::hasLightBackground() const 
Chris@224: {
Chris@287:     return ColourMapper(m_colourMap, 1.f, 255.f).hasLightBackground();
Chris@224: }
Chris@224: 
Chris@905: double
Chris@40: SpectrogramLayer::getEffectiveMinFrequency() const
Chris@40: {
Chris@907:     sv_samplerate_t sr = m_model->getSampleRate();
Chris@1087:     double minf = double(sr) / getFFTSize();
Chris@40: 
Chris@40:     if (m_minFrequency > 0.0) {
Chris@1087: 	int minbin = int((double(m_minFrequency) * getFFTSize()) / sr + 0.01);
Chris@40: 	if (minbin < 1) minbin = 1;
Chris@1087: 	minf = minbin * sr / getFFTSize();
Chris@40:     }
Chris@40: 
Chris@40:     return minf;
Chris@40: }
Chris@40: 
Chris@905: double
Chris@40: SpectrogramLayer::getEffectiveMaxFrequency() const
Chris@40: {
Chris@907:     sv_samplerate_t sr = m_model->getSampleRate();
Chris@905:     double maxf = double(sr) / 2;
Chris@40: 
Chris@40:     if (m_maxFrequency > 0.0) {
Chris@1087: 	int maxbin = int((double(m_maxFrequency) * getFFTSize()) / sr + 0.1);
Chris@1087: 	if (maxbin > getFFTSize() / 2) maxbin = getFFTSize() / 2;
Chris@1087: 	maxf = maxbin * sr / getFFTSize();
Chris@40:     }
Chris@40: 
Chris@40:     return maxf;
Chris@40: }
Chris@40: 
Chris@0: bool
Chris@918: SpectrogramLayer::getYBinRange(LayerGeometryProvider *v, int y, double &q0, double &q1) const
Chris@0: {
Chris@382:     Profiler profiler("SpectrogramLayer::getYBinRange");
Chris@382:     
Chris@918:     int h = v->getPaintHeight();
Chris@0:     if (y < 0 || y >= h) return false;
Chris@0: 
Chris@907:     sv_samplerate_t sr = m_model->getSampleRate();
Chris@905:     double minf = getEffectiveMinFrequency();
Chris@905:     double maxf = getEffectiveMaxFrequency();
Chris@0: 
Chris@1103:     bool logarithmic = (m_binScale == BinScale::Log);
Chris@38: 
Chris@44:     q0 = v->getFrequencyForY(y, minf, maxf, logarithmic);
Chris@44:     q1 = v->getFrequencyForY(y - 1, minf, maxf, logarithmic);
Chris@38: 
Chris@1086:     // Now map these on to ("proportions of") actual bins
Chris@490: 
Chris@1087:     q0 = (q0 * getFFTSize()) / sr;
Chris@1087:     q1 = (q1 * getFFTSize()) / sr;
Chris@0: 
Chris@0:     return true;
Chris@0: }
Chris@486: 
Chris@1085: double
Chris@1113: SpectrogramLayer::getYForBin(const LayerGeometryProvider *, double) const {
Chris@1115:     //!!! not implemented -- needed for non-opaque render
Chris@1115: 
Chris@1115:     //!!! this will crash, as it stands, when a spectrogram without
Chris@1115:     //!!! smoothing is zoomed in
Chris@1115:     
Chris@1115:     //!!! overlap with range methods
Chris@1085:     throw std::logic_error("not implemented");
Chris@1085: }
Chris@1085: 
Chris@1085: double
Chris@1113: SpectrogramLayer::getBinForY(const LayerGeometryProvider *v, double y) const
Chris@1085: {
Chris@1085:     //!!! overlap with range methods above (but using double arg)
Chris@1085:     //!!! tidy this
Chris@38:     
Chris@1085:     int h = v->getPaintHeight();
Chris@1085:     if (y < 0 || y >= h) return false;
Chris@1085: 
Chris@1085:     sv_samplerate_t sr = m_model->getSampleRate();
Chris@1085:     double minf = getEffectiveMinFrequency();
Chris@1085:     double maxf = getEffectiveMaxFrequency();
Chris@1085: 
Chris@1103:     bool logarithmic = (m_binScale == BinScale::Log);
Chris@1085: 
Chris@1085:     double q = v->getFrequencyForY(y, minf, maxf, logarithmic);
Chris@1085: 
Chris@1086:     // Now map on to ("proportions of") actual bins
Chris@1085: 
Chris@1087:     q = (q * getFFTSize()) / sr;
Chris@1085: 
Chris@1085:     return q;
Chris@1085: }
Chris@1085: 
Chris@0: bool
Chris@918: SpectrogramLayer::getXBinRange(LayerGeometryProvider *v, int x, double &s0, double &s1) const
Chris@0: {
Chris@907:     sv_frame_t modelStart = m_model->getStartFrame();
Chris@907:     sv_frame_t modelEnd = m_model->getEndFrame();
Chris@0: 
Chris@0:     // Each pixel column covers an exact range of sample frames:
Chris@907:     sv_frame_t f0 = v->getFrameForX(x) - modelStart;
Chris@907:     sv_frame_t f1 = v->getFrameForX(x + 1) - modelStart - 1;
Chris@20: 
Chris@41:     if (f1 < int(modelStart) || f0 > int(modelEnd)) {
Chris@41: 	return false;
Chris@41:     }
Chris@20:       
Chris@0:     // And that range may be drawn from a possibly non-integral
Chris@0:     // range of spectrogram windows:
Chris@0: 
Chris@805:     int windowIncrement = getWindowIncrement();
Chris@905:     s0 = double(f0) / windowIncrement;
Chris@905:     s1 = double(f1) / windowIncrement;
Chris@0: 
Chris@0:     return true;
Chris@0: }
Chris@0:  
Chris@0: bool
Chris@918: SpectrogramLayer::getXBinSourceRange(LayerGeometryProvider *v, int x, RealTime &min, RealTime &max) const
Chris@0: {
Chris@905:     double s0 = 0, s1 = 0;
Chris@44:     if (!getXBinRange(v, x, s0, s1)) return false;
Chris@0:     
Chris@0:     int s0i = int(s0 + 0.001);
Chris@0:     int s1i = int(s1);
Chris@0: 
Chris@0:     int windowIncrement = getWindowIncrement();
Chris@0:     int w0 = s0i * windowIncrement - (m_windowSize - windowIncrement)/2;
Chris@0:     int w1 = s1i * windowIncrement + windowIncrement +
Chris@0: 	(m_windowSize - windowIncrement)/2 - 1;
Chris@0:     
Chris@0:     min = RealTime::frame2RealTime(w0, m_model->getSampleRate());
Chris@0:     max = RealTime::frame2RealTime(w1, m_model->getSampleRate());
Chris@0:     return true;
Chris@0: }
Chris@0: 
Chris@0: bool
Chris@918: SpectrogramLayer::getYBinSourceRange(LayerGeometryProvider *v, int y, double &freqMin, double &freqMax)
Chris@0: const
Chris@0: {
Chris@905:     double q0 = 0, q1 = 0;
Chris@44:     if (!getYBinRange(v, y, q0, q1)) return false;
Chris@0: 
Chris@0:     int q0i = int(q0 + 0.001);
Chris@0:     int q1i = int(q1);
Chris@0: 
Chris@907:     sv_samplerate_t sr = m_model->getSampleRate();
Chris@0: 
Chris@0:     for (int q = q0i; q <= q1i; ++q) {
Chris@1087: 	if (q == q0i) freqMin = (sr * q) / getFFTSize();
Chris@1087: 	if (q == q1i) freqMax = (sr * (q+1)) / getFFTSize();
Chris@0:     }
Chris@0:     return true;
Chris@0: }
Chris@35: 
Chris@35: bool
Chris@918: SpectrogramLayer::getAdjustedYBinSourceRange(LayerGeometryProvider *v, int x, int y,
Chris@905: 					     double &freqMin, double &freqMax,
Chris@905: 					     double &adjFreqMin, double &adjFreqMax)
Chris@35: const
Chris@35: {
Chris@277:     if (!m_model || !m_model->isOK() || !m_model->isReady()) {
Chris@277: 	return false;
Chris@277:     }
Chris@277: 
Chris@1088:     FFTModel *fft = getFFTModel();
Chris@114:     if (!fft) return false;
Chris@110: 
Chris@905:     double s0 = 0, s1 = 0;
Chris@44:     if (!getXBinRange(v, x, s0, s1)) return false;
Chris@35: 
Chris@905:     double q0 = 0, q1 = 0;
Chris@44:     if (!getYBinRange(v, y, q0, q1)) return false;
Chris@35: 
Chris@35:     int s0i = int(s0 + 0.001);
Chris@35:     int s1i = int(s1);
Chris@35: 
Chris@35:     int q0i = int(q0 + 0.001);
Chris@35:     int q1i = int(q1);
Chris@35: 
Chris@907:     sv_samplerate_t sr = m_model->getSampleRate();
Chris@35: 
Chris@35:     bool haveAdj = false;
Chris@35: 
Chris@1103:     bool peaksOnly = (m_binDisplay == BinDisplay::PeakBins ||
Chris@1103: 		      m_binDisplay == BinDisplay::PeakFrequencies);
Chris@37: 
Chris@35:     for (int q = q0i; q <= q1i; ++q) {
Chris@35: 
Chris@35: 	for (int s = s0i; s <= s1i; ++s) {
Chris@35: 
Chris@905: 	    double binfreq = (double(sr) * q) / m_windowSize;
Chris@35: 	    if (q == q0i) freqMin = binfreq;
Chris@35: 	    if (q == q1i) freqMax = binfreq;
Chris@37: 
Chris@114: 	    if (peaksOnly && !fft->isLocalPeak(s, q)) continue;
Chris@38: 
Chris@1086: 	    if (!fft->isOverThreshold
Chris@1087:                 (s, q, float(m_threshold * double(getFFTSize())/2.0))) {
Chris@1086:                 continue;
Chris@1086:             }
Chris@907: 
Chris@907:             double freq = binfreq;
Chris@40: 	    
Chris@114: 	    if (s < int(fft->getWidth()) - 1) {
Chris@38: 
Chris@277:                 fft->estimateStableFrequency(s, q, freq);
Chris@35: 	    
Chris@38: 		if (!haveAdj || freq < adjFreqMin) adjFreqMin = freq;
Chris@38: 		if (!haveAdj || freq > adjFreqMax) adjFreqMax = freq;
Chris@35: 
Chris@35: 		haveAdj = true;
Chris@35: 	    }
Chris@35: 	}
Chris@35:     }
Chris@35: 
Chris@35:     if (!haveAdj) {
Chris@40: 	adjFreqMin = adjFreqMax = 0.0;
Chris@35:     }
Chris@35: 
Chris@35:     return haveAdj;
Chris@35: }
Chris@0:     
Chris@0: bool
Chris@918: SpectrogramLayer::getXYBinSourceRange(LayerGeometryProvider *v, int x, int y,
Chris@905: 				      double &min, double &max,
Chris@905: 				      double &phaseMin, double &phaseMax) const
Chris@0: {
Chris@277:     if (!m_model || !m_model->isOK() || !m_model->isReady()) {
Chris@277: 	return false;
Chris@277:     }
Chris@277: 
Chris@905:     double q0 = 0, q1 = 0;
Chris@44:     if (!getYBinRange(v, y, q0, q1)) return false;
Chris@0: 
Chris@905:     double s0 = 0, s1 = 0;
Chris@44:     if (!getXBinRange(v, x, s0, s1)) return false;
Chris@0:     
Chris@0:     int q0i = int(q0 + 0.001);
Chris@0:     int q1i = int(q1);
Chris@0: 
Chris@0:     int s0i = int(s0 + 0.001);
Chris@0:     int s1i = int(s1);
Chris@0: 
Chris@37:     bool rv = false;
Chris@37: 
Chris@1088:     FFTModel *fft = getFFTModel();
Chris@0: 
Chris@114:     if (fft) {
Chris@114: 
Chris@114:         int cw = fft->getWidth();
Chris@114:         int ch = fft->getHeight();
Chris@0: 
Chris@110:         min = 0.0;
Chris@110:         max = 0.0;
Chris@110:         phaseMin = 0.0;
Chris@110:         phaseMax = 0.0;
Chris@110:         bool have = false;
Chris@0: 
Chris@110:         for (int q = q0i; q <= q1i; ++q) {
Chris@110:             for (int s = s0i; s <= s1i; ++s) {
Chris@110:                 if (s >= 0 && q >= 0 && s < cw && q < ch) {
Chris@117: 
Chris@905:                     double value;
Chris@38: 
Chris@114:                     value = fft->getPhaseAt(s, q);
Chris@110:                     if (!have || value < phaseMin) { phaseMin = value; }
Chris@110:                     if (!have || value > phaseMax) { phaseMax = value; }
Chris@91: 
Chris@1087:                     value = fft->getMagnitudeAt(s, q) / (getFFTSize()/2.0);
Chris@110:                     if (!have || value < min) { min = value; }
Chris@110:                     if (!have || value > max) { max = value; }
Chris@110:                     
Chris@110:                     have = true;
Chris@110:                 }	
Chris@110:             }
Chris@110:         }
Chris@110:         
Chris@110:         if (have) {
Chris@110:             rv = true;
Chris@110:         }
Chris@0:     }
Chris@0: 
Chris@37:     return rv;
Chris@0: }
Chris@114: 	
Chris@130: FFTModel *
Chris@1088: SpectrogramLayer::getFFTModel() const
Chris@114: {
Chris@114:     if (!m_model) return 0;
Chris@114: 
Chris@1087:     int fftSize = getFFTSize();
Chris@114: 
Chris@1088:     //!!! it is now surely slower to do this on every getFFTModel()
Chris@1088:     //!!! request than it would be to recreate the model immediately
Chris@1088:     //!!! when something changes instead of just invalidating it
Chris@920:     
Chris@1088:     if (m_fftModel &&
Chris@1088:         m_fftModel->getHeight() == fftSize / 2 + 1 &&
Chris@1088:         m_fftModel->getWindowIncrement() == getWindowIncrement()) {
Chris@1088:         return m_fftModel;
Chris@114:     }
Chris@1088:     
Chris@1088:     delete m_peakCache;
Chris@1088:     m_peakCache = 0;
Chris@1088: 
Chris@1088:     delete m_fftModel;
Chris@1088:     m_fftModel = new FFTModel(m_model,
Chris@1088:                               m_channel,
Chris@1088:                               m_windowType,
Chris@1088:                               m_windowSize,
Chris@1088:                               getWindowIncrement(),
Chris@1088:                               fftSize);
Chris@1088: 
Chris@1088:     if (!m_fftModel->isOK()) {
Chris@1088:         QMessageBox::critical
Chris@1088:             (0, tr("FFT cache failed"),
Chris@1088:              tr("Failed to create the FFT model for this spectrogram.\n"
Chris@1088:                 "There may be insufficient memory or disc space to continue."));
Chris@1088:         delete m_fftModel;
Chris@1088:         m_fftModel = 0;
Chris@1088:         return 0;
Chris@114:     }
Chris@114: 
Chris@1088:     ((SpectrogramLayer *)this)->sliceableModelReplaced(0, m_fftModel);
Chris@1088: 
Chris@1088:     return m_fftModel;
Chris@114: }
Chris@114: 
Chris@484: Dense3DModelPeakCache *
Chris@1088: SpectrogramLayer::getPeakCache() const
Chris@484: {
Chris@1088:     //!!! see comment in getFFTModel
Chris@1088:     
Chris@1088:     if (!m_peakCache) {
Chris@1088:         FFTModel *f = getFFTModel();
Chris@484:         if (!f) return 0;
Chris@1088:         m_peakCache = new Dense3DModelPeakCache(f, m_peakCacheDivisor);
Chris@484:     }
Chris@1088:     return m_peakCache;
Chris@484: }
Chris@484: 
Chris@193: const Model *
Chris@193: SpectrogramLayer::getSliceableModel() const
Chris@193: {
Chris@1088:     return m_fftModel;
Chris@193: }
Chris@193: 
Chris@114: void
Chris@1088: SpectrogramLayer::invalidateFFTModel()
Chris@114: {
Chris@1044: #ifdef DEBUG_SPECTROGRAM
Chris@1088:     cerr << "SpectrogramLayer::invalidateFFTModel called" << endl;
Chris@1044: #endif
Chris@1088: 
Chris@1088:     emit sliceableModelReplaced(m_fftModel, 0);
Chris@1088: 
Chris@1088:     delete m_fftModel;
Chris@1088:     delete m_peakCache;
Chris@1088: 
Chris@1088:     m_fftModel = 0;
Chris@1088:     m_peakCache = 0;
Chris@114: }
Chris@114: 
Chris@0: void
Chris@119: SpectrogramLayer::invalidateMagnitudes()
Chris@119: {
Chris@1044: #ifdef DEBUG_SPECTROGRAM
Chris@1044:     cerr << "SpectrogramLayer::invalidateMagnitudes called" << endl;
Chris@1044: #endif
Chris@119:     m_viewMags.clear();
Chris@1025:     for (vector<MagnitudeRange>::iterator i = m_columnMags.begin();
Chris@119:          i != m_columnMags.end(); ++i) {
Chris@119:         *i = MagnitudeRange();
Chris@119:     }
Chris@119: }
Chris@119: 
Chris@119: bool
Chris@918: SpectrogramLayer::updateViewMagnitudes(LayerGeometryProvider *v) const
Chris@119: {
Chris@119:     MagnitudeRange mag;
Chris@119: 
Chris@918:     int x0 = 0, x1 = v->getPaintWidth();
Chris@905:     double s00 = 0, s01 = 0, s10 = 0, s11 = 0;
Chris@119:     
Chris@203:     if (!getXBinRange(v, x0, s00, s01)) {
Chris@907:         s00 = s01 = double(m_model->getStartFrame()) / getWindowIncrement();
Chris@203:     }
Chris@203: 
Chris@203:     if (!getXBinRange(v, x1, s10, s11)) {
Chris@907:         s10 = s11 = double(m_model->getEndFrame()) / getWindowIncrement();
Chris@203:     }
Chris@119: 
Chris@1025:     int s0 = int(min(s00, s10) + 0.0001);
Chris@1025:     int s1 = int(max(s01, s11) + 0.0001);
Chris@203: 
Chris@587: //    SVDEBUG << "SpectrogramLayer::updateViewMagnitudes: x0 = " << x0 << ", x1 = " << x1 << ", s00 = " << s00 << ", s11 = " << s11 << " s0 = " << s0 << ", s1 = " << s1 << endl;
Chris@119: 
Chris@248:     if (int(m_columnMags.size()) <= s1) {
Chris@119:         m_columnMags.resize(s1 + 1);
Chris@119:     }
Chris@119: 
Chris@119:     for (int s = s0; s <= s1; ++s) {
Chris@119:         if (m_columnMags[s].isSet()) {
Chris@119:             mag.sample(m_columnMags[s]);
Chris@119:         }
Chris@119:     }
Chris@119: 
Chris@184: #ifdef DEBUG_SPECTROGRAM_REPAINT
Chris@985:     cerr << "SpectrogramLayer::updateViewMagnitudes returning from cols "
Chris@1044:          << s0 << " -> " << s1 << " inclusive" << endl;
Chris@1044:     cerr << "SpectrogramLayer::updateViewMagnitudes: for view id " << v->getId()
Chris@1044:          << ": min is " << mag.getMin() << ", max is " << mag.getMax() << endl;
Chris@184: #endif
Chris@119: 
Chris@119:     if (!mag.isSet()) return false;
Chris@1030:     if (mag == m_viewMags[v->getId()]) return false;
Chris@1030:     m_viewMags[v->getId()] = mag;
Chris@119:     return true;
Chris@119: }
Chris@119: 
Chris@119: void
Chris@389: SpectrogramLayer::setSynchronousPainting(bool synchronous)
Chris@389: {
Chris@389:     m_synchronous = synchronous;
Chris@389: }
Chris@389: 
Chris@1089: Colour3DPlotRenderer *
Chris@1089: SpectrogramLayer::getRenderer(LayerGeometryProvider *v) const
Chris@1089: {
Chris@1089:     if (m_renderers.find(v->getId()) == m_renderers.end()) {
Chris@1089: 
Chris@1089:         Colour3DPlotRenderer::Sources sources;
Chris@1089:         sources.verticalBinLayer = this;
Chris@1090:         sources.fft = getFFTModel();
Chris@1090:         sources.source = sources.fft;
Chris@1090:         sources.peaks = getPeakCache();
Chris@1089: 
Chris@1092:         ColourScale::Parameters cparams;
Chris@1092:         cparams.colourMap = m_colourMap;
Chris@1092:         cparams.scale = m_colourScale;
Chris@1093:         cparams.threshold = m_threshold;
Chris@1093:         cparams.gain = m_gain;
Chris@1093: 
Chris@1105:         if (m_colourScale != ColourScaleType::Phase &&
Chris@1104:             m_normalization == ColumnNormalization::None) {
Chris@1093:             cparams.gain *= 2.f / float(getFFTSize());
Chris@1093:         }
Chris@1089:         
Chris@1089:         Colour3DPlotRenderer::Parameters params;
Chris@1092:         params.colourScale = ColourScale(cparams);
Chris@1089:         params.normalization = m_normalization;
Chris@1093:         params.binDisplay = m_binDisplay;
Chris@1093:         params.binScale = m_binScale;
Chris@1115:         params.alwaysOpaque = false;
Chris@1093:         params.invertVertical = false;
Chris@1112:         params.colourRotation = m_colourRotation;
Chris@1093: 
Chris@1093:         Preferences::SpectrogramSmoothing smoothing = 
Chris@1093:             Preferences::getInstance()->getSpectrogramSmoothing();
Chris@1093:         params.interpolate = 
Chris@1093:             (smoothing == Preferences::SpectrogramInterpolated ||
Chris@1093:              smoothing == Preferences::SpectrogramZeroPaddedAndInterpolated);
Chris@1089: 
Chris@1089:         m_renderers[v->getId()] = new Colour3DPlotRenderer(sources, params);
Chris@1089:     }
Chris@1089: 
Chris@1089:     return m_renderers[v->getId()];
Chris@1089: }
Chris@1089: 
Chris@1089: void
Chris@1106: SpectrogramLayer::paintWithRenderer(LayerGeometryProvider *v, QPainter &paint, QRect rect) const
Chris@1106: {
Chris@1106:     static int depth = 0;
Chris@1106:     
Chris@1106:     Colour3DPlotRenderer *renderer = getRenderer(v);
Chris@1106: 
Chris@1106:     if (m_synchronous) {
Chris@1106:         (void)renderer->render(v, paint, rect);
Chris@1106:         return;
Chris@1106:     }
Chris@1106: 
Chris@1106:     ++depth;
Chris@1106:     cerr << "paint depth " << depth << endl;
Chris@1106:     
Chris@1106:     (void)renderer->renderTimeConstrained(v, paint, rect);
Chris@1106: 
Chris@1106:     //!!! + mag range
Chris@1106: 
Chris@1106:     QRect uncached = renderer->getLargestUncachedRect();
Chris@1106:     if (uncached.width() > 0) {
Chris@1106:         cerr << "updating rect at " << uncached.x() << " width "
Chris@1106:              << uncached.width() << endl;
Chris@1106:         v->updatePaintRect(uncached);
Chris@1106:     }
Chris@1106: 
Chris@1106:     cerr << "exiting paint depth " << depth << endl;
Chris@1106:     --depth;
Chris@1106: }
Chris@1106: 
Chris@1106: void
Chris@916: SpectrogramLayer::paint(LayerGeometryProvider *v, QPainter &paint, QRect rect) const
Chris@0: {
Chris@334:     Profiler profiler("SpectrogramLayer::paint", false);
Chris@334: 
Chris@0: #ifdef DEBUG_SPECTROGRAM_REPAINT
Chris@1026:     cerr << "SpectrogramLayer::paint() entering: m_model is " << m_model << ", zoom level is " << v->getZoomLevel() << endl;
Chris@95:     
Chris@1026:     cerr << "SpectrogramLayer::paint(): rect is " << rect.x() << "," << rect.y() << " " << rect.width() << "x" << rect.height() << endl;
Chris@0: #endif
Chris@95: 
Chris@0:     if (!m_model || !m_model->isOK() || !m_model->isReady()) {
Chris@0: 	return;
Chris@0:     }
Chris@0: 
Chris@47:     if (isLayerDormant(v)) {
Chris@587: 	SVDEBUG << "SpectrogramLayer::paint(): Layer is dormant, making it undormant again" << endl;
Chris@29:     }
Chris@29: 
Chris@1106:     paintWithRenderer(v, paint, rect);
Chris@480: }
Chris@477: 
Chris@121: void
Chris@918: SpectrogramLayer::illuminateLocalFeatures(LayerGeometryProvider *v, QPainter &paint) const
Chris@121: {
Chris@382:     Profiler profiler("SpectrogramLayer::illuminateLocalFeatures");
Chris@382: 
Chris@121:     QPoint localPos;
Chris@121:     if (!v->shouldIlluminateLocalFeatures(this, localPos) || !m_model) {
Chris@121:         return;
Chris@121:     }
Chris@121: 
Chris@682: //    cerr << "SpectrogramLayer: illuminateLocalFeatures("
Chris@682: //              << localPos.x() << "," << localPos.y() << ")" << endl;
Chris@121: 
Chris@905:     double s0, s1;
Chris@905:     double f0, f1;
Chris@121: 
Chris@121:     if (getXBinRange(v, localPos.x(), s0, s1) &&
Chris@121:         getYBinSourceRange(v, localPos.y(), f0, f1)) {
Chris@121:         
Chris@121:         int s0i = int(s0 + 0.001);
Chris@121:         int s1i = int(s1);
Chris@121:         
Chris@121:         int x0 = v->getXForFrame(s0i * getWindowIncrement());
Chris@121:         int x1 = v->getXForFrame((s1i + 1) * getWindowIncrement());
Chris@121: 
Chris@248:         int y1 = int(getYForFrequency(v, f1));
Chris@248:         int y0 = int(getYForFrequency(v, f0));
Chris@121:         
Chris@682: //        cerr << "SpectrogramLayer: illuminate "
Chris@682: //                  << x0 << "," << y1 << " -> " << x1 << "," << y0 << endl;
Chris@121:         
Chris@287:         paint.setPen(v->getForeground());
Chris@133: 
Chris@133:         //!!! should we be using paintCrosshairs for this?
Chris@133: 
Chris@121:         paint.drawRect(x0, y1, x1 - x0 + 1, y0 - y1 + 1);
Chris@121:     }
Chris@121: }
Chris@121: 
Chris@905: double
Chris@918: SpectrogramLayer::getYForFrequency(const LayerGeometryProvider *v, double frequency) const
Chris@42: {
Chris@44:     return v->getYForFrequency(frequency,
Chris@44: 			       getEffectiveMinFrequency(),
Chris@44: 			       getEffectiveMaxFrequency(),
Chris@1103: 			       m_binScale == BinScale::Log);
Chris@42: }
Chris@42: 
Chris@905: double
Chris@918: SpectrogramLayer::getFrequencyForY(const LayerGeometryProvider *v, int y) const
Chris@42: {
Chris@44:     return v->getFrequencyForY(y,
Chris@44: 			       getEffectiveMinFrequency(),
Chris@44: 			       getEffectiveMaxFrequency(),
Chris@1103: 			       m_binScale == BinScale::Log);
Chris@42: }
Chris@42: 
Chris@0: int
Chris@1090: SpectrogramLayer::getCompletion(LayerGeometryProvider *) const
Chris@0: {
Chris@1088:     if (!m_fftModel) return 100;
Chris@1088:     int completion = m_fftModel->getCompletion();
Chris@224: #ifdef DEBUG_SPECTROGRAM_REPAINT
Chris@985:     cerr << "SpectrogramLayer::getCompletion: completion = " << completion << endl;
Chris@224: #endif
Chris@0:     return completion;
Chris@0: }
Chris@0: 
Chris@583: QString
Chris@1090: SpectrogramLayer::getError(LayerGeometryProvider *) const
Chris@583: {
Chris@1088:     if (!m_fftModel) return "";
Chris@1088:     return m_fftModel->getError();
Chris@583: }
Chris@583: 
Chris@28: bool
Chris@905: SpectrogramLayer::getValueExtents(double &min, double &max,
Chris@101:                                   bool &logarithmic, QString &unit) const
Chris@79: {
Chris@133:     if (!m_model) return false;
Chris@133: 
Chris@907:     sv_samplerate_t sr = m_model->getSampleRate();
Chris@1087:     min = double(sr) / getFFTSize();
Chris@905:     max = double(sr) / 2;
Chris@133:     
Chris@1103:     logarithmic = (m_binScale == BinScale::Log);
Chris@79:     unit = "Hz";
Chris@79:     return true;
Chris@79: }
Chris@79: 
Chris@79: bool
Chris@905: SpectrogramLayer::getDisplayExtents(double &min, double &max) const
Chris@101: {
Chris@101:     min = getEffectiveMinFrequency();
Chris@101:     max = getEffectiveMaxFrequency();
Chris@253: 
Chris@587: //    SVDEBUG << "SpectrogramLayer::getDisplayExtents: " << min << "->" << max << endl;
Chris@101:     return true;
Chris@101: }    
Chris@101: 
Chris@101: bool
Chris@905: SpectrogramLayer::setDisplayExtents(double min, double max)
Chris@120: {
Chris@120:     if (!m_model) return false;
Chris@187: 
Chris@587: //    SVDEBUG << "SpectrogramLayer::setDisplayExtents: " << min << "->" << max << endl;
Chris@187: 
Chris@120:     if (min < 0) min = 0;
Chris@907:     if (max > m_model->getSampleRate()/2.0) max = m_model->getSampleRate()/2.0;
Chris@120:     
Chris@907:     int minf = int(lrint(min));
Chris@907:     int maxf = int(lrint(max));
Chris@120: 
Chris@120:     if (m_minFrequency == minf && m_maxFrequency == maxf) return true;
Chris@120: 
Chris@1106:     invalidateRenderers();
Chris@120:     invalidateMagnitudes();
Chris@120: 
Chris@120:     m_minFrequency = minf;
Chris@120:     m_maxFrequency = maxf;
Chris@120:     
Chris@120:     emit layerParametersChanged();
Chris@120: 
Chris@133:     int vs = getCurrentVerticalZoomStep();
Chris@133:     if (vs != m_lastEmittedZoomStep) {
Chris@133:         emit verticalZoomChanged();
Chris@133:         m_lastEmittedZoomStep = vs;
Chris@133:     }
Chris@133: 
Chris@120:     return true;
Chris@120: }
Chris@120: 
Chris@120: bool
Chris@918: SpectrogramLayer::getYScaleValue(const LayerGeometryProvider *v, int y,
Chris@905:                                  double &value, QString &unit) const
Chris@261: {
Chris@261:     value = getFrequencyForY(v, y);
Chris@261:     unit = "Hz";
Chris@261:     return true;
Chris@261: }
Chris@261: 
Chris@261: bool
Chris@918: SpectrogramLayer::snapToFeatureFrame(LayerGeometryProvider *,
Chris@907:                                      sv_frame_t &frame,
Chris@805: 				     int &resolution,
Chris@28: 				     SnapType snap) const
Chris@13: {
Chris@13:     resolution = getWindowIncrement();
Chris@907:     sv_frame_t left = (frame / resolution) * resolution;
Chris@907:     sv_frame_t right = left + resolution;
Chris@28: 
Chris@28:     switch (snap) {
Chris@28:     case SnapLeft:  frame = left;  break;
Chris@28:     case SnapRight: frame = right; break;
Chris@28:     case SnapNearest:
Chris@28:     case SnapNeighbouring:
Chris@28: 	if (frame - left > right - frame) frame = right;
Chris@28: 	else frame = left;
Chris@28: 	break;
Chris@28:     }
Chris@28:     
Chris@28:     return true;
Chris@28: } 
Chris@13: 
Chris@283: void
Chris@1106: SpectrogramLayer::measureDoubleClick(LayerGeometryProvider *, QMouseEvent *)
Chris@283: {
Chris@1106: /*!!! replace this
Chris@920:     const View *view = v->getView();
Chris@1106: 
Chris@1030:     ScrollableImageCache &cache = getImageCacheReference(view);
Chris@1030: 
Chris@1030:     cerr << "cache width: " << cache.getSize().width() << ", height: "
Chris@1030:          << cache.getSize().height() << endl;
Chris@1030: 
Chris@1030:     QImage image = cache.getImage();
Chris@283: 
Chris@283:     ImageRegionFinder finder;
Chris@283:     QRect rect = finder.findRegionExtents(&image, e->pos());
Chris@283:     if (rect.isValid()) {
Chris@283:         MeasureRect mr;
Chris@283:         setMeasureRectFromPixrect(v, mr, rect);
Chris@283:         CommandHistory::getInstance()->addCommand
Chris@283:             (new AddMeasurementRectCommand(this, mr));
Chris@283:     }
Chris@1106: */
Chris@283: }
Chris@283: 
Chris@77: bool
Chris@918: SpectrogramLayer::getCrosshairExtents(LayerGeometryProvider *v, QPainter &paint,
Chris@77:                                       QPoint cursorPos,
Chris@1025:                                       vector<QRect> &extents) const
Chris@77: {
Chris@918:     QRect vertical(cursorPos.x() - 12, 0, 12, v->getPaintHeight());
Chris@77:     extents.push_back(vertical);
Chris@77: 
Chris@77:     QRect horizontal(0, cursorPos.y(), cursorPos.x(), 1);
Chris@77:     extents.push_back(horizontal);
Chris@77: 
Chris@608:     int sw = getVerticalScaleWidth(v, m_haveDetailedScale, paint);
Chris@264: 
Chris@280:     QRect freq(sw, cursorPos.y() - paint.fontMetrics().ascent() - 2,
Chris@280:                paint.fontMetrics().width("123456 Hz") + 2,
Chris@280:                paint.fontMetrics().height());
Chris@280:     extents.push_back(freq);
Chris@264: 
Chris@279:     QRect pitch(sw, cursorPos.y() + 2,
Chris@279:                 paint.fontMetrics().width("C#10+50c") + 2,
Chris@279:                 paint.fontMetrics().height());
Chris@279:     extents.push_back(pitch);
Chris@279: 
Chris@280:     QRect rt(cursorPos.x(),
Chris@918:              v->getPaintHeight() - paint.fontMetrics().height() - 2,
Chris@280:              paint.fontMetrics().width("1234.567 s"),
Chris@280:              paint.fontMetrics().height());
Chris@280:     extents.push_back(rt);
Chris@280: 
Chris@280:     int w(paint.fontMetrics().width("1234567890") + 2);
Chris@280:     QRect frame(cursorPos.x() - w - 2,
Chris@918:                 v->getPaintHeight() - paint.fontMetrics().height() - 2,
Chris@280:                 w,
Chris@280:                 paint.fontMetrics().height());
Chris@280:     extents.push_back(frame);
Chris@280: 
Chris@77:     return true;
Chris@77: }
Chris@77: 
Chris@77: void
Chris@918: SpectrogramLayer::paintCrosshairs(LayerGeometryProvider *v, QPainter &paint,
Chris@77:                                   QPoint cursorPos) const
Chris@77: {
Chris@77:     paint.save();
Chris@283: 
Chris@608:     int sw = getVerticalScaleWidth(v, m_haveDetailedScale, paint);
Chris@283: 
Chris@282:     QFont fn = paint.font();
Chris@282:     if (fn.pointSize() > 8) {
Chris@282:         fn.setPointSize(fn.pointSize() - 1);
Chris@282:         paint.setFont(fn);
Chris@282:     }
Chris@77:     paint.setPen(m_crosshairColour);
Chris@77: 
Chris@77:     paint.drawLine(0, cursorPos.y(), cursorPos.x() - 1, cursorPos.y());
Chris@918:     paint.drawLine(cursorPos.x(), 0, cursorPos.x(), v->getPaintHeight());
Chris@77:     
Chris@905:     double fundamental = getFrequencyForY(v, cursorPos.y());
Chris@77: 
Chris@1078:     PaintAssistant::drawVisibleText(v, paint,
Chris@278:                        sw + 2,
Chris@278:                        cursorPos.y() - 2,
Chris@278:                        QString("%1 Hz").arg(fundamental),
Chris@1078:                        PaintAssistant::OutlinedText);
Chris@278: 
Chris@279:     if (Pitch::isFrequencyInMidiRange(fundamental)) {
Chris@279:         QString pitchLabel = Pitch::getPitchLabelForFrequency(fundamental);
Chris@1078:         PaintAssistant::drawVisibleText(v, paint,
Chris@279:                            sw + 2,
Chris@279:                            cursorPos.y() + paint.fontMetrics().ascent() + 2,
Chris@279:                            pitchLabel,
Chris@1078:                            PaintAssistant::OutlinedText);
Chris@279:     }
Chris@279: 
Chris@907:     sv_frame_t frame = v->getFrameForX(cursorPos.x());
Chris@279:     RealTime rt = RealTime::frame2RealTime(frame, m_model->getSampleRate());
Chris@280:     QString rtLabel = QString("%1 s").arg(rt.toText(true).c_str());
Chris@280:     QString frameLabel = QString("%1").arg(frame);
Chris@1078:     PaintAssistant::drawVisibleText(v, paint,
Chris@280:                        cursorPos.x() - paint.fontMetrics().width(frameLabel) - 2,
Chris@918:                        v->getPaintHeight() - 2,
Chris@280:                        frameLabel,
Chris@1078:                        PaintAssistant::OutlinedText);
Chris@1078:     PaintAssistant::drawVisibleText(v, paint,
Chris@280:                        cursorPos.x() + 2,
Chris@918:                        v->getPaintHeight() - 2,
Chris@280:                        rtLabel,
Chris@1078:                        PaintAssistant::OutlinedText);
Chris@264: 
Chris@77:     int harmonic = 2;
Chris@77: 
Chris@77:     while (harmonic < 100) {
Chris@77: 
Chris@907:         int hy = int(lrint(getYForFrequency(v, fundamental * harmonic)));
Chris@918:         if (hy < 0 || hy > v->getPaintHeight()) break;
Chris@77:         
Chris@77:         int len = 7;
Chris@77: 
Chris@77:         if (harmonic % 2 == 0) {
Chris@77:             if (harmonic % 4 == 0) {
Chris@77:                 len = 12;
Chris@77:             } else {
Chris@77:                 len = 10;
Chris@77:             }
Chris@77:         }
Chris@77: 
Chris@77:         paint.drawLine(cursorPos.x() - len,
Chris@907:                        hy,
Chris@77:                        cursorPos.x(),
Chris@907:                        hy);
Chris@77: 
Chris@77:         ++harmonic;
Chris@77:     }
Chris@77: 
Chris@77:     paint.restore();
Chris@77: }
Chris@77: 
Chris@25: QString
Chris@918: SpectrogramLayer::getFeatureDescription(LayerGeometryProvider *v, QPoint &pos) const
Chris@25: {
Chris@25:     int x = pos.x();
Chris@25:     int y = pos.y();
Chris@0: 
Chris@25:     if (!m_model || !m_model->isOK()) return "";
Chris@0: 
Chris@905:     double magMin = 0, magMax = 0;
Chris@905:     double phaseMin = 0, phaseMax = 0;
Chris@905:     double freqMin = 0, freqMax = 0;
Chris@905:     double adjFreqMin = 0, adjFreqMax = 0;
Chris@25:     QString pitchMin, pitchMax;
Chris@0:     RealTime rtMin, rtMax;
Chris@0: 
Chris@38:     bool haveValues = false;
Chris@0: 
Chris@44:     if (!getXBinSourceRange(v, x, rtMin, rtMax)) {
Chris@38: 	return "";
Chris@38:     }
Chris@44:     if (getXYBinSourceRange(v, x, y, magMin, magMax, phaseMin, phaseMax)) {
Chris@38: 	haveValues = true;
Chris@38:     }
Chris@0: 
Chris@35:     QString adjFreqText = "", adjPitchText = "";
Chris@35: 
Chris@1103:     if (m_binDisplay == BinDisplay::PeakFrequencies) {
Chris@35: 
Chris@44: 	if (!getAdjustedYBinSourceRange(v, x, y, freqMin, freqMax,
Chris@38: 					adjFreqMin, adjFreqMax)) {
Chris@38: 	    return "";
Chris@38: 	}
Chris@35: 
Chris@35: 	if (adjFreqMin != adjFreqMax) {
Chris@65: 	    adjFreqText = tr("Peak Frequency:\t%1 - %2 Hz\n")
Chris@35: 		.arg(adjFreqMin).arg(adjFreqMax);
Chris@35: 	} else {
Chris@65: 	    adjFreqText = tr("Peak Frequency:\t%1 Hz\n")
Chris@35: 		.arg(adjFreqMin);
Chris@38: 	}
Chris@38: 
Chris@38: 	QString pmin = Pitch::getPitchLabelForFrequency(adjFreqMin);
Chris@38: 	QString pmax = Pitch::getPitchLabelForFrequency(adjFreqMax);
Chris@38: 
Chris@38: 	if (pmin != pmax) {
Chris@65: 	    adjPitchText = tr("Peak Pitch:\t%3 - %4\n").arg(pmin).arg(pmax);
Chris@38: 	} else {
Chris@65: 	    adjPitchText = tr("Peak Pitch:\t%2\n").arg(pmin);
Chris@35: 	}
Chris@35: 
Chris@35:     } else {
Chris@35: 	
Chris@44: 	if (!getYBinSourceRange(v, y, freqMin, freqMax)) return "";
Chris@35:     }
Chris@35: 
Chris@25:     QString text;
Chris@25: 
Chris@25:     if (rtMin != rtMax) {
Chris@25: 	text += tr("Time:\t%1 - %2\n")
Chris@25: 	    .arg(rtMin.toText(true).c_str())
Chris@25: 	    .arg(rtMax.toText(true).c_str());
Chris@25:     } else {
Chris@25: 	text += tr("Time:\t%1\n")
Chris@25: 	    .arg(rtMin.toText(true).c_str());
Chris@0:     }
Chris@0: 
Chris@25:     if (freqMin != freqMax) {
Chris@65: 	text += tr("%1Bin Frequency:\t%2 - %3 Hz\n%4Bin Pitch:\t%5 - %6\n")
Chris@65: 	    .arg(adjFreqText)
Chris@25: 	    .arg(freqMin)
Chris@25: 	    .arg(freqMax)
Chris@65: 	    .arg(adjPitchText)
Chris@65: 	    .arg(Pitch::getPitchLabelForFrequency(freqMin))
Chris@65: 	    .arg(Pitch::getPitchLabelForFrequency(freqMax));
Chris@65:     } else {
Chris@65: 	text += tr("%1Bin Frequency:\t%2 Hz\n%3Bin Pitch:\t%4\n")
Chris@35: 	    .arg(adjFreqText)
Chris@25: 	    .arg(freqMin)
Chris@65: 	    .arg(adjPitchText)
Chris@65: 	    .arg(Pitch::getPitchLabelForFrequency(freqMin));
Chris@25:     }	
Chris@25: 
Chris@38:     if (haveValues) {
Chris@905: 	double dbMin = AudioLevel::multiplier_to_dB(magMin);
Chris@905: 	double dbMax = AudioLevel::multiplier_to_dB(magMax);
Chris@43: 	QString dbMinString;
Chris@43: 	QString dbMaxString;
Chris@43: 	if (dbMin == AudioLevel::DB_FLOOR) {
Chris@43: 	    dbMinString = tr("-Inf");
Chris@43: 	} else {
Chris@907: 	    dbMinString = QString("%1").arg(lrint(dbMin));
Chris@43: 	}
Chris@43: 	if (dbMax == AudioLevel::DB_FLOOR) {
Chris@43: 	    dbMaxString = tr("-Inf");
Chris@43: 	} else {
Chris@907: 	    dbMaxString = QString("%1").arg(lrint(dbMax));
Chris@43: 	}
Chris@907: 	if (lrint(dbMin) != lrint(dbMax)) {
Chris@199: 	    text += tr("dB:\t%1 - %2").arg(dbMinString).arg(dbMaxString);
Chris@25: 	} else {
Chris@199: 	    text += tr("dB:\t%1").arg(dbMinString);
Chris@25: 	}
Chris@38: 	if (phaseMin != phaseMax) {
Chris@38: 	    text += tr("\nPhase:\t%1 - %2").arg(phaseMin).arg(phaseMax);
Chris@38: 	} else {
Chris@38: 	    text += tr("\nPhase:\t%1").arg(phaseMin);
Chris@38: 	}
Chris@25:     }
Chris@25: 
Chris@25:     return text;
Chris@0: }
Chris@25: 
Chris@0: int
Chris@40: SpectrogramLayer::getColourScaleWidth(QPainter &paint) const
Chris@40: {
Chris@40:     int cw;
Chris@40: 
Chris@119:     cw = paint.fontMetrics().width("-80dB");
Chris@119: 
Chris@40:     return cw;
Chris@40: }
Chris@40: 
Chris@40: int
Chris@918: SpectrogramLayer::getVerticalScaleWidth(LayerGeometryProvider *, bool detailed, QPainter &paint) const
Chris@0: {
Chris@0:     if (!m_model || !m_model->isOK()) return 0;
Chris@0: 
Chris@607:     int cw = 0;
Chris@607:     if (detailed) cw = getColourScaleWidth(paint);
Chris@40: 
Chris@0:     int tw = paint.fontMetrics().width(QString("%1")
Chris@0: 				     .arg(m_maxFrequency > 0 ?
Chris@0: 					  m_maxFrequency - 1 :
Chris@0: 					  m_model->getSampleRate() / 2));
Chris@0: 
Chris@234:     int fw = paint.fontMetrics().width(tr("43Hz"));
Chris@0:     if (tw < fw) tw = fw;
Chris@40: 
Chris@1103:     int tickw = (m_binScale == BinScale::Log ? 10 : 4);
Chris@0:     
Chris@40:     return cw + tickw + tw + 13;
Chris@0: }
Chris@0: 
Chris@0: void
Chris@918: SpectrogramLayer::paintVerticalScale(LayerGeometryProvider *v, bool detailed, QPainter &paint, QRect rect) const
Chris@0: {
Chris@0:     if (!m_model || !m_model->isOK()) {
Chris@0: 	return;
Chris@0:     }
Chris@0: 
Chris@382:     Profiler profiler("SpectrogramLayer::paintVerticalScale");
Chris@122: 
Chris@120:     //!!! cache this?
Chris@120: 
Chris@0:     int h = rect.height(), w = rect.width();
Chris@0: 
Chris@1103:     int tickw = (m_binScale == BinScale::Log ? 10 : 4);
Chris@1103:     int pkw = (m_binScale == BinScale::Log ? 10 : 0);
Chris@40: 
Chris@1087:     int bins = getFFTSize() / 2;
Chris@907:     sv_samplerate_t sr = m_model->getSampleRate();
Chris@0: 
Chris@0:     if (m_maxFrequency > 0) {
Chris@1087: 	bins = int((double(m_maxFrequency) * getFFTSize()) / sr + 0.1);
Chris@1087: 	if (bins > getFFTSize() / 2) bins = getFFTSize() / 2;
Chris@0:     }
Chris@0: 
Chris@607:     int cw = 0;
Chris@607: 
Chris@607:     if (detailed) cw = getColourScaleWidth(paint);
Chris@119:     int cbw = paint.fontMetrics().width("dB");
Chris@40: 
Chris@0:     int py = -1;
Chris@0:     int textHeight = paint.fontMetrics().height();
Chris@0:     int toff = -textHeight + paint.fontMetrics().ascent() + 2;
Chris@0: 
Chris@607:     if (detailed && (h > textHeight * 3 + 10)) {
Chris@119: 
Chris@119:         int topLines = 2;
Chris@1105:         if (m_colourScale == ColourScaleType::Phase) topLines = 1;
Chris@119: 
Chris@119: 	int ch = h - textHeight * (topLines + 1) - 8;
Chris@119: //	paint.drawRect(4, textHeight + 4, cw - 1, ch + 1);
Chris@119: 	paint.drawRect(4 + cw - cbw, textHeight * topLines + 4, cbw - 1, ch + 1);
Chris@40: 
Chris@40: 	QString top, bottom;
Chris@1030:         double min = m_viewMags[v->getId()].getMin();
Chris@1030:         double max = m_viewMags[v->getId()].getMax();
Chris@905: 
Chris@905:         double dBmin = AudioLevel::multiplier_to_dB(min);
Chris@905:         double dBmax = AudioLevel::multiplier_to_dB(max);
Chris@119: 
Chris@1044: #ifdef DEBUG_SPECTROGRAM_REPAINT
Chris@1044:         cerr << "paintVerticalScale: for view id " << v->getId()
Chris@1044:              << ": min = " << min << ", max = " << max
Chris@1044:              << ", dBmin = " << dBmin << ", dBmax = " << dBmax << endl;
Chris@1044: #endif
Chris@1044:         
Chris@120:         if (dBmax < -60.f) dBmax = -60.f;
Chris@907:         else top = QString("%1").arg(lrint(dBmax));
Chris@120: 
Chris@120:         if (dBmin < dBmax - 60.f) dBmin = dBmax - 60.f;
Chris@907:         bottom = QString("%1").arg(lrint(dBmin));
Chris@119: 
Chris@119:         //!!! & phase etc
Chris@119: 
Chris@1105:         if (m_colourScale != ColourScaleType::Phase) {
Chris@119:             paint.drawText((cw + 6 - paint.fontMetrics().width("dBFS")) / 2,
Chris@119:                            2 + textHeight + toff, "dBFS");
Chris@119:         }
Chris@119: 
Chris@119: //	paint.drawText((cw + 6 - paint.fontMetrics().width(top)) / 2,
Chris@119: 	paint.drawText(3 + cw - cbw - paint.fontMetrics().width(top),
Chris@119: 		       2 + textHeight * topLines + toff + textHeight/2, top);
Chris@119: 
Chris@119: 	paint.drawText(3 + cw - cbw - paint.fontMetrics().width(bottom),
Chris@119: 		       h + toff - 3 - textHeight/2, bottom);
Chris@40: 
Chris@40: 	paint.save();
Chris@40: 	paint.setBrush(Qt::NoBrush);
Chris@119: 
Chris@119:         int lasty = 0;
Chris@119:         int lastdb = 0;
Chris@119: 
Chris@40: 	for (int i = 0; i < ch; ++i) {
Chris@119: 
Chris@905:             double dBval = dBmin + (((dBmax - dBmin) * i) / (ch - 1));
Chris@119:             int idb = int(dBval);
Chris@119: 
Chris@1106: //!!! replace this
Chris@1106:             // double value = AudioLevel::dB_to_multiplier(dBval);
Chris@1106:             // int colour = getDisplayValue(v, value * m_gain);
Chris@1106: 	    // paint.setPen(m_palette.getColour((unsigned char)colour));
Chris@119: 
Chris@119:             int y = textHeight * topLines + 4 + ch - i;
Chris@119: 
Chris@119:             paint.drawLine(5 + cw - cbw, y, cw + 2, y);
Chris@119: 
Chris@119:             if (i == 0) {
Chris@119:                 lasty = y;
Chris@119:                 lastdb = idb;
Chris@119:             } else if (i < ch - paint.fontMetrics().ascent() &&
Chris@120:                        idb != lastdb &&
Chris@119:                        ((abs(y - lasty) > textHeight && 
Chris@119:                          idb % 10 == 0) ||
Chris@119:                         (abs(y - lasty) > paint.fontMetrics().ascent() && 
Chris@119:                          idb % 5 == 0))) {
Chris@287:                 paint.setPen(v->getBackground());
Chris@119:                 QString text = QString("%1").arg(idb);
Chris@119:                 paint.drawText(3 + cw - cbw - paint.fontMetrics().width(text),
Chris@119:                                y + toff + textHeight/2, text);
Chris@287:                 paint.setPen(v->getForeground());
Chris@119:                 paint.drawLine(5 + cw - cbw, y, 8 + cw - cbw, y);
Chris@119:                 lasty = y;
Chris@119:                 lastdb = idb;
Chris@119:             }
Chris@40: 	}
Chris@40: 	paint.restore();
Chris@40:     }
Chris@40: 
Chris@40:     paint.drawLine(cw + 7, 0, cw + 7, h);
Chris@40: 
Chris@0:     int bin = -1;
Chris@0: 
Chris@918:     for (int y = 0; y < v->getPaintHeight(); ++y) {
Chris@0: 
Chris@905: 	double q0, q1;
Chris@918: 	if (!getYBinRange(v, v->getPaintHeight() - y, q0, q1)) continue;
Chris@0: 
Chris@0: 	int vy;
Chris@0: 
Chris@0: 	if (int(q0) > bin) {
Chris@0: 	    vy = y;
Chris@0: 	    bin = int(q0);
Chris@0: 	} else {
Chris@0: 	    continue;
Chris@0: 	}
Chris@0: 
Chris@1087: 	int freq = int((sr * bin) / getFFTSize());
Chris@0: 
Chris@0: 	if (py >= 0 && (vy - py) < textHeight - 1) {
Chris@1103: 	    if (m_binScale == BinScale::Linear) {
Chris@40: 		paint.drawLine(w - tickw, h - vy, w, h - vy);
Chris@40: 	    }
Chris@0: 	    continue;
Chris@0: 	}
Chris@0: 
Chris@0: 	QString text = QString("%1").arg(freq);
Chris@234: 	if (bin == 1) text = tr("%1Hz").arg(freq); // bin 0 is DC
Chris@40: 	paint.drawLine(cw + 7, h - vy, w - pkw - 1, h - vy);
Chris@0: 
Chris@0: 	if (h - vy - textHeight >= -2) {
Chris@1025: 	    int tx = w - 3 - paint.fontMetrics().width(text) - max(tickw, pkw);
Chris@0: 	    paint.drawText(tx, h - vy + toff, text);
Chris@0: 	}
Chris@0: 
Chris@0: 	py = vy;
Chris@0:     }
Chris@40: 
Chris@1103:     if (m_binScale == BinScale::Log) {
Chris@40: 
Chris@277:         // piano keyboard
Chris@277: 
Chris@690:         PianoScale().paintPianoVertical
Chris@690:             (v, paint, QRect(w - pkw - 1, 0, pkw, h),
Chris@690:              getEffectiveMinFrequency(), getEffectiveMaxFrequency());
Chris@40:     }
Chris@608: 
Chris@608:     m_haveDetailedScale = detailed;
Chris@0: }
Chris@0: 
Chris@187: class SpectrogramRangeMapper : public RangeMapper
Chris@187: {
Chris@187: public:
Chris@901:     SpectrogramRangeMapper(sv_samplerate_t sr, int /* fftsize */) :
Chris@901:         m_dist(sr / 2),
Chris@901:         m_s2(sqrt(sqrt(2))) { }
Chris@187:     ~SpectrogramRangeMapper() { }
Chris@187:     
Chris@901:     virtual int getPositionForValue(double value) const {
Chris@901: 
Chris@901:         double dist = m_dist;
Chris@187:     
Chris@187:         int n = 0;
Chris@187: 
Chris@901:         while (dist > (value + 0.00001) && dist > 0.1) {
Chris@187:             dist /= m_s2;
Chris@187:             ++n;
Chris@187:         }
Chris@187: 
Chris@187:         return n;
Chris@187:     }
Chris@724:     
Chris@901:     virtual int getPositionForValueUnclamped(double value) const {
Chris@724:         // We don't really support this
Chris@724:         return getPositionForValue(value);
Chris@724:     }
Chris@187: 
Chris@901:     virtual double getValueForPosition(int position) const {
Chris@187: 
Chris@187:         // Vertical zoom step 0 shows the entire range from DC ->
Chris@187:         // Nyquist frequency.  Step 1 shows 2^(1/4) of the range of
Chris@187:         // step 0, and so on until the visible range is smaller than
Chris@187:         // the frequency step between bins at the current fft size.
Chris@187: 
Chris@901:         double dist = m_dist;
Chris@187:     
Chris@187:         int n = 0;
Chris@187:         while (n < position) {
Chris@187:             dist /= m_s2;
Chris@187:             ++n;
Chris@187:         }
Chris@187: 
Chris@187:         return dist;
Chris@187:     }
Chris@187:     
Chris@901:     virtual double getValueForPositionUnclamped(int position) const {
Chris@724:         // We don't really support this
Chris@724:         return getValueForPosition(position);
Chris@724:     }
Chris@724: 
Chris@187:     virtual QString getUnit() const { return "Hz"; }
Chris@187: 
Chris@187: protected:
Chris@901:     double m_dist;
Chris@901:     double m_s2;
Chris@187: };
Chris@187: 
Chris@133: int
Chris@133: SpectrogramLayer::getVerticalZoomSteps(int &defaultStep) const
Chris@133: {
Chris@135:     if (!m_model) return 0;
Chris@187: 
Chris@907:     sv_samplerate_t sr = m_model->getSampleRate();
Chris@187: 
Chris@1087:     SpectrogramRangeMapper mapper(sr, getFFTSize());
Chris@1087: 
Chris@1087: //    int maxStep = mapper.getPositionForValue((double(sr) / getFFTSize()) + 0.001);
Chris@187:     int maxStep = mapper.getPositionForValue(0);
Chris@905:     int minStep = mapper.getPositionForValue(double(sr) / 2);
Chris@250: 
Chris@805:     int initialMax = m_initialMaxFrequency;
Chris@907:     if (initialMax == 0) initialMax = int(sr / 2);
Chris@250: 
Chris@250:     defaultStep = mapper.getPositionForValue(initialMax) - minStep;
Chris@250: 
Chris@587: //    SVDEBUG << "SpectrogramLayer::getVerticalZoomSteps: " << maxStep - minStep << " (" << maxStep <<"-" << minStep << "), default is " << defaultStep << " (from initial max freq " << initialMax << ")" << endl;
Chris@187: 
Chris@187:     return maxStep - minStep;
Chris@133: }
Chris@133: 
Chris@133: int
Chris@133: SpectrogramLayer::getCurrentVerticalZoomStep() const
Chris@133: {
Chris@133:     if (!m_model) return 0;
Chris@133: 
Chris@905:     double dmin, dmax;
Chris@133:     getDisplayExtents(dmin, dmax);
Chris@133:     
Chris@1087:     SpectrogramRangeMapper mapper(m_model->getSampleRate(), getFFTSize());
Chris@187:     int n = mapper.getPositionForValue(dmax - dmin);
Chris@587: //    SVDEBUG << "SpectrogramLayer::getCurrentVerticalZoomStep: " << n << endl;
Chris@133:     return n;
Chris@133: }
Chris@133: 
Chris@133: void
Chris@133: SpectrogramLayer::setVerticalZoomStep(int step)
Chris@133: {
Chris@187:     if (!m_model) return;
Chris@187: 
Chris@905:     double dmin = m_minFrequency, dmax = m_maxFrequency;
Chris@253: //    getDisplayExtents(dmin, dmax);
Chris@253: 
Chris@682: //    cerr << "current range " << dmin << " -> " << dmax << ", range " << dmax-dmin << ", mid " << (dmax + dmin)/2 << endl;
Chris@133:     
Chris@907:     sv_samplerate_t sr = m_model->getSampleRate();
Chris@1087:     SpectrogramRangeMapper mapper(sr, getFFTSize());
Chris@905:     double newdist = mapper.getValueForPosition(step);
Chris@905: 
Chris@905:     double newmin, newmax;
Chris@253: 
Chris@1103:     if (m_binScale == BinScale::Log) {
Chris@253: 
Chris@253:         // need to pick newmin and newmax such that
Chris@253:         //
Chris@253:         // (log(newmin) + log(newmax)) / 2 == logmid
Chris@253:         // and
Chris@253:         // newmax - newmin = newdist
Chris@253:         //
Chris@253:         // so log(newmax - newdist) + log(newmax) == 2logmid
Chris@253:         // log(newmax(newmax - newdist)) == 2logmid
Chris@253:         // newmax.newmax - newmax.newdist == exp(2logmid)
Chris@253:         // newmax^2 + (-newdist)newmax + -exp(2logmid) == 0
Chris@253:         // quadratic with a = 1, b = -newdist, c = -exp(2logmid), all known
Chris@253:         // 
Chris@253:         // positive root
Chris@253:         // newmax = (newdist + sqrt(newdist^2 + 4exp(2logmid))) / 2
Chris@253:         //
Chris@253:         // but logmid = (log(dmin) + log(dmax)) / 2
Chris@253:         // so exp(2logmid) = exp(log(dmin) + log(dmax))
Chris@253:         // = exp(log(dmin.dmax))
Chris@253:         // = dmin.dmax
Chris@253:         // so newmax = (newdist + sqrtf(newdist^2 + 4dmin.dmax)) / 2
Chris@253: 
Chris@907:         newmax = (newdist + sqrt(newdist*newdist + 4*dmin*dmax)) / 2;
Chris@253:         newmin = newmax - newdist;
Chris@253: 
Chris@682: //        cerr << "newmin = " << newmin << ", newmax = " << newmax << endl;
Chris@253: 
Chris@253:     } else {
Chris@905:         double dmid = (dmax + dmin) / 2;
Chris@253:         newmin = dmid - newdist / 2;
Chris@253:         newmax = dmid + newdist / 2;
Chris@253:     }
Chris@187: 
Chris@905:     double mmin, mmax;
Chris@187:     mmin = 0;
Chris@905:     mmax = double(sr) / 2;
Chris@133:     
Chris@187:     if (newmin < mmin) {
Chris@187:         newmax += (mmin - newmin);
Chris@187:         newmin = mmin;
Chris@187:     }
Chris@187:     if (newmax > mmax) {
Chris@187:         newmax = mmax;
Chris@187:     }
Chris@133:     
Chris@587: //    SVDEBUG << "SpectrogramLayer::setVerticalZoomStep: " << step << ": " << newmin << " -> " << newmax << " (range " << newdist << ")" << endl;
Chris@253: 
Chris@907:     setMinFrequency(int(lrint(newmin)));
Chris@907:     setMaxFrequency(int(lrint(newmax)));
Chris@187: }
Chris@187: 
Chris@187: RangeMapper *
Chris@187: SpectrogramLayer::getNewVerticalZoomRangeMapper() const
Chris@187: {
Chris@187:     if (!m_model) return 0;
Chris@1087:     return new SpectrogramRangeMapper(m_model->getSampleRate(), getFFTSize());
Chris@133: }
Chris@133: 
Chris@273: void
Chris@918: SpectrogramLayer::updateMeasureRectYCoords(LayerGeometryProvider *v, const MeasureRect &r) const
Chris@273: {
Chris@273:     int y0 = 0;
Chris@907:     if (r.startY > 0.0) y0 = int(getYForFrequency(v, r.startY));
Chris@273:     
Chris@273:     int y1 = y0;
Chris@907:     if (r.endY > 0.0) y1 = int(getYForFrequency(v, r.endY));
Chris@273: 
Chris@587: //    SVDEBUG << "SpectrogramLayer::updateMeasureRectYCoords: start " << r.startY << " -> " << y0 << ", end " << r.endY << " -> " << y1 << endl;
Chris@273: 
Chris@273:     r.pixrect = QRect(r.pixrect.x(), y0, r.pixrect.width(), y1 - y0);
Chris@273: }
Chris@273: 
Chris@273: void
Chris@918: SpectrogramLayer::setMeasureRectYCoord(LayerGeometryProvider *v, MeasureRect &r, bool start, int y) const
Chris@273: {
Chris@273:     if (start) {
Chris@273:         r.startY = getFrequencyForY(v, y);
Chris@273:         r.endY = r.startY;
Chris@273:     } else {
Chris@273:         r.endY = getFrequencyForY(v, y);
Chris@273:     }
Chris@587: //    SVDEBUG << "SpectrogramLayer::setMeasureRectYCoord: start " << r.startY << " <- " << y << ", end " << r.endY << " <- " << y << endl;
Chris@273: 
Chris@273: }
Chris@273: 
Chris@316: void
Chris@316: SpectrogramLayer::toXml(QTextStream &stream,
Chris@316:                         QString indent, QString extraAttributes) const
Chris@6: {
Chris@6:     QString s;
Chris@6:     
Chris@6:     s += QString("channel=\"%1\" "
Chris@6: 		 "windowSize=\"%2\" "
Chris@153: 		 "windowHopLevel=\"%3\" "
Chris@153: 		 "gain=\"%4\" "
Chris@153: 		 "threshold=\"%5\" ")
Chris@6: 	.arg(m_channel)
Chris@6: 	.arg(m_windowSize)
Chris@97: 	.arg(m_windowHopLevel)
Chris@37: 	.arg(m_gain)
Chris@37: 	.arg(m_threshold);
Chris@37: 
Chris@37:     s += QString("minFrequency=\"%1\" "
Chris@37: 		 "maxFrequency=\"%2\" "
Chris@37: 		 "colourScale=\"%3\" "
Chris@37: 		 "colourScheme=\"%4\" "
Chris@37: 		 "colourRotation=\"%5\" "
Chris@37: 		 "frequencyScale=\"%6\" "
Chris@761: 		 "binDisplay=\"%7\" ")
Chris@37: 	.arg(m_minFrequency)
Chris@6: 	.arg(m_maxFrequency)
Chris@1104: 	.arg(convertFromColourScale(m_colourScale))
Chris@197: 	.arg(m_colourMap)
Chris@37: 	.arg(m_colourRotation)
Chris@1103: 	.arg(int(m_binScale))
Chris@1103: 	.arg(int(m_binDisplay));
Chris@761: 
Chris@1009:     // New-style normalization attributes, allowing for more types of
Chris@1009:     // normalization in future: write out the column normalization
Chris@1009:     // type separately, and then whether we are normalizing visible
Chris@1009:     // area as well afterwards
Chris@1009:     
Chris@1009:     s += QString("columnNormalization=\"%1\" ")
Chris@1104:         .arg(m_normalization == ColumnNormalization::Max1 ? "peak" :
Chris@1104:              m_normalization == ColumnNormalization::Hybrid ? "hybrid" : "none");
Chris@1009: 
Chris@1009:     // Old-style normalization attribute. We *don't* write out
Chris@1009:     // normalizeHybrid here because the only release that would accept
Chris@1009:     // it (Tony v1.0) has a totally different scale factor for
Chris@1009:     // it. We'll just have to accept that session files from Tony
Chris@1009:     // v2.0+ will look odd in Tony v1.0
Chris@1009:     
Chris@1009:     s += QString("normalizeColumns=\"%1\" ")
Chris@1104: 	.arg(m_normalization == ColumnNormalization::Max1 ? "true" : "false");
Chris@1009: 
Chris@1009:     // And this applies to both old- and new-style attributes
Chris@1009:     
Chris@1009:     s += QString("normalizeVisibleArea=\"%1\" ")
Chris@1104:         .arg(m_normalizeVisibleArea ? "true" : "false");
Chris@1009:     
Chris@316:     Layer::toXml(stream, indent, extraAttributes + " " + s);
Chris@6: }
Chris@6: 
Chris@11: void
Chris@11: SpectrogramLayer::setProperties(const QXmlAttributes &attributes)
Chris@11: {
Chris@11:     bool ok = false;
Chris@11: 
Chris@11:     int channel = attributes.value("channel").toInt(&ok);
Chris@11:     if (ok) setChannel(channel);
Chris@11: 
Chris@805:     int windowSize = attributes.value("windowSize").toUInt(&ok);
Chris@11:     if (ok) setWindowSize(windowSize);
Chris@11: 
Chris@805:     int windowHopLevel = attributes.value("windowHopLevel").toUInt(&ok);
Chris@97:     if (ok) setWindowHopLevel(windowHopLevel);
Chris@97:     else {
Chris@805:         int windowOverlap = attributes.value("windowOverlap").toUInt(&ok);
Chris@97:         // a percentage value
Chris@97:         if (ok) {
Chris@97:             if (windowOverlap == 0) setWindowHopLevel(0);
Chris@97:             else if (windowOverlap == 25) setWindowHopLevel(1);
Chris@97:             else if (windowOverlap == 50) setWindowHopLevel(2);
Chris@97:             else if (windowOverlap == 75) setWindowHopLevel(3);
Chris@97:             else if (windowOverlap == 90) setWindowHopLevel(4);
Chris@97:         }
Chris@97:     }
Chris@11: 
Chris@11:     float gain = attributes.value("gain").toFloat(&ok);
Chris@11:     if (ok) setGain(gain);
Chris@11: 
Chris@37:     float threshold = attributes.value("threshold").toFloat(&ok);
Chris@37:     if (ok) setThreshold(threshold);
Chris@37: 
Chris@805:     int minFrequency = attributes.value("minFrequency").toUInt(&ok);
Chris@187:     if (ok) {
Chris@587:         SVDEBUG << "SpectrogramLayer::setProperties: setting min freq to " << minFrequency << endl;
Chris@187:         setMinFrequency(minFrequency);
Chris@187:     }
Chris@37: 
Chris@805:     int maxFrequency = attributes.value("maxFrequency").toUInt(&ok);
Chris@187:     if (ok) {
Chris@587:         SVDEBUG << "SpectrogramLayer::setProperties: setting max freq to " << maxFrequency << endl;
Chris@187:         setMaxFrequency(maxFrequency);
Chris@187:     }
Chris@11: 
Chris@1105:     ColourScaleType colourScale = convertToColourScale
Chris@1092:         (attributes.value("colourScale").toInt(&ok));
Chris@11:     if (ok) setColourScale(colourScale);
Chris@11: 
Chris@197:     int colourMap = attributes.value("colourScheme").toInt(&ok);
Chris@197:     if (ok) setColourMap(colourMap);
Chris@11: 
Chris@37:     int colourRotation = attributes.value("colourRotation").toInt(&ok);
Chris@37:     if (ok) setColourRotation(colourRotation);
Chris@37: 
Chris@1103:     BinScale binScale = (BinScale)
Chris@11: 	attributes.value("frequencyScale").toInt(&ok);
Chris@1093:     if (ok) setBinScale(binScale);
Chris@1093: 
Chris@1103:     BinDisplay binDisplay = (BinDisplay)
Chris@37: 	attributes.value("binDisplay").toInt(&ok);
Chris@37:     if (ok) setBinDisplay(binDisplay);
Chris@36: 
Chris@1009:     bool haveNewStyleNormalization = false;
Chris@1009:     
Chris@1009:     QString columnNormalization = attributes.value("columnNormalization");
Chris@1009: 
Chris@1009:     if (columnNormalization != "") {
Chris@1009: 
Chris@1009:         haveNewStyleNormalization = true;
Chris@1009: 
Chris@1009:         if (columnNormalization == "peak") {
Chris@1104:             setNormalization(ColumnNormalization::Max1);
Chris@1009:         } else if (columnNormalization == "hybrid") {
Chris@1104:             setNormalization(ColumnNormalization::Hybrid);
Chris@1009:         } else if (columnNormalization == "none") {
Chris@1104:             setNormalization(ColumnNormalization::None);
Chris@1009:         } else {
Chris@1009:             cerr << "NOTE: Unknown or unsupported columnNormalization attribute \""
Chris@1009:                  << columnNormalization << "\"" << endl;
Chris@1009:         }
Chris@1009:     }
Chris@1009: 
Chris@1009:     if (!haveNewStyleNormalization) {
Chris@1009: 
Chris@1009:         bool normalizeColumns =
Chris@1009:             (attributes.value("normalizeColumns").trimmed() == "true");
Chris@1009:         if (normalizeColumns) {
Chris@1104:             setNormalization(ColumnNormalization::Max1);
Chris@1009:         }
Chris@1009: 
Chris@1009:         bool normalizeHybrid =
Chris@1009:             (attributes.value("normalizeHybrid").trimmed() == "true");
Chris@1009:         if (normalizeHybrid) {
Chris@1104:             setNormalization(ColumnNormalization::Hybrid);
Chris@1009:         }
Chris@862:     }
Chris@153: 
Chris@153:     bool normalizeVisibleArea =
Chris@1099:         (attributes.value("normalizeVisibleArea").trimmed() == "true");
Chris@1104:     setNormalizeVisibleArea(normalizeVisibleArea);
Chris@1104: 
Chris@1104:     if (!haveNewStyleNormalization && m_normalization == ColumnNormalization::Hybrid) {
Chris@1009:         // Tony v1.0 is (and hopefully will remain!) the only released
Chris@1009:         // SV-a-like to use old-style attributes when saving sessions
Chris@1009:         // that ask for hybrid normalization. It saves them with the
Chris@1009:         // wrong gain factor, so hack in a fix for that here -- this
Chris@1009:         // gives us backward but not forward compatibility.
Chris@1087:         setGain(m_gain / float(getFFTSize() / 2));
Chris@862:     }
Chris@11: }
Chris@11: