Real-time adaptive audio source separation

Methods and systems for audio source separation in real-time are described. In an embodiment, the present disclosure describes reading and decoding an audio source into PCM samples, fragmenting Pulse Code Modulation (PCM) samples into fragments, transforming fragments into spectrograms, performing audio source separation using a training database that includes a training dictionary and non-negative matrix factorization (NMF) to generate a set of component signals, and streaming the component signals to a playback engine. In an embodiment, a semantic equalizer graphical user allows for fading of individual component signals.

TECHNICAL FIELD

The present application relates to digital signal processing. More specifically, the example embodiment(s) of the present invention described below relate to real-time audio source separation.

BACKGROUND

In the field of digital signal processing, audio source separation is a problem where several audio signals have been mixed together into a combined signal and the goal is to recover the original component signals from the combined signal. Audio source separation has many practical applications including, but not limited to speech enhancement, speech recognition, speech denoising, voice recognition, audio post-production and remastering, spatial audio upmixing, and other audio functions. Denoising includes separating noise from speech, removing background music from speech, and removing bleed from other instruments. In the context of music production, audio source separation is sometimes referred to as “unmixing” or “de-mixing”.

Traditional audio equalizers have historically been used to emphasize or deemphasize certain content present in an audio signal. They work by boosting or attenuating certain frequencies bands present in the audio signal. The main limitation however is that they cannot discriminate and/or separate two signals that have frequencies in common. Thus, what is needed is the ability to separate mixed signals that have frequencies in common.

Several approaches have been proposed over the years to solve the audio source separation problem, however, these approaches are inadequate.

Beamforming is a technique for source separation that uses a microphone array to listen to a particular direction to capture desired signals while minimizing interfering ones. A shortcoming of the existing beamforming techniques is that they require multiple versions of the same recording to be captured with multiple microphones.

Adaptive signal processing is a technique for source separation that filters unwanted parts of a signal by self-adjusting parameters of the filter. A shortcoming of the adaptive signal processing technique is that requires prior knowledge about the statistics of the interfering signal. For example, in order to perform adaptive signal processing on a combined signal, it would be necessary to have prior knowledge regarding the types of expected noise in the combined signal.

Independent component analysis, also known as blind source separation (BSS), is a technique for source separation that uses a measure of statistical independence of the component signals in the combined signal to identify and separate the sources. This technique does not require prior knowledge about the sound sources, except for their mutual statistical independence. A shortcoming of this approach is that it assumes, among other things, that the sources have been linearly mixed into a combined signal. However, in the context of professionally produced music, processes like mastering apply several nonlinearities to a combined signal. Another shortcoming is that independent component analysis cannot separate or discriminate arbitrary sources (e.g., different types of musical instruments, such as drums, bass, piano) without prior knowledge about the sources or their semantics. Another shortcoming in the context of professionally produced music is that, in most cases, the sources are not fully statistically independent. Overall, in practice, this approach performs very poorly in the context of professional music recordings.

Classical denoising and enhancement is a technique for source separation that uses a Wiener filtering and spectral subtraction to separate the audio sources. A shortcoming of this approach is that it assumes prior knowledge of the spectral properties of the original signal and the noise.

Non-Negative Matrix Factorization (NMF) is a popular technique for processing audio, image and text. This technique has been widely used by itself or in combination with other techniques for audio source separation. A shortcoming of existing NMF techniques in the context of source separation is that it relies on a fixed set of spectral templates describing the spectral characteristics of the underlying sources in a given combined signal. Because commercial music has so much variety in terms of sounds, the effectiveness of this approach is highly dependent on the degree to which these default spectral templates represent the spectral characteristics of the underlying sources in a combined signal. In summary, the existing approach to source separation using NMF requires the arbitrary selection of a number of spectral templates and it does not use information from the input combined signal to adapt these spectral templates to better match the underlying sources, which is the key to accurate discrimination and reconstruction of the sources. What is needed is a way to improve existing NMF techniques to dynamically adapt to the combined signal with real-time performance.

While each of the figures illustrates a particular embodiment for purposes of illustrating a clear example, other embodiments may omit, add to, reorder, and/or modify any of the elements shown in the figures.

DESCRIPTION OF THE EXAMPLE EMBODIMENT(S)

In the following description, for the purposes of explanation, numerous specific details are set forth in order to provide a thorough understanding of the example embodiment(s) of the present invention. It will be apparent, however, that the example embodiment(s) may be practiced without these specific details. In other instances, well-known structures and devices are shown in block diagram form in order to avoid unnecessarily obscuring the example embodiment(s).1.0 GENERAL OVERVIEW2.0 EXAMPLE COMPUTER SYSTEM IMPLEMENTATION2.1 READING/DECODING LOGIC2.2 FRAGMENTATION LOGIC2.3 TRANSFORM LOGIC2.4 AUDIO SOURCE SEPARATION LOGIC2.4.1 IDENTIFICATION AND CLASSIFICATION OF CANDIDATE COMPONENT SIGNALS2.4.2 CUSTOM DICTIONARY CONSTRUCTION2.4.3 SPECTROGRAM DECOMPOSITION AND SIGNAL RECONSTRUCTION2.5 CONCATENATION, CROSSFADING & BUFFERING (CCB) LOGIC2.6 LOCAL BUFFERING LOGIC2.7 SEMANTIC EQUALIZER GUI2.8 AGGREGATION LOGIC2.9 PLAYBACK ENGINE2.10 APPLICATION AREAS3.0 EXAMPLE PROCESS AND ALGORITHM4.0 IMPLEMENTATION MECHANISMS—HARDWARE OVERVIEW5.0 IMPLEMENTATION MECHANISMS—SOFTWARE OVERVIEW6.0 OTHER ASPECTS OF DISCLOSURE

