[Рецепт] Зональная статистика в ArcGIS 10

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

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

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

Можно хитрее, сгенерировать одно покрытие на миллион точек, 1 раз извлечь данные, а потом будет просто делать в R столько, сколько надо:

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

sample(rnd, 10000, replace = FALSE)

Проблема - извлечение статистики по миллиону полигонов. Если ID последовательный от 0 до 999999, то ArcGIS после пары часов работы обрадует следующим:

zonal-error-arcgis10.gif (9.46 КБ) 6645 просмотров

Код ошибки бесполезный, а вот сообщение To many unique values - наталкивает на мысли.

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

Тогда все ок:

zonal-error-arcgis10-ok.gif (7.2 КБ) 6645 просмотров

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

[SS_Rebelious]

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

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

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

[novia]

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

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

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

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

[novia]

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

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

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

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

[novia]

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

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