8.4.3. Моделювання поверхонь та їх оцінка

Для моделювання безперервних поверхонь на основі дискретного масиву даних використовується процедура локальної інтерполяції, аналогічна до тієї, що застосовується при інтерполяції методом зваженого ковзного усереднювання, відповідно до якої розрахункове значення змінної z у деякій точці простору х₀ задається виразом:

Для визначення вагових коефіцієнтів 1., що забезпечують мінімум похибки при заданому масиві просторово-координованих даних, використовується оптимізована варіограмна модель. Процедура визначення вагових коефіцієнтів лінійної моделі (8.13) базується на теорії випадкових процесів, виходячи з якої дисперсія оцінки змінної z(x) може бути записана як функція значень напівдисперсії між всіма парами проб (вимірювань), а також між всіма пробами (вимірюваннями) і оцінюваною точкою (д:₀), та значень вагових коефіцієнтів:

Завдання оптимальної інтерполяції, таким чином, полягає в знаходженні такого набору вагових коефіцієнтів X., який би забезпечував максимальну точність оцінки, тобто мінімальну дисперсію а_е. Отже, постає завдання мінімізації функції дисперсії, розв'язком якої є ті вагові коефіцієнти, які цей мінімум забезпечують.

Відомо, що будь-який екстремум функції багатьох змінних супроводжується рівністю нулю всіх часткових похідних у точці екстремуму. У нашому випадку всі часткові похідні є лінійними функціями, і пошук екстремуму зводиться до розв'язання системи лінійних рівнянь. Позитивна напіввизначеність функції варіограми забезпечує, що розв'язання системи існуватиме, буде єдиним і відповідатиме саме мінімуму дисперсії, а не максимуму (Мальцев, 1993).

Для забезпечення однієї з головних вимог завдання оцінювання — вимоги незміщеності оцінки — у систему (8.14) необхідно ввести додаткове рівняння, що визначає умову рівності одиниці суми всіх вагових коефіцієнтів, або, що те ж саме, додати відповідний доданок у рівняння функції, що мінімізується:

Обчислюючи і прирівнюючи до нуля часткові похідні, одержуємо систему лінійних рівнянь:

Розв'язком системи (8.16) є і шукані вагові коефіцієнти, і значення множника Лагранжа, які дозволяють окрім, власне оцінки змінної z, у будь-якій точці простору або комірці растра визначити значення кригінгової дисперсії.

Для знаходження значень змінної в тих точках простору, де вимірювання не проводилися, використовується модель (8.13) зі знайденими ваговими коефіцієнтами. При використанні растрової моделі просторових даних оцінка (прогноз) проводиться для всіх комірок растра з невідомими значеннями змінної. У комірках, де значення змінної відомі, ці значення беруться як оцінні. У результаті будується (моделюється) безперервна поверхня z(x), що задовольняє сформульовані вище вимоги — мінімуму похибки і незміщенності.

Дисперсія відхилень оцінного (прогнозного) значення змінної від істинного, тобто похибка оцінювання (прогнозу), для кожної точки простору (комірки растра) обчислюється за формулою (8.14).

Описаний метод просторової інтерполяції відомий як звичайний лінійний (або ординарний) точковий кригінг.