Landsat-8, создание карты региона

[Rumato]

Добрый день! У есть задача -- регион РФ, необходимо сделать карту его из снимков Landsat-8.

Вопросы:
1. Можно ли как-то автоматически скачать последние снимки landsat из открытых баз данных (указать необходимую территорию и как-то отфильтровать по облачности и времени съёмки)?

2. Можно ли пакетно склеить данные при помощи, например GDAL, просто указать папку с данными и чтобы в итоге получился один большой общий растр. Затем планирую поделить его на тайлы и приделать веб-интерфейс для всей карты.


Заранее большое спасибо за помощь!

[bolotoved]

необходимо сделать карту

Видимо, вы имеете в виду не карту, а мозаику снимков? Карта в обычном понимании, это что-то имеющее смысловые категории, классы объектов, например.

автоматически скачать последние снимки landsat

Можно, причем разными способами:
1. Стандартный способ: Bulk Download Application (BDA) для EarthExplorer + https://lta.cr.usgs.gov/node/409
2. Нестандартный способ: Google gsutil. Оказывается, Google для своих целей зеркалирует архив съемки Landsat (Подробнее: http://fernerkundung.github.io/landsat- ... om-google/ ). Чтобы получить доступ к этому архиву нужно установить утилиту gsutil(доступно для Windows, Linux и Mac).
3. Через Amazon Web Services - там можно по отдельным каналам, что для L-8 критично.

Можно ли пакетно склеить данные при помощи, например GDAL

Можно, но чтобы добиться приемлемого результата придется повозиться: на простой (и даже не очень простой) мозаике будут видны границы между снимками. Потому что каждый снимок Landsat -- это уникальное произведение (сезон, атмосферные явления и пр.) и найти два похожих в пространстве-времени весьма затруднительно. Т.е. снимки нужно будет как-то друг к другу подгонять, а это уже даже не искусство, а скорее, магия.

[Zverik]

https://libra.developmentseed.org/
gdal-merge
Если умеете программировать на питоне, то https://github.com/mapbox/rasterio с https://github.com/mapbox/rio-color

[Rumato]

Большое спасибо за ответы! Как всё сделаю обязательно отпишусь чем закончилось.

[Rumato]

добрый день! впринципе в основном со всем разобрался, но вот когда делаю мозайку из снимков, то получается какая-то несуразица.

Привожу скрины как и что делаю




В чём может быть проблема?

[nadiopt]

значит, не разобрались, снимки надо поклеить поканально

[Игорь Белов]

Чтобы не думать о прозрачности пустых полей, всегда определяю nodata_value в самом начале:

Код: Выделить всё

#!/bin/sh
# USAGE: . get_landsat.sh <LANDSAT_SCENE_ID>
scene=${1}

for (( band=1; band<=8; band++ ))
do
    gdal_translate -a_nodata 0 -co COMPRESS=LZW ${scene}_B${band}.TIF b${band}.tif
done

Просто интересно, Вам действительно нужны монохромные мозаики?

[Rumato]

nadiopt, то есть склеиваются уже композиты? по одному каналу не склеивают?

Игорь Белов, нет, нужны композиты, я вроде понял в чём моя проблема.

[nadiopt]

судя по вашей картинке, вы просто пытаетесь поклеить два канала с соседних снимков. Соберите что-нибудь типа 654 (для Л8), выставьте nodata, как Игорь советует, и там должно уже сносно получиться, близкие по дате снимки можно особо не уравнивать.