1.0 GENERAL OVERVIEW

In one embodiment, an audio source separation system is programmed or configured to receive an audio source file that includes stereo combined audio signal and perform source separation on the combined audio signal to generate one or more component signals. The combined audio signal may represent a mix of various signals. For example, the individual component signals may correspond to individual musical instruments, voices, or other sound sources.

Audio source separation is performed by reading and decoding an audio source file into raw audio samples. The raw audio samples may be fragmented into raw audio fragments to provide for real-time processing. Each raw audio fragment may be transformed into a spectrogram. The spectrogram may be analyzed using Non-negative matrix factorization (NMF) techniques along with a training dictionary of reference spectral templates to identify specific content in the spectrogram. Specific content may include musical instruments, voices, or any other unique audio source. In an embodiment, the specific content is identified as component spectral templates. The component spectral templates can then be used to reconstruct the component signals of the original audio source into component signal fragments. The component signal fragments may be concatenated, cross-faded, and buffered for playback at a playback device.

In an embodiment of the present invention, a graphical user interface is provided that allows the end user to set the gains of the separated component signals independently. This is implemented by providing a user interface with a number of volume faders. By moving each fader from minimum gain (0% volume) to maximum gain (100% volume) the user can select the desired amount of signal from each component signal, which are recombined on real time and delivered in stereo format to the output listening device.

The present disclosure thus provides for real-time audio source separation of a combined audio signal into component signals. This real-time audio source separation provides a variety of useful application areas, such as visualization and manipulation of component signals in mobile devices, augmented reality, and virtual reality contexts. In an embodiment, “real-time” in this context means that processing of audio source separation is faster than the playback speed so as to prevent a bottleneck during the middle of playback. The present system provides an adaptive, highly parallel, multithreaded approach to audio source separation that can be used for breaking up a combined signal to component signals in real-time and provide for rapid playback.

2.0 EXAMPLE COMPUTER SYSTEM IMPLEMENTATION

FIG. 1illustrates an example audio source separation system100in which the techniques described herein may be practiced, according to some embodiments. In the example ofFIG. 1, audio source separation system100is programmed or configured to read an audio source file and perform audio source separation on the audio source file in real-time. Audio source separation system100may be implemented across one or more physical or virtual computing devices, none of which is intended as a generic computer, since it is loaded with instructions in a new ordered combination as otherwise disclosed herein to implement the functions and algorithms of this disclosure. The example components of audio source separation system100shown inFIG. 1are implemented at least partially by hardware at one or more computing devices, such as one or more hardware processors executing stored program instructions stored in one or more memories for performing the functions that are described herein. Or, one or more virtual machine instances in a shared computing facility such as a cloud computing center may be used. The functions described herein are intended to indicate operations that are performed using programming in a special-purpose computer or general-purpose computer, in various embodiments. Audio source separation system100illustrates only one of many possible arrangements of components configured to execute the programming described herein. Other arrangements may include fewer or different components, and the division of work between the components may vary depending on the arrangement.

In an embodiment, audio source separation system100comprises server102and client104, which may be communicatively coupled over a network. The network may be any communication network, such as communication over the Internet or a private network. In an embodiment, the network may include a wireless network via one or more wireless transmission protocols, such as 802.11 a/b/g/n, Bluetooth, LTE, 3G, CDMA, or any other wireless protocol.

In another embodiment, one or more components of server102and/or client104may be implemented on the same device. Thus, the separation of components of server102and/or client104is merely illustrative of one possible client-server architecture. In other embodiments, some or all of the components may be implemented on the same computing device.

Server102may optionally include reading/decoding logic110. Reading/decoding logic110is programmed or configured to read an audio source105and generate a plurality of raw audio samples based on the audio source105. For example, in one embodiment, each second of audio taken from an audio source that is read from a commercial compact disk (CD) music may include 44100 raw audio samples. The reading/decoding logic110is communicatively coupled to fragmentation logic112and can send audio samples to fragmentation logic112for further processing.

Audio source105may be any audio stream or audio file accessible to server102. For example, audio source105may be an audio stream or audio file implemented in an uncompressed audio format, a lossless compressed audio format, or a lossy compressed audio format. Examples of file formats of audio source105include, but are not limited to: .AAC, .AIFF, .AU, .DVF, .M4A, M4P, MP3, .OGG, .RAW, .WAV, and/or .WMA. Audio source105represents a combined audio signal that includes one or more component audio signals. In an embodiment, audio source105may be a stereo signal.

Server102includes reading/decoding logic110that is programmed or configured to read audio source105, decode audio source105if necessary, and generate a plurality of audio samples based on audio source105. A raw audio sample is an audio sample taken in the time domain. For example, a pulse-code modulation (PCM) sample is an example of a raw audio sample. A PCM signal is a sequence of digital audio samples containing data that provides the necessary information to reconstruction an original analog signal. Each PCM sample represents the amplitude of the PCM signal at a specific point in time, and the PCM samples are uniformly spaced in time. In one embodiment, the amplitude is the only information explicitly stored in the PCM sample, and it is typically stored as either an integer or a floating point number, and encoded as a binary number with a number of fixed number of digits which represents the sample's bit depth.

In some embodiments, the audio source105does not need to be decoded as the audio source105itself includes a raw audio sample. For example, a .WAV file format includes raw audio samples embedded directly in the file, thus, decoding is not necessary. On the other hand, the .MP3 file format is a compressed, lossy format, therefore decoding of the audio source105is necessary using an appropriate MP3 decoder in order to generate the raw audio sample(s) from the audio source105. Reading/decoding logic110can be programmed or configured to include various audio source105readers and decoders for various audio file formats and may select the appropriate readers and/or decoders for generating a raw audio sample based on the particular file format of the audio source105.

