Не работает корегистрация в SNAP для S1B

Программы и алгоритмы для обработки данных дистанционного зондирования: ERDAS, ENVI и другие.
Iggi1981
Участник
Сообщения: 93
Зарегистрирован: 06 мар 2020, 17:07
Репутация: 19
Откуда: Москва

Не работает корегистрация в SNAP для S1B

Сообщение Iggi1981 » 02 фев 2023, 17:03

Добрый день.
Скачал два архивных снимка Сентинель-1 за 2021 год.
Хочу создать RGB-композит по инструкции:
https://dev.asf.alaska.edu/how-to/data- ... el-1-data/

Дошел до корегистрации двух снимков, в результате получаю: слой Intensity у "мастера" - есть, а слой Intensity у "слейва" - нулевой. Из-за этого и когерентность не считается.
Что я делаю не так?
Вложения
slave_histo.jpg
slave_histo.jpg (195.06 КБ) 2523 просмотра
coreg_slave.jpg
coreg_slave.jpg (171.56 КБ) 2523 просмотра
coreg_master.jpg
coreg_master.jpg (433.25 КБ) 2523 просмотра

Константин Силкин
Завсегдатай
Сообщения: 420
Зарегистрирован: 21 мар 2012, 07:37
Репутация: 67
Откуда: Воронеж

Re: Не работает корегистрация в SNAP для S1B

Сообщение Константин Силкин » 02 фев 2023, 21:25

Да, давно назрела необходимость написать на сайте толковую инструкцию по обработке радарных данных. Как там переходить к высотам, до сих пор используется какая-то самопальная линуксовая утилитка или прогресс шагнул уже вперёд?

Iggi1981
Участник
Сообщения: 93
Зарегистрирован: 06 мар 2020, 17:07
Репутация: 19
Откуда: Москва

Re: Не работает корегистрация в SNAP для S1B

Сообщение Iggi1981 » 03 фев 2023, 09:43

Константин Силкин, а зачем переходить к высотам? Я несколько лет обрабатываю снимки, правда, в другой программе, и ни разу мне это не потребовалось.

Константин Силкин
Завсегдатай
Сообщения: 420
Зарегистрирован: 21 мар 2012, 07:37
Репутация: 67
Откуда: Воронеж

Re: Не работает корегистрация в SNAP для S1B

Сообщение Константин Силкин » 03 фев 2023, 09:51

Отлично, Iggi1981, поделитесь тогда своим опытом

Iggi1981
Участник
Сообщения: 93
Зарегистрирован: 06 мар 2020, 17:07
Репутация: 19
Откуда: Москва

Re: Не работает корегистрация в SNAP для S1B

Сообщение Iggi1981 » 03 фев 2023, 10:02

Константин Силкин, я не уверен, что готов претендовать на эксперта в этой области :D
Мы решаем задачи мониторинга оседаний земной поверхности вследствие техногенных либо природных процессов (откачка газа/нефти, подземная добыча полезных ископаемых, оползни, землетрясения).
SNAP я открыл вчера первый раз в жизни, так то я в Гамме работаю. Я ниже как-то спрашивал, хотел создать RGB-композит для геолога. Так красиво - не получается.

Константин Силкин
Завсегдатай
Сообщения: 420
Зарегистрирован: 21 мар 2012, 07:37
Репутация: 67
Откуда: Воронеж

Re: Не работает корегистрация в SNAP для S1B

Сообщение Константин Силкин » 03 фев 2023, 11:10

Iggi1981, я и не предлагаю вам писать библию спутниковой радарной съёмки. Однако раз у вас есть практически опыт, а здесь он востребован, то почему бы и не поделиться им? Форум ведь для этого и создан. Вы поможете кому-то, а они вам могут быть полезны. Потом подключатся другие, может даже коллективный проект возникнет. Так вместе и будем двигаться дальше.

Iggi1981
Участник
Сообщения: 93
Зарегистрирован: 06 мар 2020, 17:07
Репутация: 19
Откуда: Москва

Re: Не работает корегистрация в SNAP для S1B

Сообщение Iggi1981 » 03 фев 2023, 11:18

