Геология со спутника 🌍

Геолог Ларри Лай выложил открытый Python-пакет pyGeoMapper. По космоснимкам и рельефу он строит черновую геологическую карту, разрезы и 3D-блоковую модель участка.

➡️ Пакет выделяет характеристики рельефа по цифровой модели (ЦМР) и добавляет к ним данные Google AlphaEarth. Для каждого пикселя 10×10 м там уже посчитаны 64 признака поверхности по оптическим и радарным снимкам за 2017–2024 годы. Геолог размечает несколько эталонных участков, классификатор обучается на этих данных и переносит разметку на остальную территорию.

Будет полезен, чтобы:
1️⃣Получить черновую карту малоизученной территории, когда есть только старые отчеты.

2️⃣Отследить деформации поверхности методом радарной интерферометрии InSAR, он тоже есть в пакете. На проекте In Salah в Алжире так увидели подъем земли примерно на 5 мм в год над скважинами закачки CO₂.

3️⃣Следить за оползнями и смещениями грунта по снимкам Sentinel-2 и Landsat.

Ставится через pip, нужен Python 3.10 или новее. Лицензия AGPL-3.0: если встраивать код в свой продукт, его исходники придется открыть.

⏮Код и примеры с открытыми данными лежат на GitHub: оползни, землетрясения, ледники, пожары и паводки ⏭

Пока автор сам называет пакет учебным экспериментом. Точность карты зависит от разметки. Если обучить классификатор на трех случайных эталонных участках, он нарисует уверенную и красивую ерунду.

Видео: Ларри Лай

✈️ Telegram | 🌐 Website