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

Методика дешифрирования

Кратко.

Для участка в дехестане Соган методика ДЗЗ существует и проверена практически «на месте»: комплекс Соган прямо разбирался в работе Pournamdari, Hashim & Pour (Advances in Space Research 54: 694–709, 2014), а смежная работа тех же авторов (Resource Geology 64: 233–246, 2014) охватывает три комплекса Эсфандаге — Abdasht, Soghan, Sikhoran — и даёт рабочий рецепт: специализированное отношение бандов ASTER RGB = 4/1, 4/5, 4/7 плюс MNF и SAM с объединением признаков, цель — серпентинизированный дунит как порода-хозяин хромититов. Оманская школа (Rajendran, Nasir и др., массивы Семаил) добавляет проверенные приёмы: decorrelation stretch, композит бандов 8-4-1, композит отношений 4/7, 4/1, (2/3)×(4/3), SAM на лизардит/антигорит, картирование Мохо и переходной зоны как хромитоносного горизонта, валидацию полевым спектрометром PIMA и XRD; заявленная позиционная точность выделений 50–150 м. Физическая основа всего этого — полоса Mg-OH серпентина 2,32 мкм, которая попадает точно в банд 8 ASTER (2,295–2,365 мкм); ни Sentinel-2 (B12 центр 2202,4 нм, ширина 179 нм), ни Landsat OLI её не разрешают. Главное честное ограничение: САМ ХРОМИТИТ СО СПУТНИКА НЕ ДЕТЕКТИРУЕТСЯ — его собственные полосы (0,49; 0,59; 0,69; 0,94; 1,30; 2,00; 2,30 мкм по Cloutis 2004 и полевой спектрометрии) слишком узкие и слабые для мультиспектра, а на 30–60-метровом пикселе гиперспектра подиформное тело размером в десятки метров тонет в смеси с вмещающим дунитом; авторы иранских работ признают это прямо. Реально картируются: серпентинит, степень серпентинизации, дунит-гарцбургитовая мантийная серия, листвениты и коры выветривания — то есть поисковые индикаторы, а не руда. Дополнительные ограничения для конкретно этого полигона: ASTER SWIR пригоден только до 1 апреля 2008 г.; площадь 62 га при 30-метровом пикселе — всего ~690 пикселей (при 90-метровом TIR — ~80, тепловые индексы неприменимы); участок узкий (419 м = ~14 пикселей), нужен буфер 2–3 км; район горный, поэтому топографическая коррекция (SCS+C или Minnaert по DEM 30 м) обязательна, а не желательна; действующая добыча даёт отвалы и карьеры как источник ложных целей. Реалистичный ориентир по достигаемой точности — 72–80 % общей точности классификации (Mawat 73→79,3 %, Афганистан 72,2 %); по иранским офиолитам матриц ошибок вообще не публикуют, а численных порогов (углы SAM, отсечки по отношениям бандов) нет ни в одной проверенной работе — калибровать придётся локально.

1. Что установлено для САМОГО участка и его комплекса (Соган)

Участок попадает в комплекс Соган — один из 31 ультраосновного комплекса восточно-западного пояса Хаджиабад–Эсфандаге на юге провинции Керман. Административно комплекс Соган относится именно к шахрестану Урзуийе (Arzuiyeh County), то есть к тому же району, что и заданный полигон. Для этого комплекса существует прицельная работа по ДЗЗ: Pournamdari, Hashim & Pour (2014), «Spectral transformation of ASTER and Landsat TM bands for lithological mapping of Soghan ophiolite complex, south Iran» (Advances in Space Research, т. 54, с. 694–709). В ней применены коэффициент корреляции (CC), фактор оптимального индекса (OIF), отношения бандов (BR) и метод главных компонент (PCA) к ASTER и Landsat TM. Смежная работа тех же авторов (Resource Geology, 2014, т. 64, с. 233–246) охватывает уже три комплекса Эсфандаге — Abdasht, Soghan и Sikhoran — и формулирует главный вывод: цель картирования не сам хромитит, а серпентинизированный дунит как порода-хозяин хромититовых тел. Отдельная работа по Abdasht: Pournamdari & Hashim (2013), Arabian Journal of Geosciences, DOI 10.1007/s12517-013-0927-0. ВАЖНО: ни в одной из этих работ мы не нашли ни точных ID сцен, ни границ обработанных площадей, ни численных оценок точности (матрицы ошибок, kappa) — то есть покрытие «нашего» полигона следует считать районным, а не подтверждённым попиксельно. Локальная геология по независимому источнику (Sepidbar et al., Sikhoran–Soghan, препринт/ORA Oxford): подиформные хромититы обрамлены тонкими дунитовыми оболочками; шпинель лерцолитов–гарцбургитов Cr# = 42,0–52,6, шпинель дунитов и высокохромистых хромититов Cr# = 52,4–76,4; степень частичного плавления >17 об.%; двухпироксеновый термометр даёт 879–912 °C. Хромититы ведут себя как жёсткие тела и обрастают серпентином и карбонатами — это и есть спектрально видимый признак.

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

