Форумы GIS-Lab.info

подскажите, как сделать следующую вещь:
на входе есть растр (очень большой по экстенту, к примеру вся россия, разрешение 1km), в нем значения: 1 и nodata. Ячеек равных 1 может быть от 15 до 5000000 (то есть от очень мало, до очень много).
есть еще N-гридов (любой другой формат тоже могу сделать - хоть многоканальный тиф) такого же экстента и разрешения
Нужно получить таблицу в которой для каждой ячейки равной 1 первого растра, будут присутствовать значения всех N-гридов в этой точке (наподобие инструмента extract grid values)
что смущает:
1) combine обязательно отругает за количество гридов (более 20) и ячеек
2) если конвертировать первый растр в точечную тему и экстрагировать по точкам, то экстракция будет занимать пару недель (да и в целом - нужна только таблица, поэтому еще и компьютерные силы\время тратить на сам шейп как-то не хочется)

что бы можно было бы сделать? любым доступным инструментом (не только аркгис\аркинфо)?

Думаю арчик справится с этим легко и непринужденно.

ох, а как же арчик справится с понятием пространственной выборки?
ну-ка поподробнее о библиотеке пожалуйста:)
и вообще - какое ограничение будет на размер растра к примеру подгружаемому через библиотеку гдал? а то у меня здоровучие многоканальные данные...

выборка задается растром, где 0 и 1?

Я для начала попробовал бы просто: а потом выкинул бы в растр rgdal'ом.

data и ingrid какой формат из себя представляют? какой библиотекой читались в R и из какого формата? в виде чего хранятся при вот такой вот операции? какая библиотека подключена была для вот такой вот обработки?
поподробнее пожалуйста!:)

это была примерная схема, если опыта с растрами в R нет, то первым делом читаем и тренируемся вот на этом:
http://gis-lab.info/qa/rgdal.html

спасибо - буду пробовать:) по результатам - отпишусь