Всем привет! Я новичок на форуме, так что если я сделал что-то неверное, пожалуйста, не ругайтесь на меня
Мне нужно построить вертикальные смещения Земли на большую территорию за столько лет, за сколько это возможно. Я никогда не делал этого раньше, поэтому я прочёл множество руководств с различными советами и алгоритмами действий. Работал я в 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, я прав? Если нет, то могу ли я их получить как-то иначе?
Спасибо заранее за ваше внимание и ответы!
Правильный способ для получения вертикальных смещений Земли
-
- Новоприбывший
- Сообщения: 7
- Зарегистрирован: 30 янв 2019, 15:16
- Репутация: 0
- Откуда: Тюмень
-
- Новоприбывший
- Сообщения: 1
- Зарегистрирован: 10 окт 2020, 13:32
- Репутация: 0
- Откуда: RF
Re: Правильный способ для получения вертикальных смещений Земли
актуальная тема, тоже ее изучаю. если еще актуально - дай знать. есть кое-какие мысли
-
- Гуру
- Сообщения: 4205
- Зарегистрирован: 10 апр 2006, 22:34
- Репутация: 433
- Откуда: Париж
Re: Правильный способ для получения вертикальных смещений Земли
Я в вертикальном смещении и прочей SAR не силен, но одно заню точно - как только вы пересчитали исходные данные с целью смещения - вы внесли ошибку в пределах вашего алгоритма пересчета величин. Вы некоторым образом усреднили показания по 4-м, 9-ти или сколько там точек соседей использовал алгоритм кореегистрации для привязки снимков друг к другу.SpacePeacemaker писал(а): ↑09 апр 2019, 10:18последовательность:
- выбрал SLC-снимки в режиме IW спутника Sentinel-1B с показателями Bperp меньше 50 метров и Btemp 12 дней.
- корегистрация
Сшивать данные можно только получения конечного результата. И только ясно понимаю, что именно алгоримт сшивки с ними делает. Мы намечаете трансформацию ДО обработки, если это возможно только по исходным данным, выполняете расчеты, и из предположения, что эти ваши расчеты не меняют пространственную привязку, применяете трансформацию к полученным результатам.
Есть заявленные алгоритмы корегистрации, которые сохраняют исходные спектральные свойства, но их точность не доказана практикой.
Кто сейчас на конференции
Сейчас этот форум просматривают: нет зарегистрированных пользователей и 4 гостя