Зараз я працюю над оцінками біомаси за допомогою супутникових знімків. Я швидко визначу передумови свого запитання, а потім поясню статистичне питання, над яким я працюю.
Фон
Проблема
Я намагаюся оцінити біомасу на території Франції. Моя відповідь - об'ємна щільність парового дерева (в ), яка більш-менш пропорційна біомасі (залежно від щільності деревини ...).
Незалежні мені змінні - це індекси вегетації, отримані від вимірюваних коефіцієнтів відбиття на цій ділянці (супутник, який використовується у дослідженні, є MODIS для тих, хто це знає). Такі індекси, наприклад, NDVI, EVI тощо. У мене є карти індексів, а роздільна здатність карт - 250м.
Існують сильні кореляції між цими показниками та обсягом в одному типі лісу (біома та клімат). Тому я намагаюся регресувати густину гучності щодо цих показників (фактично їх часових рядів) на інвентарних сюжетах, де я знаю обсяг.
Лісові запаси
Обсяг на цих ділянках оцінюється наступним методом вибірки:
- Вузли інвентаризації розміщуються на звичайній сітці, що охоплює область.
- До кожного вузла додається ділянка, і на цій ділянці відбувається процес інвентаризації (типи дерев, обсяги, висота навісу тощо). Звичайно, мене цікавить лише інвентаризаційна ділянка, і значення моїх рослинних індексів - це значення пікселя, що містить ділянку.
Процес інвентаризації на ділянці такий:
- Виміряйте дерева, які мають діаметр> 37,5 см у колі радіуса 15 м
- Виміряйте дерева, діаметр яких> 22,5 см у радіусі кола 9м
- Виміряйте дерева, які мають діаметр> 7,5 см у колі радіуса 6м
Потім об'ємна щільність обчислюється за допомогою коефіцієнтів розширення.
Для кожної ділянки я маю доступ до даних для всіх виміряних дерев.
Більше того, для кожного окремого дерева я маю невизначеність щодо обсягу через використання алометричних рівнянь (скажімо 10%).
Там, де важлива статистика ...
Щоб мої регресії були більш точними, мені потрібно для кожної оцінки обсягу дисперсія / CI цього заходу. Це залежить від ІМО від кількості вибірених дерев та знайденої густоти.
Тож у мене є дві проблеми:
Як врахувати той факт, що мої показники вегетації вимірюються в пікселі 250м?
Я можу припустити, що густина гучності є постійною понад один піксель, і що я відбираю цей піксель за допомогою однієї графіки інвентаризації.
Як оцінити мінливість моєї густини густини?
Я думаю, що я міг би використовувати завантажувальну групу на популяції дерев. Але моя загальна кількість виміряних дерев може бути досить маленькою (від 7 до 20 ...). Крім того, як я можу взяти до уваги той факт, що я вимірюю дерева на різних колах залежно від їх розмірів? І як має змінюватися мінливість, якщо я дивлюся на цілий піксель?
Я також думав, що можу використовувати моделювання Монте-Карло для імітації лісу, а потім випадковим чином відібрати цей ліс з ділянками, щоб побачити, що відбувається ...
У мене немає сильного статистичного досвіду, тому я трохи загублений!