Data Science · відкриті супутникові дані · НП Дадья–Лефкімі–Суфлі, Греція

Що ми виміряли в лісах Дадьї, а що ні

Для кожної задачі ми порівняли конкуруючі підходи на лейблах з незалежних супутникових продуктів і перевірили їх на даних, яких вони не бачили. Нижче — що вийшло й що це означає для рейнджерів. Усі запуски записані локально, а датасети мають версії за вмістом.

виміряно

Правильна відповідь існує без нас: лейбл із незалежного продукту, наприклад Hansen, EFFIS, FIRMS чи розмітка еталона. Помилку моделі видно.

сценарний аналіз

Місця й події реальні, але частина відповіді це наше припущення: дальність огляду, напрямок із камери, шлях до рубки. Число каже «що було б, якщо».

симуляція

І дані, і правильна відповідь із генератора. Перевіряє лише, чи метод розплутує те, що генератор заплутав.

Усі задачі в одній таблиці

Еталон це правильна відповідь, з якою звіряємо. Простий підхід це найпростіше рішення, яке обраний має перемогти. Різниця це обраний мінус простий на тих самих тестових прикладах, з 95 % інтервалом парного бутстрепа. Якщо інтервал перетинає нуль, переваги не доведено.

Як читати сторінку

Одна система для всіх графіків

Кожна задача розказана в однаковому порядку. Кольори й позначки означають одне й те саме на всій сторінці.

Кольори та позначки

Помаранчевий — пожежа / згоріло. Позитивний клас у даних: день з пожежею, згорілий піксель чи гекс, дим у кадрі.
Сірий — ні / фон. Негативний клас.
Синій — ML-модель або алгоритм, що вчиться чи кластеризує (LogReg, LightGBM, ST-DBSCAN…).
Зелений — індекс або правило з порогом, без навчання (падіння NDVI, «сухі дні»…).
Світло-сірий — базова лінія: найпростіше рішення, яке треба перемогти.
чемпіонобраний протоколом на даних валідації — не на тесті. у продукті — модель, що працює в застосунку.
Вус довкола точки — 95 % довірчий інтервал: де опиниться результат на інших ділянках чи днях.

Метрики простими словами

Precisionз того, що позначили, яка частка правильна. «Не кричимо даремно»
Recallяку частку реальних подій знайшли. «Нічого не пропускаємо»
F1 · FD, згарищебаланс precision і recall. 1 — ідеальний збіг, 0 — жодного.
PR-AUC · FG, день пожежіякість ранжування рідкісних подій. Випадковий прогноз = частка подій (~0,03).
ROC-AUC · FG, ризик ділянокшанс, що згоріла ділянка отримає вищий ризик за незгорілу. 0,5 — випадково, 1 — ідеально.
Recall у k %частка подій, яку ловимо, стежачи лише за k % найризиковіших днів або площі.
1Дані2Підходи3Як перевіряли й чим міряли4Результат5Діагностика6Що це означає
Набори даних

Звідки взялись дані

Це сирі набори, які ми завантажили. З них зібрані датасети. Обери набір справа.

Датасети

Що лежить у кожному наборі даних

Обери датасет справа. Лейбл це правильна відповідь, він є в усіх рядках, зокрема й у тесті. Лише блок «Зв'язок із лейблом» рахуємо на 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
FD
FD · детекція пожежі

Що саме згоріло і скільки гектарів?

1

Дані

Один рядок — піксель лісу 30 м. Лейбл — чи загинув ліс за Hansen у 2023–24.

NASA FIRMS · пожежа 2023, яку міряє FD

Як пожежа 2023 рухалася парком

Кожна точка — теплова детекція супутника. Пересувайте повзунок або натисніть «відтворити».

Як читати: світлі точки — останні 12 годин, тьмяні — вже пройдені.

Лейбли й розбиття на калібрування / тест

Сітка 120 м (агрегат пікселів 30 м). Перемкніть, щоб побачити лейбл або розбиття.

Як читати: помаранчевий — ліс загинув у 2023–24 за Hansen; на вкладці «Розбиття» світлі блоки — тестова половина парку.

2

Підходи

FIRMS бачить вогонь (тепло), MODIS і Sentinel-2 бачать наслідок: зникла зелень і змінився відбиток згарища.

3

Як перевіряли й чим міряли

Протокол: просторове розбиття

