Форумы GIS-Lab.info

Поделюсь, пока не забыл, возможно кому-то пригодится.
Допустим есть задача, 100 раз повторить случайную выборку из растра, скажем по 10 тыс буферизованных точек-полигонов.

Можно попрограммировать на питоне, можно вручную сгенерировать 100 покрытий и 100 раз повторить операцию извлечения - не очень хорошее решение.

Можно хитрее, сгенерировать одно покрытие на миллион точек, 1 раз извлечь данные, а потом будет просто делать в R столько, сколько надо: Проблема - извлечение статистики по миллиону полигонов. Если ID последовательный от 0 до 999999, то ArcGIS после пары часов работы обрадует следующим:
Код ошибки бесполезный, а вот сообщение To many unique values - наталкивает на мысли.

Решается очень неочевидно, поэтому следим за пальцами:
Customize\ArcGIS options\Raster\Maximum number of unique values to render: (1000001)

Тогда все ок: Тоже самое можно сделать в плагине для зональной статистике в QGIS, проблема в том, что за 4 часа в которые уложился ArcGIS, QGIS смог продвинутся только на 31%.

А как называется плагин для зональной статистики в QGIS?

Zonal statistics, плагин ядра на C++, через установщик плагинов искать не нужно

Решается очень неочевидно, поэтому следим за пальцами:
Customize\ArcGIS options\Raster\Maximum number of unique values to render: (1000001)

Баг отправлен в Esri GeoProcessing Team

и считается похоже неправильно, увы, выборочная проверка показывает что некоторые цифры не соответствуют реальным

Поделитесь примером, если не жалко, который а) ломает работу инструмента б) подсчитан неправильно. Укажите места, где неправильно. Так будет проще и быстрее объяснить в чем баг. Можно в личку.

[Сообщение с мобильного устройства]

http://gis-lab.info/share/million.rar
если подтвердится что ломает, то попробую найти где неправильно, а то может у меня конфигурация какая-то не такая.

Попробую сломать на arcgis 10.0 sp2 @ win7 x64, результат сообщу.

[Сообщение с мобильного устройства]