def get_voxel_time_series(self, x, y, z, var=0, mode=0):
"""
Retrieve for a given voxel (x,y,z) the entire timeline.
:param x: int coordinate
:param y: int coordinate
:param z: int coordinate
:return: A complex dictionary with information about current voxel.
The main part will be a vector with all the values over time from the x,y,z coordinates.
"""
if self.region_mapping_volume is None:
raise exceptions.TVBException("Invalid method called for TS without Volume Mapping!")
volume_rm = self.region_mapping_volume
x, y, z = preprocess_space_parameters(x, y, z, volume_rm.length_1d, volume_rm.length_2d, volume_rm.length_3d)
idx_slices = slice(x, x + 1), slice(y, y + 1), slice(z, z + 1)
idx = int(volume_rm.get_data('array_data', idx_slices))
time_length = self.read_data_shape()[0]
var, mode = int(var), int(mode)
voxel_slices = prepare_time_slice(time_length), slice(var, var + 1), slice(idx, idx + 1), slice(mode, mode + 1)
label = volume_rm.connectivity.region_labels[idx]
background, back_min, back_max = None, None, None
if idx < 0:
back_min, back_max = self.get_min_max_values()
background = numpy.ones((time_length, 1)) * self.out_of_range(back_min)
label = 'background'
result = postprocess_voxel_ts(self, voxel_slices, background, back_min, back_max, label)
return result
def get_voxel_time_series(self, x, y, z, **kwargs):
"""
Retrieve for a given voxel (x,y,z) the entire timeline.
:param x: int coordinate
:param y: int coordinate
:param z: int coordinate
:return: A complex dictionary with information about current voxel.
The main part will be a vector with all the values over time from the x,y,z coordinates.
"""
overall_shape = self.read_data_shape()
x, y, z = preprocess_space_parameters(x, y, z, overall_shape[1], overall_shape[2], overall_shape[3])
slices = prepare_time_slice(overall_shape[0]), slice(x, x + 1), slice(y, y + 1), slice(z, z + 1)
result = postprocess_voxel_ts(self, slices)
return result