Парк поділено на блоки H3 (~5 км²). На половині блоків — калібрування: там обираємо поріг і гіперпараметри (82 конфігурації). Друга половина — тест: туди не дивимось до звіту.

  • Вибір підходу — лише за калібруванням.
  • 95 % CI — bootstrap по тестових блоках.
  • Стійкість: та сама метрика з лейблом «лише 2023».

Метрика: F1 (+ precision, recall, IoU)

Площу згарища оцінюють збігом з еталоном (стандарт валідації згарищ CEOS). F1 карає і за перебільшення (низька precision), і за пропуски (низький recall). Accuracy не підходить: фон домінує.

4

Результат

Тестова половина парку.

Сімейства підходів

Як читати: праворуч — краще; заповнена точка — тест, кільце — калібрування (на ньому обирали). Чим ближче вони, тим менше перенавчання.

5

Діагностика

Чому переміг саме цей підхід.

Де підхід збігається з еталоном, а де помиляється

Як читати: помаранчевий — згорів і знайдено; світлий — позначено, але ліс живий (перебільшення); сірий — згорів, але пропущено.

збігперебільшенняпропускліс без змін

Площа згарища: підхід vs еталон

Як читати: вертикальна лінія — площа за Hansen; праворуч від неї — перебільшення.

Як обирали поріг

Як читати: зі зростанням порогу precision росте, recall падає; обираємо поріг з максимальним F1.

MODIS MOD13Q1 · 250 м

Зелень до, після і відновлення 2024–2026

Як читати: «Ураження» — наскільки впала зелень; «Відновлення» — яка частка втраченої зелені повернулась (100 % — як до пожежі).

Скільки зелені повернулося

Як читати: медіана по згорілих пікселях MODIS; 100 % — повне відновлення NDVI (трава й чагарники, ще не дерева).

NDVI трьох ділянок 2021–2026

Як читати: провал у серпні 2023 — пожежа; контрольна ділянка показує звичайні сезонні коливання.

6

Що це означає

FM
FM · матчинг пожежі

Які теплові точки це одна пожежа?

1

Дані

Один рядок — сигнал FIRMS. Лейбл — згарище EFFIS, якому він належить.

2

Підходи

Від «кожна точка окремо» до зведення в часі й просторі.

3

Як перевіряли й чим міряли

Протокол: вибір на минулому

Параметри обрано на 2021–2023, тест на 2024–2025. Оцінюються лише сигнали, що мають еталон у EFFIS.

Метрика: BCubed F1

Оцінює кожен сигнал окремо: чи стоїть він разом зі своїми і не змішаний із чужими. Великі пожежі не затуляють малі. Вибірка: 279 згарищ.

4

Результат

Тест 2024–2025.

BCubed F1

Як читати: праворуч краще; вус це 95 % довірчий інтервал.

5

Діагностика

Де саме помиляється кожен підхід

Найбільші пожежі парку за 2021–2025

Як читати: коло це одна пожежа після зведення, що більше коло, то більше точок.

6

Що це означає

FG
FG · день пожежі

Чи буде завтра пожежа — вмикати посилене чергування?

1

Дані

Один рядок — доба. Ознаки — погода NASA POWER. Лейбл — чи була детекція FIRMS.

Календар пожежних днів 2021–2025

Як читати: кожна клітинка — доба; помаранчева — у зоні парку була пожежа, з білою рамкою — почалося нове займання.

день з пожежеюнове займаннябез пожежі

Погодні ознаки та пожежні дні

Як читати: лінія — обрана ознака; помаранчеві точки над графіком — дні з пожежею, кільця — нові займання.

2

Підходи

Правило з однієї ознаки проти моделі, що комбінує погоду.

3

Як перевіряли й чим міряли

Протокол: перевірка вперед у часі

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

Метрика: PR-AUC (+ recall у 10 % днів)

Пожежних днів ~3 %, тому ROC-AUC тут оманливо високий. PR-AUC чесно показує, наскільки прогноз кращий за випадковий (≈ 0,03). Recall у 10 % днів — операційна відповідь: якщо чергувати посилено в 10 % найризиковіших днів, яку частку пожеж ми накриємо.

  • CI — bootstrap 14-денними блоками (пожежі йдуть серіями).
  • Дві цілі: день з пожежею і нове займання.
4

Результат

