ПРИМЕЧАНИЕ. Сообщение выглядит длиннее, чем должно быть, из-за строк документации и массива, состоящего из 40 дат и времени.
У меня есть данные временного ряда. Для примера предположим, что у меня есть три параметра, каждый из которых состоит из 40 точек данных: даты и времени (данные dts), скорость (данные vobs) и истекший час (данные els), которые объединены ключом в словарь data_dict .
dts = np.array(['2006/01/01 02:30:04', '2006/01/01 03:30:04', '2006/01/01 03:54:04'
,'2006/01/01 05:30:04', '2006/01/01 06:30:04', '2006/01/01 07:30:04'
,'2006/01/01 08:30:04', '2006/01/01 09:30:04', '2006/01/01 10:30:04'
,'2006/01/01 11:30:04', '2006/01/01 12:30:04', '2006/01/01 13:30:04'
,'2006/01/01 14:30:04', '2006/01/01 15:30:04', '2006/01/01 16:30:04'
,'2006/01/01 17:30:04', '2006/01/01 18:30:04', '2006/01/01 19:30:04'
,'2006/01/01 20:30:04', '2006/01/01 21:30:04', '2006/01/01 21:54:05'
,'2006/01/01 23:30:04', '2006/01/02 00:30:04', '2006/01/02 01:30:04'
,'2006/01/02 02:30:04', '2006/01/02 03:30:04', '2006/01/02 04:30:04'
,'2006/01/02 05:30:04', '2006/01/02 06:30:04', '2006/01/02 07:30:04'
,'2006/01/02 08:30:04', '2006/01/02 09:30:04', '2006/01/02 10:30:04'
,'2006/01/02 11:30:04', '2006/01/02 12:30:04', '2006/01/02 13:30:04'
,'2006/01/02 14:30:04', '2006/01/02 15:30:04', '2006/01/02 16:30:04'
,'2006/01/02 17:30:04'])
vobs = np.array([158, 1, 496, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1
, 1, 1, 823, 1, 1, 1, 1, 303, 1, 1, 1, 1, 253, 1, 1, 1, 408, 1
, 1, 1, 1, 321])
els = np.array([i for i in range(len(vobs))])
data_dictionary = {'datetime' : dts, 'values' : vobs, 'elapsed' : els}
У меня есть функция, которая принимает словарь в качестве входных данных и выводит одно скалярное значение type <float> или type <int>. Приведенная ниже функция проще, чем мой реальный вариант использования, и приведена для примера.
def get_z(dictionary):
""" This function returns a scalar value. """
return np.sum(dictionary['elapsed'] / dictionary['values'])
Я хотел бы посмотреть, как изменяется вывод этой функции при изменении временного интервала. Итак, я создал функцию, которая принимает словарь в качестве входных данных и выводит новый словарь, значения массива которого нарезаются по входным индексам для каждого из ключей во входном словаре. Обратите внимание, что последовательные прошедшие часы могут служить индексами.
def subsect(dictionary, indices):
""" This function returns a dictionary, the array values
of which are sliced at the input indices. """
return {key : dictionary[key][indices] for key in list(dictionary.keys())}
Чтобы убедиться, что вышеуказанные функции работают, можно запустить цикл for, содержащий функцию read_dictionary(...) ниже.
def read_dictionary(dictionary):
""" This function prints the input dictionary as a check. """
for key in list(dictionary.keys()):
print(" .. KEY = {}\n{}\n".format(key, dictionary[key]))
print("\nORIGINAL DATA DICTIONARY\n")
read_dictionary(data_dictionary)
# for i in range(1, 38):
# mod_dictionary = subsect(data_dictionary, indices=slice(i, 39, 1))
# print("\n{}th MODIFIED DATA DICTIONARY\n".format(i))
# read_dictionary(mod_dictionary)
Моя проблема в том, что я хотел бы контурный сюжет. Ось X будет содержать нижнюю границу интервала даты и времени (первая запись mod_dictionary[i]), а ось Y будет содержать верхнюю границу интервала даты и времени (последняя запись mod_dictioary[i]). Обычно при построении контурного графика имеется массив из (x,y) значений, которые превращаются в сетку (X,Y) через numpy.meshgrid. Поскольку моя фактическая функция (не та, что в примере) не векторизована, я могу использовать X.copy().reshape(-1) и изменить форму обратно, используя (...).reshape(X.shape).
Моя точная проблема заключается в том, что я не знаю, как я могу создать сетку разных параметров, используя один словарь в качестве входных данных для функции, которая выводит одно скалярное значение. Есть ли способ сделать это?
