Landsat-8, создание карты региона
-
- Активный участник
- Сообщения: 104
- Зарегистрирован: 06 окт 2012, 15:35
- Репутация: 0
- Контактная информация:
Landsat-8, создание карты региона
Добрый день! У есть задача -- регион РФ, необходимо сделать карту его из снимков Landsat-8.
Вопросы:
1. Можно ли как-то автоматически скачать последние снимки landsat из открытых баз данных (указать необходимую территорию и как-то отфильтровать по облачности и времени съёмки)?
2. Можно ли пакетно склеить данные при помощи, например GDAL, просто указать папку с данными и чтобы в итоге получился один большой общий растр. Затем планирую поделить его на тайлы и приделать веб-интерфейс для всей карты.
Заранее большое спасибо за помощь!
Вопросы:
1. Можно ли как-то автоматически скачать последние снимки landsat из открытых баз данных (указать необходимую территорию и как-то отфильтровать по облачности и времени съёмки)?
2. Можно ли пакетно склеить данные при помощи, например GDAL, просто указать папку с данными и чтобы в итоге получился один большой общий растр. Затем планирую поделить его на тайлы и приделать веб-интерфейс для всей карты.
Заранее большое спасибо за помощь!
-
- Гуру
- Сообщения: 920
- Зарегистрирован: 30 дек 2008, 14:11
- Репутация: 236
- Откуда: Ханты-Мансийск
- Контактная информация:
Re: Landsat-8, создание карты региона
Видимо, вы имеете в виду не карту, а мозаику снимков? Карта в обычном понимании, это что-то имеющее смысловые категории, классы объектов, например.Rumato писал(а):необходимо сделать карту
Можно, причем разными способами:Rumato писал(а):автоматически скачать последние снимки 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 критично.
Можно, но чтобы добиться приемлемого результата придется повозиться: на простой (и даже не очень простой) мозаике будут видны границы между снимками. Потому что каждый снимок Landsat -- это уникальное произведение (сезон, атмосферные явления и пр.) и найти два похожих в пространстве-времени весьма затруднительно. Т.е. снимки нужно будет как-то друг к другу подгонять, а это уже даже не искусство, а скорее, магия.Rumato писал(а):Можно ли пакетно склеить данные при помощи, например GDAL
-
- Активный участник
- Сообщения: 158
- Зарегистрирован: 11 сен 2011, 14:11
- Репутация: 96
- Откуда: Москва
- Контактная информация:
Re: Landsat-8, создание карты региона
https://libra.developmentseed.org/
gdal-merge
Если умеете программировать на питоне, то https://github.com/mapbox/rasterio с https://github.com/mapbox/rio-color
gdal-merge
Если умеете программировать на питоне, то https://github.com/mapbox/rasterio с https://github.com/mapbox/rio-color
-
- Активный участник
- Сообщения: 104
- Зарегистрирован: 06 окт 2012, 15:35
- Репутация: 0
- Контактная информация:
Re: Landsat-8, создание карты региона
Большое спасибо за ответы! Как всё сделаю обязательно отпишусь чем закончилось.
-
- Активный участник
- Сообщения: 104
- Зарегистрирован: 06 окт 2012, 15:35
- Репутация: 0
- Контактная информация:
Re: Landsat-8, создание карты региона
добрый день! впринципе в основном со всем разобрался, но вот когда делаю мозайку из снимков, то получается какая-то несуразица.
Привожу скрины как и что делаю
В чём может быть проблема?
Привожу скрины как и что делаю
В чём может быть проблема?
- nadiopt
- Гуру
- Сообщения: 4744
- Зарегистрирован: 29 янв 2009, 13:27
- Репутация: 495
- Ваше звание: хрюкало
- Откуда: Санкт-Петербург
- Контактная информация:
Re: Landsat-8, создание карты региона
значит, не разобрались, снимки надо поклеить поканально
ин гроссен фамилен нихт клювен клац клац
- Игорь Белов
- Гуру
- Сообщения: 2231
- Зарегистрирован: 04 янв 2011, 22:00
- Репутация: 1504
- Откуда: Казань
Re: Landsat-8, создание карты региона
Чтобы не думать о прозрачности пустых полей, всегда определяю 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
The purpose of computing is insight, not numbers
-
- Активный участник
- Сообщения: 104
- Зарегистрирован: 06 окт 2012, 15:35
- Репутация: 0
- Контактная информация:
Re: Landsat-8, создание карты региона
nadiopt, то есть склеиваются уже композиты? по одному каналу не склеивают?
Игорь Белов, нет, нужны композиты, я вроде понял в чём моя проблема.
Игорь Белов, нет, нужны композиты, я вроде понял в чём моя проблема.
- nadiopt
- Гуру
- Сообщения: 4744
- Зарегистрирован: 29 янв 2009, 13:27
- Репутация: 495
- Ваше звание: хрюкало
- Откуда: Санкт-Петербург
- Контактная информация:
Re: Landsat-8, создание карты региона
судя по вашей картинке, вы просто пытаетесь поклеить два канала с соседних снимков. Соберите что-нибудь типа 654 (для Л8), выставьте nodata, как Игорь советует, и там должно уже сносно получиться, близкие по дате снимки можно особо не уравнивать.
ин гроссен фамилен нихт клювен клац клац
Кто сейчас на конференции
Сейчас этот форум просматривают: нет зарегистрированных пользователей и 32 гостя