Прогнози 2023–2025, яких модель не бачила.

Ціль

PR-AUC підходів

Як читати: праворуч — краще; вертикальна лінія — випадковий прогноз. * — не прогноз: потрібне вчорашнє спостереження.

5

Діагностика

Щоденний ризик чемпіона й реальні пожежі

Як читати: лінія — оцінка ризику дня; вище горизонтальної лінії — 10 % найризиковіших днів («тривога»); помаранчеві точки — дні з подією.

Що модель вважає важливим

Як читати: стандартизовані коефіцієнти чемпіона; помаранчеві підвищують ризик, сірі — знижують.

Precision–recall криві

Як читати: що вище крива, то менше хибних тривог при тому самому recall. Наведіть, щоб підсвітити.

PR-AUC за тестовими роками

Як читати: 2023 — рік великої пожежі; 2024–2025 майже без пожеж, тому там метрики шумні.

Посуха → пожежі

Як читати: частка пожежних днів у кожному децилі «днів без дощу»; Д10 — найдовші посухи.

6

Що це означає

FG
FG · ризик ділянок

Де ризик згоріти найвищий?

1

Дані

Один рядок — гекс H3 ≈ 0,1 км². Ознаки — рельєф, ліс, люди. Лейбл — чи згорів гекс.

Лейбли та ознаки гексів

Як читати: лейбли — помаранчевий «так», сірий «ні»; ознаки — від темного (мало) до світлого (багато).

2

Підходи

Довга історія багатьох подій проти однієї свіжої пожежі.

3

Як перевіряли й чим міряли

Протокол: вибір на 2022, звіт на 2023

Підхід обираємо за пожежею 2022 (валідація), звітуємо на пожежі 2023, якої не бачив жоден підхід.

  • CI — bootstrap по блоках H3 (~5 км²).
  • Стеля — що дала б модель, знаючи лейбли того ж 2023 року.
  • Стійкість — ROC-AUC з трьома визначеннями «згоріло».

Метрика: ROC-AUC (+ recall у 20 % площі)

Нам потрібне ранжування ділянок. ROC-AUC — шанс, що згоріла ділянка стоїть вище незгорілої. У 2023 році згоріло 67 % парку, тому PR-AUC майже тривіальний. Recall у 20 % площі — частка згарища в 20 % площі під посиленим наглядом.

4

Результат

Тест — пожежа 2023.

ROC-AUC підходів на 2023 році

Як читати: праворуч — краще; ліва лінія — випадково (0,5), права — стеля. Форма точки — на чому вчився підхід.

5

Діагностика

ключовий графік

Валідація 2022 vs тест 2023

Як читати: якби вибір на 2022 році переносився на 2023-й, точки лежали б уздовж діагоналі. Вони лежать майже навпаки.

Чому валідація 2022 обманює

Як читати: частка гексів, що згоріли у 2023 році, серед тих, що вже горіли у 2022-му, і серед решти.

Що важливо для моделі на історії (LightGBM)

Як читати: на скільки падає ROC-AUC, якщо перемішати ознаку (permutation importance, тест 2023).

Стійкість до визначення «згоріло»

Як читати: синє — краще за випадкове, червоне — гірше. Добрий підхід синій у всіх трьох колонках.

Прогноз, факт і ризик 2026

6

Що це означає

FG
FG · поширення пожежі

Пожежа вже йде: куди вона піде завтра й післязавтра?

1

Дані

Один рядок — (пожежа, момент t, ще не згорілий гекс H3 ≈ 0,74 км² у радіусі 12 км). Лейбл — чи згорить до t+24 / t+48 год.

28 багатоденних пожеж Греції 2021–2025

Як читати: кожне коло — пожежа з FIRMS VIIRS (розмір — кількість детекцій); помаранчева з рамкою — Еврос-2023, недоторканий тест. Решта 27 — для навчання й вибору підходу.

2

Підходи

Від «росте рівномірно» до моделі, що знає вітер, схил і паливо.

3

Як перевіряли й чим міряли

Протокол: leave-one-fire-out + недоторканий тест

Підхід порівнюємо так: вчимося на 26 пожежах, перевіряємо на 27-й — і так для кожної (параметри правил і поріг моделі підбираються лише на навчальних пожежах). Переможця обираємо за цим. Еврос-2023 не бере участі ні в чому — на ньому лише фінальний звіт.

  • Крок — доба (VIIRS пролітає ~2 рази на добу); «зріз» о 13:00 UTC.
  • Вітер на наступну добу — фактичний NASA POWER як заміна прогнозу погоди.
  • CI — bootstrap по пожежах.