2. Физическая основа: какие полосы поглощения реально работают

Ключ ко всей методике — понимание, ЧТО именно поглощает. Для офиолитов работают три группы признаков. (а) Mg-OH серпентина: узкая асимметричная полоса ~2,32 мкм (комбинационная Mg-OH), V-образная ~2,52 мкм, обертон OH ~1,38–1,39 мкм и слабая мелкая полоса ~2,12 мкм. Это главный и практически единственный надёжный дистанционный признак серпентинизированных ультрабазитов. (б) Железо: Fe3+ в видимом (переносы заряда и кристаллическое поле, ~0,48–0,55 и ~0,65–0,90 мкм), Fe2+ силикатов — широкая полоса около 1,0–1,05 мкм у оливина и ~0,9 и ~1,8–2,0 мкм у ортопироксена. (в) Карбонаты и тальк в лиственитах: CO3 около 2,30–2,35 мкм плюс кремнезём в тепловом ИК. Собственные полосы ХРОМИТА: 0,49; 0,59; 0,69; 0,94; 1,30; 1,75; 2,00 и 2,30 мкм (Spectral Evolution, полевая спектрометрия), в лабораторном исследовании Cloutis, Sunshine & Morris (2004) диагностическими названы 0,46; 0,49; 0,55; 0,59; 0,66; 0,69; 0,93; 1,3 и ~2 мкм. Критерий разделения: центр интенсивной полосы Fe2+ около 2 мкм у шпинели всегда короче 2,1 мкм, у хромита — всегда длиннее 2,1 мкм. Добавка 5 мас.% хромита к оливину лабораторно заметна. Но всё это получено на порошках <45 мкм при спектральном разрешении 5 нм и лабораторной геометрии (i=30°, e=0°) — не со спутника.

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

3. Рецепты ASTER, проверенные на офиолитах Ирана

Ниже — конкретные, воспроизводимые сочетания, каждое привязано к иранскому офиолиту. (1) ЭСФАНДАГЕ/СОГАН/АБДАШТ/СИКХОРАН (Pournamdari et al., 2014): специализированное отношение бандов ASTER RGB = 4/1, 4/5, 4/7 — выделяет серпентинизированный дунит как породу-хозяина хромититов; дополняется компонентами MNF (ASTER и Landsat TM) для регионального разделения офиолитовых блоков и SAM с последующим объединением признаков (feature-level fusion). Логика: 4/1 ловит контраст SWIR/видимого у тёмных ультрабазитов, 4/5 и 4/7 — падение отражения в 2,15–2,29 мкм у гидроксил-содержащих фаз. (2) НЕЙРИЗ (Eslami, Ghaderi, Rajendran, Pour, Hashim, 2015, Resource Geology): decorrelation stretch RGB по Landsat TM 7-5-4; PC1/PC2/PC3 по SWIR-бандам ASTER; отношения бандов RGB = 5/3, 5/1, 7/5 — уверенно разделяют большинство пород и, что важно, выделяют переходную зону (transition zone) и мантийный гарцбургит как зоны с высокой хромитоносностью; SAM по VNIR+SWIR ASTER для минералов. (3) ХОЙ, СЗ Иран (Mehdikhani & Imamalipour, 2021, Minerals 11:960): из 63 трёхбандовых комбинаций максимум OIF у бандов 3, 6, 8 → композит RGB 8-6-3; лучшие композиты отношений — RGB ((4+2)/3, (7+5)/6, (9+7)/8), RGB (4/1, 4/7, 4/5) и RGB (4/3 × 2/3, 3/4, 4/7); PCA по схеме Crosta/Amer, в т.ч. PCA(4,5,2); SAM по библиотечным спектрам, пересчитанным на банды ASTER; атмосферная коррекция — log residual; данные ASTER L1B 2002 г. (4) САБЗЕВАР, СВ Иран (Eskandari, Hosseini, Nicotra, 2023, Minerals 13:251): интеграция ASTER + Landsat + Sentinel-2 с change detection, band rationing и SAM; ультрабазиты разбиты на 4 группы по степени серпентинизации через среднюю спектральную яркость; серпентинит с предполагаемым дунитовым протолитом (опознаётся по Fe-богатой Ni-содержащей латеритной корке) — самый продуктивный на хромит, особенно на контактах с мафическими дайками; проверка орто-мозаикой с БПЛА GSD ~6 см на площади 35 км² и последующим полевым опробованием — найдены новые проявления хромита.

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

