Соган · предресерч

Покрытие открытыми данными ДЗЗ

Кратко.

Участок покрыт открытыми данными ДЗЗ очень хорошо — заметно лучше среднего, и по нескольким позициям результат оказался лучше ожидаемого.

УСТАНОВЛЕНО ДЛЯ САМОГО УЧАСТКА (проверено запросами к каталогам по его координатам). Тайл Sentinel-2 — 40RDS, орбита R020, 860 сцен L2A на AWS и архив с 21.11.2015 в CDSE, съёмка каждые 4–5 суток. Облачность именно над 62-гектарным участком, посчитанная по маске SCL для всех 158 сцен 2024–2025 гг., имеет медиану 0,0 %, и 128 сцен из 158 (81 %) практически безоблачны — то есть около 64 чистых снимков в год, а «плохая» зима по тайлу оказывается вполне рабочей над самим участком. Landsat: основная сцена WRS-2 160/040 (участок в 27,4 км от края кадра) плюс краевые 159/040 и 159/041; всего 2962 сцены Collection 2 Level-2 с 13.04.1984, в 2025 г. — 84 даты съёмки. ASTER: 116 сцен L1T, и главное — 41 дневная сцена до 01.04.2008, когда SWIR ещё работал, плюс 34 готовых продукта AST_07XT того же периода; после 23.04.2008 SWIR мёртв, это подтверждено официальными предупреждениями NASA/JPL и LP DAAC. Гиперспектр есть, и это самый неожиданный результат: EMIT — 14 сцен, каждая покрывает участок на 100 %, включая готовые карты минералов EMITL2BMIN; EnMAP — 2 съёмки, обе на 100 % площади, одна (30.12.2023) с нулевой облачностью. Sentinel-1 — 951 сцена GRD и 941 SLC по трём трекам (93 и 166 нисходящие, 57 восходящая). ЦМР — все четыре набора без пустот, высоты 1592–1764 м. Ночной TIR — 33 сцены ASTER и более тысячи гранул ECOSTRESS (тайл тот же, 40RDS), из них свыше 600 ночных.

ГЛАВНЫЕ ОГРАНИЧЕНИЯ. Первое: над участком снимает только один аппарат Sentinel-1 (S1A, с 2026 г. добавился S1D); S1B и S1C дали ноль. Значит, интерферометрический повтор — 12 суток по каждому треку, шестисуточных пар нет, и планировать InSAR надо исходя из этого. Второе: рельеф крутой — медианный уклон 16,9°, 38 % площади круче 20°, перепад 157 м на 1480 м длины; без топографической нормализации оптики и без учёта layover и радиотени в SAR результаты будут артефактными. Третье: PRISMA остаётся белым пятном — доступ бесплатный, но каталог виден только после регистрации с одобрением проекта, поэтому наличие архивных сцен над участком проверить не удалось.

ОБЩАЯ ОЦЕНКА ПРИГОДНОСТИ. Для литологического и минералогического дешифрирования ультраосновного комплекса набор редкий: одновременно есть архивный ASTER SWIR (41 сцена до 2008 г.), спутниковый гиперспектр 30 м (EnMAP, в том числе безоблачная сцена) и гиперспектр 60 м с готовыми минеральными продуктами (EMIT, 14 дат). Это позволяет строить перекрёстную проверку: одни и те же спектральные признаки серпентинитов, хромититов и вмещающих пород можно независимо проверить по трём разным сенсорам, а не полагаться на один. Практически всё, кроме PRISMA, доступно либо вообще без регистрации, либо с бесплатной учётной записью.

1. Ключи поиска для участка: MGRS, WRS-2, тайлы ЦМР

Первое, что нужно для любого каталога, — правильные «ключи» участка. Центр участка 28,35520° с.ш., 56,83704° в.д.; проверено пересчётом pyproj EPSG:4326 → EPSG:32640: X = 484 030,20 м, Y = 3 136 561,50 м (совпадает с заданным X≈484030, Y≈3136561). Прямоугольник участка в WGS84: долгота 56,82948–56,84459°, широта 28,35330–28,35710°; размеры 1480,02 × 418,89 м, площадь 0,6200 км² (62,0 га) — сходится с заданием.

Ключи, вычисленные и подтверждённые запросами к каталогам: • Sentinel-2 (MGRS): тайл 40RDS (зона 40, широтная полоса R, квадрат DS). Вычислен вручную и подтверждён фактическими сценами (S2B_40RDS_…, T40RDS в именах CDSE). • Landsat (WRS-2): участок попадает сразу в три сцены — 160/040 (основная), 159/040 и 159/041 (краевое перекрытие). Проверено геометрически: расстояние от участка до ближайшего края сцены 160/040 — 27,4 км, до края 159/040 — 9,9 км, 159/041 — 9,8 км; до центров сцен 78,5 / 112,9 / 120,0 км соответственно. • ЦМР (SRTM, Copernicus, NASADEM, ALOS, ASTER GDEM): тайл 1°×1° N28E056 во всех наборах. • ECOSTRESS (плиточные продукты L2T и L1CT): тот же MGRS-тайл 40RDS. • Sentinel-1: относительные орбиты (треки) 93 (нисходящая), 166 (нисходящая), 57 (восходящая). • Sentinel-2: относительная орбита R020 (единственная над участком).

