E NABLING F ACTORIZED P IANO M USIC M ODELING AND G ENERATION WITH THE MAESTRO D ATASET
A BSTRACT
Generating musical audio directly with neural networks is notoriously difﬁcult because it requires coherently modeling structure at many different timescales.
Fortunately, most music is also highly structured and can be represented as discrete note events played on musical instruments.
Herein, we show that by using notes as an intermediate representation, we can train a suite of models capable of transcribing, composing, and synthesizing audio waveforms with coherent musical structure on timescales spanning six orders of magnitude ( ∼ 0.1 ms to ∼s), a process we call Wave2Midi2Wave.
This large advance in the state of the art is enabled by our release of the new MAESTRO (MIDI and Audio Edited for Synchronous TRacks and Organization) dataset, composed of over 172 hours of virtuosic piano performances captured with ﬁne alignment ( ≈ 3 ms) between note labels and audio waveforms.
The networks and the dataset together present a promising approach toward creating new expressive and interpretable neural models of music.
I NTRODUCTION
Since the beginning of the recent wave of deep learning research, there have been many attempts to create generative models of expressive musical audio de novo .
These models would ideally generate audio that is both musically and sonically realistic to the point of being indistinguishable to a listener from music composed and performed by humans.
However, modeling music has proven extremely difﬁcult due to dependencies across the wide range of timescales that give rise to the characteristics of pitch and timbre (short-term) as well as those of rhythm (medium-term) and song structure (long-term).
On the other hand, much of music has a large hierarchy of discrete structure embedded in its generative process: a composer creates songs, sections, and notes, and a performer realizes those notes with discrete events on their instrument, creating sound.
The division between notes and sound is in many ways analogous to the division between symbolic language and utterances in speech.
The WaveNet model by van den Oord et al.
(2016) may be the ﬁrst breakthrough in generating musical audio directly with a neural network.
Using an autoregressive architecture, the authors trained a model on audio from piano performances that could then generate new piano audio sample-bysample.
However, as opposed to their highly convincing speech examples, which were conditioned on linguistic features, the authors lacked a conditioning signal for their piano model.
The result was audio that sounded very realistic at very short time scales (1 or 2 seconds), but that veered off into chaos beyond that.
Dieleman et al.
(2018) made great strides towards providing longer term structure to WaveNet synthesis by implicitly modeling the discrete musical structure described above.
This was achieved by training a hierarchy of VQ-VAE models at multiple time-scales, ending with a WaveNet decoder to generate piano audio as waveforms.
While the results are impressive in their ability to capture long-term structure directly from audio waveforms, the resulting sound suffers from various artifacts at the ﬁne-scale not present in the unconditional WaveNet, clearly distinguishing it from real musical audio.
Also, while the model learns a version of discrete structure from the audio, it is not
T r an sc r i p t i o n : on s e t s & f r a m e s ( s e c t i on 4 )
[Figure]
directly reﬂective of the underlying generative process and thus not interpretable or manipulable by a musician or user.
Manzelli et al.
(2018) propose a model that uses a WaveNet to generate solo cello music conditioned on MIDI notation.
This overcomes the inability to manipulate the generated sequence.
However, their model requires a large training corpus of labeled audio because they do not train a transcription model, and it is limited to monophonic sequences.
In this work, we seek to explicitly factorize the problem informed by our prior understanding of the generative process of performer and instrument:
[Equation]
which can be thought of as a generative model with a discrete latent code of musical notes.
Since the latent representation is discrete, and the scale of the problem is too large to jointly train, we split the model into three separately trained modules that are each state-of-the-art in their respective domains:
Encoder, P ( notes | audio ) :
An Onsets and Frames (Hawthorne et al., 2018) transcription model to produce a symbolic representation (MIDI) from raw audio.
Prior, P ( notes ) :
A self-attention-based music language model (Huang et al., 2018) to generate new performances in MIDI format based on those transcribed in (1).
Decoder, P ( audio | notes ) :
A WaveNet (van den Oord et al., 2016) synthesis model to generate audio of the performances conditioned on MIDI generated in (2).
We call this process Wave2Midi2Wave.
One hindrance to training such a stack of models is the lack of large-scale annotated datasets like those that exist for images.
We overcome this barrier by curating and publicly releasing alongside this work a piano performance dataset containing well-aligned audio and symbolic performances an order of magnitude larger than the previous benchmarks.
In addition to the high quality of the samples our method produces (see https://goo.gl/ magenta/maestro-examples ), training a suite of models according to the natural musician/instrument division has a number of other advantages.
First, the intermediate representation used is more suitable for human interpretation and manipulation.
Similarly, factorizing the model in this way provides better modularity: it is easy to independently swap out different performance and instrument models.
Using an explicit performance representation with modern language models also allows us to model structure at much larger time scales, up to a minute or so of music.
Finally, we can take advantage of the large amount of prior work in the areas of symbolic music generation and conditional audio generation.
And by using a state-of-the-art music transcription model, we can make use of the same wealth of unlabeled audio recordings previously only usable for training end-to-end models by transcribing unlabeled audio recordings and feeding them into the rest of our model.
C ONTRIBUTIONS OF THIS PAPER
Our contributions are as follows:
We combine a transcription model, a language model, and a MIDI-conditioned WaveNet model to produce a factorized approach to musical audio modeling capable of generating about one minute of coherent piano music.
We provide a new dataset of piano performance recordings and aligned MIDI, an order of magnitude larger than previous datasets.
Using an existing transcription model architecture trained on our new dataset, we achieve state-of-the-art results on a piano transcription benchmark.
D ATASET
We partnered with organizers of the International Piano-e-Competition 1 for the raw data used in this dataset.
During each installment of the competition, virtuoso pianists perform on Yamaha
Disklaviers which, in addition to being concert-quality acoustic grand pianos, utilize an integrated high-precision MIDI capture and playback system.
Recorded MIDI data is of sufﬁcient ﬁdelity to allow the audition stage of the competition to be judged remotely by listening to contestant performances reproduced over the wire on another Disklavier instrument.
The dataset introduced in this paper, which we name MAESTRO (“MIDI and Audio Edited for
Synchronous TRacks and Organization”), contains over a week of paired audio and MIDI recordings from nine years of International Piano-e-Competition events.
The MIDI data includes key strike velocities and sustain pedal positions.
Audio and MIDI ﬁles are aligned with ≈ 3 ms accuracy and sliced to individual musical pieces, which are annotated with composer, title, and year of performance.
Uncompressed audio is of CD quality or higher (44.1–48 kHz 16-bit PCM stereo).
A train/validation/test split conﬁguration is also proposed, so that the same composition, even if performed by multiple contestants, does not appear in multiple subsets.
Repertoire is mostly classical, including composers from the 17 th to early 20 th century.
Table 1 contains aggregate statistics of the
MAESTRO dataset.
[Table]
We make the new dataset (MIDI, audio, metadata, and train/validation/test split conﬁguration) available at https://g.co/magenta/maestro-dataset under a Creative Commons Attribution Non-Commercial Share-Alike 4.0 license.
All results in this paper are from the v1.0.0 release of the dataset.
Several datasets of paired piano audio and MIDI have been published previously and have enabled signiﬁcant advances in automatic piano transcription and related topics.
We are making MAESTRO available because we believe it provides several advantages over existing datasets.
Most signiﬁcantly, as evident from table 2, MAESTRO is around an order of magnitude larger.
Existing datasets also have different properties than MAESTRO that affect model training:
MusicNet (Thickstun et al., 2017) contains recordings of human performances, but separatelysourced scores.
As discussed in Hawthorne et al.
(2018), the alignment between audio and score is not fully accurate.
One advantage of MusicNet is that it contains instruments other than piano (not counted in table 2) and a wider variety of recording environments.
MAPS (Emiya et al., 2010) contains Disklavier recordings and synthesized audio created from
MIDI ﬁles that were originally entered via sequencer.
As such, the “performances” are not as natural as the MAESTRO performances captured from live performances.
In addition, synthesized audio makes up a large fraction of the MAPS dataset.
MAPS also contains syntheses and recordings of individual notes and chords, not counted in table 2.
Saarland Music Data (SMD) (M¨uller et al., 2011) is similar to MAESTRO in that it contains recordings and aligned MIDI of human performances on a Disklavier, but is 30 times smaller.
[Table]
3. A LIGNMENT
Our goal in processing the data from International Piano-e-Competition was to produce pairs of audio and MIDI ﬁles time-aligned to represent the same musical events.
The data we received from the organizers was a combination of MIDI ﬁles recorded by Disklaviers themselves and WAV audio captured with conventional recording equipment.
However, because the recording streams were independent, they differed widely in start times and durations, and they were also subject to jitter.
Due to the large volume of content, we developed an automated process for aligning, slicing, and time-warping provided audio and MIDI to ensure a precise match between the two.
Our approach is based on globally minimizing the distance between CQT frames from the real audio and synthesized MIDI (using FluidSynth 3 ).
Obtaining a highly accurate alignment is non-trivial, and we provide full details in the appendix.
3. D ATASET SPLITTING
For all experiments in this paper, we use a single train/validation/test split designed to satisfy the following criteria:
• No composition should appear in more than one split.
• Train/validation/test should make up roughly 80/10/10 percent of the dataset (in time), respectively.
These proportions should be true globally and also within each composer.
Maintaining these proportions is not always possible because some composers have too few compositions in the dataset.
• The validation and test splits should contain a variety of compositions.
Extremely popular compositions performed by many performers should be placed in the training split.
For comparison with our results, we recommend using the splits which we have provided.
We do not necessarily expect these splits to be suitable for all purposes; future researchers are free to use alternate experimental methodologies.
P IANO T RANSCRIPTION
The large MAESTRO dataset enables training an automatic piano music transcription model that achieves a new state of the art.
We base our model on Onsets and Frames, with several modiﬁcations informed by a coarse hyperparameter search using the validation split.
For full details of the base model architecture and training procedure, refer to Hawthorne et al.
One important modiﬁcation was adding an offset detection head, inspired by Kelz et al.
The offset head feeds into the frame detector but is not directly used during decoding.
The offset labels are deﬁned to be the 32ms following the end of each note.
We also increased the size of the bidirectional LSTM layers from 128 to 256 units, changed the number of ﬁlters in the convolutional layers from 32/32/64 to 48/48/96, and increased the units in the fully connected layer from 512 to 768.
We also stopped gradient propagation into the onset subnetwork from the frame network, disabled weighted frame loss, and switched to HTK frequency spacing (Young et al., 2006) for the mel-frequency spectrogram input.
In general, we found that the best ways to get higher performance with the larger dataset were to make the model larger and simpler.
The ﬁnal important change we made was to start using audio augmentation during training using an approach similar to the one described in McFee et al.
During training, every input sample was modiﬁed using random parameters for the SoX 4 audio tool.
The parameters, ranges, and random sampling methods are described in table 3.
We found that this was particularly important when evaluating on the MAPS dataset, likely because the audio augmentation made the model more robust to differences in recording environment and piano qualities.
The differences in training results are summarized in table 6.
When evaluating against the MAESTRO dataset, we found that audio augmentation made results slightly worse, so results in table 5 are presented without audio augmentation.
[Table]
After training on the MAESTRO training split for 670k steps, we achieved state of the art results described in table 4 for the MAPS dataset.
We also present our results on the train, validation
[Table]
and test splits of the MAESTRO dataset as a new baseline score in table 5.
Note that for calculating the scores of the train split, we use the full duration of the ﬁles without splitting them into 20-second chunks as is done during training
[Table]
In sections 5 and
[Table]
we demonstrate how using this transcription model enables training language and synthesis models on a large set of unlabeled piano data.
To do this, we transcribe the audio in the
MAESTRO training set, although in theory any large set of unlabeled piano music would work.
We
[Table]
[Table]
call this new, transcribed version of the training set MAESTRO-T. While it is true that the audio transcribed for MAESTRO-T was also used to train the transcription model, table 5 shows that the model performance is not signiﬁcantly different between the training split and the validation or test splits, and we needed the larger split to enable training the other models.
M USIC T RANSFORMER T RAINING
For our generative language model, we use the decoder portion of a Transformer (Vaswani et al., 2017) with relative self-attention, which has previously shown compelling results in generating music with longer-term coherence (Huang et al., 2018).
We trained two models, one on MIDI data from the MAESTRO dataset and another on MIDI transcriptions inferred by Onsets and Frames from audio in MAESTRO, referred to as MAESTRO-T in section 4.
For full details of the model architecture and training procedure, refer to Huang et al.
We used the same training procedure for both datasets.
We trained on random crops ofevents and employed transposition and time compression/stretching data augmentation.
The transpositions were uniformly sampled in the range of a minor third below and above the original piece.
The time stretches were at discrete amounts and uniformly sampled from the set { 0 .
We evaluated both of the models on their respective validation splits.
[Table]
Samples outputs from the Music Transformer model can be heard in the Online Supplement ( https://goo.gl/magenta/maestro-examples ).
P IANO S YNTHESIS
Most commercially available systems that are able to synthesize a MIDI sequence into a piano audio signal are concatenative: they stitch together snippets of audio from a large library of recordings of individual notes.
While this stitching process can be quite ingenious, it does not optimally capture the various interactions between notes, whether they are played simultaneously or in sequence.
An alternative but less popular strategy is to simulate a physical model of the instrument.
Constructing an accurate model constitutes a considerable engineering effort and is a ﬁeld of research by itself (Bank et al., 2010; Valimaki et al., 2012).
WaveNet (van den Oord et al., 2016) is able to synthesize realistic instrument sounds directly in the waveform domain, but it is not as adept at capturing musical structure at timescales of seconds or longer.
However, if we provide a MIDI sequence to a WaveNet model as conditioning information, we eliminate the need for capturing large scale structure, and the model can focus on local structure instead, i.e., instrument timbre and local interactions between notes.
Conditional WaveNets are also used for text-to-speech (TTS), and have been shown to excel at generating realistic speech signals conditioned on linguistic features extracted from textual data.
This indicates that the same setup could work well for music audio synthesis from MIDI sequences.
Our WaveNet model uses a similar autoregressive architecture to van den Oord et al.
(2016), but with a larger receptive ﬁeld: 6 (instead of 3) sequential stacks with 10 residual block layers each.
We found that a deeper context stack, namely 2 stacks with 6 layers each arranged in a series, worked better for this task.
We also updated the model to produce 16-bit output using a mixture of logistics as described in van den Oord et al.
The input to the context stack is an onset “piano roll” representation, a size-88 vector signaling the onset of any keys on the keyboard, with 4ms bins (250Hz).
Each element of the vector is a ﬂoat that represents the strike velocity of a piano key in the 4ms frame, scaled to the range [0, 1].
When there is no onset for a key at a given time, the value is 0.
We initially trained three models:
Unconditioned Trained only with the audio from the combined MAESTRO training/validation splits with no conditioning signal.
Ground Trained with the ground truth audio/MIDI pairs from the combined MAESTRO training/validation splits.
Transcribed Trained with ground truth audio and MIDI inferred from the audio using the Onsets and Frames method, referred to as MAESTRO-T in section 4.
The resulting losses after 1M training steps were 3.72, 3.70 and 3.84, respectively.
Due to teacher forcing, these numbers do not reﬂect the quality of conditioning, so we rely on human judgment for evaluation, which we address in the following section.
It is interesting to note that the WaveNet model recreates non-piano subtleties of the recording, including the response of the room, breathing of the player, and shufﬂing of listeners in their seats.
These results are encouraging and indicate that such methods could also capture the sound of more dynamic instruments (such as string and wind instruments) for which convincing synthesis/sampling methods lag behind piano.
Due to the heterogeneity of the ground truth audio quality in terms of microphone placement, ambient noise, etc., we sometime notice “timbral shifts” during longer outputs from these models.
We therefore additionally trained a model conditioned on a one-hot year vector at each timestep (similar to speaker conditioning in TTS), which succeeds in producing consistent timbres and ambient qualities during long outputs (see Online Supplement).
A side effect of arbitrary windowing of the training data across note boundaries is a sonic crash that often occurs at the beginning of generated outputs.
To sidestep this issue, we simply trim the ﬁrst 2 seconds of all model outputs reported in this paper, and in the Online Supplement ( https: //goo.gl/magenta/maestro-examples ).
L ISTENING T ESTS
Since our ultimate goal is to create realistic musical audio, we carried out a listening study to determine the perceived quality of our method.
To separately assess the effects of transcription, language modeling, and synthesis on the listeners’ responses, we presented users with two 20-second clipsdrawn from the following sets, each relying on an additional model from our factorization:
Ground Truth Recordings Clips randomly selected from the MAESTRO validation audio split.
WaveNet Unconditioned Clips generated by the Unconditioned WaveNet model described in section 6.
WaveNet Ground/Test Clips generated by the Ground WaveNet model described in section 6, conditioned on random 20-second MIDI subsequences from the MAESTRO test split.
WaveNet Transcribed/Test Clips generated by the Transcribed WaveNet model described in section 6, conditioned on random 20-second subsequences from the MAESTRO test split.
WaveNet Transcribed/Transformer Clips generated by the Transcribed WaveNet model described in section 6, conditioned on random 20-second subsequences from the Music Transformer model described in section 5 that was trained on MAESTRO-T.
The ﬁnal set of samples demonstrates the full end-to-end ability of taking unlabeled piano performances, inferring MIDI labels via transcription, generating new performances with a language model trained on the inferred MIDI, and rendering new audio as though it were played on a similar piano—all without any information other than raw audio recordings of piano performances.
Participants were asked which clip they thought sounded more like a recording of somebody playing a musical piece on a real piano, on a Likert scale.
640 ratings were collected, with each source involved in 128 pair-wise comparisons.
Figure 2 shows the number of comparisons in which performances from each source were selected as more realistic.
[Figure]
A Kruskal-Wallis H test of the ratings showed that there is at least one statistically signiﬁcant difference between the models: χ 2 (2)
= 67 .
63 , p < 0 .
A post-hoc analysis using the Wilcoxon signed-rank test with Bonferroni correction showed that there was not a statistically signiﬁcant difference in participant ratings between real recordings and samples from the WaveNet Ground/Test and WaveNet Transcribed/Test models with p > 0 .
Audio of some of the examples used in the listening tests is available in the Online Supplement ( https://goo.gl/magenta/maestro-examples ).
C ONCLUSION
We have demonstrated the Wave2Midi2Wave system of models for factorized piano music modeling, all enabled by the new MAESTRO dataset.
In this paper we have demonstrated all capabilities on the same dataset, but thanks to the new state-of-the-art piano transcription capabilities, any large set of piano recordings could be used, 6 which we plan to do in future work.
After transcribing the recordings, the transcriptions could be used to train a WaveNet and a Music Transformer model, and then new compositions could be generated with the Transformer and rendered with the WaveNet.
These new compositions would have similar musical characteristics to the music in the original dataset, and the audio renderings would have similar acoustical characteristics to the source piano.
The most promising future work would be to extend this approach to other instruments or even multiple simultaneous instruments.
Finding a suitable training dataset and achieving sufﬁcient transcription performance will likely be the limiting factors.
The new dataset (MIDI, audio, metadata, and train/validation/test split conﬁgurations) is available at https://g.co/magenta/maestro-dataset under a Creative Commons Attribution NonCommercial Share-Alike 4.0 license.
The Online Supplement, including audio examples, is available at https://goo.gl/magenta/maestro-examples . A CKNOWLEDGMENTS
We would like to thank Michael E. Jones and Stella Sick for their help in coordinating the release of the source data and Colin Raffel for his careful review and comments on this paper.
A PPENDIX
In this appendix, we describe in detail how the MAESTRO dataset from section 3 was aligned and segmented.
C OARSE A LIGNMENT
The key idea for the alignment process was that even an untrained human can recognize whether two performances are of the same score based on raw audio, disregarding differences in the instrument or recording equipment used.
Hence, we synthesized the provided MIDI (using FluidSynth with a
SoundFont sampled from recordings of a Disklavier 7 ) and sought to deﬁne an audio-based difference metric that could be minimized to ﬁnd the best-alignment shift for every audio/MIDI pair.
We wanted the metric to take harmonic features into account, so as a ﬁrst step we used librosa (McFee et al., 2017) to compute the Constant-Q Transform (Brown, 1991; Sch¨orkhuber & Klapuri, 2010) of both original and synthesized audio.
For the initial alignment stage we picked a hop length of 4096 samples ( ∼ 90 ms) as a trade-off between speed and accuracy, which proved reasonable for most of the repertoire.
8 Microphone setup varied between competition years and stages, resulting in varying frequency response and overall amplitude levels in recordings, especially in the lower and higher ends of the piano range.
To account for that, we limited the CQT to 48 buckets aligned with MIDI notes C2–B5, and also converted amplitude levels to dB scale with maximum absolute amplitude as a reference point and a hard cut-off at -80 dB. Original and synthesized audio also differed in sound decay rate, so we normalized the resulting CQT arrays time-wise by dividing each hop column by its minimum value (averaged over a 5-hop window).
A single MIDI ﬁle from a Disklavier typically covered several hours of material corresponding to a sequence of shorter audio ﬁles from several seconds up to an hour long.
We slid the normalized
CQT of each such original audio ﬁle against a window of synthesized MIDI CQT of the same length and used mean squared error (MSE) between the two as the difference metric.
9 Minimum error determined best alignment, after which we attempted to align the next audio ﬁle in sequence with the remaining length of the corresponding MIDI ﬁle.
Due to the length of MIDI ﬁles, it was impractical to calculate MSE at each possible shift, so instead we trimmed silence at the beginning of audio, and attempted to align it with the ﬁrst “note on” event of the MIDI ﬁle, within ± 12 minutes tolerance.
If the minimum error was still high, we attempted alignment at the next “note on” event after a 30-second silence.
This approach allowed us to skip over unusable sections of MIDI recordings that did not correspond to audio, e.g., instrument tuning and warm-ups, and also non-musical segments of audio such as applause and announcements.
Non-piano sounds also considerably increased the MSE metric for very short audio ﬁles, so we had to either concatenate those with their longer neighbors if they had any musical material or exclude them completely.
Events that were present at the beginning of audio ﬁles beyond the chosen shift tolerance which did not correspond to MIDI had to be cut off manually.
In order to recover all musically useful data we also had to manually repair several MIDI ﬁles where the clock had erroneously jumped, causing the remainder of the ﬁle to be out of sync with corresponding audio.
After tuning process parameters and addressing the misaligned audio/MIDI pairs detected by unusually high CQT MSE, we have reached the state where each competition year (i.e., different audio recording setup) has ﬁnal metric values for all pairs within a close range.
Spot-checking the pairs with the highest MSE values for each year conﬁrmed proper alignment, which allowed us to proceed to the segmentation stage.
S EGMENTATION
Since certain compositions were performed by multiple contestants, 10 we needed to segment the aligned audio/MIDI pairs further into individual musical pieces, so as to enable splitting the data into train, validation, and test sets disjoint on compositions.
While the organizers provided the list of composition metadata for each audio ﬁle, for some competition years timing information was missing.
In such cases we greedily sliced audio/MIDI pairs at the longest silences between MIDI notes up to the expected number of musical pieces.
Where expected piece duration data was available, we applied search with backtracking roughly as follows.
As an invariant, the segmentation algorithm maintained a list of intervals as start–end time offsets along with a list of expected piece durations, so that the total length of the piece durations corresponding to each interval was less than the interval duration (within a certain tolerance).
At each step we picked the next longest MIDI silence and determined which interval it belonged to.
Then we split that interval in two at the silence and attempted to split the corresponding sequence of durations as well, satisfying the invariant.
For each suitable split the algorithm continued to the next longest silence.
If multiple splits were possible, the algorithm preferred the ones that divided the piece durations more evenly according to a heuristic.
If no such split was possible, the algorithm either skipped current silence if it was short 11 and attempted to split at the next one, or backtracked otherwise.
It also backtracked if no more silences longer than 3 seconds were available.
The algorithm succeeded as soon as each interval corresponded to exactly one expected piece duration.
Once a suitable segmentation was found, we sliced each audio/MIDI pair at resulting intervals, additionally trimming short clusters of notes at the beginning or end of each segment that appeared next to long MIDI silences in order to cut off additional non-music events (e.g., tuning or contestants testing the instrument during applause), and adding an extra 1 second of padding at both ends before making the ﬁnal cut.
F INE ALIGNMENT
After the initial alignment and segmentation, we applied Dynamic Time Warping (DTW) to account for any jitter in either the audio or MIDI recordings.
DTW has seen wide use in audio-to-MIDI alignment; for an overview see M¨uller (2015).
We follow the align midi example from pretty midi (Raffel & Ellis, 2014), except that we use a custom C++ DTW implementation for improved speed and memory efﬁciency to allow for aligning long sequences.
First, in Python, we use librosa to load the audio and resample it to a 22,050Hz mono signal.
Next, we load the MIDI and synthesize it at the same sample rate, using the same FluidSynth process as above.
Then, we pad the end of the shorter of the two sample arrays
so they are the same length.
We use the same procedure as align midi to extract CQTs from both sample arrays, except that we use a hop length of 64 to achieve a resolution of ∼ 3ms.
We then pass these CQTs to our C++ DTW implementation.
To avoid calculating the full distance matrix and taking its mean to get a penalty value, we instead sample 100k random pairs and use the mean of their cosine distances.
We use the same DTW algorithm as implemented in librosa except that we calculate cosine distances only within a Sakoe-Chiba band radius (Sakoe & Chiba, 1978) of 2.5 seconds instead of calculating distances for all pairs.
Staying within this small band limits the number of calculations we need to make and the number of distances we have to store in memory.
This is possible because we know from the previous alignment pass that the sequences are already mostly aligned and we just need to account for small constant offsets due to the lower resolution of the previous process and apply small sequence warps to recover from any occasional jitter.