Распаковка QA для данных MODIS

[Максим Дубинин]

Наверняка не я один вот таким озаботился? Может кто уже делал все это?

[Максим Дубинин]

Сделал перекодировочные таблицы для MCD15A2, прикладываю.

Было бы неплохо собрать их для всех продуктов MODIS.

[Александр Мурый]

Было бы неплохо собрать их для всех продуктов MODIS.

А как оно делается, если не секрет?

[Максим Дубинин]

1. Выбираем продукт, например MCD15A2
2. Читаем про него доку в части таблиц QA (там где Layers)
3. Составляем опорную таблицу для распаковки целочисленных значений (8,16,32 бит) в битовые строки и битовые наборы (типа как у меня в примере, оно само считает).
3. Собираем из тайлов глобальную мозаику и смотрим уникальные значения
4. Распаковываем все возможные целочисленные значения и смотрим какие комбинации вообще возможны.

В общем результат такой, что можно по целочисленному значению без распаковки "прочитать" что же происходит с пикселем.

Например для MCD15A2 видим 105, а на самом деле это:
MODLAND_QC: Other Quality, Sensor: Terra, DeadDetector: Detectors OK, CloudState: Significant clouds, SCF_QC: Emprical, not geometry

[Ariki]

Более канонично через битовые операции, будет как-то так в LibreOffice:

Код: Выделить всё

=IF(BITAND(A13; BITLSHIFT(1;1)); "Aqua"; "Terra")

Может, имеет смысл приложение для QGIS сделать?

[Максим Дубинин]

это я пробовал, с вариантом 0/1 оно еще норм, а с двух-трехбитными начинаются многоступенчатые вложенные IF-ы, которые становятся совершенно нечитаемыми и хочется сразу забить на Calc и переделать все на питоне.

Про QGIS не уверен, это же десктопное приложение, тут скорее нужна какая-то библиотечка где собраны все "расшифрованные" варианты QA для всех продуктов.

Пока сам не придумал, как это лучше оформить и нужно ли это вообще.

Но то что люди боятся QA и не хотят его использовать из-за этих сверток, это факт.