Местное время съёмки (IRST = UTC+3:30, летнего перевода в Иране нет с 2022 г.): Sentinel-2 ≈ 07:01 UTC = 10:31 местного; Landsat 160/040 ≈ 06:45 UTC = 10:15 местного; EnMAP ≈ 07:22–07:29 UTC = 10:52–10:59; ASTER день ≈ 06:56–07:15 UTC = 10:26–10:45, ASTER ночь ≈ 17:10–18:26 UTC = 20:40–21:56; Sentinel-1 нисходящие ≈ 02:07–02:15 UTC = 05:37–05:45, восходящая ≈ 14:18 UTC = 17:48.

Источники (5)

2. Sentinel-2 L1C/L2A — базовый мультиспектральный слой (тайл 40RDS)

Участок целиком лежит в тайле T40RDS, орбита R020, съёмка около 07:01 UTC (10:31 местного). Архив над участком начинается 21.11.2015 (S2A_MSIL1C_20151121T065202_N0500_R020_T40RDS; в CDSE тот же снимок доступен и как L2A после переобработки baseline N0500). На AWS (Earth Search) L2A-архив мельче: самая ранняя сцена там S2A_40RDS_20170104_0_L2A от 04.01.2017, всего 860 сцен L2A над точкой.

Частота. Одна орбита R020, номинально 5 суток при двух аппаратах. Фактически за 2024–2025 гг. — 158 сцен L2A за 157 уникальных дат: интервал 5 суток встречается 128 раз, остальное — 2–3-суточные интервалы периода совместной работы трёх аппаратов (S2A 53 сцены, S2B 72, S2C 33). За 2025 г. в CDSE 86 продуктов L2A — в среднем съёмка каждые 4,2 суток.

ОБЛАЧНОСТЬ. Здесь надо различать два числа. Штатное поле eo:cloud_cover считается по всему тайлу 110×110 км и для 40RDS даёт медиану 4,5 % и сильную сезонность (медиана 32,5 % в феврале против 0,1 % в ноябре). Но участок — 62 га, и над ним облачность гораздо ниже. Проверено прямым чтением маски SCL (20 м) из COG на AWS для всех 158 сцен 2024–2025 гг. по окну участка (1554 пикселя SCL): медиана облачности над участком 0,0 %; 128 сцен из 158 (81 %) практически без облаков (менее 0,1 % пикселей классов 3, 8, 9, 10). Помесячно чистых (менее 1 % облаков) сцен: янв 7 из 12, фев 7 из 12, мар 10 из 14, апр 7 из 12, май 13 из 14, июн 11 из 13, июл 14 из 15, авг 11 из 12, сен 13 из 13, окт 13 из 14, ноя 11 из 13, дек 11 из 14. Вывод: около 64 полностью чистых снимков в год; худший сезон — январь–апрель (58–71 % чистых), лучший — май–ноябрь (85–100 %).