Метрика: F1 нової площі (+ похибка фронту)

Рахуємо лише нову площу — те, що ще не горіло. F1 карає і за перебільшення, і за пропуски; macro — кожна пожежа має однакову вагу (інакше Еврос заглушив би все). Похибка фронту — медіана відстані від реально згорілого гекса до найближчого спрогнозованого: «на скільки метрів промахнулись».

4

Результат

Вибір — 27 пожеж (LOFO), тест — Еврос-2023.

Горизонт

F1 нової площі

Як читати: праворуч — краще; кільце з вусом — середнє по 27 пожежах з 95 % CI (на цьому обирали); заповнена точка — тест на Евросі-2023.

5

Діагностика

Еврос-2023 день за днем: прогноз проти того, що сталося

Як читати: темні гекси — вже згоріло на момент прогнозу; помаранчевий — згоріло й спрогнозовано; світлий — спрогнозовано, але не згоріло; сірий — згоріло, але пропущено. Стрілка — куди дме вітер.

уже згорілозбігперебільшенняпропуск

Якість прогнозу по днях пожежі

Як читати: F1 нової площі на кожну добу Евросу-2023; у дні вибухового росту прогнозувати найважче.

На яких пожежах підхід працює

Як читати: кожна точка — одна з 27 пожеж (LOFO); праворуч — краще. Чемпіон проти базової «рівномірний ріст».

6

Що це означає

AR
AR · ареалізація появи

Де щось регулярно повторюється?

1

Дані

Реальні пожежі Греції 2021–2025: точки FIRMS, зведені правилом FM. Один рядок це пожежа. Потік демо тут не використовується.

2

Підходи

Чи дають контури ареалів більше за простіші способи показати гарячі місця.

3

Як перевіряли й чим міряли

Протокол: лише минуле

Ареали року Y будуються лише з пожеж попередніх років. Радіус і вікно обрано на 2023, тест на 2024 і 2025. Усі методи отримують рівно однакову площу: 1,5 % гексів H3.

Метрика: PAI

У скільки разів щільність нових пожеж у зоні вища за середню. 1 це випадковий вибір. Також частка нових пожеж, що потрапила в зону.

4

Як результат залежить від параметрів

Обрано на 2023, клітинка з рамкою.

PAI на валідації: радіус × вікно

Як читати: рядки це радіус DBSCAN, стовпці кількість минулих років. Що яскравіше, то краще.

5

Результат

Тест 2024–2025.

PAI на тесті, однакова площа 1,5 %

Як читати: праворуч краще, 1 це випадковий вибір; вус це 95 % довірчий інтервал.

Усі числа

6

Діагностика

Ареали за всі роки

Як читати: контур це ареал, що яскравіше, то більше пожеж у ньому. Показано 60 найбільших.

Ареали демо

Це потік демо, здебільшого симульований. Тут метрика не рахується.

7

Що це означає

LD
LD · детекція вирубки

Чи можна знайти вирубку там, де її майже немає?

1

Дані

Один рядок — піксель лісу 30 м за рік. Лейбл — втрата крон за Hansen. Згарища викинуті за EFFIS і FIRMS.

2

Підходи

Від одного знімка до моделі на всіх смугах.

3

Як перевіряли й чим міряли

Протокол: половина парку × час

Парк один, тому вибір на одній половині за 2018–2021, тест на другій за 2022–2024, і навпаки.

Метрика: F1 площі

Баланс точності й повноти за площею. У тесті лише 14 га втрат, тому інтервали широкі.

4

Результат

Правило обирає варіант за калібруванням, а не за тестом.

F1 площі на тесті

Як читати: праворуч краще; вус це 95 % довірчий інтервал; жирний рядок обрано на калібруванні.

5

Діагностика

Точність і повнота

Кандидати за 2024 рік

Як читати: зелені є й у NASA OPERA, жовті лише в Sentinel-2. Це кандидати, а не доведені рубки.

6

Що це означає

LM
LM · матчинг вирубки

Чи допомагає збіг двох супутникових джерел?

1

Дані

