def test_pitch_raw_using_creaper(self): # Test against previously generated data to make sure nothing has # broken and that there are no cross platform or REAPER version issues. # Data was generated on Manjaro Linux. for fn in wav_fns: # Estimate F0 using REAPER, use default VoiceSauce values t_raw, F0_raw = creaper_pitch(fn, reaper_path=default_reaper_path, frame_shift=1, max_pitch=500, min_pitch=40, high_pass=True, hilbert_transform=False, inter_mark=10) # Get sample time data time_sample = get_sample_data(fn, 'reaper', 'rtF0', '1ms') # Check number of entries is consistent self.assertEqual(len(t_raw), len(time_sample)) # Check that computed time and sample_data are "close enough" for # floating precision self.assertAllClose(t_raw, time_sample, rtol=1e-05, atol=1e-08) # Get sample F0 data F0_sample = get_sample_data(fn, 'reaper', 'reaperF0', '1ms') # Check number of entries is consistent self.assertEqual(len(F0_raw), len(F0_sample)) # Check that computed F0 and sample_data are "close enough" for # floating precision self.assertAllClose(F0_raw, F0_sample, rtol=1e-05, atol=1e-08, equal_nan=True)
def test_pitch_raw(self): # Test against previously generated data to make sure nothing has # broken and that there are no cross platform or Praat version issues. # Data was generated by Praat v6.0.29 on Manjaro Linux. for fn in wav_fns: # Estimate raw Praat F0 t_raw, F0_raw = praat_raw_pitch(fn, default_praat_path, frame_shift=1, method='cc', min_pitch=40, max_pitch=500, silence_threshold=0.03, voice_threshold=0.45, octave_cost=0.01, octave_jumpcost=0.35, voiced_unvoiced_cost=0.14, kill_octave_jumps=False, interpolate=False, smooth=False, smooth_bandwidth=5) # Get sample time data sample_data = get_sample_data(fn, 'praat', 'ptF0', '1ms') # Check number of entries is consistent self.assertEqual(len(t_raw), len(sample_data)) # Check that computed time and sample_data are "close enough" for # floating precision self.assertAllClose(t_raw, sample_data, rtol=1e-05, atol=1e-08) # Get sample F0 data sample_data = get_sample_data(fn, 'praat', 'pF0', '1ms') # Check number of entries is consistent self.assertEqual(len(F0_raw), len(sample_data)) # Check that computed F0 and sample_data are "close enough" for # floating precision self.assertTrue((np.isclose(F0_raw, sample_data, rtol=1e-05, atol=1e-08) | (np.isnan(F0_raw) & np.isnan(sample_data))).all())
def test_pitch_raw(self): # Test against previously generated data to make sure nothing has # broken and that there are no cross platform or snack version issues. # Data was generated by snack 2.2.10 on Manjaro Linux. for fn in wav_fns: F0_raw, V_raw = snack_raw_pitch(fn, snack_method, frame_shift=1, window_size=25, max_pitch=500, min_pitch=40, tcl_shell_cmd=tcl_cmd) # Check V data # Voice is 0 or 1, so (hopefully) no FP rounding issues. sample_data = get_sample_data(fn, 'snack', 'sV', '1ms') # Check that all voice data is either 0 or 1 self.assertTrue(np.all((V_raw == 1) | (V_raw == 0))) self.assertTrue(np.all((sample_data == 1) | (sample_data == 0))) # Check number of entries is consistent self.assertEqual(len(V_raw), len(sample_data)) # Check actual data values are "close enough", # within floating precision self.assertAllClose(V_raw, sample_data, rtol=1e-05, atol=1e-08) # Check F0 data sample_data = get_sample_data(fn, 'snack', 'sF0', '1ms') # Check number of entries is consistent self.assertEqual(len(F0_raw), len(sample_data)) # Check that F0 and sample_data are "close enough" for # floating precision if not np.allclose(F0_raw, sample_data, rtol=1e-05, atol=1e-08): # If first check fails, try lowering relative tolerance and # redoing the check idx = np.where(np.isclose(F0_raw, sample_data) == False)[0] print('\nChecking F0 data using rtol=1e-05, atol=1e-08 in {}:'. format(fn)) print( 'Out of {} array entries in F0 snack data, discrepancies in these indices' .format(len(F0_raw))) for i in idx: print('idx {}, OpenSauce F0 = {}, sample F0 = {}'.format( i, F0_raw[i], sample_data[i])) print( 'Reducing relative tolerance to rtol=3e-05 and redoing check:' ) self.assertAllClose(F0_raw, sample_data, rtol=3e-05, atol=1e-08) print('OK') else: self.assertAllClose(F0_raw, sample_data, rtol=1e-05, atol=1e-08)
def test_pitch_using_pyreaper(self): # Test against previously generated data to make sure nothing has # broken and that there are no cross platform or REAPER version issues. # Data was generated on Manjaro Linux. for fn in wav_fns: # Frame shift f_len = 1 # Need ns (number of samples) and sampling rate (Fs) from wav file # to compute data length sound_file = SoundFile(fn) data_len = np.int_(np.floor(sound_file.ns / sound_file.fs / f_len * 1000)) # Estimate F0 using REAPER, use default VoiceSauce values F0_os = reaper_pitch(sound_file, data_len, use_pyreaper=True, frame_shift=f_len, max_pitch=500, min_pitch=40, high_pass=True, hilbert_transform=False, inter_mark=10) # Get sample F0 data F0_sample_from_data = get_sample_data(fn, 'reaper', 'reaperF0', '1ms') # Fill full length F0 vector F0_sample = np.full(data_len, np.nan) F0_sample[:len(F0_sample_from_data)] = F0_sample_from_data # Check that F0 values are "close" self.assertAllClose(F0_os, F0_sample, rtol=1e-05, atol=1e-08, equal_nan=True)
def test_formants_raw(self): # Test against previously generated data to make sure nothing has # broken and that there are no cross platform or Praat version issues. # Data was generated by Praat v6.0.29 on Manjaro Linux. for fn in wav_fns: # Get raw Praat formant estimates estimates_raw = praat_raw_formants(fn, default_praat_path, frame_shift=1, window_size=25, num_formants=4, max_formant_freq=6000) for n in self.formants4_names: # Get sample data corresponding to key name sample_data = get_sample_data(fn, 'praat', n, '1ms') # Check number of entries is consistent self.assertEqual(len(estimates_raw[n]), len(sample_data)) # Check that our estimates and sample_data are "close enough" # for floating precision self.assertTrue((np.isclose( estimates_raw[n], sample_data, rtol=1e-05, atol=1e-08) | (np.isnan(estimates_raw[n]) & np.isnan(sample_data))).all())
def test_formants_raw(self): # Test against previously generated data to make sure nothing has # broken and that there are no cross platform or snack version issues. # Data was generated by snack 2.2.10 on Manjaro Linux. for fn in wav_fns: estimates_raw = snack_raw_formants(fn, snack_method, frame_shift=1, window_size=25, pre_emphasis=0.96, lpc_order=12, tcl_shell_cmd=tcl_cmd) # Get sample data sample_data = {} for n in sformant_names: sample_data[n] = get_sample_data(fn, 'snack', n, '1ms') # Check number of entries is consistent for n in sformant_names: self.assertEqual(len(estimates_raw[n]), len(sample_data[n])) # Check that estimates and sample_data are "close enough" for # floating precision for n in sformant_names: # Increase rtol from 1e-5 to 1e-3 to account for random seed # used in Snack formants self.assertAllClose(estimates_raw[n], sample_data[n], rtol=1e-03, atol=1e-08)
def test_pitch_raw(self): # Test against previously generated data to make sure nothing has # broken and that there are no cross platform or Praat version issues. # Data was generated by Praat v6.0.29 on Manjaro Linux. for fn in wav_fns: # Estimate raw Praat F0 t_raw, F0_raw = praat_raw_pitch(fn, default_praat_path, frame_shift=1, method='cc', min_pitch=40, max_pitch=500, silence_threshold=0.03, voice_threshold=0.45, octave_cost=0.01, octave_jumpcost=0.35, voiced_unvoiced_cost=0.14, kill_octave_jumps=False, interpolate=False, smooth=False, smooth_bandwidth=5) # Get sample time data sample_data = get_sample_data(fn, 'praat', 'ptF0', '1ms') # Check number of entries is consistent self.assertEqual(len(t_raw), len(sample_data)) # Check that computed time and sample_data are "close enough" for # floating precision self.assertAllClose(t_raw, sample_data, rtol=1e-05, atol=1e-08) # Get sample F0 data sample_data = get_sample_data(fn, 'praat', 'pF0', '1ms') # Check number of entries is consistent self.assertEqual(len(F0_raw), len(sample_data)) # Check that computed F0 and sample_data are "close enough" for # floating precision self.assertTrue( (np.isclose(F0_raw, sample_data, rtol=1e-05, atol=1e-08) | (np.isnan(F0_raw) & np.isnan(sample_data))).all())
def test_pitch_raw(self): # Test against previously generated data to make sure nothing has # broken and that there are no cross platform or snack version issues. # Data was generated by snack 2.2.10 on Manjaro Linux. for fn in wav_fns: F0_raw, V_raw = snack_raw_pitch(fn, snack_method, frame_shift=1, window_size=25, max_pitch=500, min_pitch=40, tcl_shell_cmd=tcl_cmd) # Check V data # Voice is 0 or 1, so (hopefully) no FP rounding issues. sample_data = get_sample_data(fn, 'snack', 'sV', '1ms') # Check that all voice data is either 0 or 1 self.assertTrue(np.all((V_raw == 1) | (V_raw == 0))) self.assertTrue(np.all((sample_data == 1) | (sample_data == 0))) # Check number of entries is consistent self.assertEqual(len(V_raw), len(sample_data)) # Check actual data values are "close enough", # within floating precision self.assertAllClose(V_raw, sample_data, rtol=1e-05, atol=1e-08) # Check F0 data sample_data = get_sample_data(fn, 'snack', 'sF0', '1ms') # Check number of entries is consistent self.assertEqual(len(F0_raw), len(sample_data)) # Check that F0 and sample_data are "close enough" for # floating precision if not np.allclose(F0_raw, sample_data, rtol=1e-05, atol=1e-08): # If first check fails, try lowering relative tolerance and # redoing the check idx = np.where(np.isclose(F0_raw, sample_data) == False)[0] print('\nChecking F0 data using rtol=1e-05, atol=1e-08 in {}:'.format(fn)) print('Out of {} array entries in F0 snack data, discrepancies in these indices'.format(len(F0_raw))) for i in idx: print('idx {}, OpenSauce F0 = {}, sample F0 = {}'.format(i, F0_raw[i], sample_data[i])) print('Reducing relative tolerance to rtol=3e-05 and redoing check:') self.assertAllClose(F0_raw, sample_data, rtol=3e-05, atol=1e-08) print('OK') else: self.assertAllClose(F0_raw, sample_data, rtol=1e-05, atol=1e-08)
def test_formants_raw(self): # Test against previously generated data to make sure nothing has # broken and that there are no cross platform or Praat version issues. # Data was generated by Praat v6.0.29 on Manjaro Linux. for fn in wav_fns: # Get raw Praat formant estimates estimates_raw = praat_raw_formants(fn, default_praat_path, frame_shift=1, window_size=25, num_formants=4, max_formant_freq=6000) for n in self.formants4_names: # Get sample data corresponding to key name sample_data = get_sample_data(fn, 'praat', n, '1ms') # Check number of entries is consistent self.assertEqual(len(estimates_raw[n]), len(sample_data)) # Check that our estimates and sample_data are "close enough" # for floating precision self.assertTrue((np.isclose(estimates_raw[n], sample_data, rtol=1e-05, atol=1e-08) | (np.isnan(estimates_raw[n]) & np.isnan(sample_data))).all())
def test_pitch_using_creaper(self): # Test against previously generated data to make sure nothing has # broken and that there are no cross platform or REAPER version issues. # Data was generated on Manjaro Linux. for fn in wav_fns: # Frame shift f_len = 1 # Need ns (number of samples) and sampling rate (Fs) from wav file # to compute data length sound_file = SoundFile(fn) data_len = np.int_( np.floor(sound_file.ns / sound_file.fs / f_len * 1000)) # Estimate F0 using REAPER, use default VoiceSauce values F0_os = reaper_pitch(sound_file, data_len, use_pyreaper=False, reaper_path=default_reaper_path, frame_shift=f_len, max_pitch=500, min_pitch=40, high_pass=True, hilbert_transform=False, inter_mark=10) # Get sample F0 data F0_sample_from_data = get_sample_data(fn, 'reaper', 'reaperF0', '1ms') # Fill full length F0 vector F0_sample = np.full(data_len, np.nan) F0_sample[:len(F0_sample_from_data)] = F0_sample_from_data # Check that F0 values are "close" self.assertAllClose(F0_os, F0_sample, rtol=1e-05, atol=1e-08, equal_nan=True)