Вопрос по работе в программе S1Toolbox

[Василий Почвовед]

Уважаемые Дамы и Господа, добрый день!

Позвольте задать следующий вопрос:

На основании данных 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) Если возможно, то как это можно сделать? Может, мы что-то недоучли в процессе работы?

Заранее огромнейшее спасибо за ответ!

С уважением,
Василий.

[bolotoved]

А какого размера эти микропонижения? Если они меньше пространственного разрешения снимка, то скорее всего их не видно.

В какой период времени эти микропонижения заполнены водой, совпадает ли это время с периодом в который был получен этот конкретный снимок?

[Василий Почвовед]

Здравствуйте! Спасибо большое за ответ!
Дело в том, что за даты 15.07.2016 (Sentinel 2) и 26.06.2016 (Landsat 8 OLI), на изображениях эти микропонижения и влага в них видны. Территория та же самая. Разрешение...у Landsat 8 OLI в пикселе 30 м. У Sentinel 2 изображение получилось более чётким.
Высылаю во вложении иллюстрацию.
Как раз на снимках (Сверху - Sentinel 2, а снизу Landsat 8 OLI) видны насыщенно-тёмные пятна и точки. Это и есть те самые микропонижения с влагой.

[bolotoved]

Дело в том, что за даты 15.07.2016 (Sentinel 2) и 26.06.2016 (Landsat 8 OLI), на изображениях эти микропонижения и влага в них видны.

Обратите, внимание, что ваш объект очень динамичен: на снимке за 26 июня видны множество мелких сильно обводненных (возможно даже с открытой водой) участков, а спустя 3 недели это уже какие-то болотины без открытой воды. Нужно уточнить за какую дату у вас радарный Sentinel-1.

Меня смущает, что NDVI вымочек на ландсате (средний кадр в нижнем ряду) очень высокий. Увеличение интенсивности зеленого -- это у вас возрастание индекса NDVI? Также смущает инвертирование NDVI на ландсате по сравнению с сентинель-2: посмотите, треугольный объект на сентинеле светлый, а окружение его темно-зеленое, а на ландсате - наоборот. Это странно.

[Василий Почвовед]

Хорошо, спасибо большое, на изображение с NDVI посмотрю ещё раз внимательнее.

Насчёт даты Sentinel 1 дата 04.07.2016. То есть, между снимками Landsat 8 OLI и Sentinel 2.

*Уточнил. NDVI на Sentinel 2 и Landsat 8 OLI отражены одинаково в плане возрастания.

[bolotoved]

Уточнил. NDVI на Sentinel 2 и Landsat 8 OLI отражены одинаково в плане возрастания.

1. В таком случае, NDVI рассчитанный по Ландсат, скорее всего, посчитан неправильно. Об этом же говорит высокий NDVI на месте облачка.

2. Вырезка снимка Sentinel-1 слишком маленькая, чтобы сориентироваться, что чему соответсвует. Сделайте побольше, чтобы было видно тот же объект, что на ландсате.

[Василий Почвовед]

Сделал скриншот того снимка, где уже проведена корректировка по местности (Range-Doppler Terrain Correction). Это то же самое поле, что и на ландсате.

[bolotoved]

Что-то моя нейронная сеть не может идентифицировать тоже самое поле.

[Василий Почвовед]

Давайте попробую ещё раз:

[artemiziy]

Напишите имя исходного снимка S-1 и приложите район интереса.
По описанным вами "шагам обработки" куча вопросов, не совсем понятно, зачем сделана часть шагов.

[Василий Почвовед]

Добрый день! Имя исходного снимка: S1A_IW_GRDH_1SDV_20160704T032943_20160704T033008_011993_01281C_A9E6

Область интереса: Тамбовская область, Кирсановский район. Сейчас приложу карты и покажу, как он выглядит...

Вот...

[bolotoved]

Давайте попробую ещё раз

Не помогло. На ландсате хорошо видно сверху треугольное острие, а снизу справа - перекресток, на sentinel-1 не могу разглядеть ничего похожего.

[Василий Почвовед]

Вот, у меня получилось такое изображение...может, я что-то не понял в алгоритме, и сделал что-то не так, поэтому всё так размыто?

[artemiziy]

scihub. подвисает, как развиснет - скачаю снимок, посмотрю.

Если вы правильно записали свои "шаги обработки", то:
1. т.к. вы использовали продукт GRD, вам точно не нужно делать:
- процедура некогерентного накопления (Multilooking);
2. возможно нужно использовать другой фильтр для фильтрации спекл-шума.
но это только после выполнения пункта №1.

и почему вы использовали поляризацию VH? в VV результат должен быть более четким.

[Василий Почвовед]

Спасибо большое за ответ!
Поляризацию VV использовал потому что она была приведена в примере, насколько я помню. Радарные снимки для меня дело новое, поэтому я и пробовал сначала сделать "как по образцу". Попробую сейчас сделать VV.
Попробую без некогерентного накопления и "поиграть" с фильтрами спекл-шума.

Ещё раз огромное спасибо!