Once the raw audio samples have been generated, reading/decoding logic110can send raw audio sample(s) to fragmentation logic112.

Server102includes fragmentation logic112which is programmed or configured to select a fixed number of raw audio samples for further processing. Given the large amount of data that can be present in a full audio source105, it is inefficient to process all raw audio samples simultaneously. Processing all such raw audio samples simultaneously would reduce the ability to provide audio source separation in real-time. Instead, by fragmenting the raw audio samples, fragmentation logic112can break up the processing into smaller chunks that will allow for real-time audio source separation and playback at client104. Thus, by fragmenting the incoming raw audio samples, fragmentation logic112allows the audio source separation system100to continue processing the incoming data as an ongoing pipeline.

Fragmentation logic112is programmed or configured to select a fixed number of raw audio samples for further processing based on a fragment size. The fragment size used by fragmentation logic112may be programmed or configured via a configuration setting. Varying the fragment size can provide various tradeoffs. A smaller fragment size will ensure that the downstream audio source separation logic116receives data for processing sooner thereby providing quicker playback to client104, but at the expense of having to process more fragments. Alternatively, a larger fragment size may delay the processing of the fragments by the audio source separation logic116thereby delaying playback by client104, however, it will ensure that audio source separation logic116has fewer fragments to process thereby reducing computing.

In one embodiment, a fragment size is selected that represents greater than or equal to 1.5 second and less than or equal to 3.0 seconds of audio. This is roughly equivalent to 44,100-132,300 samples at a sampling rate of 44100 Hz.

The fragmentation logic112continually fragments the incoming raw audio samples and sending the fragments of the raw audio samples to transform logic114for further processing.

Server102includes transform logic102that is programmed or configured to receive a raw audio fragment from fragmentation logic112and generate a spectrogram for the fragment by computing a time-frequency transform of the fragment. A spectrogram is a representation of the spectrum of frequencies in a sound as they vary with time. In an embodiment, the spectrogram is a complex-valued transform made up of complex numbers in the form of z=a+ix, where “a” represents the real part of the number and “x” represents the imaginary part of the number.

In an embodiment, a magnitude matrix may be calculated based on a spectrogram, by calculating the absolute value of the spectrogram where the y-axis represents frequency and the x-axis represents timeframes. The value in each time-frequency slot represent a specific energy level of the input signal. The magnitude matrix may be used for calculations of audio source separation and decomposition using NMF by the audio source separation logic116.

In another embodiment, a phase matrix may be calculated based on a spectrogram, which represents phase information of the audio signal. The phase matrix includes phase values which are calculated by computing the angle of each complex value as: atan(imag(z)/real(z)). In this equation, the phase value is calculated as the arctangent of the quotient between the imaginary part over the real part of each complex number. The phase matrix may be used by audio source separation logic116to combine each individual component signal's magnitude matrix with the original combined signal's phase matrix to be able to reconstruct each component signal in the combined signal.

Transform logic102may be programmed or configured to generate a spectrogram for a fragment by using any time-frequency transform. Examples of time-frequency transforms include Short-Time Fast Fourier Transforms (STFT), Constant-Q transforms, Wavelets, filter banks, or any other similar transform function.

In an embodiment, transform logic114may optionally apply non-linear interpolation to the spectrogram before sending the spectrogram to audio source separation logic116. Applying non-linear interpolation to the spectrogram is intended to provide more weight and improved resolution to frequency bands where the human auditory system has better resolution, thereby improving the ability to perform source separation for those frequencies. The human auditory system generally has a hearing range of 20 Hz to 20 kHz, although under ideal laboratory conditions the hearing range may be as wide as 12 Hz to 28 kHz. Thus, from a practical perspective for human listening, audio source separation is only relevant for frequencies in the human hearing range.

The non-linear interpolation that is applied by transform logic114may use one of many different interpolations based on the context of human auditory scales, sometimes referred to as psychoacoustic scales. Examples of such psychoacoustic scales include the mel scale, the Bark Scale, and/or the Equivalent Rectangular Bandwidth (ERB) scale. Thus, transform logic114may use one or more of these psychoacoustic scales to apply non-linear interpolation to the spectrogram.

In an embodiment, the spectrogram generated by transform logic114may be sent to audio source separation logic116for further processing. In another embodiment, the magnitude matrix and/or phase matrix generated by transform logic114may be sent to audio source separation logic116for further processing.

2.4 Audio Source Separation Logic

Server102includes audio source separation logic116. Audio source separation logic116takes a spectrogram or magnitude matrix of a spectrogram, analyzes it, and generates N output component signals representing the individual component signal fragments present in the spectrogram. A component signal can be any semantically meaningful element or pattern in the in the spectrogram, such as a musical instrument (e.g., piano, snare drum, bass drum, guitar, violin, etc.), a particular note of a musical instrument (e.g. each note of a piano may be a different instrument, each note of a violin may be a different instrument, etc.), a particular voice (e.g., for multiple human voices each voice may be a different instrument, etc.), or any other pattern that is identifiable in the spectrogram (e.g., sound of an explosion, an ambulance siren, etc.).

Audio source separation logic116identification and classification of component signals on the input spectrogram, builds a custom dictionary best suited for the input spectrogram, performs spectrogram decomposition using the custom dictionary to generate component spectrograms from the input spectrogram, and finally performs signal reconstruction on the component spectrograms to generate component signal fragments. The component signal fragments are audio samples and may be formatted as PCM sample fragments according to one embodiment. The component signal fragments may then be sent to concatenation, crossfading & buffering logic118for further processing.

2.4.1 Identification and Classification of Candidate Component Signals

