Какая температура почвы сейчас в москве

Константой отсутствия данных является “9999”.

Признаки качества:

0 – значение достоверно.

1 – значение ошибочно, забраковано методом 1.

2 – значение ошибочно, забраковано методом 2.

3 – значение ошибочно, забраковано методом 3.

4– значение ошибочно, забраковано методом 4.

5– значение сомнительно, по методу 1 (т.к. данный метод контроля применить не удалось).

6 – значение сомнительно, по методу 2.

7– значение сомнительно, по методу 3.

8–значение ошибочно: забраковано более, чем одним методом одновременно.

9 – отсутствие данных.

Методы контроля

Представленный массив данных был проконтролирован с помощью четырех методов статистического контроля. Статистический контроль был необходим, поскольку исходные данные наблюдений имеют некоторые ошибки, вызванные разными причинами.

Как известно, статистические методы дают более точные (хорошо трактуемые, понятные) результаты, если применять их к непрерывному последовательному ряду значений (Это идеальный случай). Однако на практике последовательный ряд наблюдений имеет пропуски значений, т.е. некоторые данные отсутствуют.

Вследствие этого при выполнении контроля, приходилось прибегать к некоторым ограничениям применимости того или иного метода контроля или к ограничениям на выполнение каких-либо внутренних расчетов внутри метода. Это делалось, чтобы избежать некорректного контроля, в случаях, если недостаточно данных для выполнения конкретного метода контроля. Эти ограничения часто приводили к  тому, что те или иные значения остались не проконтролированы каким-то из методов, при этом в массиве имеются пометки в виде флагов (признак качества) рядом с каждым значением, указывающие на это обстоятельство.

Метод 1. «Метод гистограмм»

Метод предназначен для обнаружения грубых ошибок, выходящих за допустимые пределы правильных экстремальных значений. Он основан на анализе гистограмм распределения температуры по градациям.

Общая идея метода:

Гистограмма отражает повторяемость температуры по отдельным градациям температуры. Повторяемость – это число случаев с температурой в градации, деленное на общее число членов ряда. При достаточном количестве данных гистограмма должна быть гладкой на краях, а ее значения на краях слева и справа должны плавно переходить в нули.

Общая идея метода заключается в определении граничных минимальной и максимальной градаций безошибочных значений температуры с целью последующего отсечения значений, которые значительно выходят за пределы этих граничных значений.

 Предполагается, что ошибочные значения превышают экстремальные правильные значения на величину более чем одна градация. Такие ошибочные значения должны выходить за пределы правильных значений температуры и отделяться от основных значений градациями с нулевой повторяемостью в левой и правой частях гистограммы. Значения принимаются ошибочными, если на краях гистограммы они попадают в градации, следующие после градации с нулевой повторяемостью.

Примечание: алгоритм применяется на каждой станции и на каждой глубине отдельно.

Метод 2. «Метод проверки на сигмы»

Общая идея метода:

Идея метода заключается в том, что при отсутствии ошибок в исходных данных среднеквадратическое отклонение (σ) всех исходных данных температуры характеризует наиболее вероятные пределы отклонений исходных данных от среднего значения (от нормы). В связи с этим, метод 2 применяется после метода 1, к данным которые признаны корректными по методу 1 – это позволяется снизить влияние грубых ошибок на вычисление среднеквадратического отклонения (σ). В интервале ±4σ (количество сигм подобрано эмпирически) от среднего значения находятся практически все правильные исходные данные (в предположении, что температура имеет нормальное распределение).

Метод 3. «Проверка на связанность соседних по времени значений»

Значения температуры в соседние дни не могут сильно отличаться из-за инерционности изменения температуры почвы. Это свойство используется для выявления скачкообразных ошибок, при анализе  последовательных дней.

Общая идея метода:

Идея метода заключается в том, чтобы проверить каждое суточное значение температуры на допустимые пределы отклонений от линейного изменения температуры между предыдущим и последующим днем.

Для проверки каждого суточного значения температуры вычисляется соответствующее ему интерполированное значение температуры между предыдущим и последующим днем. Т.е., например, для проверки температуры 2 января вычисляется интерполированное значение ti на этот день по данным за первое T i-1 и третье января T i+1 .

ti = (T i+1 + T i-1)/2,  где 

i – номер дня года

T i+1 –температура в последующий день

T i-1 –температура в предшествующий день

ti – интерполированная температура на проверяемый день

Вычисленное ti  и фактическое Ti значения на 2 января могут отличаться в некоторых пределах. Необходимо определить допустимые пределы величины ΔTi = Ti – ti , при которых Ti можно признать правильным значением. Допустимые пределы определяются по анализу среднеквадратического отклонения σо(i) значений ΔTi на проверяемый день за все годы. Если ΔTi выходит за установленные пределы ±5σо(i), то соответствующее проверяемое  значение температуры признается некорректным (в нашем примере некорректным признается исходное значение на 2 января).

Метод 3 применяется после метода 1, к данным которые признаны корректными по методу 1. Это позволяется снизить влияние грубых ошибок на вычисление среднеквадратического отклонения (σо(i)). Интервал ±5σо подобран эмпирически.

Метод 4. «Проверка соседних значений на ошибку обратного знака»

Метод 4 является дополнительным к первым трем методам. Это очень простой метод контроля, призванный найти единичные ошибки “в знаке” значения.

В исходном массиве очень редко встречаются такие ошибки, когда в ряду положительных температур вдруг попадается такое же по величине число, но с обратным знаком.

Метод 4 применяется, как и метод 1, к исходному массиву, в отличие от методов 2 и 3, которые применяются к массиву, проконтролированному методом 1.

Общая идея метода:

В каждых трех последовательных значениях температуры проверяется соответствие знака второго значения со знаком первого и третьего значения температуры.  Если знак второго значения отличается, то выполнятся сравнение значений, чтобы определить переход ли это через 0 оС (т.е. корректное значение) или же это появление одиночного резкого выброса (т.е. появление ошибки).

Конечно, возможно появление неверного знака и в других различных сочетаниях, например, появление двух последовательных чисел с неверным знаком подряд и т. д. – такие ситуации этот метод не учитывает, но такие ситуации редкие.