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) :
23 m_envelopeSize(producer->get_length()),
27 m_envelopeStdDevCalculated(false),
28 m_envelopeIsNormalized(false)
30 // make a copy of the producer to avoid audio playback issues
31 m_producer = new Mlt::Producer(*(producer->profile()), producer->get("resource"));
32 m_info = new AudioInfo(m_producer);
34 Q_ASSERT(m_offset >= 0);
36 Q_ASSERT(m_length+m_offset <= m_envelopeSize);
37 m_envelopeSize = m_length;
41 AudioEnvelope::~AudioEnvelope()
43 if (m_envelope != NULL) {
52 const int64_t *AudioEnvelope::envelope()
54 if (m_envelope == NULL) {
59 int AudioEnvelope::envelopeSize() const
61 return m_envelopeSize;
63 int64_t AudioEnvelope::maxValue() const
71 void AudioEnvelope::loadEnvelope()
73 Q_ASSERT(m_envelope == NULL);
75 std::cout << "Loading envelope ..." << std::endl;
77 int samplingRate = m_info->info(0)->samplingRate();
78 mlt_audio_format format_s16 = mlt_audio_s16;
85 m_envelope = new int64_t[m_envelopeSize];
92 m_producer->seek(m_offset);
93 m_producer->set_speed(1.0); // This is necessary, otherwise we don't get any new frames in the 2nd run.
94 for (int i = 0; i < m_envelopeSize; i++) {
96 frame = m_producer->get_frame(i);
97 position = mlt_frame_get_position(frame->get_frame());
98 samples = mlt_sample_calculator(m_producer->get_fps(), samplingRate, position);
100 int16_t *data = static_cast<int16_t*>(frame->get_audio(format_s16, samplingRate, channels, samples));
103 for (int k = 0; k < samples; k++) {
104 sum += fabs(data[k]);
108 m_envelopeMean += sum;
109 if (sum > m_envelopeMax) {
113 // std::cout << position << "|" << m_producer->get_playtime()
114 // << "-" << m_producer->get_in() << "+" << m_producer->get_out() << " ";
119 if (m_length > 0 && count > m_length) {
123 m_envelopeMean /= m_envelopeSize;
124 std::cout << "Calculating the envelope (" << m_envelopeSize << " frames) took "
125 << t.elapsed() << " ms." << std::endl;
128 int64_t AudioEnvelope::loadStdDev()
130 if (m_envelopeStdDevCalculated) {
131 std::cout << "Standard deviation already calculated, not re-calculating." << std::endl;
134 if (m_envelope == NULL) {
138 m_envelopeStdDev = 0;
139 for (int i = 0; i < m_envelopeSize; i++) {
140 m_envelopeStdDev += sqrt((m_envelope[i]-m_envelopeMean)*(m_envelope[i]-m_envelopeMean)/m_envelopeSize);
142 m_envelopeStdDevCalculated = true;
145 return m_envelopeStdDev;
148 void AudioEnvelope::normalizeEnvelope(bool clampTo0)
150 if (m_envelope == NULL) {
154 if (!m_envelopeIsNormalized) {
158 for (int i = 0; i < m_envelopeSize; i++) {
160 m_envelope[i] -= m_envelopeMean;
163 if (m_envelope[i] < 0) { m_envelope[i] = 0; }
166 if (m_envelope[i] > m_envelopeMax) {
167 m_envelopeMax = m_envelope[i];
170 newMean += m_envelope[i];
172 m_envelopeMean = newMean / m_envelopeSize;
174 m_envelopeIsNormalized = true;
179 QImage AudioEnvelope::drawEnvelope()
181 if (m_envelope == NULL) {
185 QImage img(m_envelopeSize, 400, QImage::Format_ARGB32);
186 img.fill(qRgb(255,255,255));
188 for (int x = 0; x < img.width(); x++) {
189 fy = m_envelope[x]/double(m_envelopeMax) * img.height();
190 for (int y = img.height()-1; y > img.height()-1-fy; y--) {
191 img.setPixel(x,y, qRgb(50, 50, 50));
197 void AudioEnvelope::dumpInfo() const
199 if (m_envelope == NULL) {
200 std::cout << "Envelope not generated, no information available." << std::endl;
202 std::cout << "Envelope info" << std::endl
203 << "* size = " << m_envelopeSize << std::endl
204 << "* max = " << m_envelopeMax << std::endl
205 << "* ยต = " << m_envelopeMean << std::endl
207 if (m_envelopeStdDevCalculated) {
208 std::cout << "* s = " << m_envelopeStdDev << std::endl;