Выборка значений из нескольких грид
Добавлено: 06 дек 2012, 10:58
Возник вопрос, можно ли выбрать из нескольких грид наименьшие/наибольшие значения? Например, есть 5 грид на одну территорию, надо составить из них один интегральный грид, в котором будут только наименьшие значения из пяти.
Если можно, опишите последовательность действий. Есть ArcView 3 и QGIS.
Пока алгоритм вижу следующий:
1. Строим маску по первому грид с помощью картографического калькулятора, т.е. выделяем области со значениями ниже чем во 2-5 грид;
2. Присваиваем в маске значениям с истиной - 0, с ложью - любую страшную цифру, которую легко будет в дальнейшем выявить;
3. Плюсуем грид-1 и маску для него, обнуляем ячейки со страшными цифрами, и у нас в результирующем гриде (грид-1м) остаются те значения, которые нам нужны.
4. повторяем пункт 1-3 для 2, 3, 4 и 5 грид;
5. Суммируем грид-1м, грид-2м, грид-3м, грид-4м и грид-5м, получаем то что требовалось в самом начале.
Оно работает, но жутко непродуктивно, нельзя ли как-то это упростить?
Ну и еще вопрос вдогон. В NASA LPDAAC Collection есть данные Modis с уже рассчитанными индексами NDVI и EVI, чем они хороши/плохи? Насколько я понимаю, дело упирается в недостаточно высокое разрешение (250 и 500 метров)? С другой стороны, эти индексы рассчитаны по одной методе с одинаковых сенсоров (если брать с 2000 по 2012 год). Есть ли какие-то подводные камни (помимо разрешения), делающие использование этих индексов нецелесообразным?
Если можно, опишите последовательность действий. Есть ArcView 3 и QGIS.
Пока алгоритм вижу следующий:
1. Строим маску по первому грид с помощью картографического калькулятора, т.е. выделяем области со значениями ниже чем во 2-5 грид;
2. Присваиваем в маске значениям с истиной - 0, с ложью - любую страшную цифру, которую легко будет в дальнейшем выявить;
3. Плюсуем грид-1 и маску для него, обнуляем ячейки со страшными цифрами, и у нас в результирующем гриде (грид-1м) остаются те значения, которые нам нужны.
4. повторяем пункт 1-3 для 2, 3, 4 и 5 грид;
5. Суммируем грид-1м, грид-2м, грид-3м, грид-4м и грид-5м, получаем то что требовалось в самом начале.
Оно работает, но жутко непродуктивно, нельзя ли как-то это упростить?
Ну и еще вопрос вдогон. В NASA LPDAAC Collection есть данные Modis с уже рассчитанными индексами NDVI и EVI, чем они хороши/плохи? Насколько я понимаю, дело упирается в недостаточно высокое разрешение (250 и 500 метров)? С другой стороны, эти индексы рассчитаны по одной методе с одинаковых сенсоров (если брать с 2000 по 2012 год). Есть ли какие-то подводные камни (помимо разрешения), делающие использование этих индексов нецелесообразным?