4. Рецепты ASTER, проверенные на офиолитах Омана (Семаил)

Оманская школа (Rajendran, Nasir и соавторы, Университет султана Кабуса) — самый обширный набор проверенных рецептов по офиолитам аридной зоны. (1) Rajendran et al. (2012), Ore Geology Reviews 44: 121–135, «ASTER detection of chromite bearing mineralized zones in Semail Ophiolite Massifs»: хромитовые залежи приурочены к базальному дуниту и гарцбургиту мантийной секции; применены decorrelation stretch, различные отношения бандов и PCA по VNIR и SWIR; результат — картирование серпентинизированного гарцбургита, вмещающего хромититы; полевая заверка показала пространственную точность выделения 50–150 м. Это единственная найденная нами количественная оценка позиционной точности по офиолитам Омана. (2) Rajendran & Nasir (2014), Ore Geology Reviews 62: 211–226, Вади Хиби: использованы спектральные поглощения OH и Mg-OH серпентина; цветной композит RGB по бандам ASTER 8, 4, 1 показывает выветрелые перидотиты; гидротермально изменённые зоны выделены композитом отношений RGB = 4/7, 4/1, 2/3 × 4/3 (это классический «рецепт Sultan», см. раздел 5); лизардит и антигорит разделены методом SAM качественно и количественно. (3) Rajendran & Nasir (2015), Tectonophysics 657: 63–80: картирование Мохо и переходной зоны Мохо (MTZ) в массивах Семаил по ASTER; авторы прямо указывают, что изображения способны картировать хромитоносную зону и дуниты MTZ; спектральные свойства пород контролировались полевым спектрометром PIMA, минералы подтверждены рентгенофазовым анализом (XRD). Для нашей задачи это ключевой методический ориентир: цель — не рудное тело, а MTZ/дунитовая оболочка. (4) Rajendran et al. (2013), Ore Geology Reviews: лиственитовые (listwaenite) зоны Центрального Омана по ASTER. (5) Rajendran & Nasir (2019), Ore Geology Reviews 108: 33–53 — сводный обзор возможностей ASTER по минеральным ресурсам Омана. (6) Гиперспектральный эталон: Roy R., Launeau P., Carrère V., Pinet P., Ceuleneer G., Clénet H. и др. (2009), G-Cubed, авиасъёмка HyMap над массивом Сумаил: выделены три эндмембера мантийной серии (гарцбургит, дунит, гарцбургит с интерстициальным карбонатом) и четыре коровых (массивный габбро, амфиболизированное габбро с дайками, верлит с высоким содержанием серпентина, габбронорит); характерные узкие Mg-OH-полосы дают именно оливин-богатые породы с высоким содержанием серпентина; плагиоклаз в этом диапазоне спектрально «немой» и не детектируется; авторы делят данные на мантийный и коровый поднаборы по наличию Mg-OH и работают синтетической библиотекой всех сочетаний пород в пикселе с оценкой доли покрытия.

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

5. Справочные индексы и отношения бандов ASTER (первоисточник — Geoscience Australia)