Audio source separation logic116is programmed or configured to perform identification and classification of component signals on the input spectrogram to identify one or more candidate component signals found in the input spectrogram. Audio source separation logic116may be communicatively coupled to training database115. In one embodiment, training database116may be accessible to audio source separation logic116via an exposed application programming interface (API). Training database116is a repository of information that may assist in the identification, classification, and source separation of an audio signal.

For example, in one embodiment, training database116may store information regarding expected metadata for component signal identification. Audio source105may include various metadata. For example, an audio source105in an .MP3 format may include various metadata for a musical track, including, but not limited to title, artist, album, year, and genre. Other formats of audio source105may store different types of metadata. The metadata of an audio source105may be used to identify candidate component signals in the audio source105. For example, an audio source that includes metadata indicating a “classical” musical genre is more likely to include musical instruments featuring violins and cellos. On the other hand, an audio source that includes metadata indicating a “rock” musical genre is more likely to include musical instruments featuring electric guitars and bass instruments. Training database116may include data that correlates audio source metadata to expected component signals in the audio source. Audio source separation logic116can thus use the stored data in training database116to identify candidate component signals in the audio source105.

In another embodiment, audio source separation logic116may analyze the input spectrogram to identify candidate instruments. Training database116may store a training dictionary of reference spectral templates that can be used to identify candidate instruments in a spectrogram. A spectral template is a vector containing a specific spectral pattern of energy across different frequencies. A spectral template thus contains a specific statistical distribution of magnitude or energy across frequencies and thus represents a spectral signature of a specific sound in a semantic context. For example, a spectral template for a kick drum will concentrate higher energy in the low frequency area compared to a spectral template for a hi-hat cymbal, which will concentrate most of the energy in a high frequency range. The training database116may store a training dictionary of pre-stored reference spectral templates that each represents a signature of particular candidate component signal, such as a musical instrument, a note on a particular musical instrument, a chord on a particular musical instrument, a voice, and/or some other distinctive sound effect. For example, in one embodiment, a different spectral template may be stored for every note that can be played on a piano.

Audio source separation logic116may then use the spectral templates in training database116to identify and classify the component signals present in a particular input spectrogram. For example, in one embodiment, audio source separation logic116may use non-negative matrix factorization (NMF) to identify a subset of spectral templates from the training dictionary that are most often found in the input spectrogram. NMF is an algorithm in linear algebra to provide matrix factorization. In NMF, a matrix V is approximated by two smaller matrices W and H, which, when multiplied, approximately reconstruct V. Thus, in NMF, two or more submatrices are multiplied together to get a final matrix using the following equation: W*H=V. The NMF algorithm can be applied to this equation to determine one or both of the submatrices. In the present case, the input spectrogram is represented as V. The reference spectral templates in the training dictionary are represented as basis vectors that collectively make up matrix W. The NMF algorithm performs a number of iterations to estimate values for H so that when multiplied by the fix matrix W, best approximates the input matrix V. The number of iterations can be arbitrarily chosen or alternatively, can be dynamically set based on a predefined convergence criteria. The activations matrix H is a matrix that identifies which spectral templates of all the spectral templates in the training dictionary are activated at each timeframe and with which gain. This configuration is sometimes known as supervised NMF because the system provides the matrix W in the form of reference spectral templates, versus the pure version of NMF, where both W and H are randomly initialized and the algorithm estimates the values of both W and H.

In one embodiment, audio source separation logic116may further use correlation analysis on activations matrix H to identify which spectral templates of the ones identified have a minimum correlation with spectral templates from other families of component signals. For example, in one embodiment, the Pearson correlation coefficient may be used for to measure correlation between variables in a correlation analysis. A family of component signals may represent spectral templates that are interrelated. For example, each note for a piano keyboard may be represented by its own spectral template and the collection of all such notes may be the family of component signals. Templates corresponding to different drums sounds (e.g., kick drum, snare drum, hi-hat cymbal) can be another example of a family. A goal of the correlation analysis is to load a subset of spectral templates that are used consistently but that have minimum correlation with templates for other families. Minimum correlation means that there is more spectral dissimilarity between spectral templates from different families, therefore minimizing interference between the different families when later reconstructing signals. This correlation analysis directly impacts the quality of separation and it is a key step towards adaptive source separation, which takes into account statistics from the input combined signal itself

In another embodiment, instrument identification and classification may be performed using deep learning techniques.

The final result of the component signal identification and classification analysis by the audio source separation logic116is a selection of a subset of candidate spectral templates from training dictionary that are present and most commonly found in the input spectrogram while maximizing inter-family discrimination of templates by the algorithm for a given input combined signal.

2.4.2 Custom Dictionary Construction

Based on the information determined during instrument identification and classification, audio source separation logic116is programmed or configured to generate a custom dictionary of spectral templates. The custom dictionary of spectral templates describes the special characteristics of the component signals present in the current fragment of the input spectrogram. Thus, the custom dictionary can be used for source separation and reconstruction. The custom dictionary is a subset of all the spectral templates available in the training dictionary. By focusing on this subset of spectral templates instead of all spectral templates in the training dictionary, the audio source separation logic116can improve performance and limit its subsequent source separation and reconstruction processing to instruments present in the current fragment of the input spectrogram without having to process unnecessary spectral templates in the training dictionary, while maximizing inter-family discrimination of templates by the algorithm for a given input combined signal. The use of a custom dictionary will reduce interference caused by using an incorrect spectral template to model a component signal that does not belong to the same family. For example, using a using a drum spectral template to model part of the piano can cause unwanted interference.

2.4.3 Spectrogram Decomposition and Signal Reconstruction