Ті самі пікселі, що в LD, плюс знімки стану NASA OPERA DIST-ALERT кожні 10 днів.

2

Підходи

Два джерела окремо й два способи їх поєднати.

3

Як перевіряли й чим міряли

Протокол

Sentinel-2 навчається на 2018–2021, правило OPERA обирається на 2023, тест на 2024 на іншій половині парку.

Метрика: точність

Скільки з позначеного справді втрата лісу. Це важливо для виїзду: помилковий кандидат коштує дня роботи.

4

Результат

Лише 5,5 га втрат у тесті.

Точність на тесті

Як читати: праворуч краще; показано точність, бо заради неї поєднують джерела.

6

Що це означає

LG
LG · розвиток вирубки

Де рубатимуть наступного року?

1

Дані

Один рядок — гекс H3 × рік. Лейбл — щонайменше 0,5 га нової втрати без згарищ.

2

Підходи

Від випадкового вибору до моделі на історії й ландшафті.

3

Як перевіряли й чим міряли

Протокол: лише минуле

Прогноз на рік Y вчиться на роках до Y. Вибір 2018–2020, тест 2021–2024.

Метрика: частка втрат у 11 % гексів

Якщо перевірити лише 11 % парку, скільки нових втрат туди потрапить. Випадковий вибір дає близько 11 %, на цьому тесті 6 %.

4

Результат

У тесті 39 га втрат.

Частка втрат у 11 % гексів

Як читати: праворуч краще; вус це 95 % довірчий інтервал; жирний рядок обрано на валідації.

5

Діагностика

По роках

Як читати: прогноз на рік Y вчиться лише на роках до Y. 2018–2020 валідація, далі тест.

Частка втрат залежно від площі, яку перевіряємо

Як читати: замість одного порога 11 % показано всю криву: скільки нових втрат потрапляє в найризиковіші 2, 5, 11, 20, 30 і 50 % гексів. Тест 2021–2024.

Прогноз на 2025 рік

Як читати: фіолетові гекси це 11 % найризикованіших.

6

Що це означає

SD
SD · детекція зі супутника

Чи видно машини згори, і за якої роздільності?

1

Дані

Відкритий набір VEDAI з рамками машин, розмитий до трьох роздільностей.

2

Підходи

Дві готові моделі без донавчання.

3

Як перевіряли й чим міряли

Протокол: 50/50 кадрів

Поріг впевненості обрано на половині кадрів, числа на другій. Модель не вчилась на цих даних.

Метрика: AP@0.5

Якість рамок незалежно від порога: рамка правильна, якщо перекривається з еталоном щонайменше наполовину.

4

Результат

12,5 → 50 см це симуляція зниження якості.

AP@0.5 за роздільністю

Як читати: праворуч краще; вус це 95 % довірчий інтервал.

Усі числа

Як читати: поріг обрано на половині кадрів, решта чисел на другій половині; інтервал 95 %.

Скільки пікселів займає легкове авто

Як читати: що менше пікселів, то гірше. Перехід від 25 до 50 см робить авто вдвічі меншим, а AP падає більш ніж удвічі.

5

Відкриті знімки парку

Яку роздільність дають відкриті знімки Esri парку і що з неї видно. Розмітки для парку немає, тому точності тут немає.

Чому людей не видно із супутника

Як читати: скільки пікселів займає об'єкт на відкритому знімку Esri (0,45 м/піксель); детекторам потрібно ~10 px.

6

Що це означає

CD
CD · детекція з фото й відео

Наскільки точні детектори камер і фотопасток?

1

Дані

D-Fire: дим і вогонь із рамками. ENA24: кадри фотопасток. Це фото, відео в еталоні немає.

2

Підходи

GroundingDINO за текстовим запитом, без донавчання.

3

Як перевіряли й чим міряли

Протокол: 50/50 кадрів

Поріг на половині кадрів, числа на другій. На ENA24 лише 22 кадри з машинами.

Метрики

AP@0.5 там, де є рамки. На ENA24 лише лейбл кадру, тому ROC-AUC, повнота й частка хибних кадрів.

4

Результат

AP@0.5, D-Fire

Як читати: праворуч краще; вус це 95 % довірчий інтервал.

Усі числа

5

Перевірка на відео й фото

Запис пожежі в Карсті 2022 року і 8 фото з парку. Розмітки немає, тому це не точність, а звірка з порогом, обраним на еталоні.

