Форумы GIS-Lab.info

Здравствуйте, в процессе создания модели в Modelbuilder столкнулся с проблемой - у меня есть регулярная сетка на 4.5 млн. точек по которой должна осуществляться диффузийная интерполяция. Для точек сетки производится некоторый расчет. И одним из компонентов этого расчета должны стать здания - полигоны, на месте которых (т.е. там, где точка с ними пересекается) расчета быть не должно, соответственно все расчеты в точке, по-моему, стоит домножать на 0, а в местах отсутствия здания наоборот домножать на 1, чтобы расчет оставался неизменным. Вот отсюда вопрос - как мне с точках сетки спросить о наличии под ним здания, чтобы в результирующем поле были значения 1 и 0.

Были мысли - конвертировать это дело в растр, и оттуда extract value to points поднимать значения, и все значения, что больше 0 заменять на 1, но это получается очень громоздко для модели, и нормально перевести это дело в растр у меня не получилось, получилось эдакое мракобесие
Или играться пространственной выборкой, выделять по пересечению, создавать из них новый слой, там задавать нули, оставшимся задавать единицы, а потом merge их, но опять же это мракобесие и автоматизация расчета в моделбилдер слабо совместимо по-моему.

Заранее спасибо.

каким методом делаете интерполяцию?

novia писал(а):каким методом делаете интерполяцию?

Diffusion Interpolation With Barriers (Geostatisical Analyst) - но при 12000 барьеров - зданий он не очень хочет все это дело считать

И если методически - на месте зданий получаются "дырки", но расчетное поле не корректируется, словно расчет произвелся без барьеров, а они потом просто вырезались оттуда, хотя в ресурсном центре ЭСРИ в описании метода показаны расчеты, без барьера и с барьером, где максимумы не выходят за пределы барьера, у меня так не выходит.

сделайте интерполяцию без зданий, а потом как сами написали - сделайте умножение на 1 и 0 (сделайте растр со зданиями)

Решил так: Конвертировал полигональные здания в растр с маленькой ячейкой (1) с полем высоты сооружения, потом извлек значения по сетке, все что больше нуля приравнял к единице. Однако стало ясно что 4.5 млн точек маловато, пришлось вдвое увеличивать плотность сетки...Вот только где теперь процессор соответствующий изыскать. Но в целом спасибо больше за помощь