получение значений растров в таблицу

[Анна]

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

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

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

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

[Анна]

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

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

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

Я для начала попробовал бы просто:

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

vals = vector()
for (i in 1:length(data)) {
 if (ingrid != 0) {
    val = outgrid[i]
    vals = c(vals,val)
 }
}

а потом выкинул бы в растр rgdal'ом.

[Анна]

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

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

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

[Анна]

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