YOLO-World на фото з парку

Як читати: рамки — знайдені об'єкти з впевненістю; зелений — людина, світлий — авто.

6

Що це означає

TP
TP · де ставити вежі

Чи побачили б вежі майбутні пожежі?

1

Дані

Точки NASA FIRMS 2021–2025 у парку, зведені в пожежі правилом FM (2 км, 5 діб). Рельєф SRTM дає, що видно з кожного місця. Вежа: огляд 360°, щогла 15 м, 10 км, дим 20 м над землею.

2

Підходи

Три бази й жадібне покриття з різною вагою: кожна наступна вежа додає найбільше нового. Параметр: скільки найвищих місць розглядаємо, 30, 100 чи 300.

3

Як перевіряли й чим міряли

Протокол: лише минуле

Вежі ставимо за тим, що відомо до року Y. Валідація: дані до 2022, пожежі 2022. Тест: дані до 2023, пожежі 2023–2025. Вибір лише на валідації, а за рівності береться простіший варіант. Вага з вкладки «Ризик» вчилась на пожежі 2022, тому на валідації її немає: вона лише для порівняння на тесті.

Метрики

Перша точка: частка пожеж, чия перша точка FIRMS у полі зору веж. Це головна метрика. Хоч одна точка: пожежу бачили бодай колись. Згорілі гекси: частка гексів пожежі 2023 у полі зору. Інтервал це бутстреп за пожежами.

4

Результат

Усі рядки цього графіка на 300 місцях-кандидатах, тест 2023–2025.

Яку частку пожеж бачать вежі

Як читати: праворуч краще; вус це 95 % інтервал. Вага ризику не проходила валідацію, її показано для порівняння. Усі варіанти «покриття» дали ті самі три місця, тому їхні стовпчики однакові.

Підібраний параметр: скільки місць розглядати

Як читати: кожна клітинка це частка пожеж за першою точкою на валідації й тесті. Вибір робився лише за числом до риски.

Усі числа, тест

5

Діагностика: скільки веж і за яких умов

Місця ті самі, змінюємо лише умови. Це не вибір, а перевірка стійкості.

По роках

Як читати: прогресивний протокол: кожен рік вежі стоять лише за даними до нього. 2022 це валідація, далі тест.

Скільки веж потрібно

Як читати: кожна нова вежа додана жадібно, тест 2023–2025.

Якщо умови інші

Як читати: перший рядок це базові умови, далі змінюється одне число.

Обрані місця

Як читати: остаточні вежі за всіма даними до 2026 року. Кольорові гекси видно з веж, колір це рівні ризику. Кружок це перша точка пожежі 2021–2025: світлий у полі зору, сірий ні.

6

Що це означає

CG
CG · камера на координати

Де на мапі дим, який камера бачить у кадрі?

1

Дані

Реальні гекси парку, де NASA FIRMS бачив вогонь у 2022 і 2023 роках, і рельєф SRTM. Реальних детекцій диму з відомими координатами немає, тому напрямок із камери симульовано.

2

Підходи

Геометрія, без навчання. Спільний параметр лише один: ширина клина.

3

Як перевіряли й чим міряли

Протокол: валідація 2022, тест 2023

Для кожного гекса з пожежею беремо справжній напрямок і додаємо похибку з відхиленням 1, 3, 5 або 10°. Ширину клина (1, 2 чи 3 відхилення) обрано на 2022: найменша зона, що накриває місце щонайменше в 90 % випадків.

Метрики

Влучання: справжній гекс у зоні. Площа зони в км². Відстань від місця до центру зони.

4

Результат

Це перевірка геометрії на симуляції, а не точність детектора.

Площа зони пошуку на тесті 2023

Як читати: ліворуч краще: менша зона, менше їхати перевіряти. Перемикач змінює похибку напрямку.

Усі числа

Яку частку гексів із пожежею 2023 бачить камера

Як читати: три колонки міняють висоту щогли за дальності 10 км, дві міняють дальність за щогли 15 м. Вежі стоять там, де їх поставила гілка TP.

5

Як це виглядає: відео й мапа синхронізовані

Запустіть відео.

Камера CAM-02: детекція в кадрі → клин на мапі

