Правильный способ для получения вертикальных смещений Земли

Не знаете, где задать вопрос? Задавайте здесь.
Ответить
SpacePeacemaker
Новоприбывший
Сообщения: 7
Зарегистрирован: 30 янв 2019, 15:16
Репутация: 0
Откуда: Тюмень

Правильный способ для получения вертикальных смещений Земли

Сообщение SpacePeacemaker » 09 апр 2019, 10:18

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

Мне нужно построить вертикальные смещения Земли на большую территорию за столько лет, за сколько это возможно. Я никогда не делал этого раньше, поэтому я прочёл множество руководств с различными советами и алгоритмами действий. Работал я в 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, я прав? Если нет, то могу ли я их получить как-то иначе?

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

el_gyrgeno
Новоприбывший
Сообщения: 1
Зарегистрирован: 10 окт 2020, 13:32
Репутация: 0
Откуда: RF

Re: Правильный способ для получения вертикальных смещений Земли

Сообщение el_gyrgeno » 10 окт 2020, 13:40

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

Boris
Гуру
Сообщения: 4205
Зарегистрирован: 10 апр 2006, 22:34
Репутация: 433
Откуда: Париж

Re: Правильный способ для получения вертикальных смещений Земли

Сообщение Boris » 16 окт 2020, 19:23

SpacePeacemaker писал(а):
09 апр 2019, 10:18
последовательность:

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

Ответить

Вернуться в «Я новичок!»

Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей и 4 гостя