Що ми виміряли в лісах Дадьї, а що ні
Для кожної задачі ми порівняли конкуруючі підходи на лейблах з незалежних супутникових продуктів і перевірили їх на даних, яких вони не бачили. Нижче — що вийшло й що це означає для рейнджерів. Усі запуски записані локально, а датасети мають версії за вмістом.
виміряно
Правильна відповідь існує без нас: лейбл із незалежного продукту, наприклад Hansen, EFFIS, FIRMS чи розмітка еталона. Помилку моделі видно.
сценарний аналіз
Місця й події реальні, але частина відповіді це наше припущення: дальність огляду, напрямок із камери, шлях до рубки. Число каже «що було б, якщо».
симуляція
І дані, і правильна відповідь із генератора. Перевіряє лише, чи метод розплутує те, що генератор заплутав.
Усі задачі в одній таблиці
Еталон це правильна відповідь, з якою звіряємо. Простий підхід це найпростіше рішення, яке обраний має перемогти. Різниця це обраний мінус простий на тих самих тестових прикладах, з 95 % інтервалом парного бутстрепа. Якщо інтервал перетинає нуль, переваги не доведено.
Одна система для всіх графіків
Кожна задача розказана в однаковому порядку. Кольори й позначки означають одне й те саме на всій сторінці.
Кольори та позначки
Метрики простими словами
Звідки взялись дані
Це сирі набори, які ми завантажили. З них зібрані датасети. Обери набір справа.
Що лежить у кожному наборі даних
Обери датасет справа. Лейбл це правильна відповідь, він є в усіх рядках, зокрема й у тесті. Лише блок «Зв'язок із лейблом» рахуємо на train, щоб не підглядати в тест. Решта описує всі дані.
Як ми знаходимо лідера кожної гілки
16 аналізів, з яких складається перевірка. Праворуч гілки, де аналіз є.
- Аналіз даних (EDA)Кожен сирий набір і кожен датасет: що в ньому, розподіли, лейбли, мапи, яка ознака відділяє лейбл, кореляціїусі
- Вибір лідера на валідаціїУсі підходи гілки порівнюються на валідаційних даних, лідер обирається там, а не на тестіусі
- Тест на нових данихЛідер перевіряється на даних, яких не бачив: інша половина парку, наступні роки, інша пожежаусі
- Довірчий інтервалБутстреп за просторовими чи часовими блоками: наскільки число стійкеFD FM FG LD LM LG AR TP CP SD CD
- Різниця з базовою лінієюЧи кращий лідер за просте правило на тих самих блоках: FIRMS, кліматологія, випадковий вибір, «де вже горіло»FD FG AR LG TP CP
- Стійкість до еталонаТой самий результат проти іншого лейбла: EFFIS, лише 2023, Hansen, падіння NDVIFD FG
- Лише минулеКожен рік прогнозується тільки з даних до нього, як у реальній роботіFG LG TP CP
- Без однієї пожежіМодель вчиться на всіх пожежах, крім однієї, і перевіряється на ній; так для кожноїFG
- Розрізи за рокамиЯк якість змінюється від року до рокуFG LG TP CP
- Криві ROC і PRЯк лідер відділяє позитивні випадки від решти на всіх порогахFG
- Верхня межаСкільки дала б модель, якби лейбли того ж року були відоміFG
- Важливість ознакЯкі ознаки найбільше впливають на рішення моделіFG
- Мапа помилок лідераДе лідер влучив, де помилився, де пропустивFD
- ЧутливістьСкільки веж чи камер треба, інша висота щогли, інші під'їздиTP CP
- Крива часткиЯку частку втрат ловить прогноз залежно від площі, яку перевіряютьLG
- Роздільність знімкаЯк точність падає з 12 до 50 см на піксельSD
Що саме згоріло і скільки гектарів?
Дані
Один рядок — піксель лісу 30 м. Лейбл — чи загинув ліс за Hansen у 2023–24.
Як пожежа 2023 рухалася парком
Кожна точка — теплова детекція супутника. Пересувайте повзунок або натисніть «відтворити».
Як читати: світлі точки — останні 12 годин, тьмяні — вже пройдені.
Лейбли й розбиття на калібрування / тест
Сітка 120 м (агрегат пікселів 30 м). Перемкніть, щоб побачити лейбл або розбиття.
Як читати: помаранчевий — ліс загинув у 2023–24 за Hansen; на вкладці «Розбиття» світлі блоки — тестова половина парку.
Підходи
FIRMS бачить вогонь (тепло), MODIS і Sentinel-2 бачать наслідок: зникла зелень і змінився відбиток згарища.
Як перевіряли й чим міряли
Протокол: просторове розбиття
Парк поділено на блоки H3 (~5 км²). На половині блоків — калібрування: там обираємо поріг і гіперпараметри (82 конфігурації). Друга половина — тест: туди не дивимось до звіту.
- Вибір підходу — лише за калібруванням.
- 95 % CI — bootstrap по тестових блоках.
- Стійкість: та сама метрика з лейблом «лише 2023».
Метрика: F1 (+ precision, recall, IoU)
Площу згарища оцінюють збігом з еталоном (стандарт валідації згарищ CEOS). F1 карає і за перебільшення (низька precision), і за пропуски (низький recall). Accuracy не підходить: фон домінує.
Результат
Тестова половина парку.
Сімейства підходів
Як читати: праворуч — краще; заповнена точка — тест, кільце — калібрування (на ньому обирали). Чим ближче вони, тим менше перенавчання.
Діагностика
Чому переміг саме цей підхід.
Де підхід збігається з еталоном, а де помиляється
Як читати: помаранчевий — згорів і знайдено; світлий — позначено, але ліс живий (перебільшення); сірий — згорів, але пропущено.
Площа згарища: підхід vs еталон
Як читати: вертикальна лінія — площа за Hansen; праворуч від неї — перебільшення.
Як обирали поріг
Як читати: зі зростанням порогу precision росте, recall падає; обираємо поріг з максимальним F1.
Зелень до, після і відновлення 2024–2026
Як читати: «Ураження» — наскільки впала зелень; «Відновлення» — яка частка втраченої зелені повернулась (100 % — як до пожежі).
Скільки зелені повернулося
Як читати: медіана по згорілих пікселях MODIS; 100 % — повне відновлення NDVI (трава й чагарники, ще не дерева).
NDVI трьох ділянок 2021–2026
Як читати: провал у серпні 2023 — пожежа; контрольна ділянка показує звичайні сезонні коливання.
Що це означає
Які теплові точки це одна пожежа?
Дані
Один рядок — сигнал FIRMS. Лейбл — згарище EFFIS, якому він належить.
Підходи
Від «кожна точка окремо» до зведення в часі й просторі.
Як перевіряли й чим міряли
Протокол: вибір на минулому
Параметри обрано на 2021–2023, тест на 2024–2025. Оцінюються лише сигнали, що мають еталон у EFFIS.
Метрика: BCubed F1
Оцінює кожен сигнал окремо: чи стоїть він разом зі своїми і не змішаний із чужими. Великі пожежі не затуляють малі. Вибірка: 279 згарищ.
Результат
Тест 2024–2025.
BCubed F1
Як читати: праворуч краще; вус це 95 % довірчий інтервал.
Діагностика
Де саме помиляється кожен підхід
Найбільші пожежі парку за 2021–2025
Як читати: коло це одна пожежа після зведення, що більше коло, то більше точок.
Що це означає
Чи буде завтра пожежа — вмикати посилене чергування?
Дані
Один рядок — доба. Ознаки — погода NASA POWER. Лейбл — чи була детекція FIRMS.
Календар пожежних днів 2021–2025
Як читати: кожна клітинка — доба; помаранчева — у зоні парку була пожежа, з білою рамкою — почалося нове займання.
Погодні ознаки та пожежні дні
Як читати: лінія — обрана ознака; помаранчеві точки над графіком — дні з пожежею, кільця — нові займання.
Підходи
Правило з однієї ознаки проти моделі, що комбінує погоду.
Як перевіряли й чим міряли
Протокол: перевірка вперед у часі
Модель ніколи не бачить майбутнього: навчаємось на минулих роках, перевіряємо на наступному.
Метрика: PR-AUC (+ recall у 10 % днів)
Пожежних днів ~3 %, тому ROC-AUC тут оманливо високий. PR-AUC чесно показує, наскільки прогноз кращий за випадковий (≈ 0,03). Recall у 10 % днів — операційна відповідь: якщо чергувати посилено в 10 % найризиковіших днів, яку частку пожеж ми накриємо.
- CI — bootstrap 14-денними блоками (пожежі йдуть серіями).
- Дві цілі: день з пожежею і нове займання.
Результат
Прогнози 2023–2025, яких модель не бачила.
PR-AUC підходів
Як читати: праворуч — краще; вертикальна лінія — випадковий прогноз. * — не прогноз: потрібне вчорашнє спостереження.
Діагностика
Щоденний ризик чемпіона й реальні пожежі
Як читати: лінія — оцінка ризику дня; вище горизонтальної лінії — 10 % найризиковіших днів («тривога»); помаранчеві точки — дні з подією.
Що модель вважає важливим
Як читати: стандартизовані коефіцієнти чемпіона; помаранчеві підвищують ризик, сірі — знижують.
Precision–recall криві
Як читати: що вище крива, то менше хибних тривог при тому самому recall. Наведіть, щоб підсвітити.
PR-AUC за тестовими роками
Як читати: 2023 — рік великої пожежі; 2024–2025 майже без пожеж, тому там метрики шумні.
Посуха → пожежі
Як читати: частка пожежних днів у кожному децилі «днів без дощу»; Д10 — найдовші посухи.
Що це означає
Де ризик згоріти найвищий?
Дані
Один рядок — гекс H3 ≈ 0,1 км². Ознаки — рельєф, ліс, люди. Лейбл — чи згорів гекс.
Лейбли та ознаки гексів
Як читати: лейбли — помаранчевий «так», сірий «ні»; ознаки — від темного (мало) до світлого (багато).
Підходи
Довга історія багатьох подій проти однієї свіжої пожежі.
Як перевіряли й чим міряли
Протокол: вибір на 2022, звіт на 2023
Підхід обираємо за пожежею 2022 (валідація), звітуємо на пожежі 2023, якої не бачив жоден підхід.
- CI — bootstrap по блоках H3 (~5 км²).
- Стеля — що дала б модель, знаючи лейбли того ж 2023 року.
- Стійкість — ROC-AUC з трьома визначеннями «згоріло».
Метрика: ROC-AUC (+ recall у 20 % площі)
Нам потрібне ранжування ділянок. ROC-AUC — шанс, що згоріла ділянка стоїть вище незгорілої. У 2023 році згоріло 67 % парку, тому PR-AUC майже тривіальний. Recall у 20 % площі — частка згарища в 20 % площі під посиленим наглядом.
Результат
Тест — пожежа 2023.
ROC-AUC підходів на 2023 році
Як читати: праворуч — краще; ліва лінія — випадково (0,5), права — стеля. Форма точки — на чому вчився підхід.
Діагностика
Валідація 2022 vs тест 2023
Як читати: якби вибір на 2022 році переносився на 2023-й, точки лежали б уздовж діагоналі. Вони лежать майже навпаки.
Чому валідація 2022 обманює
Як читати: частка гексів, що згоріли у 2023 році, серед тих, що вже горіли у 2022-му, і серед решти.
Що важливо для моделі на історії (LightGBM)
Як читати: на скільки падає ROC-AUC, якщо перемішати ознаку (permutation importance, тест 2023).
Стійкість до визначення «згоріло»
Як читати: синє — краще за випадкове, червоне — гірше. Добрий підхід синій у всіх трьох колонках.
Прогноз, факт і ризик 2026
Що це означає
Пожежа вже йде: куди вона піде завтра й післязавтра?
Дані
Один рядок — (пожежа, момент t, ще не згорілий гекс H3 ≈ 0,74 км² у радіусі 12 км). Лейбл — чи згорить до t+24 / t+48 год.
28 багатоденних пожеж Греції 2021–2025
Як читати: кожне коло — пожежа з FIRMS VIIRS (розмір — кількість детекцій); помаранчева з рамкою — Еврос-2023, недоторканий тест. Решта 27 — для навчання й вибору підходу.
Підходи
Від «росте рівномірно» до моделі, що знає вітер, схил і паливо.
Як перевіряли й чим міряли
Протокол: leave-one-fire-out + недоторканий тест
Підхід порівнюємо так: вчимося на 26 пожежах, перевіряємо на 27-й — і так для кожної (параметри правил і поріг моделі підбираються лише на навчальних пожежах). Переможця обираємо за цим. Еврос-2023 не бере участі ні в чому — на ньому лише фінальний звіт.
- Крок — доба (VIIRS пролітає ~2 рази на добу); «зріз» о 13:00 UTC.
- Вітер на наступну добу — фактичний NASA POWER як заміна прогнозу погоди.
- CI — bootstrap по пожежах.
Метрика: F1 нової площі (+ похибка фронту)
Рахуємо лише нову площу — те, що ще не горіло. F1 карає і за перебільшення, і за пропуски; macro — кожна пожежа має однакову вагу (інакше Еврос заглушив би все). Похибка фронту — медіана відстані від реально згорілого гекса до найближчого спрогнозованого: «на скільки метрів промахнулись».
Результат
Вибір — 27 пожеж (LOFO), тест — Еврос-2023.
F1 нової площі
Як читати: праворуч — краще; кільце з вусом — середнє по 27 пожежах з 95 % CI (на цьому обирали); заповнена точка — тест на Евросі-2023.
Діагностика
Еврос-2023 день за днем: прогноз проти того, що сталося
Як читати: темні гекси — вже згоріло на момент прогнозу; помаранчевий — згоріло й спрогнозовано; світлий — спрогнозовано, але не згоріло; сірий — згоріло, але пропущено. Стрілка — куди дме вітер.
Якість прогнозу по днях пожежі
Як читати: F1 нової площі на кожну добу Евросу-2023; у дні вибухового росту прогнозувати найважче.
На яких пожежах підхід працює
Як читати: кожна точка — одна з 27 пожеж (LOFO); праворуч — краще. Чемпіон проти базової «рівномірний ріст».
Що це означає
Де щось регулярно повторюється?
Дані
Реальні пожежі Греції 2021–2025: точки FIRMS, зведені правилом FM. Один рядок це пожежа. Потік демо тут не використовується.
Підходи
Чи дають контури ареалів більше за простіші способи показати гарячі місця.
Як перевіряли й чим міряли
Протокол: лише минуле
Ареали року Y будуються лише з пожеж попередніх років. Радіус і вікно обрано на 2023, тест на 2024 і 2025. Усі методи отримують рівно однакову площу: 1,5 % гексів H3.
Метрика: PAI
У скільки разів щільність нових пожеж у зоні вища за середню. 1 це випадковий вибір. Також частка нових пожеж, що потрапила в зону.
Як результат залежить від параметрів
Обрано на 2023, клітинка з рамкою.
PAI на валідації: радіус × вікно
Як читати: рядки це радіус DBSCAN, стовпці кількість минулих років. Що яскравіше, то краще.
Результат
Тест 2024–2025.
PAI на тесті, однакова площа 1,5 %
Як читати: праворуч краще, 1 це випадковий вибір; вус це 95 % довірчий інтервал.
Усі числа
Діагностика
Ареали за всі роки
Як читати: контур це ареал, що яскравіше, то більше пожеж у ньому. Показано 60 найбільших.
Ареали демо
Це потік демо, здебільшого симульований. Тут метрика не рахується.
Що це означає
Чи можна знайти вирубку там, де її майже немає?
Дані
Один рядок — піксель лісу 30 м за рік. Лейбл — втрата крон за Hansen. Згарища викинуті за EFFIS і FIRMS.
Підходи
Від одного знімка до моделі на всіх смугах.
Як перевіряли й чим міряли
Протокол: половина парку × час
Парк один, тому вибір на одній половині за 2018–2021, тест на другій за 2022–2024, і навпаки.
Метрика: F1 площі
Баланс точності й повноти за площею. У тесті лише 14 га втрат, тому інтервали широкі.
Результат
Правило обирає варіант за калібруванням, а не за тестом.
F1 площі на тесті
Як читати: праворуч краще; вус це 95 % довірчий інтервал; жирний рядок обрано на калібруванні.
Діагностика
Точність і повнота
Кандидати за 2024 рік
Як читати: зелені є й у NASA OPERA, жовті лише в Sentinel-2. Це кандидати, а не доведені рубки.
Що це означає
Чи допомагає збіг двох супутникових джерел?
Дані
Ті самі пікселі, що в LD, плюс знімки стану NASA OPERA DIST-ALERT кожні 10 днів.
Підходи
Два джерела окремо й два способи їх поєднати.
Як перевіряли й чим міряли
Протокол
Sentinel-2 навчається на 2018–2021, правило OPERA обирається на 2023, тест на 2024 на іншій половині парку.
Метрика: точність
Скільки з позначеного справді втрата лісу. Це важливо для виїзду: помилковий кандидат коштує дня роботи.
Результат
Лише 5,5 га втрат у тесті.
Точність на тесті
Як читати: праворуч краще; показано точність, бо заради неї поєднують джерела.
Що це означає
Де рубатимуть наступного року?
Дані
Один рядок — гекс H3 × рік. Лейбл — щонайменше 0,5 га нової втрати без згарищ.
Підходи
Від випадкового вибору до моделі на історії й ландшафті.
Як перевіряли й чим міряли
Протокол: лише минуле
Прогноз на рік Y вчиться на роках до Y. Вибір 2018–2020, тест 2021–2024.
Метрика: частка втрат у 11 % гексів
Якщо перевірити лише 11 % парку, скільки нових втрат туди потрапить. Випадковий вибір дає близько 11 %, на цьому тесті 6 %.
Результат
У тесті 39 га втрат.
Частка втрат у 11 % гексів
Як читати: праворуч краще; вус це 95 % довірчий інтервал; жирний рядок обрано на валідації.
Діагностика
По роках
Як читати: прогноз на рік Y вчиться лише на роках до Y. 2018–2020 валідація, далі тест.
Частка втрат залежно від площі, яку перевіряємо
Як читати: замість одного порога 11 % показано всю криву: скільки нових втрат потрапляє в найризиковіші 2, 5, 11, 20, 30 і 50 % гексів. Тест 2021–2024.
Прогноз на 2025 рік
Як читати: фіолетові гекси це 11 % найризикованіших.
Що це означає
Чи видно машини згори, і за якої роздільності?
Дані
Відкритий набір VEDAI з рамками машин, розмитий до трьох роздільностей.
Підходи
Дві готові моделі без донавчання.
Як перевіряли й чим міряли
Протокол: 50/50 кадрів
Поріг впевненості обрано на половині кадрів, числа на другій. Модель не вчилась на цих даних.
Метрика: AP@0.5
Якість рамок незалежно від порога: рамка правильна, якщо перекривається з еталоном щонайменше наполовину.
Результат
12,5 → 50 см це симуляція зниження якості.
AP@0.5 за роздільністю
Як читати: праворуч краще; вус це 95 % довірчий інтервал.
Усі числа
Як читати: поріг обрано на половині кадрів, решта чисел на другій половині; інтервал 95 %.
Скільки пікселів займає легкове авто
Як читати: що менше пікселів, то гірше. Перехід від 25 до 50 см робить авто вдвічі меншим, а AP падає більш ніж удвічі.
Відкриті знімки парку
Яку роздільність дають відкриті знімки Esri парку і що з неї видно. Розмітки для парку немає, тому точності тут немає.
Чому людей не видно із супутника
Як читати: скільки пікселів займає об'єкт на відкритому знімку Esri (0,45 м/піксель); детекторам потрібно ~10 px.
Що це означає
Наскільки точні детектори камер і фотопасток?
Дані
D-Fire: дим і вогонь із рамками. ENA24: кадри фотопасток. Це фото, відео в еталоні немає.
Підходи
GroundingDINO за текстовим запитом, без донавчання.
Як перевіряли й чим міряли
Протокол: 50/50 кадрів
Поріг на половині кадрів, числа на другій. На ENA24 лише 22 кадри з машинами.
Метрики
AP@0.5 там, де є рамки. На ENA24 лише лейбл кадру, тому ROC-AUC, повнота й частка хибних кадрів.
Результат
AP@0.5, D-Fire
Як читати: праворуч краще; вус це 95 % довірчий інтервал.
Усі числа
Перевірка на відео й фото
Запис пожежі в Карсті 2022 року і 8 фото з парку. Розмітки немає, тому це не точність, а звірка з порогом, обраним на еталоні.
YOLO-World на фото з парку
Як читати: рамки — знайдені об'єкти з впевненістю; зелений — людина, світлий — авто.
Що це означає
Чи побачили б вежі майбутні пожежі?
Дані
Точки NASA FIRMS 2021–2025 у парку, зведені в пожежі правилом FM (2 км, 5 діб). Рельєф SRTM дає, що видно з кожного місця. Вежа: огляд 360°, щогла 15 м, 10 км, дим 20 м над землею.
Підходи
Три бази й жадібне покриття з різною вагою: кожна наступна вежа додає найбільше нового. Параметр: скільки найвищих місць розглядаємо, 30, 100 чи 300.
Як перевіряли й чим міряли
Протокол: лише минуле
Вежі ставимо за тим, що відомо до року Y. Валідація: дані до 2022, пожежі 2022. Тест: дані до 2023, пожежі 2023–2025. Вибір лише на валідації, а за рівності береться простіший варіант. Вага з вкладки «Ризик» вчилась на пожежі 2022, тому на валідації її немає: вона лише для порівняння на тесті.
Метрики
Перша точка: частка пожеж, чия перша точка FIRMS у полі зору веж. Це головна метрика. Хоч одна точка: пожежу бачили бодай колись. Згорілі гекси: частка гексів пожежі 2023 у полі зору. Інтервал це бутстреп за пожежами.
Результат
Усі рядки цього графіка на 300 місцях-кандидатах, тест 2023–2025.
Яку частку пожеж бачать вежі
Як читати: праворуч краще; вус це 95 % інтервал. Вага ризику не проходила валідацію, її показано для порівняння. Усі варіанти «покриття» дали ті самі три місця, тому їхні стовпчики однакові.
Підібраний параметр: скільки місць розглядати
Як читати: кожна клітинка це частка пожеж за першою точкою на валідації й тесті. Вибір робився лише за числом до риски.
Усі числа, тест
Діагностика: скільки веж і за яких умов
Місця ті самі, змінюємо лише умови. Це не вибір, а перевірка стійкості.
По роках
Як читати: прогресивний протокол: кожен рік вежі стоять лише за даними до нього. 2022 це валідація, далі тест.
Скільки веж потрібно
Як читати: кожна нова вежа додана жадібно, тест 2023–2025.
Якщо умови інші
Як читати: перший рядок це базові умови, далі змінюється одне число.
Обрані місця
Як читати: остаточні вежі за всіма даними до 2026 року. Кольорові гекси видно з веж, колір це рівні ризику. Кружок це перша точка пожежі 2021–2025: світлий у полі зору, сірий ні.
Що це означає
Де на мапі дим, який камера бачить у кадрі?
Дані
Реальні гекси парку, де NASA FIRMS бачив вогонь у 2022 і 2023 роках, і рельєф SRTM. Реальних детекцій диму з відомими координатами немає, тому напрямок із камери симульовано.
Підходи
Геометрія, без навчання. Спільний параметр лише один: ширина клина.
Як перевіряли й чим міряли
Протокол: валідація 2022, тест 2023
Для кожного гекса з пожежею беремо справжній напрямок і додаємо похибку з відхиленням 1, 3, 5 або 10°. Ширину клина (1, 2 чи 3 відхилення) обрано на 2022: найменша зона, що накриває місце щонайменше в 90 % випадків.
Метрики
Влучання: справжній гекс у зоні. Площа зони в км². Відстань від місця до центру зони.
Результат
Це перевірка геометрії на симуляції, а не точність детектора.
Площа зони пошуку на тесті 2023
Як читати: ліворуч краще: менша зона, менше їхати перевіряти. Перемикач змінює похибку напрямку.
Усі числа
Яку частку гексів із пожежею 2023 бачить камера
Як читати: три колонки міняють висоту щогли за дальності 10 км, дві міняють дальність за щогли 15 м. Вежі стоять там, де їх поставила гілка TP.
Як це виглядає: відео й мапа синхронізовані
Запустіть відео.
Камера CAM-02: детекція в кадрі → клин на мапі
Як читати: положення рамки по горизонталі дає азимут; дальність камера не знає, тому підсвічуємо весь клин — але лише ті гекси, які видно з вишки з урахуванням рельєфу. Помаранчевий — дим, зелений — люди. Відео це запис пожежі в Карсті 2022 року, показаний як демо: він ілюструє механіку, а не Дадью.
Що це означає
Які місця на дорогах перехоплять під'їзд до рубок?
Дані
Дороги й стежки OpenStreetMap, нові втрати лісу Hansen поза згарищами, пожежі FIRMS. Наземна камера чи фотопастка бачить лише свій відрізок дороги, тому важить не рельєф, а те, через що проходить під'їзд.
Підходи
Від кожного в'їзду в парк будуємо дерево найкоротших під'їздів. Камера перехоплює все, що лежить за нею на цьому дереві. Параметр: дозволяти всі дороги чи лише лісові.
Як перевіряли й чим міряли
Протокол: лише минуле
Місця обираємо за тим, що відомо до року Y. Валідація: дані до 2018, рубки 2018–2020, пожежі 2021. Тест: дані до 2021, рубки 2021–2024, пожежі 2022–2025. Вибір лише на валідації. Перший варіант тесту, 2023–2024, мав лише 2 рубки і дав 0 %. Тому тест розширено, а той результат лишається в звіті нижче.
Метрики
Рубки: частка рубок, найкоротший під'їзд до яких проходить через камеру. Це головна метрика. Пожежі: те саме для перших точок пожеж. Ліс: частка лісових гексів за камерами. Інтервал це бутстреп за подіями.
Результат
Тест 2021–2024, 10 камер.
Що перехоплюють камери
Як читати: праворуч краще; вус це 95 % інтервал для рубок.
Усі числа
Як читати: вибір робився лише за валідацією.
Діагностика: скільки камер і чи тримається
Місця обрані, тепер перевіряємо стійкість. Це не вибір.
По роках
Як читати: прогресивний протокол: кожен рік місця обрано лише за даними до нього. 2018–2020 валідація, далі тест.
Скільки камер потрібно
Перевірка припущень
Як читати: перший рядок це основний тест, далі те, що може зменшити результат.
Обрані місця
Як читати: остаточні 10 місць за всіма даними. Жовті лінії це під'їзди, які проходять через камери. Кружок це рубка: зелений перехоплено, сірий ні.
Що це означає
Які сигнали з усіх джерел це одна подія?
Дані
Потік демо за 90 днів: один рядок це сигнал. Більшість сигналів симульовані.
Сигнали за джерелом
Підходи
Три способи зведення й сітка параметрів для найскладнішого.
Як перевіряли й чим міряли
Протокол: вибір на перших 60 днях
Для кожного типу події окремо: пожежа, вирубка, машини. Параметри обрано на перших 60 днях, тест на останніх 30. Перебрано 36 пар радіус × час.
Еталон і метрика
Еталон це задум симуляції: сигнали одного сценарію й одного дня вважаються однією подією. Це не спостереження. Метрика BCubed F1, як у FM.
Як результат залежить від параметрів
Обрано на валідації, клітинка з рамкою.
BCubed F1 на валідації: радіус × час
Як читати: рядки це радіус, стовпці час. Що темніше, то гірше.
Результат
Тест, останні 30 днів.
BCubed F1 на тесті: обраний варіант проти простіших
Як читати: BCubed F1 від 0 до 1, праворуч краще. Він показує, чи стоїть кожен сигнал разом зі своєю подією і не змішаний із чужими. Тест це останні 30 днів, еталон це симуляція.
Що це означає
Які зони інтересу є в парку, якого типу, з якими рисами в часі й просторі і з яких джерел?
Ареали застосунку кажуть, що відбувається в кожній зоні: тип, чи вона активна, згасла чи відновилась, чи розширюється, як часто там щось стається й хто це бачить. Перевіряємо, наскільки ці риси збігаються із задумом сценарію.
Постановка
AS · Зони інтересу на симуляції: чи правильно ареали застосунку читають, що відбувається на місці.
Задача ареали застосунку (scripts/03_signals.py, superareals) після зведення сигналів (XM) кажуть про кожну зону інтересу: тип
(пожежа, вирубка, незаконний проїзд), життєвий цикл (активний, згас, відновився, новий), динаміку (розширюється,
звужується, стабільний), регулярність і джерела. Тут перевіряємо, наскільки ці риси збігаються з тим, що закладено в сценарій.
Дані генератор потоку застосунку, прогнаний багато разів з різною випадковістю (місця на треках OSM, які джерела побачили подію,
розкид точок) окремо для кожного парку. Знімки Esri не качаються: супутник «бачить» у рамці ~0,9 × 0,9 км навколо місця.
Еталон задум генератора: сценарій = зона інтересу з відомими типом, життєвим циклом, динамікою, регулярністю й джерелами.
Це еталон симуляції, а не спостереження.
Протокол прогони 0…N_VAL−1 — вибір підходу, решта — тест. Тест не бере участі у виборі. CI — bootstrap по прогонах.
Метрики на кожен прогін п'ять оцінок: зведення сигналів у зони (BCubed F1: сигнали одного сценарію в одному ареалі);
точність життєвого циклу, динаміки (лише де ареал має ≥ 4 подій) і регулярності на рівні сценарію; збіг джерел (Jaccard).
Головна метрика — середнє п'яти.
Підходи радіус ареалу (множник до AREAL_EPS_KM застосунку) × поріг динаміки (зміна розкиду, %); «застосунок» — ×1 і 30 %.
Базова лінія «прості правила»: згас, якщо остання подія давніше 21 дня, інакше активний; динаміка завжди «стабільний».
Обмеження результат показує, як правила ареалів поводяться на нашій симуляції; генератор сам задає розкид і частоту подій,
тож висока точність тут не доводить роботу на реальному потоці.
Як обрано лідера
найкраща середня точність на прогонах вибору; потік симульований. Тест на вибір не впливає; негативні результати лишаються в таблиці.
Лідер
| група | підхід | вибір | тест: середня точність рис зони | 95 % CI |
|---|---|---|---|---|
| — | радіус ×1, поріг 50 % | 0.785 | 0.785 | 0.78–0.79 |
Усі підходи
Помаранчевий це лідер; риска це метрика вибору; вуса це 95 % довірчий інтервал на тесті.
| група | підхід | вибір | тест | CI від | CI до | лідер | гіпотеза |
|---|---|---|---|---|---|---|---|
| — | радіус ×1, поріг 50 % | 0.785 | 0.785 | 0.78 | 0.79 | так | |
| — | радіус ×1.5, поріг 50 % | 0.784 | 0.784 | 0.779 | 0.789 | ||
| — | радіус ×0.5, поріг 50 % | 0.776 | 0.776 | 0.771 | 0.781 | ||
| — | застосунок: радіус ×1, поріг 30 % | 0.774 | 0.772 | 0.766 | 0.777 | ||
| — | радіус ×1.5, поріг 30 % | 0.771 | 0.77 | 0.765 | 0.776 | ||
| — | радіус ×0.5, поріг 30 % | 0.768 | 0.767 | 0.761 | 0.772 | ||
| — | прості правила | 0.752 | 0.753 | 0.751 | 0.755 | ||
| — | радіус ×1, поріг 15 % | 0.736 | 0.737 | 0.732 | 0.743 | ||
| — | радіус ×0.5, поріг 15 % | 0.734 | 0.735 | 0.73 | 0.741 | ||
| — | радіус ×1.5, поріг 15 % | 0.734 | 0.736 | 0.73 | 0.742 |
Усі метрики алгоритмів
Кожен рядок це алгоритм або його налаштування, кожна колонка окрема метрика.
| approach | eps | thr | val_score | test_score | ci_lo | ci_hi | n_val | n_test | areals | test_grouping | test_lifecycle | test_dynamics | test_regularity | test_sources | test_score_dadia | test_score_evros_delta | test_score_emt |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| радіус ×0.5, поріг 15 % | 0.5 | 15 | 0.734 | 0.735 | 0.73 | 0.741 | 150 | 150 | 15.2 | 0.988 | 0.728 | 0.528 | 0.648 | 0.785 | 0.736 | 0.732 | 0.738 |
| радіус ×0.5, поріг 30 % | 0.5 | 30 | 0.768 | 0.767 | 0.761 | 0.772 | 150 | 150 | 15.2 | 0.988 | 0.728 | 0.684 | 0.648 | 0.785 | 0.77 | 0.766 | 0.765 |
| радіус ×0.5, поріг 50 % | 0.5 | 50 | 0.776 | 0.776 | 0.771 | 0.781 | 150 | 150 | 15.2 | 0.988 | 0.728 | 0.731 | 0.648 | 0.785 | 0.779 | 0.773 | 0.776 |
| радіус ×1, поріг 15 % | 1 | 15 | 0.736 | 0.737 | 0.732 | 0.743 | 150 | 150 | 12.2 | 0.998 | 0.728 | 0.53 | 0.648 | 0.785 | 0.739 | 0.734 | 0.74 |
| застосунок: радіус ×1, поріг 30 % | 1 | 30 | 0.774 | 0.772 | 0.766 | 0.777 | 150 | 150 | 12.2 | 0.998 | 0.728 | 0.7 | 0.648 | 0.785 | 0.775 | 0.769 | 0.771 |
| радіус ×1, поріг 50 % | 1 | 50 | 0.785 | 0.785 | 0.78 | 0.79 | 150 | 150 | 12.2 | 0.998 | 0.728 | 0.766 | 0.648 | 0.785 | 0.792 | 0.779 | 0.784 |
| радіус ×1.5, поріг 15 % | 1.5 | 15 | 0.734 | 0.736 | 0.73 | 0.742 | 150 | 150 | 11.9 | 0.995 | 0.729 | 0.524 | 0.647 | 0.783 | 0.739 | 0.73 | 0.739 |
| радіус ×1.5, поріг 30 % | 1.5 | 30 | 0.771 | 0.77 | 0.765 | 0.776 | 150 | 150 | 11.9 | 0.995 | 0.729 | 0.698 | 0.647 | 0.783 | 0.776 | 0.765 | 0.77 |
| радіус ×1.5, поріг 50 % | 1.5 | 50 | 0.784 | 0.784 | 0.779 | 0.789 | 150 | 150 | 11.9 | 0.995 | 0.729 | 0.766 | 0.647 | 0.783 | 0.793 | 0.775 | 0.784 |
| прості правила | 1 | 30 | 0.752 | 0.753 | 0.751 | 0.755 | 150 | 150 | 12.2 | 0.998 | 0.583 | 0.75 | 0.648 | 0.785 | 0.753 | 0.752 | 0.753 |
Що ми знаємо — і чого ще ні
Наступні кроки
- Більше подій. Пожежних днів і рубок у тестах замало: розширити дані на сусідні парки й більше років.
- Валідація на кількох сезонах і регіонах для ризику ділянок замість однієї пожежі.
- Розмітка кадрів реальної вежі й супутникових знімків лісу, щоб виміряти детектори в умовах парку.
- Прив'язка реальної вежі. Зіставити кадри з відомими точками на мапі й виміряти справжню похибку напрямку: зараз вона симульована.
- Журнали проїздів. Перевірити місця камер на реальних GPS-треках або записах: зараз шлях припущений.
- Реальний потік сигналів. Перевірити зведення на справжніх камерах і дзвінках, а надійність джерел навчити на підтвердженнях рейнджерів: зараз вона задана вручну.
Загрози валідності
- Одна велика пожежа. Еврос-2023 домінує в згарищі, ризику ділянок, поширенні й камері на координати. Результати можуть не переноситися на дрібні пожежі.
- Лейбли з продуктів, а не з поля. Hansen бачить загибель крон, але не її причину. EFFIS пропускає частину дрібних пожеж. FIRMS бачить вогонь лише під час прольоту супутника.
- Мало подій у тестах. 30 пожежних днів, кілька десятків рубок і пожеж у задачах розміщення, кілька гектарів у детекції вирубки. Тому інтервали широкі, а різниці з базою часто не доведені.
- Погода в одній точці. NASA POWER у центрі парку не бачить різниці між ділянками.
- Сценарні припущення й симуляція. Дальність огляду й висота диму у вежах, похибка напрямку в камері, найкоротший шлях до рубки, генератор потоку сигналів. Ці задачі позначено окремо.
- Зміни протоколу. У «Де ставити камери» перший тест мав 2 рубки, після цього всі три задачі розміщення перейшли на прогресивний протокол. У «Дні пожежі» підхід раніше обирали на тестових роках, тепер на 2022. Обидві зміни описано в блоках.
- Детектори не донавчались. Їх перевірено на відкритих еталонах, а не на кадрах із парку.
Джерела
- Hansen M.C. et al. (2013). High-resolution global maps of 21st-century forest cover change. Science 342, 850–853. Global Forest Change v1.12.
- Schroeder W. et al. (2014). The New VIIRS 375 m active fire detection data product. Remote Sensing of Environment 143, 85–96. NASA FIRMS.
- Giglio L. et al. (2016). The collection 6 MODIS active fire detection algorithm and fire products. Remote Sensing of Environment 178, 31–41.
- Didan K. (2021). MODIS/Terra Vegetation Indices 16-Day L3 Global 250 m, MOD13Q1 v061. NASA LP DAAC; вирізки ORNL DAAC.
- NASA POWER Project, Prediction Of Worldwide Energy Resources, NASA Langley Research Center.
- Copernicus Sentinel-2 L2A, ESA; каталог Earth Search, Element 84.
- OPERA Land Surface Disturbance Alert (DIST-ALERT-HLS), NASA JPL; шар NASA GIBS.
- San-Miguel-Ayanz J. et al. (2012). Comprehensive monitoring of wildfires in Europe: the European Forest Fire Information System (EFFIS).
- Van Wagner C.E. (1987). Development and structure of the Canadian Forest Fire Weather Index System. Canadian Forestry Service, Technical Report 35.
- Key C.H., Benson N.C. (2006). Landscape Assessment: ground measure of severity, the Composite Burn Index, and remote sensing of severity, the Normalized Burn Ratio. FIREMON, USDA Forest Service.
- Anderson H.E. (1983). Predicting wind-driven wild land fire size and shape. USDA Forest Service, Research Paper INT-305.
- Birant D., Kut A. (2007). ST-DBSCAN: an algorithm for clustering spatial–temporal data. Data & Knowledge Engineering 60, 208–221.
- Amigó E. et al. (2009). A comparison of extrinsic clustering evaluation metrics based on formal constraints. Information Retrieval 12, 461–486. BCubed.
- Razakarivony S., Jurie F. (2016). Vehicle detection in aerial imagery: a small target detection benchmark. JVCIR 34, 187–203. VEDAI.
- Venâncio P.V.A.B. et al. (2022). An automatic fire detection system based on deep convolutional neural networks for low-power, resource-constrained devices. Neural Computing and Applications. D-Fire.
- Yousif H. et al. (2019). Dynamic programming selection of object proposals for sequence-level animal species classification in the wild. ENA24, LILA BC.
- Liu S. et al. (2023). Grounding DINO. arXiv:2303.05499. Cheng T. et al. (2024). YOLO-World. CVPR. Ultralytics YOLO11.
- AWS Terrain Tiles (SRTM). OpenStreetMap contributors, ODbL. H3, Uber Technologies.
Використання ШІ
Код, аналіз і тексти цієї сторінки підготовлено з допомогою ШІ-асистента Claude (Anthropic). Протоколи, перевірку чисел і висновки визначали автори.