Audio source separation logic116is programmed or configured analyze the spectrogram using an audio source separation algorithm based on NMF, decompose the spectrogram, and perform signal reconstruction of the component signals. As described earlier, in NMF, a matrix V is approximated by two smaller matrices W and H, which, when multiplied, approximately reconstruct V: W*H=V. In the context of spectrogram decomposition, the matrix V is the magnitude matrix of the spectrogram of the input fragment received from transform logic114. Many other implementations of NMF may iteratively update both W and H to obtain the final approximation of V, however, in one embodiment of the present disclosure, a semi-supervised NMF implementation is applied instead. In a semi-supervised NMF implementation, the matrix W is dynamically constructed based on the custom dictionary of spectral templates that was previously calculated by audio source separation logic116and is left unchanged by the separation algorithm. W is a concatenation of the 1-dimensional spectral templates in the custom dictionary. The separation algorithm using NMF techniques thus only updates and iterates the matrix H. Once the separation algorithm converges for the NMF equation, the result is an activations matrix H.

Audio source separation logic116can then using the activations matrix H to reconstruct the signals corresponding to the component signals. In one embodiment, audio source separation logic116uses the activations matrix H to perform a direct reconstruction of the signals. In direct reconstruction, the reconstruction for each family of component signals is achieved by multiplying a subset of the columns of the W matrix from the separation algorithm by a corresponding subset of the rows of the H matrix. For example, in one embodiment, the W matrix includes 30 columns where columns 1-10 represent a drum family of spectral templates, column 11-20 represent a guitar family spectral templates, and columns 21-30 represent a piano family of spectral templates. Thus, in order to reconstruct the spectrogram corresponding to the drums, columns 1-10 of the W matrix and multiplied against rows 1-10 of the H matrix. This reconstructed spectrogram thus represents the drums musical instrument. The reconstructed magnitude matrix of the spectrogram can then then be used as an input to an Inverse Short Time Fourier Transform (ISTFT) to generate a reconstructed signal that represent the drums musical instrument, by combining the magnitude matrix with the original phase matrix of the spectrogram. Similar processing can be performed for the guitar musical instrument and the piano musical instrument. Direct reconstruction gives better isolation of signals and less interference between signals at the expense of some distortion and coloring in the sound.

In another embodiment, audio source separation logic116uses the activations matrix H to perform a source filtering using a soft mask to reconstruct the signals corresponding to musical instruments. A soft mask is a filter matrix M that can be calculated as follows:

In this equation, wiand hiare the spectral templates and activations, respectively, corresponding to the ithsource, where K denotes the number of component signals to be reconstructed. The filter matrix M can be applied to the original mixed signal's spectrogram to obtain the reconstructed component signals. Source filtering provides better quality of sound but may introduce greater interference between the component signals.

In one embodiment, audio source separation logic116is programmed or configured to use direct reconstruction for reconstructing individual component signals. In another embodiment, audio source separation logic116is programmed or configured to use source filtering for reconstructing the individual component signals. In another embodiment, audio source separation logic116is programmed or configured to combine both direct reconstruction source filtering to reconstruct individual component signals. In the combined approach, audio source separation logic116performs both direct reconstruction and source filtering in parallel and takes the weighted average of the reconstructed signals.

The final output of the audio source separation logic116is a series of 1 to N component signal fragments, each one corresponding to one musical instrument or other source identifier.

Server102may include concatenation, crossfading & buffering (CCB) logic118that is programmed or configured to concatenate subsequent component signal fragments for a particular component signal, crossfade between component signal fragments and buffer each individual component signal. CCB logic118receives each fragment for a given component signal and concatenates them into a dedicated buffer, each corresponding to a particular component signal, such as a musical instrument. Since the separation is done on independent fragments of the component signal, a direct concatenation of the separated signal fragments might create audible discontinuities in the separated signal during playback. In order to avoid that undesired effect, a cross-fading mechanism is applied to these fragments prior to the concatenation, which avoids the audible discontinuity.FIG. 2illustrates an example signal200for a particular musical component signal where a crossfading effect is applied between fragments. In signal200, frame n corresponds to a component signal fragment for a particular musical instrument, and frame n−1 represents the next subsequent component signal fragment for the particular musical instrument. A crossfading algorithm is applied to smooth the overlapping region between frame n and frame n−1, thereby improving signal quality. The resulting output of the CCB logic118is a series of N audio signals where each signal corresponds to a component signal of the audio source105.

2.6 Local Buffering Logic

In an embodiment, client104includes local buffering logic120that receives the N component audio signals from CCB logic118and buffers them locally in preparation for audio playback.

Client104may be communicatively coupled to semantic equalizer graphical user interface (GUI)122. Semantic equalizer GUI122allows a user to manipulate the gain(s) applied to each component audio signal.FIG. 3illustrates an example of the semantic equalizer GUI300, according to one embodiment in the musical context. Each volume fader in semantic equalizer GUI300corresponds to one component signal. For example, in this example, each volume fader corresponds to a musical instrument, in this case “Drums”, “Bass guitar”, “Strings” and “Vocals”. In this example, semantic equalizer GUI122can receive one or more user inputs, such as touch gesture, mouse gesture, or keyboard input to modify a gain level applied to a particular component audio signal such as a signal that corresponds to a musical instrument. By modifying the component audio signal, the fader's current value can be used to multiply the corresponding component audio signals stored in the local buffers. Moreover, since each instrument is itself stereo, meaning it has a right channel and a left channel and the addition is done channel-wise, this leads to a stereo mix.

Client104includes aggregation logic124that is programmed or configured to take the component audio signals modified by gains, and aggregate them into a new stereo audio mixture in which the amount or presence of each instrument is regulated by the values of the faders in the semantic equalizer GUI122.