Чтобы не пересказывать рецепты по цитатам, приводим формулы из первоисточника — «ASTER Mineral Index Processing Manual», составители Aleks Kalinowski и Simon Oliver, Remote Sensing Applications, Geoscience Australia, октябрь 2004 (PDF GA7833, выкачан и разобран). Железо: Fe3+ = 2/1 (Rowan; CSIRO); Fe2+ = 5/3 + 1/2 (Rowan); латерит = 4/5 (Bierwith); госсан = 4/2 (Volesky); железистые силикаты (биотит, хлорит, амфибол) = 5/4 (CSIRO); окислы Fe3+ = 4/3 (CSIRO, «может быть неоднозначно»). Карбонаты и мафические минералы: карбонат/хлорит/эпидот = (7+9)/8 (Rowan); эпидот/хлорит/амфибол = (6+9)/(7+8) (CSIRO); амфибол/MgOH = (6+9)/8 (Hewson, с прямой оговоркой «может быть либо MgOH, либо карбонат»); амфибол = 6/8 (Bierwith); доломит = (6+8)/7 (Rowan, USGS); карбонат в TIR = 13/14. Кремнезём: кварцевые породы = 14/12; кремнезём = (11×11)/(10×12) (Nimoyima); индекс основности (гранат, клинопироксен, эпидот, хлорит) = 12/13; SiO2 = 13/12. Готовые RGB-композиты: «Sultan» = 4/7, 4/1, (2/3)×(4/3); «Abrams (USGS)» = 4/7, 4/3, 2/1; «госсан/изменение/вмещающая порода» = 4/2, 4/5, 5/6 (Volesky). Тепловые индексы Ninomiya (2004; Ninomiya, Fu, Cudahy, 2005): индекс кварца QI = (b11 × b11)/(b10 × b12); карбонатный индекс CI = b13/b14; мафический индекс MI при разделении карбонатов n = 3: MI = (b12 × b14³)/b13⁴; эмиссивность растёт на коротких волнах (b12) у ультрамафитов и на длинных (b13) — у фельзических пород. Индексы QI и MI устойчивы к вариациям температуры поверхности при хороших атмосферных условиях; CI напрямую зависит от температуры и требует нормировки радиояркости. Индексы специально «под хромитит» (Дхарвар, Индия, Advances in Space Research 56: 1453, 2015): Chromitite Index = (band2 × band2)/(band1 + band3), Ultramafic Index = band1/band3 — основаны на том, что у хромитита отражение в банде 2 ASTER выше, чем в бандах 1 и 3. ВАЖНО: последние два индекса опираются только на VNIR-контраст (по сути на форму низкоальбедной кривой), проверены в Индии, на офиолитах Ирана/Омана нам подтверждения не встретилось.

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

6. Sentinel-2, Landsat и WorldView-3: что можно и чего нельзя

Главный ограничитель — ширина и положение SWIR-полос. У Sentinel-2 MSI всего две SWIR-полосы: B11 (S2A: центр 1613,7 нм, ширина 88 нм) и B12 (S2A: центр 2202,4 нм, ширина 179 нм), обе 20 м. B12 покрывает примерно 2113–2292 нм, то есть диагностическая полоса Mg-OH серпентина 2,32 мкм лежит ЗА краем полосы. Поэтому Sentinel-2 в принципе не различает Al-OH (2,20) и Mg-OH (2,32) — он видит только «есть какое-то поглощение в 2,1–2,3 мкм». То же у Landsat 8/9 OLI (SWIR2 ~2,11–2,29 мкм). Отсюда: Sentinel-2 и Landsat годятся для морфоструктуры, контрастов альбедо, оксидов железа и растительной маски, но не для минералогического разделения серпентинита от карбонатов и глин. Практика: в китайском офиолите Шибанжин применяли совмещённое отношение (S3 + A9)/(S12 + A8), где S — банды Sentinel-2A, A — банды ASTER, то есть Sentinel-2 использовали как «добавку разрешения» к ASTER, а не самостоятельно. ASTER принципиально лучше: банд 8 (2,295–2,365 мкм) точно попадает на Mg-OH 2,32 мкм, банды 5–7 — на Al-OH 2,20 мкм. Отсюда все рабочие рецепты (6+9)/8, (6+9)/(7+8), 4/7, 4/5. Но ASTER SWIR мёртв: с 1 апреля 2008 г. аномалии в SWIR сделали данные непригодными, и SWIR не включается в продукт AST_L1T; данные после апреля 2008 г. показывают аномальное насыщение и полосатость. То есть для нашего участка пригодны только архивные сцены ASTER 2000–начала 2008 гг. Радикальное решение по разрешению — WorldView-3: 8 SWIR-полос с GSD 3,7 м, из них SWIR-7 2228–2292 нм (Al-OH) и SWIR-8 2285–2373 нм (Mg-OH) — прямой аналог пары ASTER b7/b8, но в 8 раз детальнее. На нашем участке 1480 × 419 м это даёт ~400 × 113 пикселей SWIR вместо ~49 × 14 у ASTER. Подтверждённый пример: Yin et al. (2024), Economic Geology 119(3): 665, картирование мелких (<100 м) мафит-ультрамафитовых тел по WV-3 с композитами трёх отношений бандов, новым SWIR-индексом и глубоким обучением (spatial-spectral transformer); в поле подтверждены 7 перспективных участков, в 4 из них найдены халькопирит, пентландит, пирит и ХРОМИТ. WorldView-3 — коммерческие данные (Maxar), бесплатного архива нет.

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