Как получить. Без регистрации — STAC-поиск Earth Search (POST https://earth-search.aws.element84.com/v1/search, collections=["sentinel-2-l2a"], intersects=точка) и прямое чтение COG из открытого бакета s3://sentinel-cogs (регион us-west-2, AWS_NO_SIGN_REQUEST=YES) — именно так и получена статистика выше. С бесплатной регистрацией — Copernicus Browser, CDSE OData или S3 (архивы SAFE, полная глубина с 2015 г.). Формат: JPEG2000 в SAFE (CDSE) или COG по полосам (AWS). Разрешение 10 м (B2, B3, B4, B8), 20 м (B5–B7, B8A, B11, B12), 60 м (B1, B9, B10). Лицензия — открытые данные Copernicus.

Источники (4)

3. Landsat 8/9 OLI-TIRS — WRS-2 160/040 (плюс 159/040 и 159/041)

Участок попадает в перекрытие трёх сцен WRS-2, но «своя» сцена — path 160 / row 040: участок в 27,4 км от ближайшего края кадра, тогда как в 159/040 и 159/041 он в 9,8–9,9 км от края (краевая зона, где выше искажения и риск попасть на границу кадра). Съёмка 160/040 ≈ 06:45 UTC (10:15 местного), высота Солнца в июле около 66°.

Наличие сцен (Planetary Computer, коллекция landsat-c2-l2, поиск по точке участка): за 2023–2025 гг. Landsat 8 и 9 дают 100 сцен на выборку (52 L8 и 48 L9). За 2025 г. над участком 125 сцен по трём path/row, 84 уникальные даты съёмки — благодаря перекрытию трёх кадров интервал чередуется: 38 интервалов по 7 суток и 38 по 1 суткам за год, то есть де-факто снимок Landsat примерно каждые 3,5–4 суток вместо номинальных 8 (L8+L9) или 16 (один аппарат).

Облачность (уровень сцены 185×180 км, path/row 160/040, весь архив): 817 из 1143 сцен (71,5 %) имеют менее 10 % облачности. Для 159/041 в 2025 г. с марта по ноябрь большинство сцен имеет 0–13 % облачности.

Как получить. (а) Без регистрации — STAC Planetary Computer (collections=["landsat-c2-l2"]) с подписью ссылок через /api/sas/v1/sign. (б) USGS EarthExplorer (бесплатная регистрация ERS) — набор данных «Landsat Collection 2 Level-2», фильтр Path=160, Row=040. (в) Программно — USGS M2M API (https://m2m.cr.usgs.gov/api/api/json/stable/): требуется учётная запись ERS и отдельно одобряемый запрос на доступ к M2M (документация закрыта логином — проверено, запрос к /api/docs/json/ перенаправляет на ers.cr.usgs.gov/login). (г) AWS Open Data: s3://usgs-landsat (requester-pays). Формат: Cloud-Optimized GeoTIFF, Collection 2 Level-2 (отражение поверхности 30 м и температура поверхности, снятая радиометром 100 м и пересэмплированная в 30 м), плюс маска QA_PIXEL. Категория T1 у 2892 из 2962 сцен архива.

Источники (3)

4. Исторический архив Landsat 4/5 TM и 7 ETM+ — глубина с 1984 года

Проверено поиском по точке участка в коллекции landsat-c2-l2 Planetary Computer (это Collection 2 Level-2, то есть уже атмосферно скорректированное отражение поверхности): всего 2962 сцены над участком за период 13.04.1984 – 07.07.2026. Разбивка по аппаратам: Landsat 4 — 16, Landsat 5 (TM) — 873, Landsat 7 (ETM+) — 913, Landsat 8 — 874, Landsat 9 — 286. По кадрам: 160/040 — 1143, 159/040 — 1121, 159/041 — 698.

Ограничения, о которых надо помнить при работе с историей: • Landsat 7 ETM+ с 31 мая 2003 г. работает в режиме SLC-off — в каждом кадре клиновидные полосы пропусков, расширяющиеся к краям; в центре кадра потери меньше. Участок в 160/040 находится далеко от края, но полосы всё равно возможны — проверять маску QA_PIXEL перед использованием. • TM и ETM+ дают только две полосы SWIR (1,55–1,75 и 2,08–2,35 мкм) — для различения минералов гидроксильной группы этого мало; серьёзная минералогия по этим данным не делается, они полезны как временной ряд и для структурного дешифрирования. • Тепловой канал TM/ETM+ — 120 и 60 м соответственно (в продукте пересэмплирован в 30 м); съёмка дневная.

Как получить: те же маршруты, что и для L8/L9 (EarthExplorer, Planetary Computer, s3://usgs-landsat). В EarthExplorer это наборы «Landsat 4-5 TM C2 L2» и «Landsat 7 ETM+ C2 L2».

Источники (2)

5. ASTER — ключевой архив SWIR до 23 апреля 2008 года

ASTER — самый важный из «старых» источников для ультраосновного комплекса: 6 полос SWIR (2,145–2,430 мкм, 30 м) и 5 полос TIR (8,125–11,65 мкм, 90 м). Но SWIR мёртв: 23 апреля 2008 г. температура детекторов SWIR резко выросла, полосы 5–9 ушли в насыщение; процедура перезапуска 7 мая 2008 г. не помогла, в январе 2009 г. NASA официально признала SWIR неработающим. Все SWIR-данные, снятые после 1 апреля 2008 г., помечены как непригодные (насыщение и сильный полосатый шум). Поэтому для минералогии нужны только сцены до апреля 2008 г.

ФАКТИЧЕСКОЕ НАЛИЧИЕ (проверено через NASA CMR, поиск гранул по точке участка): • AST_L1T (ортотрансформированная радиояркость) — 116 сцен, первая 29.06.2000. Из них 83 дневных и 33 ночных. • До 01.04.2008 — 43 сцены, из них 41 дневная. Это и есть рабочий пул для SWIR. Полный список дат дневных сцен с исправным SWIR: 2000-06-29, 2000-08-16, 2001-02-15, 2001-05-22, 2001-06-23, 2001-07-18, 2001-08-03, 2001-09-20, 2001-10-29, 2001-11-07, 2002-01-26, 2002-04-07, 2002-07-05, 2002-08-22, 2003-01-20, 2003-03-02, 2003-06-06, 2003-07-24, 2003-08-25, 2003-09-01, 2003-10-03, 2003-11-13, 2004-02-08, 2004-03-27, 2004-07-01, 2004-07-26, 2004-10-05, 2005-03-14, 2005-03-30, 2005-05-17, 2005-05-26, 2005-12-27, 2006-03-17, 2006-04-02, 2006-04-18, 2006-06-30, 2006-07-23, 2006-08-24, 2007-04-30, 2007-06-01, 2007-07-03. • AST_07XT (отражение поверхности VNIR и SWIR с коррекцией перекрёстных наводок, crosstalk) — 70 гранул, из них 34 дневных до 01.04.2008. Это самый удобный готовый продукт: атмосферная и crosstalk-коррекция уже выполнены. • AST_09XT (радиояркость поверхности с crosstalk-коррекцией) — 70 гранул. • AST_L1A — 88 гранул, AST_L1B — 100 гранул. • AST_05 (излучательная способность TIR) и AST_08 (кинетическая температура) — по 98 гранул, из них по 28 ночных.

Пример идентификатора: AST_L1T_00406292000071204_20250125084644 (формат: MMDDYYYYhhmmss съёмки плюс дата обработки). Файлы раздаются по полосам как GeoTIFF, например …SWIR_B04.tif … _SWIR_B09.tif, …_VNIR*.tif, …_TIR.tif — SWIR-полосы можно качать отдельно, не забирая всю сцену.

Как искать. (а) EarthExplorer: раздел «NASA LPDAAC Collections» → «ASTER Level 1T», координаты участка, обязательно ограничить дату «до 2008-04-01» и Day/Night = Day, фильтр по облачности (в EarthExplorer поле Scene Cloud Cover есть; в CMR его нет — проверено, у гранул AST_L1T среди AdditionalAttributes облачности нет). (б) Earthdata Search (https://search.earthdata.nasa.gov) — набор AST_L1T или AST_07XT, поиск по точке. (в) Программно — CMR: https://cmr.earthdata.nasa.gov/search/granules.umm_json?short_name=AST_L1T&point=56.83704,28.35520&page_size=200 (без регистрации). Сама загрузка — с бесплатной учётной записью NASA Earthdata Login по ссылкам вида https://data.lpdaac.earthdatacloud.nasa.gov/lp-prod-protected/AST_L1T.004/…

Отдельно: ASTER GED (глобальный набор излучательной способности) — гранулы AG100.v003.29.056.0001 (100 м) и AG1km.v003.29.056.0010 (1 км) покрывают участок; ASTER GDEM v3 — гранула ASTGTMV003_N28E056.

Источники (6)

6. EMIT (гиперспектрометр на МКС) — 14 сцен, участок покрыт на 100 %

Одна из главных находок предисследования. EMIT — визуально-инфракрасный изображающий спектрометр NASA/JPL на МКС: диапазон 380–2500 нм, около 285 полос с шагом примерно 7,4 нм, пространственное разрешение около 60 м, полоса захвата около 75 км. Миссия изначально нацелена именно на минералогию засушливых пылеобразующих районов и снимает освещённую сушу между 52° с.ш. и 52° ю.ш. С 2024 г. миссия в фазе продления и снимает всю доступную МКС сушу.

ПРОВЕРЕНО через NASA CMR: над участком 14 гранул EMITL2ARFL (отражение поверхности), и все 14 покрывают участок на 100 % — проверено пересечением полигонов гранул с прямоугольником участка. Те же 14 дат есть в EMITL1BRAD (радиояркость), EMITL2AMASK (маски) и, что особенно ценно, в EMITL2BMIN — готовых картах минералов (два «наиболее подходящих» минерала на пиксель по библиотеке USGS и алгоритму Tetracorder).

Даты и идентификаторы (UTC): 2022-08-17T11:00 EMIT_L2A_RFL_001_20220817T110034_2222907_013; 2023-02-05T06:21 EMIT_L2A_RFL_001_20230205T062127_2303605_022; 2023-04-08T05:50 EMIT_L2A_RFL_001_20230408T055003_2309804_026; 2023-04-19T10:01 EMIT_L2A_RFL_001_20230419T100117_2310907_011; 2023-06-26T07:08 EMIT_L2A_RFL_001_20230626T070810_2317705_012; 2023-10-04T06:56 EMIT_L2A_RFL_001_20231004T065608_2327705_026; 2024-02-26T06:08 EMIT_L2A_RFL_001_20240226T060843_2405704_016; 2024-06-10T12:28 EMIT_L2A_RFL_001_20240610T122837_2416208_008; 2024-10-04T06:00 EMIT_L2A_RFL_001_20241004T060044_2427804_025; 2025-02-16T09:05 EMIT_L2A_RFL_001_20250216T090521_2504706_010; 2026-02-01T05:58 EMIT_L2A_RFL_001_20260201T055833_2603203_025; 2026-04-18T08:29 EMIT_L2A_RFL_001_20260418T082952_2610805_043; 2026-04-22T06:57 EMIT_L2A_RFL_001_20260422T065738_2611204_030; 2026-04-26T05:24 EMIT_L2A_RFL_001_20260426T052448_2611603_008.

Отдельная ценность: орбита МКС не солнечно-синхронная, поэтому времена съёмки разбросаны от 05:24 до 12:28 UTC (08:54–15:58 местного) — получается набор с разными углами освещения на одном и том же участке, что помогает отделять спектральные эффекты от эффектов рельефа и теней.

Как проверить и скачать. Веб: Earthdata Search (https://search.earthdata.nasa.gov), ввести «EMIT L2A Reflectance», задать точку 28.35520 с.ш., 56.83704 в.д. Программно без регистрации: https://cmr.earthdata.nasa.gov/search/granules.umm_json?short_name=EMITL2ARFL&point=56.83704,28.35520&page_size=50. Загрузка — с бесплатным NASA Earthdata Login по ссылкам вида https://data.lpdaac.earthdatacloud.nasa.gov/lp-prod-protected/EMITL2ARFL.001//.nc. Формат — NetCDF-4 (несетчатый, с массивами геолокации; для наложения на проектную СК нужен ортоперенос, у JPL для этого есть пакет emit-utils). Данные полностью открыты.

Источники (5)

7. EnMAP (DLR) — 2 съёмки участка, одна с нулевой облачностью

Вторая важная находка. EnMAP — немецкий гиперспектральный спутник (запуск 1 апреля 2022 г., штатная работа с ноября 2022 г.): 224 полосы в диапазоне 420–2450 нм, средний спектральный шаг 6,5 нм в VNIR и 10 нм в SWIR, разрешение на местности 30×30 м, полоса захвата 30 км, повтор в надире 27 суток, при отклонении на 30° — до 4 суток. Данные бесплатны и открыты с ноября 2022 г.

ПРОВЕРЕНО в STAC-каталоге DLR EOC Geoservice (коллекция ENMAP_HSI_L2A, поиск по прямоугольнику участка): найдены две даты съёмки, обе покрывают участок на 100 %: • 30.12.2023, 07:22:48 UTC, облачность 0 %, высота Солнца 36,9° — ENMAP01-_L2A-DT000005585920231230T072248Z_002_V010502_… (доступны две версии обработки: V010502 от 2025-03-20 и переработанная от 2026-01-21, у последней proj:epsg = 32640); • 12.11.2025, 07:29:29 UTC, облачность 16 %, высота Солнца 43,3° — ENMAP01-__L2A-DT0000164112_20251112T072929Z_020_V010505_… В окрестности (квадрат около 1°×1° вокруг участка, bbox 56,3–57,4 в.д. / 27,9–28,9 с.ш.) поиск вернул 50 продуктов L2A, включая безоблачные 01.01.2026 и серию 08.10.2025 — то есть район снимается регулярно, и есть смысл подать заявку на новую съёмку именно нашего прямоугольника.

Как получить. Каталог открыт без регистрации: https://geoservice.dlr.de/eoc/ogc/stac/v1/search?collections=ENMAP_HSI_L2A&bbox=56.82948,28.35330,56.84459,28.35710. Ассеты (SPECTRAL_IMAGE_COG.TIF — весь куб в Cloud-Optimized GeoTIFF, плюс QL_VNIR и QL_SWIR, маски облаков, теней, дымки, перистых, снега, дефектных пикселей, METADATA.XML) лежат на https://download.geoservice.dlr.de/ENMAP/files/L2A/… — но проверка HTTP-запросом показала, что прямая загрузка отдаёт HTML-страницу входа, то есть нужна бесплатная регистрация в DLR EOC Geoservice. Альтернативы: EOWEB GeoPortal (https://eoweb.dlr.de) — весь архив с обработкой по требованию до L1B/L1C/L2A; EO-Lab (https://eo-lab.org). Новая съёмка заказывается через Instrument Planning Portal (https://planning.enmap.org) — для этого нужна соответствующая роль пользователя и предложение (заявки Cat-1 проходят научное рецензирование).

Источники (5)

8. PRISMA (ASI) — доступ есть, но каталог закрыт логином: наличие сцен не проверено

PRISMA — итальянский гиперспектральный спутник (запуск 22 марта 2019 г., миссия продлена): VNIR 400–1010 нм (66 каналов), SWIR 920–2505 нм (173 канала), всего 239 гиперспектральных полос плюс панхроматический канал; спектральное разрешение не хуже 12 нм; гиперспектр 30 м, панхром 5 м; полоса захвата 30 км; повторный обзор за счёт кренов — 7 суток при цикле повторения орбиты 29 суток. Основная зона интересов — вся суша между 70° с.ш. и 70° ю.ш., до 200 000 км² в сутки, так что провинция Керман принципиально в зоне съёмки.

ЧЕСТНО О ГЛАВНОМ: наличие сцен PRISMA конкретно над участком проверить не удалось — открытого API или общедоступного каталога у ASI нет, поиск ведётся только из личного кабинета после регистрации и одобрения проектной формы. Публикаций по применению PRISMA в Иране достаточно (например, картирование гидротермальных изменений на медно-порфировом месторождении Daralu на юго-востоке Ирана; совместный анализ EMIT и PRISMA в металлогенической зоне Саве–Йезд), но это говорит только о том, что снимки над Ираном существуют, — про наш прямоугольник из этого ничего не следует.

Порядок получения. 1) Регистрация на https://prismauserregistration.asi.it/ (выбрать Institutional/Research или Commercial), указать аффилиацию, заполнить Project Form и подписать лицензию (LICENCE_TO_USE_PRISMA_DATA — запрет перераспределения, обязательная ссылка на ASI в публикациях). 2) После одобрения (обычно несколько дней) открывается доступ к порталу https://prisma.asi.it/ и системе CAS — там и каталог архива, и подача заявок на новую съёмку по заданной области. 3) Продукты уровней L0, L1 и L2 (L2B/L2C/L2D — геокодированное отражение поверхности). Доступ бесплатный для зарегистрированных пользователей.

Рекомендация: подать заявку на регистрацию сразу — даже если архива над участком нет, PRISMA принимает запросы на новую съёмку, а 30-километровая полоса при 62-гектарном участке гарантированно накроет его целиком.

Источники (6)

9. Sentinel-1 (SAR) — три трека, но повтор 12 суток, а не 6

Покрытие плотное. По Planetary Computer (sentinel-1-grd, поиск по точке): 951 сцена IW GRD за 12.10.2014 – 24.07.2026; поляризации VV+VH у 863 сцен, только VV у 88 (ранние годы). Три относительные орбиты: 93 (нисходящая, 324 сцены), 166 (нисходящая, 303), 57 (восходящая, 324). По CDSE над точкой 1891 продукт GRD и 941 продукт SLC — SLC и есть то, что нужно для InSAR.

ВАЖНОЕ ОГРАНИЧЕНИЕ ДЛЯ InSAR. Проверка по CDSE с разбивкой по аппаратам дала: S1A — 933 продукта SLC, S1B — 0, S1C — 0, S1D — 8 (с 2026 г.). Контрольный запрос показал, что S1C в каталоге присутствует (163 продукта над югом Германии за 2026 г.), то есть нуль над участком — это не сбой запроса, а реальное отсутствие съёмки. Проверка второй точкой (центральный Иран, 53° в.д., 32° с.ш.) тоже дала 0 продуктов S1C. Практический вывод: над участком исторически работал только один аппарат, поэтому повтор по каждому треку — 12 суток, шестисуточных пар (как в Европе в 2016–2021 гг.) здесь нет. Подтверждено эмпирически: интервалы между датами за 2022–2026 гг. по треку 93 — 128 раз ровно 12 суток, по треку 166 — 122 раза, по треку 57 — 114 раз. Совмещая три трека, снимок SAR получается каждые 3–5 суток, но для интерферометрии треки не складываются — геометрия разная.

Времена: нисходящие пролёты около 02:07 и 02:15 UTC (05:37 и 05:45 местного, раннее утро), восходящий около 14:18 UTC (17:48 местного).

Рельеф участка крутой (медиана уклона 16,9°, 38 % площади круче 20° — см. раздел о ЦМР), поэтому при выборе трека надо считать зоны наложения (layover) и радиотени; наличие и восходящего (57), и нисходящих (93, 166) треков позволяет закрыть оба склона.

Как получить: CDSE (бесплатная регистрация) — OData, STAC или S3, продукты SLC и GRD в формате SAFE; Alaska Satellite Facility (ASF) Vertex (https://search.asf.alaska.edu, бесплатная регистрация NASA Earthdata) — там же готовые RTC-продукты и сервис обработки HyP3 (интерферограммы «под ключ»); Planetary Computer — GRD и RTC как COG, без регистрации, через SAS-подпись.

Источники (3)

10. Цифровые модели рельефа — все основные наборы покрывают участок

Тайл во всех наборах — N28E056. Проверено не только наличие, но и фактическое содержимое: по каждой ЦМР прочитано окно участка и посчитаны высоты.

• Copernicus DEM GLO-30 — Copernicus_DSM_COG_10_N28_00_E056_00_DEM. Высоты над участком: 1597,5–1754,5 м, средняя 1664,5 м, перепад 157,0 м. Источник — интерферометрия TanDEM-X (X-диапазон), съёмка 2011–2015 гг.; заявленная абсолютная вертикальная точность менее 4 м (90 % LE), относительная — менее 2 м на склонах до 20° и менее 4 м на склонах круче 20°. Лицензия свободная (30-метровый набор открыт с 2023 г.). Это лучший бесплатный вариант для района: новее SRTM, без пустот, поверхность типа DSM. • NASADEM (переобработанный SRTM) — NASADEM_HGT_n28e056. Высоты 1595,0–1742,0 м, средняя 1663,7 м, перепад 147,0 м. Пустот в окне участка нет. Съёмка февраля 2000 г. • SRTM GL1 (1 угловая секунда) — гранула N28E056.SRTMGL1.hgt, есть в CMR. • ALOS AW3D30 — ALPSMLC30_N028E056_DSM. Высоты 1592,0–1764,0 м, средняя 1667,4 м, перепад 172,0 м. Формируется из стереопар PRISM/ALOS-1 2006–2011 гг. • ASTER GDEM v3 — ASTGTMV003_N28E056. Существует, но для инженерных задач уступает трём предыдущим (шум и артефакты стереообработки).

ХАРАКТЕР РЕЛЬЕФА (расчёт по GLO-30 в окне 56,825–56,850 в.д. / 28,348–28,362 с.ш., сетка 50×90 пикселей, размер пикселя 30,87 м по широте и 27,17 м по долготе): высоты 1526–1803 м, медианный уклон 16,9°, средний 17,6°, максимальный 42,7°; 38 % площади круче 20°, 10 % круче 30°. Это существенно: понадобятся топографическая нормализация оптических снимков (иначе тени и разная освещённость склонов замаскируют спектральные различия), поправка на рельеф при геометрической коррекции и учёт layover и радиотени при работе с SAR.

Как получить: Copernicus GLO-30 — Planetary Computer (cop-dem-glo-30, поиск STAC по точке, подпись ссылок), CDSE, AWS Open Data (s3://copernicus-dem-30m, открыт); NASADEM, SRTM и ASTER GDEM — NASA Earthdata (бесплатная учётная запись) либо Planetary Computer; ALOS AW3D30 — JAXA (бесплатная регистрация) либо Planetary Computer (alos-dem).

Источники (5)

11. Ночной тепловой канал: ASTER TIR и ECOSTRESS

Ночной TIR интересен тем, что убирает влияние текущей освещённости и подсвечивает различия по тепловой инерции — плотные ультраосновные породы, серпентиниты и рыхлые отложения остывают по-разному.

ASTER TIR (90 м, 5 полос 8,125–11,65 мкм). Над участком 33 ночные сцены AST_L1T и по 28 ночных гранул AST_05 (излучательная способность) и AST_08 (кинетическая температура). Ночные даты AST_L1T: 2005-05-13, 2005-08-17, 2008-11-22, 2015-04-07, 2015-07-12, 2015-07-21, 2016-05-27, 2016-08-15, 2016-11-19, 2017-01-22, 2017-10-21, 2017-11-22, 2017-12-08, 2017-12-17, 2018-08-14 (две сцены), 2019-01-05, 2019-03-01, 2019-07-16, 2020-05-22, 2021-05-09, 2022-05-12, 2022-10-31, 2022-12-08, 2023-01-01, 2023-07-02, 2023-12-27, 2024-02-06 (две сцены), 2024-05-08, 2024-05-25, 2025-10-15, 2025-10-30. Важно: TIR ASTER продолжает работать и после 2008 г. (в отличие от SWIR), так что весь ночной ряд пригоден. Время ночных пролётов ASTER: 17:10–18:26 UTC = 20:40–21:56 местного.

ECOSTRESS (МКС, 70 м, температура и излучательная способность поверхности). Плиточные продукты используют тот же MGRS-тайл 40RDS. Над участком: ECO_L2T_LSTE v002 — 864 гранулы (522 ночных, 342 дневных), первая 11.01.2019; ECO_L2T_LSTE v003 — 150 гранул (87 ночных, 63 дневных) с 26.10.2025 по 25.07.2026; ECO_L1CT_RAD v002 — 838 гранул (469 ночных, 369 дневных). Пример идентификатора: ECOv003_L2T_LSTE_41435_004_40RDS_20251026T202628_02. Благодаря нессинхронной орбите МКС ночные съёмки ECOSTRESS приходятся на самые разные часы (в выборке встречались 20:26, 00:38 и 02:14 UTC, то есть 23:56, 04:08 и 05:44 местного), что даёт возможность построить суточный ход температуры и оценить тепловую инерцию — по сути бесплатный аналог тепловой съёмки «перед рассветом».

Дополнительно: ASTER GED (AG100 v003 — 100 м, AG1km v003 — 1 км) даёт усреднённую по времени излучательную способность в 5 полосах TIR — удобная фоновая подложка для литологии; гранулы AG100.v003.29.056.0001 и AG1km.v003.29.056.0010 покрывают участок.

Как получить: всё через NASA Earthdata (CMR или Earthdata Search) с бесплатной регистрацией; плиточные продукты ECOSTRESS — Cloud-Optimized GeoTIFF по тайлам MGRS, ASTER — GeoTIFF по полосам.

Источники (3)

12. Что ещё нашлось, чего нет и что с регистрацией

ALOS-1 PALSAR (L-диапазон, JAXA и ASF). Поиск в ASF по точке участка дал 65 продуктов, среди них L1.1 (SLC, пригодны для InSAR) и L1.5, режимы FBS и FBD, а также готовые RTC низкого и высокого разрешения. L-диапазон глубже проникает сквозь растительность и щебень и лучше держит когерентность на длинных интервалах, чем C-диапазон Sentinel-1 — полезно для медленных смещений и структурного дешифрирования. Доступ — ASF Vertex, бесплатная регистрация NASA Earthdata.

DESIS (гиперспектрометр на МКС, DLR). Коллекция DESIS_HSI_L2A в STAC DLR проверена: над участком и в окрестности 1°×1° — 0 сцен.

Sentinel-1C. Над участком и над центральным Ираном — 0 продуктов (при заведомо рабочем каталоге, см. раздел 9).

СВОДКА ПО РЕГИСТРАЦИИ: • Без регистрации: STAC Earth Search (Sentinel-2 COG на AWS), STAC Planetary Computer (Landsat C2 L2, Sentinel-1 GRD и RTC, ЦМР — с подписью SAS), STAC DLR (метаданные EnMAP), CMR (метаданные всех наборов NASA), бакеты s3://sentinel-cogs и s3://copernicus-dem-30m. • Бесплатная регистрация: Copernicus Data Space Ecosystem (полный архив Sentinel-1 и -2, SAFE и SLC), NASA Earthdata Login (ASTER, EMIT, ECOSTRESS, NASADEM, SRTM, ASF), USGS ERS и EarthExplorer (Landsat, ASTER), DLR EOC Geoservice (файлы EnMAP), JAXA (AW3D30). • Регистрация с одобрением проекта: PRISMA (ASI), заявки на новую съёмку EnMAP (Instrument Planning Portal), USGS M2M API.

ФОРМАТЫ: COG и GeoTIFF (Landsat, EnMAP L2A, ECOSTRESS L2T, ЦМР, ASTER L1T по полосам), JPEG2000 в SAFE (Sentinel-2 из CDSE), SAFE (Sentinel-1), NetCDF-4 (EMIT), HDF-EOS и HDF5 (часть продуктов ASTER и PRISMA).

ПРОЕКТНАЯ СК. Все перечисленные источники приводятся к EPSG:32640 (UTM 40N) — той же зоне, в которой задан участок. Обратите внимание: в конфигурации ядра проекта strata (config/default.yaml) целевая СК сейчас EPSG:32642 (UTM 42N) — для этого участка её нужно переключить на EPSG:32640, иначе снимки будут перепроецироваться в чужую зону с лишними искажениями.

Источники (5)

13. Готовые запросы для воспроизведения проверки

Все числа выше получены запросами, которые можно повторить. Точка участка: lon = 56.83704, lat = 28.35520. Прямоугольник: bbox = 56.82948, 28.35330, 56.84459, 28.35710. Окно в UTM 40N: 483289.9961, 3136351.7277, 484770.0135, 3136770.6136.

1) Sentinel-2 L2A (без регистрации): POST на https://earth-search.aws.element84.com/v1/search с телом {"collections":["sentinel-2-l2a"],"intersects":{"type":"Point","coordinates":[56.83704,28.35520]},"datetime":"2024-01-01T00:00:00Z/2025-12-31T23:59:59Z","limit":100}

2) Облачность именно над участком: взять ассет scl каждой сцены (COG в s3://sentinel-cogs) и посчитать долю пикселей классов 3, 8, 9, 10 в окне UTM выше (rasterio.windows.from_bounds, переменная окружения AWS_NO_SIGN_REQUEST=YES).

3) Landsat: POST на https://planetarycomputer.microsoft.com/api/stac/v1/search с collections=["landsat-c2-l2"] и той же точкой; ссылки подписывать через https://planetarycomputer.microsoft.com/api/sas/v1/sign?href=…

4) ASTER до апреля 2008 (без регистрации, только метаданные): GET https://cmr.earthdata.nasa.gov/search/granules.umm_json?short_name=AST_L1T&point=56.83704,28.35520&page_size=200&sort_key=start_date — затем отфильтровать по TemporalExtent < 2008-04-01 и DataGranule.DayNightFlag = Day. Для готового отражения — short_name=AST_07XT.

5) EMIT: GET https://cmr.earthdata.nasa.gov/search/granules.umm_json?short_name=EMITL2ARFL&point=56.83704,28.35520&page_size=50 (и short_name=EMITL2BMIN для карт минералов).

6) EnMAP: GET https://geoservice.dlr.de/eoc/ogc/stac/v1/search?collections=ENMAP_HSI_L2A&bbox=56.82948,28.35330,56.84459,28.35710&limit=50

7) Sentinel-1 SLC (CDSE, метаданные без регистрации): GET https://catalogue.dataspace.copernicus.eu/odata/v1/Products?$filter=Collection/Name eq 'SENTINEL-1' and OData.CSC.Intersects(area=geography'SRID=4326;POINT(56.83704 28.35520)') and contains(Name,'SLC')&$count=true

