Карта продуктивности (NDVI)

[athe]

Доброго времени суток уважаемые форумчане.
Появилась необходимость делать так называемые "карты продуктивности" на основе индекса NDVI. В итоге из растра необходимо получить мультиполигон, который разбит на зоны соответственно определённому диапазону значений NDVI растра (0-0,2;0,21-0,3...). Пример к сожалению не очень удачный, но в общем передаёт суть вопроса:

Я пробовал это делать через Saga -> Raster values to points, после чего по правилу выделять строки атрибутики и объединять, но это занимает чудовищную уйму времени.
Подскажите пожалуйста каким альтернативным образом возможно выполнить эту задачу?
Заранее благодарен.
С уважением, athe.

[_taras_]

Если правильно понял задачу, то можно сделать так

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

(("NDVI"<0)*1)+
(("NDVI"=0)*2)+
((("NDVI">0) AND ("NDVI"<0.2))*1)+
((("NDVI">=0.2) AND ("NDVI"<0.3))*3)+
((("NDVI">=0.3) AND ("NDVI"<0.6))*4)+
((("NDVI">=0.6) AND ("NDVI"<0.8))*5)+
(("NDVI">=0.8)*6)

Результат

[athe]

_taras_, добрый день. Спасибо за ответ.
Как я понимаю указанный вами код необходимо использовать в "field calculator"? Я попробовал через "обновление существующего поля" и в итоге NDVI пересчитался в значения от 1 до 5, а 300к+ строк в атрибутике как было- так и осталось. Я явно делаю что-то не так).

[_taras_]

Нет, я использовал калькулятор растров (растр --- калькулятор растров).

[bim2010]

Module Reclassify Grid Values
http://www.saga-gis.org/saga_tool_doc/2 ... ls_15.html
https://gis.stackexchange.com/questions ... es-in-qgis

[antonv]

athe, идея _taras_ была в том, чтобы сначала переклассифицировать значения в самом растре, а уже затем его векторизовать.

Если вы уже получили векторизованные полигончики с нужными значениями, в меню Vector > Geoprocessing > Dissolve их можно объединить по нужному атрибуту. Тогда они станут мультиполигонами, по одному на каждое уникальное значение.

[athe]

_taras_,antonv, спасибо, это сработало. Значения в растре пересчитались в целые, что дало возможность использовать "Polygonize". Вышло вот что:

Ещё раз спасибо.

[airmadina-1]

good afternoon. Thanks for the answer.
As I understand it, the code you specified should be used in the "field calculator"? I tried through "updating an existing field" and as a result, NDVI was converted to values ​​from 1 to 5, and 300k + lines in the attributes as it was, it remains. I'm obviously doing something wrong).

[Ethel11-]

If you have already received vectorized polygons with the desired values, in the Vector> Geoprocessing> Dissolve menu you can combine them by the desired attribute. Then they will become multipolygons, one for each unique value.