annotate plugins/SegmenterPlugin.cpp @ 42:0f85778f1b53

* some more fixes, + add chroma option to similarity plugin
author Chris Cannam <c.cannam@qmul.ac.uk>
date Mon, 14 Jan 2008 18:14:55 +0000
parents 77e394a5f3c9
children 5d7ce1d87301
rev   line source
c@38 1 /* -*- c-basic-offset: 4 indent-tabs-mode: nil -*- vi:set ts=8 sts=4 sw=4: */
c@38 2
c@36 3 /*
c@38 4 * SegmenterPlugin.cpp
c@36 5 *
c@38 6 * Created by Mark Levy on 24/03/2006.
c@38 7 * Copyright 2006 Centre for Digital Music, Queen Mary, University of London.
c@38 8 * All rights reserved.
c@36 9 */
c@36 10
c@36 11 #include <iostream>
c@36 12 #include <sstream>
c@36 13
c@36 14 #include "SegmenterPlugin.h"
c@37 15 #include "dsp/segmentation/ClusterMeltSegmenter.h"
c@36 16
c@36 17 using std::string;
c@36 18 using std::vector;
c@36 19 using std::cerr;
c@36 20 using std::endl;
c@36 21 using std::ostringstream;
c@36 22
c@36 23 SegmenterPlugin::SegmenterPlugin(float inputSampleRate) :
c@38 24 Plugin(inputSampleRate),
c@38 25 segmenter(0),
c@38 26 nSegmentTypes(10),
c@38 27 featureType(feature_types(1))
c@36 28 {
c@36 29
c@36 30 }
c@36 31
c@36 32 SegmenterPlugin::~SegmenterPlugin()
c@36 33 {
c@38 34 delete segmenter;
c@36 35 }
c@36 36
c@36 37 string
c@36 38 SegmenterPlugin::getMaker() const
c@36 39 {
c@36 40 return "Mark Levy, Queen Mary, University of London";
c@36 41 }
c@36 42
c@36 43 int
c@36 44 SegmenterPlugin::getPluginVersion() const
c@36 45 {
c@36 46 return 2;
c@36 47 }
c@36 48
c@36 49 string
c@36 50 SegmenterPlugin::getCopyright() const
c@36 51 {
c@37 52 return "Copyright (c) 2006-2008 - All Rights Reserved";
c@36 53 }
c@36 54
c@36 55 bool
c@36 56 SegmenterPlugin::initialise(size_t channels, size_t stepSize, size_t blockSize)
c@36 57 {
c@36 58 if (channels < getMinChannelCount() ||
c@36 59 channels > getMaxChannelCount()) return false;
c@36 60
c@36 61 if (!segmenter) makeSegmenter();
c@36 62
c@39 63 if (stepSize != hopsize) {
c@39 64 std::cerr << "SegmenterPlugin::initialise: supplied step size "
c@39 65 << stepSize << " differs from required step size " << hopsize
c@39 66 << std::endl;
c@39 67 return false;
c@39 68 }
c@39 69
c@39 70 if (blockSize != windowsize) {
c@39 71 std::cerr << "SegmenterPlugin::initialise: supplied block size "
c@39 72 << blockSize << " differs from required block size " << windowsize
c@39 73 << std::endl;
c@39 74 return false;
c@39 75 }
c@36 76
c@36 77 return true;
c@36 78 }
c@36 79
c@36 80 void
c@36 81 SegmenterPlugin::reset()
c@36 82 {
c@38 83 //!!!
c@36 84 }
c@36 85
c@36 86 size_t
c@36 87 SegmenterPlugin::getPreferredStepSize() const
c@36 88 {
c@38 89 if (!segmenter) makeSegmenter();
c@38 90 return hopsize;
c@36 91 }
c@36 92
c@36 93 size_t
c@36 94 SegmenterPlugin::getPreferredBlockSize() const
c@36 95 {
c@38 96 if (!segmenter) makeSegmenter();
c@38 97 return windowsize;
c@36 98 }
c@36 99
c@36 100 SegmenterPlugin::ParameterList SegmenterPlugin::getParameterDescriptors() const
c@36 101 {
c@36 102 ParameterList list;
c@36 103
c@36 104 ParameterDescriptor desc;
c@36 105 desc.identifier = "nSegmentTypes";
c@36 106 desc.name = "Number of segment-types";
c@36 107 desc.description = "Maximum number of different kinds of segment to find";
c@36 108 desc.unit = "";
c@36 109 desc.minValue = 2;
c@36 110 desc.maxValue = 12;
c@36 111 desc.defaultValue = 10;
c@36 112 desc.isQuantized = true;
c@36 113 desc.quantizeStep = 1;
c@36 114 list.push_back(desc);
c@36 115
c@36 116 ParameterDescriptor desc2;
c@36 117 desc2.identifier = "featureType";
c@36 118 desc2.name = "Feature Type";
c@39 119 desc2.description = "Try Chromatic for acoustic or pre-1980 recordings, otherwise use Hybrid";
c@36 120 desc2.unit = "";
c@36 121 desc2.minValue = 1;
c@39 122 desc2.maxValue = 3;
c@36 123 desc2.defaultValue = 1;
c@36 124 desc2.isQuantized = true;
c@36 125 desc2.quantizeStep = 1;
c@39 126 desc2.valueNames.push_back("Hybrid (Constant-Q)");
c@39 127 desc2.valueNames.push_back("Chromatic (Chroma)");
c@39 128 desc2.valueNames.push_back("Timbral (MFCC)");
c@36 129 list.push_back(desc2);
c@36 130
c@36 131 return list;
c@36 132 }
c@36 133
c@36 134 float
c@36 135 SegmenterPlugin::getParameter(std::string param) const
c@36 136 {
c@36 137 if (param == "nSegmentTypes") {
c@36 138 return nSegmentTypes;
c@36 139 }
c@36 140
c@38 141 if (param == "featureType") {
c@38 142 return featureType;
c@38 143 }
c@36 144
c@38 145 std::cerr << "WARNING: SegmenterPlugin::getParameter: unknown parameter \""
c@38 146 << param << "\"" << std::endl;
c@36 147 return 0.0;
c@36 148 }
c@36 149
c@36 150 void
c@36 151 SegmenterPlugin::setParameter(std::string param, float value)
c@36 152 {
c@36 153 if (param == "nSegmentTypes") {
c@38 154
c@42 155 nSegmentTypes = int(value + 0.0001);
c@38 156
c@38 157 } else {
c@38 158
c@38 159 if (param == "featureType") {
c@38 160 if (featureType != feature_types(value)) // feature type changed, create a new segmenter
c@38 161 {
c@38 162 featureType = feature_types(value);
c@38 163 makeSegmenter();
c@38 164 }
c@38 165 }
c@38 166 else
c@38 167 {
c@38 168 std::cerr << "WARNING: SegmenterPlugin::setParameter: unknown parameter \""
c@38 169 << param << "\"" << std::endl;
c@38 170 }
c@38 171 }
c@36 172 }
c@36 173
c@36 174 void
c@36 175 SegmenterPlugin::makeSegmenter() const
c@36 176 {
c@38 177 ClusterMeltSegmenterParams params = ClusterMeltSegmenterParams();
c@38 178 params.featureType = (feature_types) featureType;
c@38 179
c@38 180 if (params.featureType == FEATURE_TYPE_CONSTQ)
c@38 181 {
c@38 182 params.ncomponents = 20;
c@38 183 params.neighbourhoodLimit = 30;
c@38 184 }
c@38 185 if (params.featureType == FEATURE_TYPE_CHROMA)
c@38 186 {
c@38 187 params.hopSize = 0.1;
c@38 188 params.windowSize = 0.372;
c@38 189 params.nbins = 12;
c@38 190 params.histogramLength = 20;
c@38 191 params.neighbourhoodLimit = 40;
c@38 192 }
c@39 193 if (params.featureType == FEATURE_TYPE_MFCC)
c@39 194 {
c@39 195 params.ncomponents = 20;
c@39 196 params.neighbourhoodLimit = 30;
c@39 197 }
c@38 198 delete segmenter;
c@38 199
c@38 200 segmenter = new ClusterMeltSegmenter(params);
c@38 201 segmenter->initialise(m_inputSampleRate);
c@38 202 hopsize = segmenter->getHopsize();
c@38 203 windowsize = segmenter->getWindowsize();
c@38 204
c@38 205 std::cerr << "segmenter window size: " << segmenter->getWindowsize()
c@38 206 << std::endl;
c@36 207 }
c@36 208
c@36 209 SegmenterPlugin::OutputList
c@36 210 SegmenterPlugin::getOutputDescriptors() const
c@36 211 {
c@36 212 OutputList list;
c@36 213
c@38 214 OutputDescriptor segmentation;
c@38 215 segmentation.identifier = "segmentation";
c@36 216 segmentation.name = "Segmentation";
c@36 217 segmentation.description = "Segmentation";
c@36 218 segmentation.unit = "segment-type";
c@36 219 segmentation.hasFixedBinCount = true;
c@36 220 segmentation.binCount = 1;
c@40 221 segmentation.hasKnownExtents = true;
c@38 222 segmentation.minValue = 1;
c@38 223 segmentation.maxValue = nSegmentTypes;
c@38 224 segmentation.isQuantized = true;
c@38 225 segmentation.quantizeStep = 1;
c@36 226 segmentation.sampleType = OutputDescriptor::VariableSampleRate;
c@36 227 segmentation.sampleRate = m_inputSampleRate / getPreferredStepSize();
c@36 228
c@36 229 list.push_back(segmentation);
c@36 230
c@38 231 return list;
c@36 232 }
c@36 233
c@36 234 SegmenterPlugin::FeatureSet
c@36 235 SegmenterPlugin::process(const float *const *inputBuffers, Vamp::RealTime /* timestamp */)
c@36 236 {
c@36 237 // convert float* to double*
c@36 238 double *tempBuffer = new double[windowsize];
c@36 239 for (size_t i = 0; i < windowsize; ++i) {
c@38 240 tempBuffer[i] = inputBuffers[0][i];
c@36 241 }
c@38 242
c@38 243 segmenter->extractFeatures(tempBuffer, segmenter->getWindowsize());
c@38 244
c@38 245 delete [] tempBuffer;
c@36 246
c@38 247 return FeatureSet();
c@36 248 }
c@36 249
c@36 250 SegmenterPlugin::FeatureSet
c@36 251 SegmenterPlugin::getRemainingFeatures()
c@36 252 {
c@36 253 segmenter->segment(nSegmentTypes);
c@38 254 Segmentation segm = segmenter->getSegmentation();
c@36 255
c@38 256 FeatureSet returnFeatures;
c@36 257
c@36 258 for (int i = 0; i < segm.segments.size(); ++i) {
c@36 259
c@38 260 Segment s = segm.segments[i];
c@36 261
c@38 262 Feature feature;
c@38 263 feature.hasTimestamp = true;
c@38 264 feature.timestamp = Vamp::RealTime::frame2RealTime(s.start, static_cast<unsigned int>(m_inputSampleRate));
c@36 265
c@38 266 vector<float> floatval;
c@38 267 floatval.push_back(s.type);
c@38 268 feature.values = floatval;
c@36 269
c@38 270 ostringstream oss;
c@38 271 oss << s.type;
c@38 272 feature.label = oss.str();
c@36 273
c@38 274 returnFeatures[0].push_back(feature);
c@36 275 }
c@36 276
c@36 277 return returnFeatures;
c@36 278 }