Константин Силкин, ну, я же участник форума, если возникают вопросы, на которые я могу ответить - я отвечаю.
Проблема в том, что Гамма - это коммерческий продукт, а бесплатный SNAP - я в нем ни разу не эксперт. Если нужна последовательность обработки - так это не секрет, она стандартная.
Спойлер
1. Распаковка данных и вырезание области счета (кроп).
2. Интерполяция (оверсамплинг) по дальности, в 2 раза, мультилукинг проводится с учетом оверсамплинга.
3. Подготовка ЦМР (цифровой модели рельефа) для области счета и таблицы пересчета радарных координат в наземные и обратно.
4. Корегистрация второго (slave) снимка на первый (master).
5. Расчет фазы от рельефа дифференциальной интерферограммы.
6. Вычисление дифференциальной интерферограммы.
7. Мультилукинг корегистрированного снимка.
8. Расчет когерентности по интерферограмме и двум снимкам после мультилукинга.
9. Используя полученную когерентность, преобразуем ее в растровый формат.
10. Фильтрация интерферограммы для уменьшения шума.
11. Выполняется развертка интерферограммы, используя маску растра, методом «Minimum cost flow».
12. Устранение возможного тренда и перевод развернутой фазы в смещения земной поверхности (в LOS или в вертикальные при условии "малости" горизонтальных смещений).
13. Перевод поля смещений из радарных координат в наземные, создание файла geotiff для работы в графических пакетах.

gamm
Гуру
Сообщения: 4046
Зарегистрирован: 15 окт 2010, 08:33
Репутация: 1050
Ваше звание: программист
Откуда: Казань

Re: Не работает корегистрация в SNAP для S1B

Сообщение gamm » 03 фев 2023, 11:52

пожалуйста покажите (скопируйте под спойлер) из метаданных обоих снимков mbr (bounding box), такое ощущение, что снимки не перекрываются.

Iggi1981
Участник
Сообщения: 93
Зарегистрирован: 06 мар 2020, 17:07
Репутация: 19
Откуда: Москва

Re: Не работает корегистрация в SNAP для S1B

Сообщение Iggi1981 » 03 фев 2023, 12:05