2.9 Playback Engine

Client104includes a playback engine126that is programmed or configured to provide playback of the final stereo audio mixture. In one embodiment, playback engine126may be implemented, for example, via a web client, a native mobile application, a virtual reality audio system, or any other audio output device.

2.10 Application Areas

The present invention can separate an audio signal incrementally while it is being captured, which makes the technology suitable to be used with audio that is being streamed over the internet or signals that are being captured in real time. The invention achieves this by automatically dividing the input signal into smaller pieces, separating each piece independently and then recombining the individual outputs. This parallelized implementation provides the ability to perform audio source separation in real-time with no delays and minimal computing power.

3.0 EXAMPLE PROCESS AND ALGORITHM

FIG. 4illustrates a flow diagram of an example process400for real-time audio source separation using audio source separation system100, according to one embodiment.FIG. 4is intended to disclose an algorithm or functional description that may be used as a basis of writing computer programs to implement the functions that are described herein and which cause a computer to operate in the new manner that is disclosed herein. Further,FIG. 4is provided to communicate such an algorithm at the same level of detail that is normally used, by persons of skill in the art to which this disclosure is directed, to communicate among themselves about plans, designs, specifications and algorithms for other computer programs of a similar level of complexity.

The process400may start at step402. At step402, reading/decoding logic110is programmed or configured to read an audio source105that represents a combined audio signal comprising a plurality of component audio signals. The audio source105may be formatted in any available audio format. In an embodiment, reading/decoding logic110is programmed or configured to generate a raw audio sample based on audio source105. For example, in one embodiment, reading/decoding logic110may use one or more decoders to decode audio source105into raw audio samples. A raw audio sample is an audio sample taken in the time domain, such as a PCM sample. The reading/decoding logic110sends the generated raw audio samples to fragmentation logic112. The process400may then proceed to step404.

In step404, fragmentation logic112is programmed or configured to fragment the raw audio samples received from reading/decoding logic110. Fragmentation logic112may use a preconfigured fragment size to fragment the raw audio samples into a plurality of audio fragments. In one embodiment, each audio fragment may represent greater than or equal to 1.5 seconds of audio and less than or equal to 3.0 seconds of audio. The fragmentation of the raw audio samples into audio fragments is necessary to allow the audio source separation system100to process the audio source105in real-time. The fragmentation logic112sends each audio fragment to transform logic114for further process. The process400may then proceed to step406.

In step406, transform logic114is programmed or configured to transform audio fragments into a spectrogram. For each audio fragment of raw audio samples received from fragmentation logic112, transform logic114can apply a transform function to the audio fragment to generate a spectrogram. For example, in one embodiment, transform logic114can apply a Short-Time Fast Fourier Transform function to a fragment to generate a spectrogram. In an embodiment, transform logic114may optionally be further programmed or configured to apply non-linear interpolation to the spectrogram to provide more weight and improved resolution to the auditory frequencies in the human auditory range. For example, in one embodiment, transform logic114may perform non-linear interpolation to the spectrogram using the mel scale, which is one well known psychoacoustic scale for the human auditory range. The resulting spectrogram is sent to audio source separation logic116. In an embodiment, a magnitude matrix and/or phase matrix of the spectrogram is sent to the audio source separation logic116. The process400may then proceed to step408.

In step408, audio source separation logic116is programmed to identify candidate component signals in a spectrogram. In one embodiment, audio source separation logic116may compare metadata from audio source105to a set of metadata in training database115to identify one or more candidate component signals likely to exist in the spectrogram. In another embodiment, audio source separation logic116is programmed or configured to use a training dictionary comprising a plurality of reference spectral templates for known component signal in a NMF analysis to identify the candidate component signals present in the input spectrogram. In one embodiment, audio source separation logic116may programmed or configured to use correlation analysis to identify the component signals in the spectrogram with minimum correlation with spectral templates from other families of component signals. The final result of this analysis is a selection of candidate spectral templates from the training dictionary that are present in the spectrogram and most commonly found in the spectrogram, while minimizing inter-family correlation between templates, thereby minimizing the resulting interference after reconstruction of the sources. The process400may then proceed to step410.

In step410, the audio source separation logic116is programmed or configured to build a custom dictionary of spectral templates for the spectrogram. The custom dictionary thus includes the spectral templates for the candidate component signals identified in the prior step. The process400may then proceed to step412.

In step412, the audio source separation logic116is programmed or configured to use NMF to decompose the magnitude matrix of the spectrogram and reconstruct the component signals. The source separation algorithm uses the custom dictionary and the spectrogram along with NMF to decompose the spectrogram into an activations matrix H. In an embodiment, activations matrix H is a matrix of concatenated spectral templates that represent the individual component signals of the input spectrogram. Audio source separation logic116then can use this activations matrix H to reconstruct the component spectrograms of the input spectrogram. In one embodiment, audio source separation logic116uses direct reconstruction to generate the component spectrograms, and then uses Inverse Short Time Fourier Transform (ISTFT) to generate the component signal fragments for each component spectrogram. In another embodiment, audio source separation logic116uses source filtering using a soft mask to generate the component signal fragments. In another embodiment, both direct reconstruction and source filtering are performed in parallel to generate intermediate component signal fragments, which are then used to generate a weighted average for the final component signal fragments. The reconstructed component signals fragments may then be sent to CCB logic118for further processing. The process F00 may then proceed to step414.

In step414, CCB logic118is programmed or configured to concatenate the input component signal fragments for each component signal. For example, component signal fragment stream1may include component signal fragments for a piano. Each fragment in the component signal fragment stream1may be concatenated with one another to generate a full stream for the piano. In one embodiment, CCB logic118is programmed to perform crossfading to smooth out the fragments of a reconstructed component signal. In one embodiment, CCB logic118may buffer the individual components signals1through N. CCB logic118may be programmed or configured to send the component signals to client104. The process400may then proceed to step416.

