Новости космоса (МКС, Марса, Луны), спутники и их запуски

Убираем черные полосы в снимках Landsat 7 +ETM

Узнаем, как убрать черные полосы (дефект CLS-off) у снимках Landsat 7 +ETM с помощью ENVI, Erdas.

Недавно меня попросили сделать мозаику тернопольской области территории Украины за август 2008 года по космическим снимкам с КА ДЗЗ Landsat 7 +ETM.

В принципе ничего сложно здесь нет. Выкачал бесплатно снимки за указанный период из сервиса «EarthExplorer» и как я думал в любом ПО (а я привык к ArcGis, ENVI, Erdas) без проблем справлюсь с созданием мозаики. Но я ошибался…

(Между прочим, если дорогие читатели нужно сделать инструкцию по скачиванию космических снимков с КА ДЗЗ Landsat 4, 5, 7, 8 и Sentanal 2 с сервиса «EarthExplorer» пишите в комментариях)

Дело в том, что все космические снимки, сделанные с КА Landsat 7 +ETM после 31 мая 2003 года имеют известный дефект CLS-off (черные полосы), снимки имеют расширяющиеся от центра к краям черные полосы, где изображение отсутствует. Эти "пропуски" (gap-ы) существенно затрудняют использование снимков в автоматическом дешифрировании, и несколько мешают при ручном.

Благодаря тому, что при повторном пролете gap-ы не пересекаются, возможно изготовление композитов из 2-3 снимков одной территории или с помощью методов интерполяции заполнять пропуски соседними пикселями. С одной стороны, при классификации территории или расчете вегетационных индексов мы вносим погрешности, но для создания мозаики конкретной территории смотрится даже очень ничего.

Убираем черные полосы в снимках Landsat 7 +ETM

Убираем черные полосы (SLC correction) в снимках Landsat 7 с помощью ERDAS

В основу этого алгоритма заложено циклическое повторение операции Focal analysis для устранения артефактов и заполнение пропущенных данных (значений равно 0). Операция должна повторятся до тех пор, пока, черные полосы не будут устранены.

По факту, на месте черных полос, данные отсутствуют, значение пикселей равно 0. Для устранения этих пробелов мы грубо говоря проведем операцию интерполяции с учетом соседних пикселей, которые содержат реальные значения.

Для каждой ячейки, значение которой равно 0, вычисляется значение на основе оценки соседних вокруг нее ячеек. В стандартном случае соседние ячейки ищутся в окне 3х3.

В моем варианте, создан инструмент, который позволяет автоматизировать этот процесс, а также после второго прохождения меняет окно поиска (kernel) с 3х3 до 5х5, что придает более плавные переходи.

Закачать инструмент.

Этот инструмент создавался для версий ERDAS 10 и выше.

Отблагодари автора, поставив лайк или оставив положительный коммент.

Убираем черные полосы в снимках Landsat 7 +ETM с помощью ENVI

Для этого нам понадобится специальный скрипт landsat_gapfill.sav. Скопируйте содержимое архива в папку "C:\Program Files\Exelis\ENVI5[minor version]\extensions", для версии 5.1 путь будет следующим: «C:\Program Files\Exelis\ENVI51\extensions».

Описание работы инструмента

Комментарии:

Новое на сайте

Biarri Cubesat - является четырехсторонним проектом с участием Австралии, США, Великобритании и Канады для...
CXBN 1, 2 (Cosmic X-Ray Background NanoSat) – небольшие КА предназначенные для наблюдения за...
CSUNSat 1 – небольшой КА созданный в Калифорнийском государственном университете в партнерстве с Лабораторией...
IceCube (Earth 1) – небольшой КА НАСА для измерения частиц льда в облаках.
Altair Pathfinder (Altair 1) – небольшой КА типа CubeSat 6U частной компании Millennium Space...
Tianzhou 1, 2 (TZ 1, 2) - первый китайский грузовой космический корабль, обслуживающий один Tiangong...
Все спутники

Популярные материалы

Архив новостей

Апрель 2017
ПнВтСрЧтПтСбВс
272829303112
3 4 5 6
  • Космические аппараты (спутники) Индии Oceansat 3, 3A
  • Космические аппараты (спутники) Индии CartoSat 3, 3A, 3B
789
10
  • Космические аппараты (спутники) Индии RISAT 2A
  • Космические аппараты (спутники) Индии RISAT 2
11 12 13 14 15 16
1718 19 20 21 2223
24 252627282930

Ваш email:

Присоединяйся к нам в