gamm, не очень понял, те ли метаданные нужны.
Да, пробовал корегистрировать и первый на второй, и второй на первый, результат один, слейв - пустой.
Спойлер
Orbit_State_Vectors
SRGR_Coefficients
Doppler_Centroid_Coefficients
Band_IW1_VH
Band_IW1_VV
Band_IW2_VH
Band_IW2_VV
Band_IW3_VH
Band_IW3_VV
PRODUCT S1B_IW_SLC__1SDV_20210707T205535_20210707T205602_027693_034E11_57A8 ascii Product name
PRODUCT_TYPE SLC ascii Product type
SPH_DESCRIPTOR Sentinel-1 IW Level-1 SLC Product ascii Description
MISSION SENTINEL-1B ascii Satellite mission
ACQUISITION_MODE IW ascii Acquisition mode
antenna_pointing right ascii Right or left facing
BEAMS - ascii Beams used
SWATH - ascii Swath name
PROC_TIME 08-JUL-2021 04:10:48.036186 uint32 utc Processed time
Processing_system_identifier ESA Sentinel-1 IPF 003.31 ascii Processing system identifier
orbit_cycle 165 int32 Cycle
REL_ORBIT 17 int32 Track
ABS_ORBIT 27693 int32 Orbit
STATE_VECTOR_TIME 07-JUL-2021 20:54:37.000000 uint32 utc Time of orbit state vector
VECTOR_SOURCE - ascii State vector source
incidence_near 99999.0 float64 deg
incidence_far 99999.0 float64 deg
slice_num 4 int32 Slice number
data_take_id 216593 int32 Data take identifier
first_line_time 07-JUL-2021 20:55:35.976569 uint32 utc First zero doppler azimuth time
last_line_time 07-JUL-2021 20:56:02.935190 uint32 utc Last zero doppler azimuth time
first_near_lat 49.547658423351145 float64 deg
first_near_long 140.72378615889232 float64 deg
first_far_lat 49.99187674928981 float64 deg
first_far_long 137.26800822147013 float64 deg
last_near_lat 48.043412424264396 float64 deg
last_near_long 140.26552747449003 float64 deg
last_far_lat 48.484625071494555 float64 deg
last_far_long 136.91040565568375 float64 deg
PASS DESCENDING ascii ASCENDING or DESCENDING
SAMPLE_TYPE COMPLEX ascii DETECTED or COMPLEX
mds1_tx_rx_polar VH ascii Polarization
mds2_tx_rx_polar VV ascii Polarization
mds3_tx_rx_polar - ascii Polarization
mds4_tx_rx_polar - ascii Polarization
polsar_data 0 uint8 flag Polarimetric Matrix
algorithm - ascii Processing algorithm
azimuth_looks 1.0 float64
range_looks 1.0 float64
range_spacing 2.3295619999999997 float64 m Range sample spacing
azimuth_spacing 13.901110000000001 float64 m Azimuth sample spacing
range_window_type Hamming ascii Range window type
range_window_coefficient 0.75 float64 Range window coefficient
pulse_repetition_frequency 1717.128973878037 float64 Hz PRF
radar_frequency 5405.000454334349 float64 MHz Radar frequency
line_time_interval 0.002055556299999998 float64 s
total_size 7259 uint32 MB Total product size
num_output_lines 13653 uint32 lines Raster height
num_samples_per_line 69965 uint32 samples Raster width
subset_offset_x 0 uint32 samples X coordinate of UL corner of subset in original image
subset_offset_y 0 uint32 samples Y coordinate of UL corner of subset in original image
srgr_flag 0 uint8 flag SRGR applied
avg_scene_height 475.47396462255193 float64 m Average scene height ellipsoid
map_projection - ascii Map projection applied
is_terrain_corrected 0 uint8 flag orthorectification applied
DEM - ascii Digital Elevation Model used
geo_ref_system - ascii geographic reference system
lat_pixel_res 99999.0 float64 deg pixel resolution in geocoded image
lon_pixel_res 99999.0 float64 deg pixel resolution in geocoded image
slant_range_to_first_pixel 801931.2520425358 float64 m Slant range to 1st data sample
ant_elev_corr_flag 0 uint8 flag Antenna elevation applied
range_spread_comp_flag 0 uint8 flag range spread compensation applied
replica_power_corr_flag 0 uint8 flag Replica pulse power correction applied
abs_calibration_flag 0 uint8 flag Product calibrated
calibration_factor 99999.0 float64 dB Calibration constant
chirp_power 99999.0 float64 Chirp power
inc_angle_comp_flag 0 uint8 flag incidence angle compensation applied
ref_inc_angle 99999.0 float64 Reference incidence angle
ref_slant_range 99999.0 float64 Reference slant range
ref_slant_range_exp 99999.0 float64 Reference slant range exponent
rescaling_factor 99999.0 float64 Rescaling factor
bistatic_correction_applied 1 uint8 flag
range_sampling_rate 64.34523812571427 float64 MHz Range Sampling Rate
range_bandwidth 56.5 float64 MHz Bandwidth total in range
azimuth_bandwidth 327.0 float64 Hz Bandwidth total in azimuth
multilook_flag 0 uint8 flag Multilook applied
coregistered_stack 0 uint8 flag Coregistration applied
external_calibration_file - ascii External calibration file used
orbit_state_vector_file - ascii Orbit file used
metadata_version 7.0 ascii AbsMetadata version
centre_lat 49.016998300610105 float64
centre_lon 138.7920057311667 float64
centre_heading 190.0771739121946 float64
centre_heading2 10.074918539903244 float64
Спойлер
Orbit_State_Vectors
SRGR_Coefficients
Doppler_Centroid_Coefficients
Band_IW1_VH
Band_IW1_VV
Band_IW2_VH
Band_IW2_VV
Band_IW3_VH
Band_IW3_VV
PRODUCT S1B_IW_SLC__1SDV_20210625T205535_20210625T205602_027518_0348ED_F8B3 ascii Product name
PRODUCT_TYPE SLC ascii Product type
SPH_DESCRIPTOR Sentinel-1 IW Level-1 SLC Product ascii Description
MISSION SENTINEL-1B ascii Satellite mission
ACQUISITION_MODE IW ascii Acquisition mode
antenna_pointing right ascii Right or left facing
BEAMS - ascii Beams used
SWATH - ascii Swath name
PROC_TIME 26-JUN-2021 05:15:04.689393 uint32 utc Processed time
Processing_system_identifier ESA Sentinel-1 IPF 003.31 ascii Processing system identifier
orbit_cycle 164 int32 Cycle
REL_ORBIT 17 int32 Track
ABS_ORBIT 27518 int32 Orbit
STATE_VECTOR_TIME 07-JUL-2021 20:54:37.000000 uint32 utc Time of orbit state vector
VECTOR_SOURCE - ascii State vector source
incidence_near 99999.0 float64 deg
incidence_far 99999.0 float64 deg
slice_num 4 int32 Slice number
data_take_id 215277 int32 Data take identifier
first_line_time 25-JUN-2021 20:55:35.491515 uint32 utc First zero doppler azimuth time
last_line_time 25-JUN-2021 20:56:02.452191 uint32 utc Last zero doppler azimuth time
first_near_lat 49.547898749088745 float64 deg
first_near_long 140.72330078710957 float64 deg
first_far_lat 49.99224677492301 float64 deg
first_far_long 137.26742815906346 float64 deg
last_near_lat 48.04353067969139 float64 deg
last_near_long 140.26500844135842 float64 deg
last_far_lat 48.48487683676428 float64 deg
last_far_long 136.9097944777493 float64 deg
PASS DESCENDING ascii ASCENDING or DESCENDING
SAMPLE_TYPE COMPLEX ascii DETECTED or COMPLEX
mds1_tx_rx_polar VH ascii Polarization
mds2_tx_rx_polar VV ascii Polarization
mds3_tx_rx_polar - ascii Polarization
mds4_tx_rx_polar - ascii Polarization
polsar_data 0 uint8 flag Polarimetric Matrix
algorithm - ascii Processing algorithm
azimuth_looks 1.0 float64
range_looks 1.0 float64
range_spacing 2.3295619999999997 float64 m Range sample spacing
azimuth_spacing 13.901069999999999 float64 m Azimuth sample spacing
range_window_type Hamming ascii Range window type
range_window_coefficient 0.75 float64 Range window coefficient
pulse_repetition_frequency 1717.128973878037 float64 Hz PRF
radar_frequency 5405.000454334349 float64 MHz Radar frequency
line_time_interval 0.002055556299999998 float64 s
total_size 7259 uint32 MB Total product size
num_output_lines 13653 uint32 lines Raster height
num_samples_per_line 69964 uint32 samples Raster width
subset_offset_x 0 uint32 samples X coordinate of UL corner of subset in original image
subset_offset_y 0 uint32 samples Y coordinate of UL corner of subset in original image
srgr_flag 0 uint8 flag SRGR applied
avg_scene_height 473.4238802428765 float64 m Average scene height ellipsoid
map_projection - ascii Map projection applied
is_terrain_corrected 0 uint8 flag orthorectification applied
DEM - ascii Digital Elevation Model used
geo_ref_system - ascii geographic reference system
lat_pixel_res 99999.0 float64 deg pixel resolution in geocoded image
lon_pixel_res 99999.0 float64 deg pixel resolution in geocoded image
slant_range_to_first_pixel 801931.2520425358 float64 m Slant range to 1st data sample
ant_elev_corr_flag 0 uint8 flag Antenna elevation applied
range_spread_comp_flag 0 uint8 flag range spread compensation applied
replica_power_corr_flag 0 uint8 flag Replica pulse power correction applied
abs_calibration_flag 0 uint8 flag Product calibrated
calibration_factor 99999.0 float64 dB Calibration constant
chirp_power 99999.0 float64 Chirp power
inc_angle_comp_flag 0 uint8 flag incidence angle compensation applied
ref_inc_angle 99999.0 float64 Reference incidence angle
ref_slant_range 99999.0 float64 Reference slant range
ref_slant_range_exp 99999.0 float64 Reference slant range exponent
rescaling_factor 99999.0 float64 Rescaling factor
bistatic_correction_applied 1 uint8 flag
range_sampling_rate 64.34523812571427 float64 MHz Range Sampling Rate
range_bandwidth 56.5 float64 MHz Bandwidth total in range
azimuth_bandwidth 327.0 float64 Hz Bandwidth total in azimuth
multilook_flag 0 uint8 flag Multilook applied
coregistered_stack 0 uint8 flag Coregistration applied
external_calibration_file - ascii External calibration file used
orbit_state_vector_file - ascii Orbit file used
metadata_version 7.0 ascii AbsMetadata version
centre_lat 49.01724561195156 float64
centre_lon 138.7914404643422 float64
centre_heading 190.0769570588702 float64
centre_heading2 10.074701562215607 float64
Вложения
dates.jpg
dates.jpg (47.16 КБ) 2382 просмотра