7. Гиперспектр 2019–2026: PRISMA, EnMAP, EMIT над ультрамафитами

Гиперспектр — единственный класс спутниковых данных, который в принципе может подойти к собственным полосам хромита и к разделению дунит/гарцбургит. Характеристики: PRISMA (ASI, запуск 22.03.2019) — VNIR 400–1010 нм (66 каналов) + SWIR 920–2505 нм (173 канала), всего 240 полос, спектральное разрешение ≤12 нм, GSD 30 м (панхром 5 м), полоса 30 км, повтор 29 сут (перенацеливание 7 сут), данные бесплатны для зарегистрированных пользователей. EnMAP (DLR, запуск 01.04.2022) — 420–2450 нм, до 228 полос (99 VNIR + 163 SWIR), шаг 6,5 нм в VNIR и 10 нм в SWIR, GSD 30 × 30 м, полоса 30 км, SNR 175:1 на 2200 нм, повтор в надире 27 сут, при отклонении 30° — 4 сут. EMIT (NASA JPL, на борту МКС) — 380–2500 нм, спектральный шаг 7,4 нм. Что доказано на нашем типе объектов: (а) антигорит устойчиво выделяется в данных EMIT и PRISMA по Mg-OH — широкое сложное поглощение в области 1050–1100 нм плюс сильная полоса 2320–2330 нм; (б) офиолитовый серпентинит успешно разделён по данным EnMAP (первые применения EnMAP для геологического картирования Арабо-Нубийского щита); (в) по Ирану есть прямые прецеденты работы PRISMA: Kuhpanj (юг Ирана, порфировая медь) — общая точность выделения зон изменений 71,87 % (SPCA), 87,50 % (отношения бандов), 81,25 % (MTMF) и 71,87 % (LSU) при каппа 58,85 / 80,54 / 70,55 / 57,52 %; Daralu (ЮВ Иран, провинция Керман) — субпиксельное картирование минералов изменения по PRISMA; (г) EMIT + PRISMA совместно применены в металлогенической зоне Саве–Йезд (центральный Иран). Практическое следствие для участка: гиперспектр 30 м даёт всего ~49 × 14 пикселей на полигон, но зато перекрывает провал ASTER между 0,86 и 1,60 мкм, где лежит диагностическая полоса оливина ~1,05 мкм — именно она отличает дунит от гарцбургита. Ни ASTER, ни Sentinel-2, ни Landsat этой полосы не видят вообще.

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

8. Методы классификации и достигнутые точности на офиолитах

Порядок применения методов в проверенных работах устойчив: сначала контрастирование (band ratio, decorrelation stretch, OIF-подбор композита), затем ортогонализация (PCA по схеме Crosta или MNF), затем спектральное сопоставление (SAM), затем при необходимости субпиксельные методы (MTMF, LSU, CEM) и обучаемая классификация (SVM, RF, XGBoost, глубокое обучение). Достигнутые точности (все — по офиолитам или мафит-ультрамафитам): Mawat, Иракский Курдистан (Othman & Gloaguen, 2014, Remote Sensing 6: 6867) — SVM по ASTER: общая точность 73 % только по спектральным признакам и 79,3 % при добавлении геоморфометрических признаков и текстур; тем же исследованием методом SAM найдены два новых участка с высокой концентрацией хромита, один заверен в поле — тело ~0,3 км² с крупнокристаллическим хромитом в дунитовой вмещающей породе, содержание Cr2O3 ~8,46 %. Афганский Центральный блок (Ahmadi et al., 2023, Minerals 13: 1453) — FCC + собственные отношения бандов + PCA + SAM, общая точность 72,2 %; выделено 4 офиолитовых комплекса. Иран, Kuhpanj по PRISMA (не офиолит, но тот же класс задач и та же страна) — от 71,87 % до 87,50 % общей точности в зависимости от алгоритма, лучший результат у отношений бандов. Оман, Семаил (Rajendran et al., 2012) — численной точности классификации нет, но полевая заверка дала позиционную ошибку выделения 50–150 м. Иранские работы по Эсфандаге, Соган, Абдашт, Нейриз и Хой численных матриц ошибок НЕ приводят вовсе; их «валидация» — согласие с существующей геологической картой и маршрутные наблюдения. Это системная слабость всей литературы по ДЗЗ офиолитов: точность 72–80 % — реалистичный ориентир, а публикации без матрицы ошибок надо читать как качественные. Отдельная ветвь — сочетание со спутниковой геофизикой и БПЛА: подиформный хромит практически немагнитен (восприимчивость хромита 2·10⁻⁷…1,9·10⁻³ SI, оманские образцы 3,1·10⁻⁴…3,7·10⁻⁴ SI), тогда как серпентинизированные перидотиты Омана 5·10⁻⁴…7,8·10⁻³ SI — на порядок выше, поэтому магнитометрия тоже картирует не руду, а серпентинизацию (карточка проекта № 111, Parvar et al., БПЛА-магнитометрия в Самаиле: аномалия 146 нТл на высоте 20 м над известным выходом хромита исчезает уже на 60 м).

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

