выборка значений грида гридом

[Анна]

Добрый вечер,
подскажите пожалуйста как можно сделать следующую операцию:
есть два грида - 90м и 1км. Тот который 90м имеет либо data, либо nodata. Нужно выбрать все ячейки 1км грида, в которые попадает хотя бы 1 (и более) ячейки 90м грида со значением data. Те где только nodata тоже должны превратиться в nodata ( в 1км гриде)
спасибо!

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

[KolesovDmitry]

подскажите пожалуйста как можно сделать следующую операцию:
есть два грида - 90м и 1км. Тот который 90м имеет либо data, либо nodata. Нужно выбрать все ячейки 1км грида, в которые попадает хотя бы 1 (и более) ячейки 90м грида со значением data. Те где только nodata тоже должны превратиться в nodata ( в 1км гриде)

А какое программное обеспечение будете использовать? В GRASS, например, это делается одной командой, а qGIS имеет доступ к GRASSовским базам.
Если не привязываться к ПО, то общая идея, наверное, такая:
* установить разрешение километрового грида равное 90 м. (Во многом ПО этот шаг можно опустить - он будет проделан автоматически)
* Создать грид-маску по 90-метровому гриду: там, где было nodata установить значение ячейки в 0, инача в 1.
* Умножить грид-маску на исследуемый грид.
* Вернуться в разрешение 1 км.

[Анна]

Дмитрий, спасибо!
да, именно такой логике я и следовала - только в arcinfo
и еще много вариантов (типа combine, docell, con и т.д.)
но что мне остается абсолютно непонятным это результат (см картинку) - почему-то часть ячеек того грида который 1 км так и остается не выбранной - притом любым методом
итак делала я следующим методом:
precip5 - 1км дем
sp1 - 90м слой-маска
1) resample 1 км в 90 м (precip90)
2) isnull = isnull (precip90 * sp1) /* это чтобы получить маску на всю территорию, так как экстент sp1 много меньше, чем precip5. Команда дает на выходе 1 для nodata, 0 - для data
3) cond1 = con (isnull < 1 ,1,0) /* переклассификация, чтобы там где дата - стал 1. где nodata - 0
4) mult = multiple (precip90 * cond1) - дает перемножение precip90 на 0 и 1 маску шагу 3
5) resample (mult, 1 km) - обратно ресамплинг к 1 км

[KolesovDmitry]

Возможно, дело в методе resample. Посмотрите результирующий растр, который получается после преобразования километровой сетки к 90-метровой. Подозреваю, что ячейки "изчезают" уже на precip90.

[Анна]

в том то и дело что исчезают на последнем шаге - из 90 м в 1 км

[KolesovDmitry]

Странно, но все равно похоже, что дело в конкретной реализации resample (в том, как он обрабатывает ячейки, заполненные не полностью)... Попробуйте эксперимента ради тогда 90-метровый растр привести к километровому, а потом умножать.

[Анна]

Воющем в результате решение придумала такое:
взять тот грид который 1 км по нему установить setwindow snapgrid
потом аггретировать aggregate'ом 90 м грид в 1 км
и дальше уже вести обсчеты на основе аггрегированного грида. Работает:) первая команда позволила сохранять сетку грида единообразной.
Если кому-то интересно - могу все команды прописать

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