gamm
Гуру
Сообщения: 4046
Зарегистрирован: 15 окт 2010, 08:33
Репутация: 1050
Ваше звание: программист
Откуда: Казань

Re: Не работает корегистрация в SNAP для S1B

Сообщение gamm » 03 фев 2023, 12:35

метаданные нужные, снимки правтически совпадают (разница тысячные градуса, это парак илометров максимум). Значит нужно попробовать на тех же снимках, что в инструкции - может, что-то поменялось в данных или программе.

Константин Силкин
Завсегдатай
Сообщения: 420
Зарегистрирован: 21 мар 2012, 07:37
Репутация: 67
Откуда: Воронеж

Re: Не работает корегистрация в SNAP для S1B

Сообщение Константин Силкин » 03 фев 2023, 15:43

Iggi1981 , похоже волей-неволей вам всё же придётся стать экспертом по SNAP. Вам личное сообщение нельзя написать?

Iggi1981
Участник
Сообщения: 93
Зарегистрирован: 06 мар 2020, 17:07
Репутация: 19
Откуда: Москва

Re: Не работает корегистрация в SNAP для S1B

Сообщение Iggi1981 » 03 фев 2023, 15:50

Константин Силкин, не хотелось бы, мне надо решить одну частную задачу, ради этого осваивать целый пакет - слишком расточительно.
Спойлер
можно написать мне в почту pa81@ya.ru

