Форумы GIS-Lab.info

Господа.
У меня есть полигональный и точечный слой.
В каждом полигоне некоторое количество точек.
В точечном слое у меня есть данные по кислотности.
Я делаю интерполяцию.
Делаю вырезание по форме полигогального слоя.
Всё прекрасно.
Но данные точек находящиеся в отдельных полигонах влияют друг на друга.
А мне необходима интерполяция внутри полигонов по отдельности,
чтобы данные точек находящихся в других полигонах не влияли.

Возможно ли такое?

Логично было бы сначала разрезать точечный слой на отдельные по полигонам, а потом выполнять интерполяцию.

В принципе идея понятна.
Мне подобное приходило в голову.
Достаточно трудоёмко получится.
У меня порядка сотни полей на каждое хозяйство.
Делаю сотню точечных шейпов через выборку по местоположению.
Сотню квадратных интерполяций по каждому полю,
затем делаю вырезание по форме полигонального слоя?
Может получиться.

- конвертируем полигоны в линии
- интерполируем методом ОВР (IDW), указываем линейные границы полигонов в качестве "Входных линейных объектов барьеров"(in barrier polyline features)

_lenik писал(а): - интерполируем методом ОВР (IDW), указываем линейные границы полигонов в качестве "Входных линейных объектов барьеров"(in barrier polyline features)

Что-то у меня не получилось.
Полигоны сделал линейными объектами.
Интерполяцию я делаю через Geostatistical Analyst,
а там нет инструмента или такого поля как "Входные линейные объекты барьеры"(in barrier polyline features).
Я работаю в 9.3.

Воронин С.В. писал(а):Я работаю в 9.3.

это зря ...
не занимайтесь мазохизмом. Есть простой путь под названием брутфорс:
1) берем что-то типа R или Python (тут дело вкуса)
2) разбиваем исходную выборку на подмножества по полигонам.
3) строим или грузим сетку (от нее нам нужны только координаты узлов), разбиваем сетку на подмножества по полигонам
4) запускаем цикл, снаружи по переменным ("кислотность"), внутри по группам (полигонам), интерполируем группы из (2) на группы из (3), если точек много используем MBA, если мало - tension splines. Если хотим получить ерунду, используем IDW (если не хотим - не используем). По завершении внутреннего цикла все у нас готово для данной переменной, сохраняем результат.
5) если таки хочется мазохизма, делаем 2-4 в Арке ...

Воронин С.В. писал(а):Что-то у меня не получилось.
Полигоны сделал линейными объектами.
Интерполяцию я делаю через Geostatistical Analyst,
а там нет инструмента или такого поля как "Входные линейные объекты барьеры"(in barrier polyline features).
Я работаю в 9.3.

Это я имею ввиду Spatial Analist -> Interpolation -> IDW

gamm писал(а): если таки хочется мазохизма, делаем 2-4 в Арке ...

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

_lenik писал(а):Это я имею ввиду Spatial Analist -> Interpolation -> IDW

Что-то на 9.3 Spatial Analist не срабатывает, пишет ошибку какую-то.
Поставил 10.2, получилось построить, но пока не могу поймать необходимого качества.

Пойду-ка я спать, а то уже и голова не работает, да и руки устали.

_lenik писал(а):Это я имею ввиду Spatial Analist -> Interpolation -> IDW

Получается.
Но скорость подает значительно.
Очень значительно. На порядок, а то и больше.
Приходится даже принудительно задачу снимать.

Воронин С.В. писал(а):Но скорость подает значительно.
Очень значительно. На порядок, а то и больше.

это свойство метода. Если точек мало, терпите (но смотрите, что получается). Если много, то переходите к (2)-(4)