def plot_individual_derivative_heatmaps(inputfile, outputfolder, group='Group', cutoff=None, id_col='TrackID'):
"""
Plot the derivatives of spectral components of different groups over time as a heatmap.
Parameters
----------
inputfile : str
Path to the file with spectral data.
outputfolder : str
Directory to save the plotted heat maps.
group : str, optional
Column in the input data sheet to use for grouping.
Default is 'Group'.
cutoff : int, optional
The number of degrees to display.
If None, all degrees will be displayed.
id_col : str
Column in the input data sheet to group connected time points.
Default is 'TrackID'
"""
filelib.make_folders([os.path.dirname(outputfolder)])
stat = pd.read_csv(inputfile, sep='\t', index_col=0)
stat['Time'] = np.int_(np.round_(stat['Time'] / 10.)) * 10
if cutoff is not None:
stat = stat[stat['degree'] <= cutoff]
for gr in stat[group].unique():
substat = stat[stat[group] == gr]
for id in substat[id_col].unique():
subsubstat = substat[substat[id_col] == id]
time_spectrum = TimeSpectrum()
for t in subsubstat['Time'].unique():
sp = Spectrum()
sp.harmonics_csv = subsubstat[subsubstat['Time'] == t]
time_spectrum.add_spectrum(sp, timepoint=t)
time_spectrum.compute_derivative()
pl = time_spectrum.derivative_heatmap()
if pl is not None:
pl.savefig(outputfolder + '_' + gr + '_' + 'track_' + str(id) + '.png')
def test_add_spectrum_and_plotting(self):
tsp = TimeSpectrum()
sp = Spectrum()
sp.from_surface(surface=np.ones([10, 10]))
sp.convert_to_csv()
for i in range(3):
tsp.add_spectrum(sp, timepoint=i*20)
sp = Spectrum()
surf = np.ones([10, 10])
surf[3:4, 4:8] = 10
sp.from_surface(surface=surf)
sp.convert_to_csv()
for i in range(2):
tsp.add_spectrum(sp, timepoint=i*20+60)
sp = Spectrum()
sp.from_surface(surface=np.ones([10, 10]))
sp.convert_to_csv()
for i in range(3):
tsp.add_spectrum(sp, timepoint=i*20+100)
self.assertEqual(len(tsp.spectra), 8)
tsp.save_to_csv('data/test_data/time_spectrum.csv')
self.assertEqual(os.path.exists('data/test_data/time_spectrum.csv'), True)
pl = tsp.time_heatmap()
pl.savefig('data/test_data/time_heatmap.png')
self.assertEqual(os.path.exists('data/test_data/time_heatmap.png'), True)
tsp.compute_derivative()
pl = tsp.derivative_heatmap()
pl.savefig('data/test_data/derivative_heatmap.png')
self.assertEqual(os.path.exists('data/test_data/derivative_heatmap.png'), True)
pl = tsp.plot_mean_abs_derivative()
pl.savefig('data/test_data/mean_abs_derivative.png')
self.assertEqual(os.path.exists('data/test_data/mean_abs_derivative.png'), True)
shutil.rmtree('data/test_data/')