In step416, the local buffering loci120receives the individual component signal streams from CCB logic118and buffers them locally at the client104for local playback. The process400may then proceed to step418.

In step418, semantic equalizer GUI122may be used to receive user input on a volume fader to modify the gains to a particular component signal. By modifying the fader's value for a component signal on the semantic equalizer GUI122, the gain level of the corresponding component signal is modified to multiply the corresponding component signal received from local buffing logic120. The modified component signals are then received at aggregation logic124. The process400may then proceed to step420.

In step420, aggregation logic124is programmed or configured to receive the component signals modified by gains and aggregate them into a new stereo audio mixture, where the presence of each component signal is regulated by the values of the faders in the semantic equalizer GUI122. The process400may then proceed to step422.

In step422, the playback engine126is programmed or configured to play back the final stereo audio mixture provided by aggregation logic124. The process400may then end.

Referring now toFIG. 5, it is a block diagram that illustrates a basic computing device500in which the example embodiment(s) of the present invention may be embodied. Computing device500and its components, including their connections, relationships, and functions, is meant to be exemplary only, and not meant to limit implementations of the example embodiment(s). Other computing devices suitable for implementing the example embodiment(s) may have different components, including components with different connections, relationships, and functions.

Computing device500may include a bus502or other communication mechanism for addressing main memory506and for transferring data between and among the various components of device500.

Computing device500may also include one or more hardware processors504coupled with bus502for processing information. A hardware processor504may be a general purpose microprocessor, a system on a chip (SoC), or other processor.

Main memory506, such as a random access memory (RAM) or other dynamic storage device, also may be coupled to bus502for storing information and software instructions to be executed by processor(s)504. Main memory506also may be used for storing temporary variables or other intermediate information during execution of software instructions to be executed by processor(s)504.

Software instructions, when stored in storage media accessible to processor(s)504, render computing device500into a special-purpose computing device that is customized to perform the operations specified in the software instructions. The terms “software”, “software instructions”, “computer program”, “computer-executable instructions”, and “processor-executable instructions” are to be broadly construed to cover any machine-readable information, whether or not human-readable, for instructing a computing device to perform specific operations, and including, but not limited to, application software, desktop applications, scripts, binaries, operating systems, device drivers, boot loaders, shells, utilities, system software, JAVASCRIPT, web pages, web applications, plugins, embedded software, microcode, compilers, debuggers, interpreters, virtual machines, linkers, and text editors.

Computing device500also may include read only memory (ROM)508or other static storage device coupled to bus502for storing static information and software instructions for processor(s)504.

One or more mass storage devices510may be coupled to bus502for persistently storing information and software instructions on fixed or removable media, such as magnetic, optical, solid-state, magnetic-optical, flash memory, or any other available mass storage technology. The mass storage may be shared on a network, or it may be dedicated mass storage. Typically, at least one of the mass storage devices510(e.g., the main hard disk for the device) stores a body of program and data for directing operation of the computing device, including an operating system, user application programs, driver and other support files, as well as other data files of all sorts.

Computing device500may be coupled via bus502to display512, such as a liquid crystal display (LCD) or other electronic visual display, for displaying information to a computer user. In some configurations, a touch sensitive surface incorporating touch detection technology (e.g., resistive, capacitive, etc.) may be overlaid on display512to form a touch sensitive display for communicating touch gesture (e.g., finger or stylus) input to processor(s)504.

An input device514, including alphanumeric and other keys, may be coupled to bus502for communicating information and command selections to processor504. In addition to or instead of alphanumeric and other keys, input device514may include one or more physical buttons or switches such as, for example, a power (on/off) button, a “home” button, volume control buttons, or the like.

Another type of user input device may be a cursor control516, such as a mouse, a trackball, or cursor direction keys for communicating direction information and command selections to processor504and for controlling cursor movement on display512. This input device typically has two degrees of freedom in two axes, a first axis (e.g., x) and a second axis (e.g., y), that allows the device to specify positions in a plane.

While in some configurations, such as the configuration depicted inFIG. 5, one or more of display512, input device514, and cursor control516are external components (i.e., peripheral devices) of computing device500, some or all of display512, input device514, and cursor control516are integrated as part of the form factor of computing device500in other configurations.

Functions of the disclosed systems, methods, and modules may be performed by computing device500in response to processor(s)504executing one or more programs of software instructions contained in main memory506. Such software instructions may be read into main memory506from another storage medium, such as storage device(s)510. Execution of the software instructions contained in main memory506cause processor(s)504to perform the functions of the example embodiment(s).

While functions and operations of the example embodiment(s) may be implemented entirely with software instructions, hard-wired or programmable circuitry of computing device500(e.g., an ASIC, a FPGA, or the like) may be used in other embodiments in place of or in combination with software instructions to perform the functions, according to the requirements of the particular implementation at hand.

The term “storage media” as used herein refers to any non-transitory media that store data and/or software instructions that cause a computing device to operate in a specific fashion. Such storage media may comprise non-volatile media and/or volatile media. Non-volatile media includes, for example, non-volatile random access memory (NVRAM), flash memory, optical disks, magnetic disks, or solid-state drives, such as storage device510. Volatile media includes dynamic memory, such as main memory506. Common forms of storage media include, for example, a floppy disk, a flexible disk, hard disk, solid-state drive, magnetic tape, or any other magnetic data storage medium, a CD-ROM, any other optical data storage medium, any physical medium with patterns of holes, a RAM, a PROM, and EPROM, a FLASH-EPROM, NVRAM, flash memory, any other memory chip or cartridge.