9. Порядок предобработки (в горном рельефе — критичен)

Рекомендуемая последовательность, собранная из проверенных работ. ШАГ 1 — выбор данных. Для ASTER брать только сцены до 1 апреля 2008 г. (после — SWIR непригоден). Оптимально — готовый продукт AST_07XT (ASTER L2 Surface Reflectance VNIR and Crosstalk Corrected SWIR): в нём уже применены коррекция кроссталка и атмосферная коррекция; VNIR 15 м, SWIR 30 м. ШАГ 2 — коррекция кроссталка, если работаете с L1B/L1T. Кроссталк — паразитная засветка от детектора банда 4 в соседние детекторы бандов 5 и 9; она создаёт ложные полосы поглощения и искажает диагностические сигнатуры, что ведёт к прямым ошибкам определения минералов. Алгоритм — Iwasaki & Tonooka (2005); программа коррекции распространялась ERSDAC/GDS. ШАГ 3 — атмосферная коррекция. В иранских работах по офиолитам применялся log residual (Эсфандаге/Соган, Хой) — метод простой, но параметров авторы не приводят. Современная практика — FLAASH (в работе по офиолитам Ласбелы, Пакистан, применена и к ASTER, и к Hyperion: DN → радиояркость → FLAASH → отражение) или готовый L2-продукт. ВАЖНО: log residual и другие «внутренние» методы (IARR, flat field) убирают и часть полезного сигнала; для сопоставления с библиотечными спектрами (USGS, RELAB, ECOSTRESS) нужна физическая коррекция. ШАГ 4 — топографическая коррекция. Для нашего участка это не опция, а обязательный шаг: район горный, перепады рельефа создают ложные контрасты, которые band ratio подавляют лишь частично (отношения бандов частично сокращают мультипликативный эффект освещённости, но не тени и не эффект неламбертовости). По сравнительным испытаниям 10 методов топокоррекции на снимках ASTER в сложном горном рельефе (Teillet regression, VECA, Cosine, Cosine-B/C/T, C, SCS, SCS+C, Minnaert) неламбертовские методы устойчиво лучше ламбертовских; в одних испытаниях лучшим оказался Minnaert-E, затем Minnaert и SCS+C, в других — SCS+C, а Minnaert худшим. Практический вывод: применять SCS+C или Minnaert с эмпирическим k, обязательно сравнивать «до/после» по корреляции отражения с cos(i), и проверять, что коррекция не переисправила теневые склоны. DEM: Copernicus DEM 30 м или ALOS AW3D30; при GSD 15 м ASTER VNIR желателен DEM не грубее 30 м. ШАГ 5 — кросс-калибровка при совместном использовании сенсоров (ASTER + Sentinel-2 + Landsat + гиперспектр): приведение к общей сетке (для проекта — EPSG:32640 на этом участке, при том что конфиг проекта сейчас настроен на EPSG:32642 — расхождение надо снять до обработки), совместная регистрация с точностью лучше половины пикселя, ресемплинг библиотечных спектров под функции спектрального отклика каждого сенсора (шаг 13 руководства Geoscience Australia). ШАГ 6 — маски: растительность (NDVI = (b3−b2)/(b3+b2) для ASTER), вода, тени, облака, а также техногенные объекты (карьеры, отвалы, дороги) — на действующем хромитовом поле отвалы дают спектр «свежей» породы и создают ложные цели. ШАГ 7 — уменьшение размерности (MNF) до расчёта SAM/MTMF, отбор компонент по кривой собственных значений. ШАГ 8 — эталонные спектры: снимать полевым спектрорадиометром (ASD FieldSpec / PIMA, как в оманских работах) по обнажениям и керну, подтверждать минералогию XRD, только потом строить SAM/MTMF.

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

10. Честное разделение: что реально детектируется, а что нет

