Страница 1 из 2
Вопрос по работе в программе S1Toolbox
Добавлено: 10 янв 2017, 11:48
Василий Почвовед
Уважаемые Дамы и Господа, добрый день!
Позвольте задать следующий вопрос:
На основании данных Sentinel 1A (спутник оснащён радиолокационной системой с синтезированной апертурой) пытаемся получить изображения, на которых были бы видны заполненные водой микропонижения на полях.
В этих целях был использован следующий алгоритм (
http://wiki.gis-lab.info/w/Sentinel-1_T ... _Tutorial)).
Используемое ПО: S1Toolbox.
Вкратце вот он:
1) после открытия исходных изображений (Intensity VH), была произведена Калибровка данных (Radar --- Radiometric --- Calibrate);
2) Была выполнена процедура некогерентного накопления (Multilooking);
3) Затем, была произведена фильтрация спекл-шума (Speckle-filtering);
4) Была произведена корректировка по местности (Terrain Correction --- Range-Doppler Terrain Correction).
В итоге получились следующие изображения (файл во вложении).
Получилось, что микропонижения, заполненные водой не видны.
Вопросы состоят в следующем:
1) Возможно ли на радарном снимке чётко увидеть эти микропонижения рельефа, заполненные водой?
2) Если возможно, то как это можно сделать? Может, мы что-то недоучли в процессе работы?
Заранее огромнейшее спасибо за ответ!
С уважением,
Василий.
Re: Вопрос по работе в программе S1Toolbox
Добавлено: 10 янв 2017, 12:44
bolotoved
А какого размера эти микропонижения? Если они меньше пространственного разрешения снимка, то скорее всего их не видно.
В какой период времени эти микропонижения заполнены водой, совпадает ли это время с периодом в который был получен этот конкретный снимок?
Re: Вопрос по работе в программе S1Toolbox
Добавлено: 10 янв 2017, 13:14
Василий Почвовед
Здравствуйте! Спасибо большое за ответ!
Дело в том, что за даты 15.07.2016 (Sentinel 2) и 26.06.2016 (Landsat 8 OLI), на изображениях эти микропонижения и влага в них видны. Территория та же самая. Разрешение...у Landsat 8 OLI в пикселе 30 м. У Sentinel 2 изображение получилось более чётким.
Высылаю во вложении иллюстрацию.
Как раз на снимках (Сверху - Sentinel 2, а снизу Landsat 8 OLI) видны насыщенно-тёмные пятна и точки. Это и есть те самые микропонижения с влагой.
Re: Вопрос по работе в программе S1Toolbox
Добавлено: 10 янв 2017, 14:21
bolotoved
Дело в том, что за даты 15.07.2016 (Sentinel 2) и 26.06.2016 (Landsat 8 OLI), на изображениях эти микропонижения и влага в них видны.
Обратите, внимание, что ваш объект очень динамичен: на снимке за 26 июня видны множество мелких сильно обводненных (возможно даже с открытой водой) участков, а спустя 3 недели это уже какие-то болотины без открытой воды.
Нужно уточнить за какую дату у вас радарный Sentinel-1.
Меня смущает, что NDVI вымочек на ландсате (средний кадр в нижнем ряду) очень высокий. Увеличение интенсивности зеленого -- это у вас возрастание индекса NDVI? Также смущает инвертирование NDVI на ландсате по сравнению с сентинель-2: посмотите, треугольный объект на сентинеле светлый, а окружение его темно-зеленое, а на ландсате - наоборот. Это странно.
Re: Вопрос по работе в программе S1Toolbox
Добавлено: 10 янв 2017, 14:35
Василий Почвовед
Хорошо, спасибо большое, на изображение с NDVI посмотрю ещё раз внимательнее.
Насчёт даты Sentinel 1 дата 04.07.2016. То есть, между снимками Landsat 8 OLI и Sentinel 2.
*Уточнил. NDVI на Sentinel 2 и Landsat 8 OLI отражены одинаково в плане возрастания.
Re: Вопрос по работе в программе S1Toolbox
Добавлено: 10 янв 2017, 15:11
bolotoved
Уточнил. NDVI на Sentinel 2 и Landsat 8 OLI отражены одинаково в плане возрастания.
1. В таком случае, NDVI рассчитанный по Ландсат, скорее всего, посчитан неправильно. Об этом же говорит высокий NDVI на месте облачка.
2. Вырезка снимка Sentinel-1 слишком маленькая, чтобы сориентироваться, что чему соответсвует. Сделайте побольше, чтобы было видно тот же объект, что на ландсате.
Re: Вопрос по работе в программе S1Toolbox
Добавлено: 10 янв 2017, 15:27
Василий Почвовед
Сделал скриншот того снимка, где уже проведена корректировка по местности (Range-Doppler Terrain Correction). Это то же самое поле, что и на ландсате.
Re: Вопрос по работе в программе S1Toolbox
Добавлено: 10 янв 2017, 15:45
bolotoved
Что-то моя нейронная сеть не может идентифицировать тоже самое поле.
Re: Вопрос по работе в программе S1Toolbox
Добавлено: 10 янв 2017, 15:50
Василий Почвовед
Давайте попробую ещё раз:
Re: Вопрос по работе в программе S1Toolbox
Добавлено: 10 янв 2017, 15:57
artemiziy
Напишите имя исходного снимка S-1 и приложите район интереса.
По описанным вами "шагам обработки" куча вопросов, не совсем понятно, зачем сделана часть шагов.
Re: Вопрос по работе в программе S1Toolbox
Добавлено: 10 янв 2017, 16:10
Василий Почвовед
Добрый день! Имя исходного снимка: S1A_IW_GRDH_1SDV_20160704T032943_20160704T033008_011993_01281C_A9E6
Область интереса: Тамбовская область, Кирсановский район. Сейчас приложу карты и покажу, как он выглядит...
Вот...
Re: Вопрос по работе в программе S1Toolbox
Добавлено: 10 янв 2017, 16:16
bolotoved
Давайте попробую ещё раз
Не помогло. На ландсате хорошо видно сверху треугольное острие, а снизу справа - перекресток, на sentinel-1 не могу разглядеть ничего похожего.
Re: Вопрос по работе в программе S1Toolbox
Добавлено: 10 янв 2017, 16:19
Василий Почвовед
Вот, у меня получилось такое изображение...может, я что-то не понял в алгоритме, и сделал что-то не так, поэтому всё так размыто?
Re: Вопрос по работе в программе S1Toolbox
Добавлено: 10 янв 2017, 16:33
artemiziy
scihub. подвисает, как развиснет - скачаю снимок, посмотрю.
Если вы правильно записали свои "шаги обработки", то:
1. т.к. вы использовали продукт GRD, вам точно не нужно делать:
- процедура некогерентного накопления (Multilooking);
2. возможно нужно использовать другой фильтр для фильтрации спекл-шума.
но это только после выполнения пункта №1.
и почему вы использовали поляризацию VH? в VV результат должен быть более четким.
Re: Вопрос по работе в программе S1Toolbox
Добавлено: 10 янв 2017, 16:39
Василий Почвовед
Спасибо большое за ответ!
Поляризацию VV использовал потому что она была приведена в примере, насколько я помню. Радарные снимки для меня дело новое, поэтому я и пробовал сначала сделать "как по образцу". Попробую сейчас сделать VV.
Попробую без некогерентного накопления и "поиграть" с фильтрами спекл-шума.
Ещё раз огромное спасибо!