Various forms of media may be involved in carrying one or more sequences of one or more software instructions to processor(s)504for execution. For example, the software instructions may initially be carried on a magnetic disk or solid-state drive of a remote computer. The remote computer can load the software instructions into its dynamic memory and send the software instructions over a telephone line using a modem. A modem local to computing device500can receive the data on the telephone line and use an infrared transmitter to convert the data to an infrared signal. An infrared detector can receive the data carried in the infrared signal and appropriate circuitry can place the data on bus502. Bus502carries the data to main memory506, from which processor(s)504retrieves and executes the software instructions. The software instructions received by main memory506may optionally be stored on storage device(s)510either before or after execution by processor(s)504.

Computing device500also may include one or more communication interface(s)518coupled to bus502. A communication interface518provides a two-way data communication coupling to a wired or wireless network link520that is connected to a local network522(e.g., Ethernet network, Wireless Local Area Network, cellular phone network, Bluetooth wireless network, or the like). Communication interface518sends and receives electrical, electromagnetic, or optical signals that carry digital data streams representing various types of information. For example, communication interface518may be a wired network interface card, a wireless network interface card with an integrated radio antenna, or a modem (e.g., ISDN, DSL, or cable modem).

Network link(s)520typically provide data communication through one or more networks to other data devices. For example, a network link520may provide a connection through a local network522to a host computer524or to data equipment operated by an Internet Service Provider (ISP)526. ISP526in turn provides data communication services through the world wide packet data communication network now commonly referred to as the “Internet”528. Local network(s)522and Internet528use electrical, electromagnetic or optical signals that carry digital data streams. The signals through the various networks and the signals on network link(s)520and through communication interface(s)518, which carry the digital data to and from computing device500, are example forms of transmission media.

Computing device500can send messages and receive data, including program code, through the network(s), network link(s)520and communication interface(s)518. In the Internet example, a server530might transmit a requested code for an application program through Internet528, ISP526, local network(s)522and communication interface(s)518.

FIG. 6is a block diagram of a basic software system600that may be employed for controlling the operation of computing device500. Software system600and its components, including their connections, relationships, and functions, is meant to be exemplary only, and not meant to limit implementations of the example embodiment(s). Other software systems suitable for implementing the example embodiment(s) may have different components, including components with different connections, relationships, and functions.

Software system600is provided for directing the operation of computing device500. Software system600, which may be stored in system memory (RAM)506and on fixed storage (e.g., hard disk or flash memory)510, includes a kernel or operating system (OS)610.

The OS610manages low-level aspects of computer operation, including managing execution of processes, memory allocation, file input and output (I/O), and device I/O. One or more application programs, represented as602A,602B,602C . . .602N, may be “loaded” (e.g., transferred from fixed storage510into memory506) for execution by the system600. The applications or other software intended for use on device600may also be stored as a set of downloadable computer-executable instructions, for example, for downloading and installation from an Internet location (e.g., a Web server, an app store, or other online service).

Software system600includes a graphical user interface (GUI)615, for receiving user commands and data in a graphical (e.g., “point-and-click” or “touch gesture”) fashion. These inputs, in turn, may be acted upon by the system600in accordance with instructions from operating system610and/or application(s)602. The GUI615also serves to display the results of operation from the OS610and application(s)602, whereupon the user may supply additional inputs or terminate the session (e.g., log off).

OS610can execute directly on the bare hardware620(e.g., processor(s)504) of device500. Alternatively, a hypervisor or virtual machine monitor (VMM)630may be interposed between the bare hardware620and the OS610. In this configuration, VMM630acts as a software “cushion” or virtualization layer between the OS610and the bare hardware620of the device500.

VMM630instantiates and runs one or more virtual machine instances (“guest machines”). Each guest machine comprises a “guest” operating system, such as OS610, and one or more applications, such as application(s)602, designed to execute on the guest operating system. The VMM630presents the guest operating systems with a virtual operating platform and manages the execution of the guest operating systems.

In some instances, the VMM630may allow a guest operating system to run as if it is running on the bare hardware620of device500directly. In these instances, the same version of the guest operating system configured to execute on the bare hardware620directly may also execute on VMM630without modification or reconfiguration. In other words, VMM630may provide full hardware and CPU virtualization to a guest operating system in some instances.

In other instances, a guest operating system may be specially designed or configured to execute on VMM630for efficiency. In these instances, the guest operating system is “aware” that it executes on a virtual machine monitor. In other words, VMM630may provide para-virtualization to a guest operating system in some instances.

6.0 OTHER ASPECTS OF DISCLOSURE

Although some of the figures described in the foregoing specification include flow diagrams with steps that are shown in an order, the steps may be performed in any order, and are not limited to the order shown in those flowcharts. Additionally, some steps may be optional, may be performed multiple times, and/or may be performed by different components. All steps, operations and functions of a flow diagram that are described herein are intended to indicate operations that are performed using programming in a special-purpose computer or general-purpose computer, in various embodiments. In other words, each flow diagram in this disclosure, in combination with the related text herein, is a guide, plan or specification of all or part of an algorithm for programming a computer to execute the functions that are described. The level of skill in the field associated with this disclosure is known to be high, and therefore the flow diagrams and related text in this disclosure have been prepared to convey information at a level of sufficiency and detail that is normally expected in the field when skilled persons communicate among themselves with respect to programs, algorithms and their implementation.

In the foregoing specification, the example embodiment(s) of the present invention have been described with reference to numerous specific details. However, the details may vary from implementation to implementation according to the requirements of the particular implement at hand. The example embodiment(s) are, accordingly, to be regarded in an illustrative rather than a restrictive sense.