MBG
Активный участник
Сообщения: 123
Зарегистрирован: 02 июн 2012, 12:02
Репутация: 35

Re: Не работает корегистрация в SNAP для S1B

Сообщение MBG » 11 фев 2023, 09:04

В открытом пакете PyGMTSAR вся обработка делается в браузере без установки какого-либо софта (используя бесплатный Google Colab), см. примеры на гитхаб странице проекта: https://github.com/mobigroup/gmtsar Также можно докер образы использовать для локальной обработки под любой операционной системой.
Изображение

MBG
Активный участник
Сообщения: 123
Зарегистрирован: 02 июн 2012, 12:02
Репутация: 35

Re: Не работает корегистрация в SNAP для S1B

Сообщение MBG » 11 фев 2023, 09:17

Константин Силкин писал(а):
02 фев 2023, 21:25
Да, давно назрела необходимость написать на сайте толковую инструкцию по обработке радарных данных
Есть же интерактивные примеры, куда уж проще. И даже на русском легко найти информацию, хоть на хабре: PyGMTSAR (Python GMTSAR) — Спутниковая интерферометрия для всех или Спутниковая интерферометрия для танцующих гор Ирана на Google Colab или PyGMTSAR, или спутниковая интерферометрия для всех с примерами Jupyter Python ноутбуков на Google Colab и так далее, включая изложение физмат основ. PyGMTSAR доступен и на облачных серверах НАСА (нужен платный аккаунт).

Константин Силкин
Завсегдатай
Сообщения: 420
Зарегистрирован: 21 мар 2012, 07:37
Репутация: 67
Откуда: Воронеж

Re: Не работает корегистрация в SNAP для S1B

Сообщение Константин Силкин » 11 фев 2023, 09:49

Вы правы, MBG, но мы-то тут что ушами хлопаем?

Ответить

Вернуться в «Обработка ДДЗ»

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

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