CC BY
0
СПУТНИКОВАЯ ГИДРОФИЗИКА
Научная статья
УДК 551.46.09 EDN: GYAIRZ
Применение спутниковых радиолокационных измерений ЗвпИпвЬЛА/Б для оценки дрейфа морского льда в прикромочной зоне Арктики
Е. В. Плотников 1 И. Е. Козлов \ Е. В. Жук \ А. В. Марченко 2
1 Морской гидрофизический институт РАН, Севастополь, Россия
2 Свальбардский международный университет, Лонгйир, Норвегия
н ev.plotnikov@ya.ru
Поступила в редакцию 06.09.2023; одобрена после рецензирования 29.12.2023; принята к публикации 18.01.2024.
Аннотация
Цель. Целью работы является построение автоматизированной системы расчета полей скорости дрейфа морского льда по данным радиолокационных (РЛ) измерений спутников 8епЫпе1-1А/Е на основе метода нормализованной максимальной кросс-корреляции. Изложены условия и результаты численного эксперимента, направленного на оценку эффективности указанной методики для 63 пар изображений района пролива Фрама в летне-осенние периоды 2017 и 2018 гг. Подробно описан алгоритм проведения расчетов, приведены качественные и количественные характеристики результатов. Изложены соображения об эффективности указанного подхода для регулярного мониторинга ледяного дрейфа.
Методы и результаты. Для расчетов использовался метод максимальной кросс-корреляции (МКК), основанный на автоматизированном поиске фотографически сходных фрагментов на парах изображений, для которых известен временной интервал съемки. При этом в качестве метрики близости использовался коэффициент корреляции Пирсона. В результате построены 63 поля скорости дрейфа морского льда в районе пролива Фрама, каждое из которых имеет пространственный масштаб порядка нескольких сотен тысяч квадратных километров. Предложен метод фильтрации ложных корреляций.
Выводы. Используемый в работе подход позволяет в автоматическом режиме восстанавливать динамику дрейфа морского льда по спутниковым снимкам с высоким пространственным разрешением (40 м). Восстановленные поля скорости охватывают значительные по площади фрагменты поверхности океана. Предложенный метод фильтрации ложных корреляций позволяет эффективно выделять фрагменты результатов расчетов с искажениями, обусловленными ограничениями алгоритма МКК.
Ключевые слова: динамика дрейфа морского льда, морской лед, оптический поток, метод максимальной кросс-корреляции, РЛ-снимки 8епйпе1-1А/Е, пролив Фрама, Арктика
Благодарности: Исследование выполнено при финансовой поддержке гранта Российского научного фонда № 21-17-00278 (анализ, валидация и создание веб-сервиса для хранения полей дрейфа льда), а также в рамках темы государственного задания FNNN-2024-0017 (создание методики расчета полей скорости дрейфа льда по спутниковым радиолокационным данным).
Для цитирования: Применение спутниковых радиолокационных измерений 5епйпе1-1А/Е для оценки дрейфа морского льда в прикромочной зоне Арктики / Е. В. Плотников [и др.] // Морской гидрофизический журнал. 2024. Т. 40, № 2. С. 312-324. EDN GYAIRZ.
© Плотников Е. В., Козлов И. Е., Жук Е. В., Марченко А. В., 2024
Original article
Evaluation of Sea Ice Drift in the Arctic Marginal Ice Zone based on Sentinel-1A/B Satellite Radar Measurements
E. V. Plotnikov 1, I. E. Kozlov 1, E. V. Zhuk 1, A. V. Marchenko 2
1 Marine Hydrophysical Institute, Russian Academy of Sciences, Sevastopol, Russia 2 University Centre in Svalbard, Longyearbyen, Norway H ev.plotnikov@ya.ru
Abstract
Purpose. The object of the work is to construct an automated system for calculating sea ice drift velocity fields using Sentinel-1A/B radar measurements based on the normalized maximum cross-correlation approach. The conditions and results of a numerical experiment aimed at evaluating the effectiveness of this technique for 63 pairs of radar images of the Fram Strait region for the summer-autumn periods in 2017 and 2018 are presented. Both the calculation algorithm and the qualitative and quantitative characteristics of the results are described in details. The effectiveness of the approach being applied to regular monitoring of ice drift is considered.
Methods and Results. The maximum cross-correlation (MCC) method was used for calculations. It is based on an automated finding of photographically similar fragments in the pairs of images, for which the sensing time interval is known. The Pearson correlation coefficient was applied as a proximity metric. As a result, 63 sea ice drift velocity fields were constructed for the Fram Strait region, each of which has a spatial scale of approximately several hundred thousand square kilometers. The method for filtering false correlations is proposed.
Conclusions. The approach applied in the study makes it possible to obtain automatically the sea ice drift velocity fields from the satellite data with high spatial resolution (40 m). The reconstructed velocity fields cover significant areas of the ocean surface. The method proposed for filtering false correlations permits to extract effectively the fragments with distortions resulting from the MCC algorithm limitations, from the calculation results.
Keywords: sea ice drift dynamics, sea ice, optical flow, maximum cross-correlation approach, Sentinel-1A/B images, Fram Strait, Arctic Ocean
Acknowledgements: The study was carried out with financial support of the Russian Science Foundation grant No. 21-17-00278 (analysis, validation and development of a web-service to store sea ice drift fields) and within the framework of the theme of state assignment FNNN-2024-0017 (development of the method for calculating sea ice drift velocity based on satellite radar data).
For citation: Plotnikov, E.V., Kozlov, I.E., Zhuk, E.V. and Marchenko, A.V., 2024. Evaluation of Sea Ice Drift in the Arctic Marginal Ice Zone Based on the Sentinel-1A/B Satellite Radar Measurements. Physical Oceanography, 31(2), pp. 284-294.
Введение
Мониторинг дрейфа морского льда является важной составляющей исследований геофизических процессов в Арктике. Получаемые при этом данные можно использовать при решении широкого круга теоретических и практических задач. На сегодняшний день одним из наиболее перспективных подходов для решения этой задачи является анализ спутниковых снимков в видимом, инфракрасном и микроволновом диапазонах. Последовательности спутниковых изображений, относящиеся к одному и тому же фрагменту морской поверхности и узкому временному интервалу, могут быть использованы для восстановления динамики морской поверхности, в частности, скорости и направления дрейфа морского льда. При этом для решения описанной задачи часто
применяются методы расчета оптического потока. По данному подходу существует обширная литература, в которой описаны как сами алгоритмы, так и особенности их использования на практике в зависимости от специфики исходных данных 1 2 [1-7]. В силу объема имеющихся на сегодняшний день архивов спутниковых снимков актуальной является проблема автоматизации расчетов. Успешное решение этой задачи открывает возможность получения регулярных и продолжительных наблюдений за динамикой дрейфа морского льда.
Данное направление получило широкое развитие в последнее десятилетие, результаты исследований изложены, напр., в работах [8-12]. При этом используются различные алгоритмы трекинга применительно к данным MODIS, Sen-tinel-1, Sentinel-2, Landsat-8 и ряду других преимущественно в оптическом диапазоне. Важными факторами являются разрешение исходных снимков, возможность эффективной оптимизации вычислений с целью обработки больших массивов данных, автоматизация расчетов. В настоящей работе описан автоматизированный подход, позволяющий оценивать динамику морского льда по радиолокакционным (РЛ) изображениям, покрывающим площадь порядка нескольких сотен тысяч квадратных километров, с разрешением 40 м. По информации авторов, практические результаты по использованию автоматической системы обработки данных с такими параметрами публикуются впервые.
После запуска спутниковых радиолокаторов с синтезированной апертурой (РСА) Sentinel-1A/B появилась возможность получать пары РЛ-изображений (РЛИ) арктического региона с небольшой временной задержкой, что открыло возможность оценивать кинематические характеристики различных мезо-и субмезомасштабных процессов верхнего слоя арктических морей [13-17]. В приложении к оценке полей скорости дрейфа льда авторами работы [14] была показана возможность эффективно рассчитывать такие поля в областях интенсивной вихревой динамики в прикромочной ледовой зоне (ПЛЗ) с пространственным разрешением 0(100 м) на основе метода максимальной кросс-корреляции. Дальнейшее развитие это направление получило в работах [18, 19], где авторы применили методы машинного обучения для идентификации вихрей ПЛЗ в данных спутниковых РСА.
Целью настоящей работы является описание общей методики и демонстрация эффективности расчетной схемы автоматической обработки спутниковых данных Sentinel-1A/B с целью восстановления полей дрейфа морского льда в ПЛЗ. Предложенная методика была применена к трехмесячному массиву данных, относящихся к району пролива Фрама. Использованный подход позволил в автоматическом режиме обработать статистически значимый объем данных и получить качественные результаты, которые подтверждают эффективность используемого подхода для оперативного мониторинга динамики дрейфа ледяных полей в прикромочной зоне Арктики.
1 Barron J. L., Thacker N. A. Tutorial: Computing 2D and 3D Optical Flow. Manchester, 2005. 12 p. (Tina Memo No. 2004-012).
2 StarkM. Optical Flow PIV: Improving the Accuracy and Applicability of Particle Image Veloc-imetry : Master Thesis. Zürich : ETH, Department of Mechanical and Process Engineering (D-MAVT), 2013. 57 p. https://doi.org/10.3929/ethz-a-009767070
314 МОРЖОЙ ГИДРОФИЗИЧЕСКИЙ ЖУРНАЛ том 40 № 2 2024
Данные и методы
В качестве исходных данных использованы РЛИ Sentinel-1A/B, полученные для района пролива Фрама в 2017 и 2018 гг. В частности, использовались продукты L1 Ground Range Detected с размером пикселя 40 х 40 м и пространственным разрешением 93 х 87 м в режиме сьемки Extra Wide Swath (сверхширокозахватный). Спутниковые РЛИ были получены с сайта https://scihub.copernicus.eu из архивов системы Европейского центра морских прогнозов Copernicus Open Access Hub.
Всего было отобрано 63 пары изображений, полученных в сентябре 2017 г., а также в июле и сентябре 2018 г. Для каждой пары РЛИ временной диапазон между последовательными измерениями составил от 48 до 52 мин. Радиолокационные изображения выбирались исходя из соображений информативности, определяемой выраженностью поверхностных структур в поле дрейфующего льда. Данные интерполировались на регулярную сетку с разрешением 40 х 40 м, после чего выделялись пересечения сканируемых областей на каждом РЛИ. Дополнительно изображения сглаживались при помощи медианного фильтра со скользящим окном размером 20 х 20 пикселей. На рис. 1 показана пара РЛИ Sen-tinel-1A/B за 26 июля 2018 г., полученных таким образом.
и ь
Р и с. 1. Пара РЛ-изображений Sentinel-1A (a) и Sentinel-1B (b) за 15 сентября 2018 г. с наличием дрейфующего льда в прикромочной ледовой зоне пролива Фрама
F i g. 1. A pair ofSentinel-1A (a) and Sentinel-1B (b) radar images for September 15, 2018 with the drifting ice present in the marginal ice zone of Fram Strait
Основным алгоритмом, использовавшимся при обработке данных, является алгоритм нормализированной максимальной кросс-корреляции (МКК) 1 [1, 3, 5-7]. Суть метода заключается в нахождении максимально коррелирующих фрагментов на анализируемых изображениях, смещение которых друг от друга не превышает теоретически возможного. Для каждого пикселя на обоих снимках выделяются окрестности, для которых строится нормализированная кросс-корреляционная функция. Отклонение положения максимума этой функции от центральной точки принимается за искомый сдвиг, т. е. смещение фрагмента текстуры от снимка к снимку. Формально подход можно описать
следующим образом: обозначим матрицы изображений как I и /2. Рассмотрим некоторый пиксель с индексами / иПусть окрестность имеет квадратную форму и размер N х N, где N - нечетно;
к = |-(н-1)/2+(н-1)/2, I = |;-(н-1)/2+(н-1)/2Ц, М X I (к, I) и М2 (и, = X Ь (к + и, I + V).
Тогда нормализированная кросс-корреляционная функция будет иметь
вид
F (u, v, i, j ) = -
£(/ (k,l) -M )(/2 (k + u,l + v) -M2 (u, v))
£(/(k,l)-M)2 I! X(/2(k + u,l + v)-M2)2
l ,k
l, k
Здесь — u < u < u , - v < v < v , где u и v - заданные максимальные
^ max ma^' max max' ^ max max ^
возможные сдвиги по обеим координатам. Искомые значения u0 и v0 для пикселя c индексами i и j таковы, что F ((u0, v0, i, j) = max (F (u, v, i, j))), где и и v
принимают значения из указанного выше промежутка. Найдя описанные величины для всех (N-1)/2 < i < S -(N -1)/2 +1 и (N -1)/2 < j <-(N -1)/2 +1, где (Sj, S2) - размер изображений, получим матрицы U и V размером (S - N - umax + l, S2 - N - vmax + l) .
На рис. 2 продемонстрированы обе компоненты (u0 и v0) векторного поля, построенного таким образом с параметрами N = 80 umx = vmax = 50 по данным, показанным на рис. 1.
Р и с. 2. Пиксельные сдвиги, соответствующие компонентам полей скорости u0 (a), v0 (b)
и полученные 26 июля 2018 г. на основе алгоритма МКК с использованием РЛ-изображений Sentinel-1A/B
F i g. 2. Pixel shifts corresponding to the velocity field components u0 (a) and v0 (b), and obtained
using the MCC algorithm for the Sentinel-IA/B radar images for July 26, 2018
Здесь обращает на себя внимание соответствие фрагментов, содержащих полезные данные, фрагментам, где явно просматривается наличие ледяных скоплений на исходных изображениях. Остальная часть сцены либо не информативна, либо содержит шум, являющийся проявлением эффекта, который называют «ложными корреляциями». Основным общепринятым подходом для фильтрации соответствующих фрагментов является ограничение снизу на минимальное значение кросс-корреляционной функции и отсечение пикселей, не удовлетворяющих данному условию. Кроме этого, очевидно, эффективным подходом здесь мог бы стать морфологический анализ полученного поля. Однако авторы решили пойти другим путем, приняв во внимание то, что оператор (1,12) ^ (и, V) не является симметричным.
Рассмотрим некоторый пиксель с индексами ]). Обозначим индексы его окрест
Пусть (и0, у0 ) - рассчитанный сдвиг и для любых смещений фрагмент 1Х (к, /) максимально коррелирует с фрагментом 12 (к + и0, / + у0 ) . Тогда корреляция 1 (к — и0, / + V ) и /2 (к, /) также должна быть высокой. Следовательно, если провести описанный расчет, взяв 12 в качестве первого изображения, а 1 -в качестве второго, результаты должны быть близкими при изменении знаков у обеих компонент. Назовем такой расчет обратным. Для оценки близости можно использовать косинус угла между векторами двух таких полей. На рис. 3, а показаны значения этой величины. Здесь явно просматривается соответствие фрагментов с высокими значениями косинуса фрагментам с информативными данными, что подтверждает эффективность описанного подхода. Установив оптимальное пороговое значение, можно выделить маску, для которой рассчитанные поля скорости являются релевантными. Для расчетов использовался порог, равный 0,92. Кроме этого, для удаления остаточных шумов можно применить градиентный фильтр, выделяющий области высокой пространственной неоднородности. Для каждого пикселя выберем максимальное значение вариации компонент u и v по окрестности размером 3 х 3 пикселя (рис. 3, Ь). Все фрагменты с превышением порогового значения 2,51 отсекаются.
После осуществления описанной процедуры фильтрации в получаемой маске остается ряд небольших по площади отдельных фрагментов, информативность которых представляется сомнительной. Для их удаления маска разбивается на связные компоненты, для каждой из которых рассчитывается число содержимых пикселей. Если это число меньше заданного порога, компонента вырезается. В качестве оценочного значения выбрано отношение площади компоненты к площади всей сцены, в качестве порога - значение 0,25%.
Суммируя все вышеизложенное, можно кратко описать алгоритм обработки пары РЛ-изображений в виде следующих шагов:
- распаковка и репроекция исходных данных на регулярную сетку;
- прямой и обратный расчеты по методике нормализированной максимальной кросс-корреляции;
ности как
- фильтрация нерелевантных фрагментов;
- сжатие сцены с учетом границ полученной на предыдущем шаге маски;
- удаление мелких отдельных компонент маски;
- преобразование сдвигов в поле скорости (м/с).
Р и с. 3. Пространственное распределение косинуса угла между векторами рассчитанных сдвигов, полученных при прямом и обратном расчетах (a); оценка однородности компонент сдвигов по окну размером 3 х 3 пикселя (b)
F i g. 3. Spatial distribution of the cosine of angle between the vectors of calculated shifts resulted from direct and inverse calculations (a); assessment of the homogeneity of shift components over the 3 х 3 pixel window (b)
Р и с. 4. Конечный результат расчета поля скорости (м/с) дрейфа льда в проливе Фрама за 26 июля 2018 г.
F i g. 4. Final result of calculation of the ice drift velocity (m/s) field in Fram Strait for July 26, 2018
Результат всего процесса обработки данных показан на рис. 4. Видно, что диапазон наблюдаемых значений скорости дрейфа льда составляет 0,1-0,8 м/с, среднее значение скорости - 0,2-0,3 м/с, доминирующее направление дрейфа льда - южное, юго-восточное. Максимальные значения скорости дрейфа наблюдаются в северо-западном секторе, а также на юго-восточной периферии показанной области.
Результаты
По итогам расчетов получено 63 векторных поля скорости дрейфа льда, описывающих горизонтальное перемещение скоплений морского льда в ПЛЗ пролива Фрама. Для пары изображений размером 10000 х 10000 пикселей время работы использованной реализации алгоритма МКК с размером скользящего окна 80 х 80 и максимальными сдвигами — 50, V „ = 50 при рас-
шах ^ шах
чете в два потока на машине с процессором Шв1 15 12500Н под управлением операционной системы ВвЫап ОЫи/Ыпих 12 составляет около двух часов.
Объективная оценка эффективности описанного метода является затруднительной в силу отсутствия эталонных данных о поверхностной морской динамике. Наиболее перспективным подходом здесь представляется сопоставление восстановленных по спутниковым данным полей скорости с натурными данными о направлении и скорости дрейфа льда, полученными на основе установленных на льдинах стационарных геолокационных датчиков (ОР^-треке-ров). На рис. 5 показана траектория дрейфа такого датчика на фоне рассчитанного по спутниковым данным поля скорости.
76,22е
С.Ш.
76.20
76.18
76,16
ft 11 ff t f f* / / // JLts t t / / / / t f f » f f ' / /
11 Ш t\f ........Ж......m sf f 11 fJb ......A.......Л.. /Д /.....T"J / * * S s
f f t / / / 4 1 NL ' ■Я -it гууг f / S y* S S У/*]? S / It if lit 7Г Я у* ^ * yf S 7Г -Я v
25,78
25,83
25,89
25.94
В.Д.
Р и с. 5. Сопоставление векторов скорости дрейфа льда по спутниковым изображениям с натурными измерениями, полученными по данным установленных на лед GPS-трекеров. Синей кривой показана траектория движения льдин по натурным данным. Красной стрелкой выделен фрагмент натурных измерений, совпадающий по времени со спутниковой съемкой. Моменты съемки Sentinel-1: 26.04.2019 в 05:02 и в 05:51. Траектория датчика - с 26.04.2019 (00:10) по 27.04.2019 (11:30)
F i g. 5. Comparison of the ice drift velocity vectors derived from satellite data with the field measurements obtained from the data of GPS trackers installed on ice. Blue curve shows the trajectory of ice floes based on the field data. Red arrow highlights a fragment of field measurements coinciding in time with satellite imagery. Time of Sentinel-1 acquisitions: at 05:02 UTC and at 5:51 UTC on 26.04.2019. GPS tracker is traced from 26.04.2019 (00:10) to 27.04.2019 (11:30)
К сожалению, авторы не обладают статистически значимым массивом натурных данных для проведения более детального сопоставления. Тем не менее из рис. 5 видно, что измерения по двум независимым источникам хорошо совпадают по направлению и по абсолютному значению. Кроме этого, определенное представление о релевантности подхода можно получить на основе экспертной оценки структуры восстанавливаемых полей и соответствия видимых смещений контрастных структур в текстурах изображений рассчитанным сдвигам. С этой точки зрения результаты расчетов демонстрируют высокое качество. Фрагменты с ярко выделяющимися рассогласованностями эффективно отсекаются с помощью процедуры автоматической фильтрации, при этом количество выброшенных пикселей пренебрежимо мало по сравнению с их общим числом в зонах с хорошо выраженными контрастами. Существенной избыточности или недостаточности этого алгоритма замечено не было. Мезомас-штабные и субмезомасштабные динамические структуры (такие как отдельные вихри, вихревые диполи или ледовые филаменты), различимые на снимках, хорошо просматриваются также и в получаемых полях скорости.
Визуализация спутниковых продуктов с полями дрейфа льда
На следующем этапе продукты с восстановленными по квазисихронным спутниковым данным полями скорости дрейфа морского льда размещаются на веб-странице лаборатории морских полярных исследований (ЛМПИ) Морского гидрофизического института РАН по адресу http://polar-space.ru/arctic_currents. Для этого рассчитанные поля сначала интерполируются на сетку с шагом 40 м и записываются в виде «еС^^-файлов. Для размещения в сети Интернет была разработана система визуализации данных на базе клиент-серверной архитектуры. На рис. 6 показана общая структура системы визуализации данного спутникового продукта.
Р и с. 6. Общая структура системы визуализации продуктов с полями дрейфа льда, представленная на сайте http://polar-space.ru/arctic_currents
F i g. 6. General structure of the system for visualizing the products with ice drift fields represented on the site http://polar-space.ru/arctic_currents
В рамках этой системы данные с полями дрейфа льда в формате netCDF хранятся в виде файлового архива на сервере. Для отображения скоростей дрейфа льда на географической карте был разработан специальный программный модуль на языке Python, преобразующий исходные данные в изображения с единой шкалой скорости для всех файлов. Название каждого файла рисунка содержит дату, что позволяет установить однозначное соответствие между запросом пользователя и отображаемым полем дрейфа льда.
Пользовательский интерфейс (ПИ) для выбора и отображения интересующего продукта был разработан с использованием javascript библиотек jQuery. Функции картографического сервиса реализованы с помощью библиотеки mapBox GL; ПИ позволяет выбирать поле дрейфа льда за интересующую дату, пролистывая массив дат или выбирая нужную дату из списка. На рис. 7 показан пример ПИ на сайте http://polar-space.ru/arctic_currents за выбранную дату 26 июля 2018 г.
Р и с. 7. Пользовательский интерфейс системы визуализации продуктов с полями дрейфа льда в прикромочной ледовой зоне пролива Фрама за 26 июля 2018 г. Слева - пример статического отображения полей дрейфа льда, справа - динамическое поле с использованием библиотеки Wind-JS
F i g. 7. User interface of the system for visualizing the products with ice drift fields in the marginal ice zone of Fram Strait for July 26, 2018. On the left is an example of a static display of ice drift fields, on the right - a dynamic field using the Wind-JS library
В рамках разработанного ПИ существует выбор отображения данных с полями дрейфа льда в виде не только статической картинки (рис. 7, слева), но и динамического поля (рис. 7, справа). Для реализации этой задачи используется библиотека Wind-JS с сайта https://github.com/Esri/wind-js, которая была адаптирована для отображения векторов дрейфа льда. Входными данными являются файлы в формате json, содержащие значения зональной и меридиональной компонент дрейфа льда. Для генерации этих файлов в автоматическом режиме был разработан специальный программный модуль на языке Python.
К настоящему времени на веб-сайте ЛМПИ МГИ размещены все обработанные за 2017 и 2018 гг. данные. Кроме того, существующий интерфейс дополнен возможностью скачивания интересующих полей дрейфа льда в формате netCDF.
Заключение
В работе дано описание общей методологии обработки квазисихронных измерений спутниковых РСА Sentinel-1 A/B, работавших в тандеме в течение 5 лет с 2016 по 2021 гг. с целью получения регулярных по времени полей скорости дрейфа льда в арктической прикромочной зоне. Временное разрешение получаемых продуктов с полями скорости дрейфа льда составляет, как правило, одни сутки в европейском секторе Арктики, и от 1 до 5 сут в других арктических районах.
Ограниченная по набору исходных данных валидация получаемых полей скорости с натурными измерениями показала хорошее соответствие между ними. Тем не менее работа в этом направлении будет продолжена для определения объективных статистик качества рассчитываемых полей скорости дрейфа льда.
На фоне существующих продуктов по дрейфу льда в Арктике представленный подход и восстанавливаемые поля скорости отличаются значительно более высоким пространственным разрешением (40 м), относительной простотой используемых методов и нетребовательностью к вычислительным ресурсам.
Полученные продукты с полями скорости дрейфа льда в прикромочной зоне Арктики доступны на сайте http://polar-space.ru/arctic_currents и могут быть использованы специалистами для решения ряда практических задач, в том числе для оценки динамики и эволюции ПЛЗ и протекающих в ней процессов, а также для сравнения с иными спутниковыми продуктами и валида-ции моделей общей циркуляции океана более низкого разрешения.
Спутниковые данные Sentinel-1A/B были получены из архивов системы Европейского центра морских прогнозов Copernicus Open Access Hub с сайта URL: https://scihub.copernicus.eu.
СПИСОК ЛИТЕРАТУРЫ
1. Алексанин А. И., Алексанина М. Г., Карнацкий А. Ю. Автоматический расчет скоростей поверхностных течений океана по последовательности спутниковых изображений // Современные проблемы дистанционного зондирования Земли из космоса. 2013. Т. 10, № 2. С. 131-142. EDN QZNBRH.
2. Beauchemin S. S., Barron J. L. The computation of optical flow // ACM Computing Surveys. 1995. Vol. 27, iss. 3. P. 433-466. https://doi.org/10.1145/212094.212141
3. An objective method for computing advective surface velocities from sequential infrared satellite images / W. J. Emery [et al.] // Journal of Geophysical Research: Oceans. 1986. Vol. 91, iss. C11. P. 12865-12878. https://doi.org/10.1029/JC091iC11p12865
4. Fortun D., Bouthemy P., Kervrann C. Optical flow modeling and computation: A survey // Computer Vision and Image Understanding. 2015. Vol. 134. P. 1-21. https://doi.org/10.1016/j.cviu.2015.02.008
5. Luo J., Konofagou E. E. A fast normalized cross-correlation calculation method for motion estimation // IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control. 2010. Vol. 57, iss. 6. P. 1347-1357. https://doi.org/10.1109/TUFFC.2010.1554
6. Lawrence P. The derivation of sea surface velocities from satellite imagery using maximum cross correlation (MCC) // The Plymouth Student Scientist. 2016. Vol. 9, iss. 1. P. 145-161.
7. Learning Optical Flow / D. Sun [et al.] // Computer Vision - ECCV 2008. ECCV 2008. Lecture Notes in Computer Science / Eds. D. Forsyth, P. Torr, A. Zisserman. Berlin ; Heidelberg : Springer, 2008. P. 83-97. (Lecture Notes in Computer Science ; vol. 5304). https://doi.org/10.1007/978-3-540-88690-7_7
8. Sea ice motion from satellite passive microwave imagery assessed with ERS SAR and buoy motions / R. Kwok [et al.] // Journal of Geophysical Research: Oceans. 1998. Vol. 103, iss. C4. P. 8191-8214. https://doi.org/10.1029/97JC03334
9. Samardzija I. Two Applications of a Cross-Correlation Based Ice Drift Tracking Algorithm; Ship-Based Marine Radar Images and Camera Images from a Fixed Structure // Proceedings of the 24th IAHR International Symposium on Ice (Vladivostok 2018) / Eds. T. Uvarova, G. Ar-batskaia. Vladivostok : IAHR, 2018. P. 141-151.
10. Lopez-AcostaR., SchodlokM. P., WilhelmusM. M. Ice Floe Tracker: An algorithm to automatically retrieve Lagrangian trajectories via feature matching from moderate-resolution visual imagery // Remote Sensing of Environment. 2019. Vol. 234. 111406. https://doi.org/10.1016/j.rse.2019.111406
11. Wang M., König M., Oppelt N. Partial Shape Recognition for Sea Ice Motion Retrieval in the Marginal Ice Zone from Sentinel-1 and Sentinel-2 // Remote Sensing. 2021. Vol. 13, iss. 21. 4473. https://doi.org/10.3390/rs13214473
12. Howell S. E. L., Brady M., Komarov A. S. Generating large-scale sea ice motion from Sentinel-1 and the RADARSAT Constellation Mission using the Environment and Climate Change Canada automated sea ice tracking system // The Cryosphere. 2022. Vol. 16, iss. 3. P. 1125-1139. https://doi.org/10.5194/tc-16-1125-2022
13. Козлов И. Е., Михайличенко Т. В. Оценка фазовой скорости внутренних волн в Арктике по данным последовательных спутниковых РСА-измерений // Современные проблемы дистанционного зондирования Земли из космоса. 2021. Т. 18, № 5. С. 181-192. EDN WZMNHK. doi:10.21046/2070-7401-2021-18-5-181-192
14. Kozlov I. E., Plotnikov E. V., Manucharyan G. E. Brief Communication: Mesoscale and sub-mesoscale dynamics in the marginal ice zone from sequential synthetic aperture radar observations // The Cryosphere. 2020. Vol. 14, iss. 9. P. 2941-2947. https://doi.org/10.5194/tc-14-2941-2020
15. Kozlov I. E., Atadzhanova O. A. Eddies in the Marginal Ice Zone of Fram Strait and Svalbard from Spaceborne SAR Observations in Winter // Remote Sensing. 2022. Vol. 14, iss. 1. 134. https://doi.org/10.3390/rs14010134
16. Артамонова А. В., Козлов И. Е. Вихри в Норвежском и Гренландском морях по данным спутниковых РСА в летний период 2007 года // Морской гидрофизический журнал. 2023. Т. 39, № 1. С. 120-133. EDN OFTTND. https://doi.org/10.29039/0233-7584-2023-1-120-133
17. High-amplitude internal waves southeast of Spitsbergen / A. V. Marchenko [et al.] // Continental Shelf Research. 2021. Vol. 227. 104523. https://doi.org/10.1016/j.csr.2021.104523
18. Khachatrian E., Sandalyuk N. V. On the Exploitation of Multimodal Remote Sensing Data Combination for Mesoscale/Submesoscale Eddy Detection in the Marginal Ice Zone // IEEE Geoscience and Remote Sensing Letters. 2022. Vol. 19. 3513805. https://doi.org/10.1109/LGRS.2022.3215202
19. Khachatrian E., SandalyukN., Lozou P. Eddy Detection in the Marginal Ice Zone with Sen-tinel-1 Data Using YOLOv5 // Remote Sensing. 2023. Vol. 15, iss. 9. 2244. https://doi.org/10.3390/rs15092244
Об авторах:
Плотников Евгений Владимирович, младший научный сотрудник, лаборатория морских полярных исследований, отдел дистанционных методов исследования, ФГБУН ФИЦ МГИ (299011, Россия, г. Севастополь, ул. Капитанская, д. 2), ORCID ID: 0000-0003-4365-9369, Scopus Author ID: 57190382092, SciProfiles: 2325823, SPIN-код: 6601-1958, ev.plotnikov@ya.ru
Козлов Игорь Евгеньевич, ведущий научный сотрудник, заведующий лабораторией морских полярных исследований, ФГБУН ФИЦ МГИ (299011, Россия, г. Севастополь, ул. Капитанская, д. 2), кандидат физико-математических наук, ORCID ID: 0000-0001-6378-8956, Re-searcherID: G-1103-2014, Scopus Author ID: 49963767500, SPIN-код: 3267-6693, ik@mhi-ras.ru
Жук Елена Владимировна, младший научный сотрудник, отдел океанографии, ФГБУН ФИЦ МГИ (299011, Россия, г. Севастополь, ул. Капитанская, д. 2), ORCID ID: 0000-0002-42637734, Scopus Author ID: 57191412660, SPIN-код: 3814-6300, alenixx@gmail.com
Марченко Алексей Валерьевич, доктор физико-математических наук, профессор, департамент арктических технологий, Свальбардский международный университет (P.O. Box 156 N-9171, Лонгйир, Норвегия), ORCID ID: 0000-0003-4169-0063, Scopus Author ID: 7101880290, ResearcherID: GSD-3516-2022, SPIN-код: 9213-8763, aleksey.marchenko@unis.no
