Распаковка информации о качестве данных MODIS: различия между версиями
(→Теория) |
|||
| Строка 7: | Строка 7: | ||
==Расчеты "на коленке"== | ==Расчеты "на коленке"== | ||
В качестве коленки будет использовать LibreOffice Calc. Вероятно те же самые формулы заработают и в MS Excel, но автор этой статьи это не проверял. | ===План=== | ||
В качестве "коленки" будет использовать LibreOffice Calc. Вероятно те же самые формулы заработают и в MS Excel, но автор этой статьи это не проверял. | |||
Наша задача - распаковать целочисленные коды из QA в человекочитаемый формат. Действовать будем по следующему плану: | |||
#Получим список уникальных значений из всех растров чтобы не пропустить чего-то полезного, но редкого. | |||
#Переведем | |||
===Реализация=== | |||
Проблема в том, что по умолчанию фукнция DEC2BIN в Calc не работает с числами вне диапазона -512..512 (2^10 - один бит на четность). Поэтому, если битов качества больше чем 10, формулу придется модифицировать следующим образом, например для 16 bit: | Проблема в том, что по умолчанию фукнция DEC2BIN в Calc не работает с числами вне диапазона -512..512 (2^10 - один бит на четность). Поэтому, если битов качества больше чем 10, формулу придется модифицировать следующим образом, например для 16 bit: | ||
<pre>=DEC2BIN((MOD(A2,65536)/512),7)&DEC2BIN(MOD(A2,512),9)</pre> | <pre>=DEC2BIN((MOD(A2,65536)/512),7)&DEC2BIN(MOD(A2,512),9)</pre> | ||
==Результаты== | |||
Здесь можно скачать готовые перекодировочные таблицы для следующих продуктов MODIS: | |||
*MCD15A2 - Leaf Area Index - Fraction of Photosynthetically Active Radiation 8-Day L4 Global 1km | |||
*MOD13A2 - Vegetation Indices 16-Day L3 Global 1km | |||
*MOD17A2 - Gross Primary Productivity 8-Day L4 Global 1km | |||
Версия от 18:00, 8 декабря 2015
заготовка
Теория
MODIS - камера дистанционного зондирования на борту спутников Terra и Aqua, снимающая каждую точку на Земле два раза в день и производящая десятки террабайт данных ежедневно. Данные MODIS проходят интенсивную обработку, включающую контроль качества. Помимо всего прочего, большое количество проблем доставляют облака, которые необходимо определенным образом маскировать. Поэтому все продукты на базе данных MODIS включают слой качества (он же QA=Quality Assessment, QC=Quality Control).
Это статья посвящена тому, как распаковать данные MODIS QA в человекочитаемый вид.
Расчеты "на коленке"
План
В качестве "коленки" будет использовать LibreOffice Calc. Вероятно те же самые формулы заработают и в MS Excel, но автор этой статьи это не проверял.
Наша задача - распаковать целочисленные коды из QA в человекочитаемый формат. Действовать будем по следующему плану:
- Получим список уникальных значений из всех растров чтобы не пропустить чего-то полезного, но редкого.
- Переведем
Реализация
Проблема в том, что по умолчанию фукнция DEC2BIN в Calc не работает с числами вне диапазона -512..512 (2^10 - один бит на четность). Поэтому, если битов качества больше чем 10, формулу придется модифицировать следующим образом, например для 16 bit:
=DEC2BIN((MOD(A2,65536)/512),7)&DEC2BIN(MOD(A2,512),9)
Результаты
Здесь можно скачать готовые перекодировочные таблицы для следующих продуктов MODIS:
- MCD15A2 - Leaf Area Index - Fraction of Photosynthetically Active Radiation 8-Day L4 Global 1km
- MOD13A2 - Vegetation Indices 16-Day L3 Global 1km
- MOD17A2 - Gross Primary Productivity 8-Day L4 Global 1km