ДЕТЕКТИРУЕТСЯ УВЕРЕННО (мультиспектр ASTER, при исправной предобработке): серпентинизированные ультрабазиты — по поглощению Mg-OH 2,32 мкм в банде 8, через отношения (6+9)/8, (6+9)/(7+8), 4/5, 4/7 и композиты 4/1, 4/5, 4/7; окислы и гидроксиды железа (кора выветривания, «железные шляпы») — 2/1, 4/2, 4/3; листвениты (карбонатизированные ультрабазиты) — по карбонатным полосам 2,30–2,35 мкм и по TIR; разделение «ультрамафит — мафит — фельзит» по тепловому мафическому индексу Ninomiya MI (но 90 м пиксель). ДЕТЕКТИРУЕТСЯ ЧАСТИЧНО/КОСВЕННО: дунит против гарцбургита — по мультиспектру только косвенно (через степень серпентинизации, альбедо, морфологию и цвет коры выветривания); по гиперспектру — реально, через полосу оливина ~1,05 мкм и полосы ортопироксена ~0,9 и 1,8–2,0 мкм; переходная зона (MTZ) как хромитоносный горизонт — картируется по совокупности признаков (Rajendran & Nasir 2015, Eslami et al. 2015). НЕ ДЕТЕКТИРУЕТСЯ СО СПУТНИКА: сам хромитит по своим диагностическим полосам. Причины: (1) собственные полосы хромита 0,49; 0,59; 0,69 мкм узкие и слабые, требуют спектрального шага порядка 5–10 нм; у ASTER в этой области всего три широких банда (0,52–0,60; 0,63–0,69; 0,78–0,86 мкм), у Sentinel-2 — несколько узких, но не в тех местах; (2) хромитит очень тёмный, отражение низкое, отношение сигнал/шум в глубоких полосах мало; (3) даже гиперспектр (PRISMA ≤12 нм, EnMAP 6,5/10 нм, EMIT 7,4 нм) на 30–60 м пикселе смешивает подиформное тело размером десятки метров с вмещающим дунитом — нужна субпиксельная обработка (MTMF, LSU, CEM), которая даёт долю эндмембера, а не факт руды; (4) все лабораторные калибровки Cloutis получены на порошках <45 мкм при разрешении 5 нм и лабораторной геометрии — перенос на природную шероховатую поверхность с патиной не обоснован. Прямо это признают и авторы иранских работ: Mehdikhani & Imamalipour пишут, что хромит спектрально напрямую не детектируется и картируется индикатор — дунит/серпентинит. ДОПОЛНИТЕЛЬНЫЕ ЛОВУШКИ для нашего участка: (а) размер полигона 62 га — при ASTER SWIR 30 м это ~690 пикселей, при TIR 90 м всего ~80 пикселей: тепловые индексы здесь почти бессмысленны; (б) действующий горнодобывающий район — отвалы, карьеры и дороги дают «свежие» спектры ультрабазитов и создают ложные цели; (в) участок вытянутый (1480 × 419 м), при 30-метровом пикселе ширина всего ~14 пикселей — краевые эффекты фильтров и MNF будут заметны, обрабатывать надо с буфером не менее 2–3 км вокруг; (г) полоса Mg-OH может давать и хлорит, и тальк, и амфибол, и карбонат — отношение (6+9)/8 официально помечено в руководстве Geoscience Australia как «может быть либо MgOH, либо карбонат».

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

11. Готовый рабочий конвейер для участка (сводка)

