Форумы GIS-Lab.info

Здравствуйте! Такая проблема: постоил кригинг содержания (клопов на кв. км.) в модуле Geostat. analyst.
Теперь нужно вытащить полученные данные в Statistica. Модуль действительно отдает таблицу рассчитанных содержаний, но только для тех точек, по которым расчет и строился. А мне нужно получить каким-то макаром, массив значений с мелким шагом, шажком, можно сказать, по всей области оценивания. Как же это сделать?

дык а наверное построить сетку и присвоить узлам атрибуты областей? фишнет, грид аналиста в шейп и интерсект, неизящно, но я бы так сделала

мм.. сорри, я работаю в GA, но с остальным аркгис я не кутлху.. мм.. ктулху, нет , с остальным аркгис я не иван сусанин ) мне нужно поточнее, а точнее:
1. построить сетку (могу найти сам, но все же..)
1-а. присвоить узлам атрибуты областей / вы не поверите как долго я уже думал, что-то подобное сотворить, но не знал как )
2. путем выполнения 1. и 1-а = получаем, грид ?
2-а. фишнет и фишинг это что-то близкое?
3. грид аналиста в шейп и интерсект - ура, это понятно!

ваще все не так, как у Соколова в Уральских пельменях
1) Строим сетку - это инструмент create fishnet
2) преобразуем грид, который получился в результате кригинга, в шейп (если он еще не в шейпе, я никогда не делала кригинг)
3) делаем интерсект фишнета этим шейпом. по-моему, узлам фишнета должны присвоиться атрибуты шейпа, но надо бы проверить на практике, так с ходу не соображу

При крикинге по моему получается грид.
в сетку данные можно извлечь и из грида - arctoolbox\spatial anlyst \extraction (для этого вам нужны будут узлы сетки в виде точек) или зональной статистикой zonal (это если у вас шейп полигональный)

про создать фишнет не знала

ну я ж говорю, неизящно. Я когда вижу стенку, разбегаюсь и в нее со всей дури.

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

В GA в ArcToolbox есть функция GA Layer To Points.
Fishnet создаётся с лэйблами (центроиды ячеек сетки). А потом этой функцией присваиваете значения в точках центроидов.
Если нужно считать среднее значение в ячейках сетки, то нужно сначала перевести геостатистический слой в растр (.tif, например, только задавайте значение пикселя действительно малое), а потом используйте инструмент Zonal Statistics. Только теперь вы получите значения не в центроидах, а в самих ячейках сетки.

Геостатистика вообще то грид не создаёт.
В получаемом слое есть пункт Data-Export Data.
И загоняйте в растр с нужным размером ячейки.

jerry-maori писал(а):Геостатистика вообще то грид не создаёт.

"геостатистика" (если не брать блочную) всегда выполняет оценку в точках, обычно на регулярной сетке. Т.е. таки создает грид. Если нужны значения на густой сетке, то нужно ее соответственно и задавать, поскольку (в отличие от сплайнов и полиномов) непрерывной функции здесь нет, оценка выполняется индивидуально в заданных точках, и интерполировать ее нежелательно (зачем тогда вообще кригинг использовать).

P.S. чем гонять данные туда-сюда, проще обойтись и без Арки, и без Статистики, а сделать все в R (взяв модель из Арки, если таковая вообще строилась), в R с кригингом все хорошо, как и со статистикой. Кроме того, в R есть симуляция (случайные поля), которых в Арке вроде как нет, что позволяет оценить не только средние, но и дисперсии, а это бывает важно для оценки рисков и т.д.

Всем спасибо! В итоге сделал следущим образом: Создал Fishnet с точками узлов, вырезал его по нужной мне рамке, затем использовал инструмент GA Layer To Points из тулбокса (он "берет" обычный геостатистический слой и слой с точками узлов). Ок.

to gamm: пока просто нет времени разбираться с R, но в дальнейшем планирую всю эту работу автоматизировать на его основе.