Як читати: положення рамки по горизонталі дає азимут; дальність камера не знає, тому підсвічуємо весь клин — але лише ті гекси, які видно з вишки з урахуванням рельєфу. Помаранчевий — дим, зелений — люди. Відео це запис пожежі в Карсті 2022 року, показаний як демо: він ілюструє механіку, а не Дадью.

клин димуклин людейвидно з вишкиу секторі, закрито рельєфом
6

Що це означає

CP
CP · де ставити камери

Які місця на дорогах перехоплять під'їзд до рубок?

1

Дані

Дороги й стежки OpenStreetMap, нові втрати лісу Hansen поза згарищами, пожежі FIRMS. Наземна камера чи фотопастка бачить лише свій відрізок дороги, тому важить не рельєф, а те, через що проходить під'їзд.

2

Підходи

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

3

Як перевіряли й чим міряли

Протокол: лише минуле

Місця обираємо за тим, що відомо до року Y. Валідація: дані до 2018, рубки 2018–2020, пожежі 2021. Тест: дані до 2021, рубки 2021–2024, пожежі 2022–2025. Вибір лише на валідації. Перший варіант тесту, 2023–2024, мав лише 2 рубки і дав 0 %. Тому тест розширено, а той результат лишається в звіті нижче.

Метрики

Рубки: частка рубок, найкоротший під'їзд до яких проходить через камеру. Це головна метрика. Пожежі: те саме для перших точок пожеж. Ліс: частка лісових гексів за камерами. Інтервал це бутстреп за подіями.

4

Результат

Тест 2021–2024, 10 камер.

Що перехоплюють камери

Як читати: праворуч краще; вус це 95 % інтервал для рубок.

Усі числа

Як читати: вибір робився лише за валідацією.

5

Діагностика: скільки камер і чи тримається

Місця обрані, тепер перевіряємо стійкість. Це не вибір.

По роках

Як читати: прогресивний протокол: кожен рік місця обрано лише за даними до нього. 2018–2020 валідація, далі тест.

Скільки камер потрібно

Перевірка припущень

Як читати: перший рядок це основний тест, далі те, що може зменшити результат.

Обрані місця

Як читати: остаточні 10 місць за всіма даними. Жовті лінії це під'їзди, які проходять через камери. Кружок це рубка: зелений перехоплено, сірий ні.

6

Що це означає

XM
XM · крос-джерельний матчинг

Які сигнали з усіх джерел це одна подія?

1

Дані

Потік демо за 90 днів: один рядок це сигнал. Більшість сигналів симульовані.

Сигнали за джерелом

2

Підходи

Три способи зведення й сітка параметрів для найскладнішого.

3

Як перевіряли й чим міряли

Протокол: вибір на перших 60 днях

Для кожного типу події окремо: пожежа, вирубка, машини. Параметри обрано на перших 60 днях, тест на останніх 30. Перебрано 36 пар радіус × час.

Еталон і метрика

Еталон це задум симуляції: сигнали одного сценарію й одного дня вважаються однією подією. Це не спостереження. Метрика BCubed F1, як у FM.

4

Як результат залежить від параметрів

Обрано на валідації, клітинка з рамкою.

Тип події

BCubed F1 на валідації: радіус × час

Як читати: рядки це радіус, стовпці час. Що темніше, то гірше.

5

Результат

Тест, останні 30 днів.

BCubed F1 на тесті: обраний варіант проти простіших

Як читати: BCubed F1 від 0 до 1, праворуч краще. Він показує, чи стоїть кожен сигнал разом зі своєю подією і не змішаний із чужими. Тест це останні 30 днів, еталон це симуляція.

6

Що це означає

AS
AS · Які зони інтересу є в парку, якого типу, з якими рисами в часі й просторі і з яких джерел?

Які зони інтересу є в парку, якого типу, з якими рисами в часі й просторі і з яких джерел?

Ареали застосунку кажуть, що відбувається в кожній зоні: тип, чи вона активна, згасла чи відновилась, чи розширюється, як часто там щось стається й хто це бачить. Перевіряємо, наскільки ці риси збігаються із задумом сценарію.

Постановка

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
Підсумок

Що ми знаємо — і чого ще ні

