Форумы GIS-Lab.info

Всем привет! Я новичок на форуме, так что если я сделал что-то неверное, пожалуйста, не ругайтесь на меня

Мне нужно построить вертикальные смещения Земли на большую территорию за столько лет, за сколько это возможно. Я никогда не делал этого раньше, поэтому я прочёл множество руководств с различными советами и алгоритмами действий. Работал я в SNAP, так как доступного SARscape я не нашёл. Вкратце, я сделал следующую последовательность:

- выбрал SLC-снимки в режиме IW спутника Sentinel-1B с показателями Bperp меньше 50 метров и Btemp 12 дней.
- корегистрация
- формирование интерферограммы
- deburst
- topophase removal
- фильтр Гольдштейна
- развёртка фазы
- перевод фазы в смещения
- terrain correction
- экспортировал в tiff

Во всех случаях я использовал стандартные настройки, кроме этапа экспорта в развёртке - я выбрал в строках и колонках значение 20 вместо 10. Также в шаге "перевода фазы в смещения" я использовал два решения. 1) стандартная функция в SNAP; 2) я создал новый бэнд с таким выражением:

(Unw_Phase * 5.55) / (-4 * PI * cos(rad(incident_angle))).

Полученные картинки получились похожими, но с разными значениями:

Бэнд с формулой и его статистика:





Бэнд после процедуры и его статистика:





Поэтому у меня много вопросов:

1) В целом я получил изображения верно? Или я допустил ошибку? Если допустил - то где? Каков правильный способ построить вертикальные смещения?
2) Нужно ли было сделать процедуру мультилукинга перед фильтрацией?
3) Когда я соединил части (swaths) снимка через merge я получил странную картинку. Означает ли это, что я не могу соединить части изображения в одну большую карту смещений? Это вообще возможно? Или не нужно этого делать?

Бэнды с формулой и их статистика:





Бэнды после процедуры и их статистика:





4) Могу ли я посчитать горизонтальные смещения? Всё что нужно - это сделать вышеописанные шаги с поляризацией HH, я прав? Если нет, то могу ли я их получить как-то иначе?

Спасибо заранее за ваше внимание и ответы!

актуальная тема, тоже ее изучаю. если еще актуально - дай знать. есть кое-какие мысли

SpacePeacemaker писал(а): ↑

09 апр 2019, 10:18

последовательность:

- выбрал SLC-снимки в режиме IW спутника Sentinel-1B с показателями Bperp меньше 50 метров и Btemp 12 дней.
- корегистрация

Я в вертикальном смещении и прочей SAR не силен, но одно заню точно - как только вы пересчитали исходные данные с целью смещения - вы внесли ошибку в пределах вашего алгоритма пересчета величин. Вы некоторым образом усреднили показания по 4-м, 9-ти или сколько там точек соседей использовал алгоритм кореегистрации для привязки снимков друг к другу.
Сшивать данные можно только получения конечного результата. И только ясно понимаю, что именно алгоримт сшивки с ними делает. Мы намечаете трансформацию ДО обработки, если это возможно только по исходным данным, выполняете расчеты, и из предположения, что эти ваши расчеты не меняют пространственную привязку, применяете трансформацию к полученным результатам.
Есть заявленные алгоритмы корегистрации, которые сохраняют исходные спектральные свойства, но их точность не доказана практикой.