8) ECOSTRESS с разделением день и ночь: GET https://cmr.earthdata.nasa.gov/search/granules.json?short_name=ECO_L2T_LSTE&version=002&point=56.83704,28.35520&day_night_flag=night&page_size=1 — число берётся из HTTP-заголовка ответа CMR-Hits.

9) ЦМР: STAC Planetary Computer, collections=["cop-dem-glo-30"], ["nasadem"] или ["alos-dem"], intersects=точка.

Источники (5)

Что дальше по этому направлению

  1. Переключить crs.target в config/default.yaml с EPSG:32642 на EPSG:32640 — участок лежит в UTM 40N; иначе весь конвейер модуля 1 будет перепроецировать снимки в чужую зону.
  2. Скачать безоблачную сцену EnMAP L2A от 30.12.2023 (ENMAP01-____L2A-DT0000055859_20231230T072248Z_002_V010502_20260121T012840Z) — завести бесплатную учётную запись DLR EOC Geoservice; это лучший из доступных гиперспектральных снимков на участок.
  3. Забрать все 14 гранул EMIT (EMITL2ARFL и EMITL2BMIN) через NASA Earthdata Login и сопоставить готовые минеральные карты EMIT со сценой EnMAP — независимая перекрёстная проверка спектральных признаков.
  4. Отобрать из 41 дневной сцены ASTER до 01.04.2008 те, у которых в EarthExplorer облачность менее 10 % (в CMR поля облачности нет, фильтровать нужно именно в EarthExplorer), и скачать соответствующие AST_07XT — crosstalk-коррекция и атмосферная коррекция там уже выполнены.
  5. Собрать безоблачный композит Sentinel-2 по маске SCL, отобрав сцены с нулевой облачностью над участком (их 128 за 2024–2025) — база для отношений полос и для приведения остальных источников к общей геометрии.
  6. Подать заявку на регистрацию PRISMA (prismauserregistration.asi.it), проверить архив по координатам участка и, если сцен нет, оформить запрос на новую съёмку — полоса 30 км гарантированно накроет 62 га.
  7. Для InSAR выбрать один трек Sentinel-1 (93 или 166 нисходящие либо 57 восходящий), выкачать стек SLC из CDSE или собрать интерферограммы через ASF HyP3; рассчитывать на 12-суточный интервал, не на 6-суточный.
  8. Посчитать по Copernicus GLO-30 карты layover и радиотени для каждого из трёх треков Sentinel-1 и решить, какой склон каким треком снимать; там же подготовить маску уклонов круче 30° (10 % площади) для топографической нормализации оптики.
  9. Построить суточный ход температуры по ночным гранулам ECOSTRESS (522 ночных в v002 плюс 87 в v003) и ночным AST_08 (28 гранул) — оценка тепловой инерции для разделения плотных ультраосновных пород и рыхлых отложений.
  10. Оформить полученные идентификаторы сцен (Sentinel-2, Landsat, ASTER до 2008 г., EMIT, EnMAP, Sentinel-1) как реестр в стиле strata/module1_open_sources/sources.py, чтобы конвейер приёма работал по конкретным сценам, а не по общим правилам.