annotate loudness.js @ 752:fab881bdeeac

Offline Context no longer promise type for Safari OSX support.
author Nicholas Jillings <nicholas.jillings@eecs.qmul.ac.uk>
date Mon, 21 Dec 2015 13:54:49 +0000
parents 07c996307cbd
children 56ee920d5d85
rev   line source
n@749 1 /**
n@749 2 * loundess.js
n@749 3 * Loudness module for the Web Audio Evaluation Toolbox
n@749 4 * Allows for automatic calculation of loudness of Web Audio API Buffer objects,
n@749 5 * return gain values to correct for a target loudness or match loudness between
n@749 6 * multiple objects
n@749 7 */
n@749 8
n@749 9 var interval_cal_loudness_event = null;
n@749 10
nicholas@752 11 if (typeof OfflineAudioContext == "undefined"){
nicholas@752 12 var OfflineAudioContext = webkitOfflineAudioContext;
nicholas@752 13 }
nicholas@752 14
n@749 15 function calculateLoudness(buffer, timescale, target, offlineContext)
n@749 16 {
n@749 17 // This function returns the EBU R 128 specification loudness model and sets the linear gain required to match -23 LUFS
n@749 18 // buffer -> Web Audio API Buffer object
n@749 19 // timescale -> M or Momentary (returns Array), S or Short (returns Array),
n@749 20 // I or Integrated (default, returns number)
n@749 21 // target -> default is -23 LUFS but can be any LUFS measurement.
n@749 22
n@749 23 if (buffer == undefined)
n@749 24 {
n@749 25 return 0;
n@749 26 }
n@749 27 if (timescale == undefined)
n@749 28 {
n@749 29 timescale = "I";
n@749 30 }
n@749 31 if (target == undefined)
n@749 32 {
n@749 33 target = -23;
n@749 34 }
n@749 35 if (offlineContext == undefined)
n@749 36 {
n@749 37 offlineContext = new OfflineAudioContext(buffer.numberOfChannels, buffer.length, buffer.sampleRate);
n@749 38 }
n@749 39 // Create the required filters
n@749 40 var KFilter = offlineContext.createBiquadFilter();
n@749 41 KFilter.type = "highshelf";
n@749 42 KFilter.gain.value = 4;
n@749 43 KFilter.frequency.value = 1480;
n@749 44
n@749 45 var HPFilter = offlineContext.createBiquadFilter();
n@749 46 HPFilter.type = "highpass";
n@749 47 HPFilter.Q.value = 0.707;
n@749 48 HPFilter.frequency.value = 60;
n@749 49 // copy Data into the process buffer
n@749 50 var processSource = offlineContext.createBufferSource();
n@749 51 processSource.buffer = buffer;
n@749 52
n@749 53 processSource.connect(KFilter);
n@749 54 KFilter.connect(HPFilter);
n@749 55 HPFilter.connect(offlineContext.destination);
n@749 56 processSource.start();
nicholas@752 57 offlineContext.oncomplete = function(renderedBuffer) {
n@749 58 // Have the renderedBuffer information, now continue processing
nicholas@752 59 if (typeof renderedBuffer.renderedBuffer == 'object') {
nicholas@752 60 renderedBuffer = renderedBuffer.renderedBuffer;
nicholas@752 61 }
n@749 62 switch(timescale)
n@749 63 {
n@749 64 case "I":
n@749 65 var blockEnergy = calculateProcessedLoudness(renderedBuffer, 400, 0.75);
n@749 66 // Apply the absolute gate
n@749 67 var loudness = calculateLoudnessFromChannelBlocks(blockEnergy);
n@749 68 var absgatedEnergy = new Array(blockEnergy.length);
n@749 69 for (var c=0; c<blockEnergy.length; c++)
n@749 70 {
n@749 71 absgatedEnergy[c] = [];
n@749 72 }
n@749 73 for (var i=0; i<loudness.length; i++)
n@749 74 {
n@749 75 if (loudness[i] >= -70)
n@749 76 {
n@749 77 for (var c=0; c<blockEnergy.length; c++)
n@749 78 {
n@749 79 absgatedEnergy[c].push(blockEnergy[c][i]);
n@749 80 }
n@749 81 }
n@749 82 }
n@749 83 var overallAbsLoudness = calculateOverallLoudnessFromChannelBlocks(absgatedEnergy);
n@749 84
n@749 85 //applying the relative gate 8 dB down from overallAbsLoudness
n@749 86 var relGateLevel = overallAbsLoudness - 8;
n@749 87 var relgateEnergy = new Array(blockEnergy.length);
n@749 88 for (var c=0; c<blockEnergy.length; c++)
n@749 89 {
n@749 90 relgateEnergy[c] = [];
n@749 91 }
n@749 92 for (var i=0; i<loudness.length; i++)
n@749 93 {
n@749 94 if (loudness[i] >= relGateLevel)
n@749 95 {
n@749 96 for (var c=0; c<blockEnergy.length; c++)
n@749 97 {
n@749 98 relgateEnergy[c].push(blockEnergy[c][i]);
n@749 99 }
n@749 100 }
n@749 101 }
n@749 102 var overallRelLoudness = calculateOverallLoudnessFromChannelBlocks(relgateEnergy);
n@749 103 buffer.lufs = overallRelLoudness;
n@749 104 }
nicholas@752 105 };
nicholas@752 106 offlineContext.startRendering();
n@749 107 }
n@749 108
n@749 109 function calculateProcessedLoudness(buffer, winDur, overlap)
n@749 110 {
n@749 111 // Buffer Web Audio buffer node
n@749 112 // winDur Window Duration in milliseconds
n@749 113 // overlap Window overlap as normalised (0.5 = 50% overlap);
n@749 114 if (buffer == undefined)
n@749 115 {
n@749 116 return 0;
n@749 117 }
n@749 118 if (winDur == undefined)
n@749 119 {
n@749 120 winDur = 400;
n@749 121 }
n@749 122 if (overlap == undefined)
n@749 123 {
n@749 124 overlap = 0.5;
n@749 125 }
n@749 126 var winSize = buffer.sampleRate*winDur/1000;
n@749 127 var olapSize = (1-overlap)*winSize;
n@749 128 var numberOfFrames = Math.floor(buffer.length/olapSize - winSize/olapSize + 1);
n@749 129 var blockEnergy = new Array(buffer.numberOfChannels);
n@749 130 for (var channel = 0; channel < buffer.numberOfChannels; channel++)
n@749 131 {
n@749 132 blockEnergy[channel] = new Float32Array(numberOfFrames);
n@749 133 var data = buffer.getChannelData(channel);
n@749 134 for (var i=0; i<numberOfFrames; i++)
n@749 135 {
n@749 136 var sigma = 0;
n@749 137 for (var n=i*olapSize; n < i*olapSize+winSize; n++)
n@749 138 {
n@749 139 sigma += Math.pow(data[n],2);
n@749 140 }
n@749 141 blockEnergy[channel][i] = sigma/winSize;
n@749 142 }
n@749 143 }
n@749 144 return blockEnergy;
n@749 145 }
n@749 146 function calculateLoudnessFromChannelBlocks(blockEnergy)
n@749 147 {
n@749 148 // Loudness
n@749 149 var loudness = new Float32Array(blockEnergy[0].length);
n@749 150 for (var i=0; i<blockEnergy[0].length; i++)
n@749 151 {
n@749 152 var sigma = 0;
n@749 153 for (var channel = 0; channel < blockEnergy.length; channel++)
n@749 154 {
n@749 155 var G = 1.0;
n@749 156 if (channel >= 4) {G = 1.41;}
n@749 157 sigma += blockEnergy[channel][i]*G;
n@749 158 }
n@749 159 loudness[i] = -0.691 + 10*Math.log10(sigma);
n@749 160 }
n@749 161 return loudness;
n@749 162 }
n@749 163 function calculateOverallLoudnessFromChannelBlocks(blockEnergy)
n@749 164 {
n@749 165 // Loudness
n@749 166 var summation = 0;
n@749 167 for (var channel = 0; channel < blockEnergy.length; channel++)
n@749 168 {
n@749 169 var G = 1.0;
n@749 170 if (channel >= 4) {G = 1.41;}
n@749 171 var sigma = 0;
n@749 172 for (var i=0; i<blockEnergy[0].length; i++)
n@749 173 {
n@749 174 blockEnergy[channel][i] *= G;
n@749 175 sigma += blockEnergy[channel][i];
n@749 176 }
n@749 177 sigma /= blockEnergy.length;
n@749 178 summation+= sigma;
n@749 179 }
n@749 180 return -0.691 + 10*Math.log10(summation);;
n@749 181 }