Продуктите от ниво 1- открит обхват( GRD), се състоят от фокусирани данни за SAR, които са проектирани до земната повърхност, като се използва елипсоиден модел на Земята, като WGS84. Елипсоидната проекция на GRD продуктите се коригира, като се използва височинен модел на терена( DEM) – ACE30, предоставен и за територията, на която се намира ледник Перуника [ 10 ].
Обработката се извършва от минимум два продукта GRD Sentinel-1 IW в среда на софтуерния продукт SNAP [ 21 ]. За целта са използвани данни в периода 24.07-05.08.2019 г. и 08.01-20.01.2020 г., с обхват над ледника Перуника.
За да се генерира карта на скоростта на хоризонталните премествания на ледника, с помощта на инструмента за проследяване – Offset tracking, входните продукти е необходимо да са минимум два, с покритие върху една и съща зона, придобити за различен времеви период. За целта са използвани следните GDR продукти, регистрирани през 12 дни:
S1A _ IW _ GRDH _ 1SSH _ 20190724T081741 _ 20190724T081806 _ 028256 _ 033122 _ 2C13 S1A _ IW _ GRDH _ 1SSH _ 20190805T081741 _ 20190805T081806 _ 028431 _ 03367C _ 2F03 S1A _ IW _ GRDH _ 1SSH _ 20200108T081742 _ 20200108T081807 _ 030706 _ 03851C _ 63B5 S1A _ IW _ GRDH _ 1SSH _ 20200120T081742 _ 20200120T081807 _ 030881 _ 038B41 _ 9A9B
Последователността на прилагания метод на работа предоставя възможност за оценка на настъпилите хоризонтални премествания за определен времеви диапазон.
За обработка на отместването двете изображения е необходимо да бъдат регистрирани в стек( пакетирани вертикално). Изображението, което е придобито от по-ранен период, се избира за главно, а другото изображение е подчинено. Пикселите в подчинено изображение се преизчисляват, за да се приведат в съответствие с главното изображение, с помощта на орбитални данни и референтен височинен модел на терена.
Корегистрацията гарантира, че всеки наземен обект от неподвижната сцена допринася за един и същ пиксел( по обхват и азимут) както в главното, така и в подчиненото изображение.
За прилагане на“ offset tracking” се използва DEM Assisted Coregistration, като се създават продуктите строго въз основа на геометрията, DEM, орбитални данни и време. Това избягва възможно изкривяване на изображението поради движението.
Преди корегистрация на изображението, орбитални вектори в оригиналния продукт на GRD трябва да бъдат актуализирани с по-точните, предоставени в наличния орбитален файл.
Услугата за точно определяне на орбитата за SENTINEL-1 предоставя файлове за възстановени орбити и орбитални файлове с точни орбитални ефемериди. Тези файлове обхващат приблизително 28 часа и съдържат вектори на състоянието на орбитата на фиксирани стъпки от интервали от 10 секунди. Файловете се генерират по един файл на ден и се доставят в рамките на 20 дни след получаване на данните [ 10 ].
Операторът за проследяване на преместването извършва оценка на движението на повърхностите на ледника между главни и подчинени изображения в посока наклона и по азимут. Извършва се кръстосана корелация на избрана контролна точка на земната повърхност( GCP) в главни и подчинени изображения. След което скоростите на ледника върху избраните GCP се изчисляват въз основа на компенсациите, изчислени от кръстосаната корелация и се генерира картата на скоростта на движение на ледника чрез интерполация на изчислените скорости на GCP мрежата.
Проследяването на отместването се извършва в следната последователност:
• За всяка точка от зададената от потребителя GCP мрежа в главното изображение се изчислява съответното ѝ положение в подчинено изображение, използвайки нормализирана кръстосана корелация;
• Ако изчисленото изместване между главните и подчинените GCP позиции надвишава максималното изместване( изчислено от максималната скорост, определена от потребителя), тогава GCP точката се маркира като по-външна;
• Изчислява се средна стойност за компенсиране на валидни GCP точки;
• Запълват се празните участъци, причинени от външни елементи. Отместването в точката на празни участъци се заменя с ново отместване, изчислено от средно претеглена средна стойност;
• Изчисляват се скоростите за всички точки на GCP мрежата от техните компенсации;
• Изчисляват се скоростите за всички пиксели в главното изображение от скоростите в GCP мрежата чрез интерполация [ 10 ].
281