A. ДАННЫЕ. Обязательные: ASTER AST_07XT (сцена до 01.04.2008) — единственный источник, где банд 8 накрывает Mg-OH 2,32 мкм в открытом доступе; Sentinel-2 L2A (10/20 м, свежие даты, для морфоструктуры, маски растительности, оксидов железа и мультивременного контроля); Copernicus DEM 30 м (топокоррекция, морфометрия, уклоны — они же признаки для SVM по образцу Mawat). Желательные: PRISMA и/или EnMAP L2A (заявка на съёмку участка; перекрывают провал 0,86–1,60 мкм и дают полосу оливина 1,05 мкм); EMIT (проверить наличие гранул над участком). Опционально при бюджете: WorldView-3 SWIR 3,7 м — единственный способ получить Mg-OH-канал на масштабе тел десятки метров. B. ПРЕДОБРАБОТКА: см. раздел 9 (кроссталк → атмосфера → топокоррекция SCS+C/Minnaert по DEM → общая сетка EPSG:32640 → маски). C. РАСЧЁТ ПРИЗНАКОВ по ASTER: композит-«якорь» RGB (4/1, 4/5, 4/7) — рецепт, проверенный именно на Эсфандаге/Соган; контрольные композиты RGB (5/3, 5/1, 7/5) — по Нейризу, и RGB ((4+2)/3, (7+5)/6, (9+7)/8) — по Хою; индексы: MgOH/амфибол (6+9)/8, эпидот-хлорит-амфибол (6+9)/(7+8), карбонат/хлорит/эпидот (7+9)/8, Fe3+ 2/1, Fe2+ 5/3+1/2, госсан 4/2; композит Sultan 4/7, 4/1, (2/3)×(4/3) — на листвениты, по оманскому опыту; TIR — QI, CI, MI(n=3), но только как региональный фон. D. КЛАССИФИКАЦИЯ: MNF → отбор компонент → SAM по эталонным спектрам (серпентин/антигорит/лизардит, тальк, магнезит, доломит, хлорит, гётит; ресемплировать под функции отклика сенсора) → MTMF/CEM для субпиксельных долей → SVM или RF с добавлением морфометрии (по образцу Mawat, где это дало +6,3 процентных пункта к общей точности). E. ВАЛИДАЦИЯ: полевой спектрорадиометр по обнажениям и керну + XRD (оманский стандарт), матрица ошибок и каппа по независимой выборке точек — то, чего нет в иранских публикациях и что надо сделать самим. F. ПОРОГИ: готовых числовых порогов ни в одной из проверенных работ по офиолитам Ирана и Омана не опубликовано — ни углов SAM, ни отсечек по отношениям бандов. Пороги придётся калибровать локально по обучающим полигонам (известные выходы серпентинита, дунита, гарцбургита, габбро, лиственита и отвалов), а не переносить из литературы.

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

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

  1. Устранить расхождение систем координат: участок задан в EPSG:32640 (UTM 40N), а crs.target в config/default.yaml проекта — EPSG:32642; выбрать одну и зафиксировать до запуска конвейера
  2. Проверить архив ASTER (USGS EarthExplorer / NASA Earthdata) на наличие безоблачных сцен AST_07XT над 28,355° с.ш., 56,837° в.д. за период 2000 — 31.03.2008 и выписать granule ID; без такой сцены весь SWIR-блок методики неприменим
  3. Проверить наличие гранул EMIT и сцен PRISMA/EnMAP над участком; при отсутствии — подать заявку на съёмку PRISMA (бесплатно для зарегистрированных) и EnMAP (научная лицензия)
  4. Реализовать в strata/module1_open_sources расчёт набора индексов ASTER как отдельный модуль: (6+9)/8, (6+9)/(7+8), (7+9)/8, (6+8)/7, 2/1, 5/3+1/2, 4/2, 4/3, 4/5, 4/7 и композиты RGB (4/1, 4/5, 4/7), (5/3, 5/1, 7/5), ((4+2)/3, (7+5)/6, (9+7)/8), 4/7, 4/1, (2/3)×(4/3)
  5. Добавить в конвейер шаг топографической коррекции (SCS+C и Minnaert с эмпирическим k) по Copernicus DEM 30 м с контролем корреляции отражения и cos(i) до/после
  6. Собрать спектральную библиотеку эталонов под функции отклика ASTER/Sentinel-2/PRISMA: серпентин (лизардит, антигорит, хризотил), тальк, магнезит, доломит, хлорит, актинолит, гётит, гематит, оливин, ортопироксен — из USGS Spectral Library, ECOSTRESS и RELAB
  7. Попытаться получить полные тексты закрытых ключевых работ легальным путём (Pournamdari 2014 Adv Space Res и Res Geol, Rajendran 2012 и 2015, Rajendran & Nasir 2014 и 2019) — через межбиблиотечный абонемент или запрос авторам; именно там могут быть неопубликованные в аннотациях численные детали
  8. Завести в проекте карточки по недостающим темам: спектральные признаки лиственитов (Rajendran et al. 2013), картирование MTZ (Rajendran & Nasir 2015), WorldView-3 для мелких мафит-ультрамафитовых тел (Yin et al. 2024), PRISMA в Кермане (Kahnouj & Tangestani 2024)
  9. Заложить в план полевых работ измерение спектров обнажений и керна полевым спектрорадиометром (350–2500 нм) с обязательным подтверждением минералогии XRD — без этого SAM/MTMF настраивать не на чем
  10. Построить собственную матрицу ошибок и каппа по независимой выборке точек — этого нет ни в одной иранской публикации по офиолитам, и это будет реальным приращением к состоянию вопроса