1 /***************************************************************************
2 * Copyright (C) 2012 by Simon Andreas Eugster (simon.eu@gmail.com) *
3 * This file is part of kdenlive. See www.kdenlive.org. *
5 * This program is free software; you can redistribute it and/or modify *
6 * it under the terms of the GNU General Public License as published by *
7 * the Free Software Foundation; either version 2 of the License, or *
8 * (at your option) any later version. *
9 ***************************************************************************/
11 #include "audioEnvelope.h"
13 #include "audioStreamInfo.h"
19 AudioEnvelope::AudioEnvelope(Mlt::Producer *producer, int offset, int length) :
24 m_envelopeSize(producer->get_length()),
25 m_envelopeStdDevCalculated(false),
26 m_envelopeIsNormalized(false)
28 m_info = new AudioInfo(m_producer);
30 Q_ASSERT(m_offset >= 0);
32 Q_ASSERT(m_length+m_offset <= m_envelopeSize);
33 m_envelopeSize = m_length;
37 AudioEnvelope::~AudioEnvelope()
39 if (m_envelope != NULL) {
47 const int64_t *AudioEnvelope::envelope()
49 if (m_envelope == NULL) {
54 int AudioEnvelope::envelopeSize() const
56 return m_envelopeSize;
58 int64_t AudioEnvelope::maxValue() const
66 void AudioEnvelope::loadEnvelope()
68 Q_ASSERT(m_envelope == NULL);
70 std::cout << "Loading envelope ..." << std::endl;
72 int samplingRate = m_info->info(0)->samplingRate();
73 mlt_audio_format format_s16 = mlt_audio_s16;
80 m_envelope = new int64_t[m_envelopeSize];
87 m_producer->seek(m_offset);
88 m_producer->set_speed(1.0); // This is necessary, otherwise we don't get any new frames in the 2nd run.
89 for (int i = 0; i < m_envelopeSize; i++) {
91 frame = m_producer->get_frame(i);
92 position = mlt_frame_get_position(frame->get_frame());
93 samples = mlt_sample_calculator(m_producer->get_fps(), samplingRate, position);
95 int16_t *data = static_cast<int16_t*>(frame->get_audio(format_s16, samplingRate, channels, samples));
98 for (int k = 0; k < samples; k++) {
103 m_envelopeMean += sum;
104 if (sum > m_envelopeMax) {
108 // std::cout << position << "|" << m_producer->get_playtime()
109 // << "-" << m_producer->get_in() << "+" << m_producer->get_out() << " ";
114 if (m_length > 0 && count > m_length) {
118 m_envelopeMean /= m_envelopeSize;
119 std::cout << "Calculating the envelope (" << m_envelopeSize << " frames) took "
120 << t.elapsed() << " ms." << std::endl;
123 int64_t AudioEnvelope::loadStdDev()
125 if (m_envelopeStdDevCalculated) {
126 std::cout << "Standard deviation already calculated, not re-calculating." << std::endl;
129 if (m_envelope == NULL) {
133 m_envelopeStdDev = 0;
134 for (int i = 0; i < m_envelopeSize; i++) {
135 m_envelopeStdDev += sqrt((m_envelope[i]-m_envelopeMean)*(m_envelope[i]-m_envelopeMean)/m_envelopeSize);
137 m_envelopeStdDevCalculated = true;
140 return m_envelopeStdDev;
143 void AudioEnvelope::normalizeEnvelope(bool clampTo0)
145 if (m_envelope == NULL) {
149 if (!m_envelopeIsNormalized) {
153 for (int i = 0; i < m_envelopeSize; i++) {
155 m_envelope[i] -= m_envelopeMean;
158 if (m_envelope[i] < 0) { m_envelope[i] = 0; }
161 if (m_envelope[i] > m_envelopeMax) {
162 m_envelopeMax = m_envelope[i];
165 newMean += m_envelope[i];
167 m_envelopeMean = newMean / m_envelopeSize;
169 m_envelopeIsNormalized = true;
174 QImage AudioEnvelope::drawEnvelope()
176 if (m_envelope == NULL) {
180 QImage img(m_envelopeSize, 400, QImage::Format_ARGB32);
181 img.fill(qRgb(255,255,255));
183 for (int x = 0; x < img.width(); x++) {
184 fy = m_envelope[x]/double(m_envelopeMax) * img.height();
185 for (int y = img.height()-1; y > img.height()-1-fy; y--) {
186 img.setPixel(x,y, qRgb(50, 50, 50));
192 void AudioEnvelope::dumpInfo() const
194 if (m_envelope == NULL) {
195 std::cout << "Envelope not generated, no information available." << std::endl;
197 std::cout << "Envelope info" << std::endl
198 << "* size = " << m_envelopeSize << std::endl
199 << "* max = " << m_envelopeMax << std::endl
200 << "* ยต = " << m_envelopeMean << std::endl
202 if (m_envelopeStdDevCalculated) {
203 std::cout << "* s = " << m_envelopeStdDev << std::endl;