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;
62 void AudioEnvelope::loadEnvelope()
64 Q_ASSERT(m_envelope == NULL);
66 std::cout << "Loading envelope ..." << std::endl;
68 int samplingRate = m_info->info(0)->samplingRate();
69 mlt_audio_format format_s16 = mlt_audio_s16;
76 m_envelope = new int64_t[m_envelopeSize];
83 m_producer->seek(m_offset);
84 m_producer->set_speed(1.0); // This is necessary, otherwise we don't get any new frames in the 2nd run.
85 for (int i = 0; i < m_envelopeSize; i++) {
87 frame = m_producer->get_frame(i);
88 position = mlt_frame_get_position(frame->get_frame());
89 samples = mlt_sample_calculator(m_producer->get_fps(), samplingRate, position);
91 int16_t *data = static_cast<int16_t*>(frame->get_audio(format_s16, samplingRate, channels, samples));
94 for (int k = 0; k < samples; k++) {
99 m_envelopeMean += sum;
100 if (sum > m_envelopeMax) {
104 std::cout << position << "|" << m_producer->get_playtime()
105 << "-" << m_producer->get_in() << "+" << m_producer->get_out() << " ";
110 if (m_length > 0 && count > m_length) {
114 m_envelopeMean /= m_envelopeSize;
115 std::cout << "Calculating the envelope (" << m_envelopeSize << " frames) took "
116 << t.elapsed() << " ms." << std::endl;
119 int64_t AudioEnvelope::loadStdDev()
121 if (m_envelopeStdDevCalculated) {
122 std::cout << "Standard deviation already calculated, not re-calculating." << std::endl;
125 if (m_envelope == NULL) {
129 m_envelopeStdDev = 0;
130 for (int i = 0; i < m_envelopeSize; i++) {
131 m_envelopeStdDev += sqrt((m_envelope[i]-m_envelopeMean)*(m_envelope[i]-m_envelopeMean)/m_envelopeSize);
133 m_envelopeStdDevCalculated = true;
136 return m_envelopeStdDev;
139 void AudioEnvelope::normalizeEnvelope(bool clampTo0)
141 if (m_envelope == NULL) {
145 if (!m_envelopeIsNormalized) {
149 for (int i = 0; i < m_envelopeSize; i++) {
151 m_envelope[i] -= m_envelopeMean;
154 if (m_envelope[i] < 0) { m_envelope[i] = 0; }
157 if (m_envelope[i] > m_envelopeMax) {
158 m_envelopeMax = m_envelope[i];
161 newMean += m_envelope[i];
163 m_envelopeMean = newMean / m_envelopeSize;
165 m_envelopeIsNormalized = true;
170 QImage AudioEnvelope::drawEnvelope()
172 if (m_envelope == NULL) {
176 QImage img(m_envelopeSize, 400, QImage::Format_ARGB32);
177 img.fill(qRgb(255,255,255));
179 for (int x = 0; x < img.width(); x++) {
180 fy = m_envelope[x]/double(m_envelopeMax) * img.height();
181 for (int y = img.height()-1; y > img.height()-1-fy; y--) {
182 img.setPixel(x,y, qRgb(50, 50, 50));
188 void AudioEnvelope::dumpInfo() const
190 if (m_envelope == NULL) {
191 std::cout << "Envelope not generated, no information available." << std::endl;
193 std::cout << "Envelope info" << std::endl
194 << "* size = " << m_envelopeSize << std::endl
195 << "* max = " << m_envelopeMax << std::endl
196 << "* ยต = " << m_envelopeMean << std::endl
198 if (m_envelopeStdDevCalculated) {
199 std::cout << "* s = " << m_envelopeStdDev << std::endl;