Наступні кроки

  1. Більше подій. Пожежних днів і рубок у тестах замало: розширити дані на сусідні парки й більше років.
  2. Валідація на кількох сезонах і регіонах для ризику ділянок замість однієї пожежі.
  3. Розмітка кадрів реальної вежі й супутникових знімків лісу, щоб виміряти детектори в умовах парку.
  4. Прив'язка реальної вежі. Зіставити кадри з відомими точками на мапі й виміряти справжню похибку напрямку: зараз вона симульована.
  5. Журнали проїздів. Перевірити місця камер на реальних GPS-треках або записах: зараз шлях припущений.
  6. Реальний потік сигналів. Перевірити зведення на справжніх камерах і дзвінках, а надійність джерел навчити на підтвердженнях рейнджерів: зараз вона задана вручну.

Загрози валідності

  • Одна велика пожежа. Еврос-2023 домінує в згарищі, ризику ділянок, поширенні й камері на координати. Результати можуть не переноситися на дрібні пожежі.
  • Лейбли з продуктів, а не з поля. Hansen бачить загибель крон, але не її причину. EFFIS пропускає частину дрібних пожеж. FIRMS бачить вогонь лише під час прольоту супутника.
  • Мало подій у тестах. 30 пожежних днів, кілька десятків рубок і пожеж у задачах розміщення, кілька гектарів у детекції вирубки. Тому інтервали широкі, а різниці з базою часто не доведені.
  • Погода в одній точці. NASA POWER у центрі парку не бачить різниці між ділянками.
  • Сценарні припущення й симуляція. Дальність огляду й висота диму у вежах, похибка напрямку в камері, найкоротший шлях до рубки, генератор потоку сигналів. Ці задачі позначено окремо.
  • Зміни протоколу. У «Де ставити камери» перший тест мав 2 рубки, після цього всі три задачі розміщення перейшли на прогресивний протокол. У «Дні пожежі» підхід раніше обирали на тестових роках, тепер на 2022. Обидві зміни описано в блоках.
  • Детектори не донавчались. Їх перевірено на відкритих еталонах, а не на кадрах із парку.

Джерела

  1. 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.
  2. Schroeder W. et al. (2014). The New VIIRS 375 m active fire detection data product. Remote Sensing of Environment 143, 85–96. NASA FIRMS.
  3. Giglio L. et al. (2016). The collection 6 MODIS active fire detection algorithm and fire products. Remote Sensing of Environment 178, 31–41.
  4. Didan K. (2021). MODIS/Terra Vegetation Indices 16-Day L3 Global 250 m, MOD13Q1 v061. NASA LP DAAC; вирізки ORNL DAAC.
  5. NASA POWER Project, Prediction Of Worldwide Energy Resources, NASA Langley Research Center.
  6. Copernicus Sentinel-2 L2A, ESA; каталог Earth Search, Element 84.
  7. OPERA Land Surface Disturbance Alert (DIST-ALERT-HLS), NASA JPL; шар NASA GIBS.
  8. San-Miguel-Ayanz J. et al. (2012). Comprehensive monitoring of wildfires in Europe: the European Forest Fire Information System (EFFIS).
  9. Van Wagner C.E. (1987). Development and structure of the Canadian Forest Fire Weather Index System. Canadian Forestry Service, Technical Report 35.
  10. 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.
  11. Anderson H.E. (1983). Predicting wind-driven wild land fire size and shape. USDA Forest Service, Research Paper INT-305.
  12. Birant D., Kut A. (2007). ST-DBSCAN: an algorithm for clustering spatial–temporal data. Data & Knowledge Engineering 60, 208–221.
  13. Amigó E. et al. (2009). A comparison of extrinsic clustering evaluation metrics based on formal constraints. Information Retrieval 12, 461–486. BCubed.
  14. Razakarivony S., Jurie F. (2016). Vehicle detection in aerial imagery: a small target detection benchmark. JVCIR 34, 187–203. VEDAI.
  15. 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.
  16. Yousif H. et al. (2019). Dynamic programming selection of object proposals for sequence-level animal species classification in the wild. ENA24, LILA BC.
  17. Liu S. et al. (2023). Grounding DINO. arXiv:2303.05499. Cheng T. et al. (2024). YOLO-World. CVPR. Ultralytics YOLO11.
  18. AWS Terrain Tiles (SRTM). OpenStreetMap contributors, ODbL. H3, Uber Technologies.

Використання ШІ

Код, аналіз і тексти цієї сторінки підготовлено з допомогою ШІ-асистента Claude (Anthropic). Протоколи, перевірку чисел і висновки визначали автори.