search

RS-Lab / SAT-Viewer

Веб-приложение каталога спутниковых снимков
0 stars 0 forks source link

Комбинации каналов #4

Open derevnya opened 6 years ago

derevnya commented 6 years ago

MODIS 1-4-3 Естественные цвета 7-2-1 Искусственные цвета 3-6-7 Искусственные цвета NDVI Вегетационный индекс (2-1)/(2+1) T Температура поверхности по каналу 31

derevnya commented 6 years ago
  1. Что делать с ночным modis? Там у нас только 21,31,32.
  2. Нужно сделать список по viirs. Пока добавил I3,I2,I1, как аналог для 7,2,1. По естественным у меня записана комбинация I1,M1,M2. M2 нет.
Vtirko commented 6 years ago

  1. По ночному MODIS-у есть только температура поверхности. Пожарную 21-31-32, розово-белую с элементами голубизны отменили. Она и правда невзрачная. Осталась одна температура — чёрно-белая и раскрашенная в палитру.

  2. VIIRS I1-M4-M3 Естественные цвета (наиболее близкие к MODIS). M11-I2-I1 Искусственные цвета (наиболее близкие к MODIS). M3-I3-M11 Искусственные цвета (наиболее близкие к MODIS). NDVI Вегетационный индекс (I2-I1)/(I2+I1) T Температура поверхности по каналу I5.

Июнь пересчитывается с согласованным набором каналов — май или остальное пересчитывать?

derevnya commented 6 years ago

Пока пересчитывать не надо. Для теста пусть останутся пока в старом варианте.

derevnya commented 6 years ago

Какая у нас комбинация используется для ночных снимков в "rgb"? Появилась мысль сделать что-то типа канала естественных цветов для ночных снимков.

Vtirko commented 6 years ago

Для ночных снимков используется для MODIS: 21, 31, 32, для VIIRS: i4, i5, i5. Контраст всех каналов линейный (гамма == 1) в диапазоне 235-300 К радиационной температуры. Ночью же все кошки серыe? Какие такие "естественные цвета" можно для них сделать?

yakubailik commented 6 years ago

Перепутаны названия комбинаций каналов VIIRS: 'Пожары'<-->'Снег и лед'. Должно быть так: Пожары (M11-I2-I1) Снег и лед (M3-I3-M11)

Vtirko commented 6 years ago

В соответствие с совещанием от 30.08.2018 — ночной чёрно-белый вариант должен иметь негатив: 220К - белое, 300К - чёрное, и гамму 0.5.

Vtirko commented 6 years ago

В соответствие с совещанием от 30.08.2018 вариант rgb квиклуков идёт с разрешением 1500-2000 м и чёрно-белым ночным представлением. Пока пересчитывается с 21.08.2018.

derevnya commented 6 years ago

Считаю LUT. Пиксель у нас считается так: (((pixel/65536)^(1/gamma))*65536)/256. Получил что-то такое: LUT=1:0,7500:77,15000:113,30000:166,65535:255

Правильно ли я понимаю, что для ночных формула остается и все что меньше 22000 = 255 (65535), все что больше 30000 = 0. Т.е. считаем гамму, инвертируем и потом обрезаем? Или в обратном порядке?

Vtirko commented 6 years ago

Да, правильно. Гамма больше 1 единицы растягивает тёмное, а меньше — сжимает. В данном случае тёмное в исходном — это облака. Мы их сжимаем гаммой меньше единицы и инвертируем — получаются сжатые яркие облака и растянутая тёмная поверхность.

derevnya commented 6 years ago

Для ночных правильно оказалось так: растягиваем 220000-300000 в 0-255 (65535), считаем гамму и инвертируем.

День из rgb: 20180820050635-aqu-500-rgb День из multi: 20180820050635-aqu-1000 Ночь из rgb: 20180820183249-aqu-500-rgb Ночь из multi: 20180820183249-aqu-1000

Vtirko commented 6 years ago

Всё отлично! Тогда я 500-750 м RGB удаляю и RGB схемы вообще выключаю? 250-375 уже удалил.

derevnya commented 6 years ago

Общий скрипт генерации png переделаю, можно будет отключить. Напишу тут.

derevnya commented 6 years ago

Скрипт переделал, все обзорные изображения пересоздал.

Vtirko commented 6 years ago

Отключаю генерацию RGB-схемы.

derevnya commented 5 years ago

Вопрос по Sentinel и набору каналов. Аналогично с modis и viirs. Пытаюсь сделать комбинации. Скорректируйте таблицу: 4-3-2 Естественные цвета (у modis 1-4-3) ?-?-? Искусственные цвета (у modis 7-2-1) ?-?-? Искусственные цвета (у modis 3-6-7) NDVI Вегетационный индекс (8-4)/(8+4) (у modis (2-1)/(2+1)) Температуры нет.

derevnya commented 5 years ago

Гамму такую же для sentinel используем? И для rgb используем комбинацию 4-3-2?

Vtirko commented 5 years ago

4-3-2 Естественные цвета (у modis 1-4-3) 12-8-4 Искусственные цвета (у modis 7-2-1) 2-11-12 Искусственные цвета (у modis 3-6-7) NDVI Вегетационный индекс (8-4)/(8+4) (у modis (2-1)/(2+1)) Температуры нет.

11 канал добавил в конфигурацию (его не было). Данные сегодня пересчитаются вместе с исправлением ошибки в TIFFах.

Vtirko commented 5 years ago

Да, гамма не зависит от прибора, если каналы на тех же или близких длинах волн.

derevnya commented 5 years ago

Для rgb наверное правильнее будет эта "12-8-4 Искусственные цвета (у modis 7-2-1)".

Vtirko commented 5 years ago

Да, конечно — 12-8-4 — она самая цветастая круглый год.

derevnya commented 5 years ago

Статья про синтез: art13_3_8.pdf

Vtirko commented 5 years ago

Аналоги каналов Landsat-8 (OLI и TIRS): 1-OLI-1 == MODIS-9 2-OLI-2 == MODIS-3 3-OLI-3 == MODIS-4 4-OLI-4 == MODIS-1 5-OLI-5 == MODIS-2 6-OLI-6 == MODIS-6 7-OLI-7 == MODIS-7 8-OLI-8 == MODIS-1 + MODIS-4 - панхром 9-OLI-9 == MODIS-26 11-TIRS-1 == MODIS-31 12-TIRS-2 == MODIS-32

Аналоги каналов Канопус МСС: МСС-1 == MODIS-3 МСС-2 == MODIS-4 МСС-3 == MODIS-1 МСС-4 == MODIS-2

derevnya commented 5 years ago

В метаданных у ландсата 6 каналов OLI и 1 TIRS. Если делать по аналогии с модисом , то не хватает 7го канала.

derevnya commented 5 years ago

Для канопусов названия каналов надо сделать в метаданных в общем виде. Вместо MSS1, MSS2 и т.д. просто 1,2 и т.д. Для PSS аналогично. Плюс проверить проекцию.

derevnya commented 5 years ago

С Канопусам теперь все хорошо. По ландсатам вопрос остался. Вариант, который в Модисе "1-4-3 Естественные цвета" выглядит темноватым. У Ландсата это 4-3-2. Может что-то другое взять? Плюс с "11-TIRS-1 == MODIS-31" не работает. У Ландсата этот канал отличается по значениям?

derevnya commented 5 years ago

Добавил еще комбинаций, естественные цвета можно оставить как есть. Остался вопрос с температурой.

Vtirko commented 5 years ago

Нет, не отличается. Я проверил цифровые значения — на мой взгляд всё в порядке. image image На этих двух снимках TIRS Landsat-8 и MODIS Terra за 07.08.2019 температура самого южного залива Красноярского моря имеет разницу в 2 градуса 28900/29100.

derevnya commented 5 years ago

Да, я нашел ошибку и поправил. Теперь все хорошо. Остался только ресурс.

Vtirko commented 5 years ago

Вроде как закончил обработку Ресурса-П.

Несколько итогов:

Соответствие каналов с MODIS довольно странное:

derevnya commented 5 years ago

Что касается 8 и 16 бит, есть подозрение, что 8 бит это только 2017 год. 2017 год это растры из Железногорска. 2017 год это было разовое явление. Остальные пункты для всего годов встречаются?

Vtirko commented 5 years ago
derevnya commented 5 years ago

Картинки сильно светлые получаются при комбинации 4-3-1...