Геоинформационные системы (ГИС) и Дистанционное зондирование Земли
Цель интерполяции получить/заполнить значениями растр так чтобы значения были свяаны с основной линией (растр 1) и повторяли геометрию линий. На пересечении должно быть либо большее из значений либо сумма. Поскольку дальнейшие рассчеты предполагают использование значений растра из линий, не только в координатах линий, но и по нормали к ним, где в исходном растре NULL, то есть использовать значения из полученной интерполяцией/буферизацией области вне линий. Линии имеют разные значения растра. Это что-то похоже к обратной задаче выделения водотоков по ЦМР. Есть линии, надо получить типа водосборные бассейны, но при этом нет ЦМР. Идеально было бы, если это была бы не интерполяция а просто что-то типа буферизация с заполнением значениями из линии. Но, без потери/обьединения соседних мелких кусков других линейных обьектов.KolesovDmitry писал(а):А какое требование к результату интерполяции? Чем второй рисунок лучше третьего? Тем, что линии прямые или тем, что на перекрестках значения растра просуммированы или еще что?
Другими словами, какова цель интерполяции?
А мне кажется, что именно на них и похоже. В создании растра №2 могли поучаствовать модули r.grow.distance и r.buffer.SBoris писал(а):На буферные зоны - не похоже.
Посмотрел. Я импортировал tiff.amuriy писал(а):SBoris, если вы делали это в GRASS, не осталось ли у вас области/набора? Тогда можно было бы просто посмотреть использованные команды...
r.grow позволяет выростить поли область и присваивать вырощенной зоне только вручную заданное одно значение. НапримерKolesovDmitry писал(а):Судя по описанию задачи, вам должен подойти модуль r.grow (без distance). Нужно только предварительно убедиться, что в растре 1 в ячейках на перекрестках хранится максимальное значение из пересекающихся категорий.
Да и результаты, приведенные на растре 2, похожи на то, что выходит из r.grow. Хотя есть и отличия: похоже, что при его построении ячейки перекрестков были просуммированы.
Код: Выделить всё
new=6550Нет-нет, сейчас специально проверил, r.grow выращивает области и присваивает ячейкам те категории, из которых эти ячейки выросли.SBoris писал(а): r.grow позволяет выростить поли область и присваивать вырощенной зоне только вручную заданное одно значение. Например, а мне надо количество регионов = количеству линий и задать каждому региону значение, которое соответствует значению на линии. Вот в чём проблемма.Код: Выделить всё
new=6550
если я правильно понял задачу, то мы делали (в похожих обстоятельствах) примерно так:SBoris писал(а):Пол года назад сделал из растра или вектора (непомню, но для наглядности показываю растр)
Да, работает в терминале. При больших радиусах получается очень похожее на растр3.KolesovDmitry писал(а):Нет-нет, сейчас специально проверил, r.grow выращивает области и присваивает ячейкам те категории, из которых эти ячейки выросли.SBoris писал(а): r.grow позволяет выростить поли область и присваивать вырощенной зоне только вручную заданное одно значение. Например, а мне надо количество регионов = количеству линий и задать каждому региону значение, которое соответствует значению на линии. Вот в чём проблемма.Код: Выделить всё
new=6550
Действительно, если указать параметр new, то будут созданы ячейки с указанными категориями, но так вы просто не указывайте этот параметр
Считали по ОНД-86?gamm писал(а): Будет очень похожий результат; считали антропогенную нагрузку от транспорта.
не, экологи дали эмпирику. Типа на какое расстояние от дороги может забрести сборщик краснокнижных растенийSBoris писал(а):Считали по ОНД-86?
Спасибо в копилку знаний не положишь.SBoris писал(а):Спасибо всем!
Лучше расскажите, чем дело-то кончилось? Помог <r.grow> или нашёлся другой способ?r.grow - работает, но не совсем так как хочется (просто иначе и не может быть с точки зрения алгоритма буферных зон). Либо делать цикл как советовалamuriy писал(а):Спасибо в копилку знаний не положишь.SBoris писал(а):Спасибо всем!
а потом все растры сложить (пока что не хочу питонить - интуитивно лень почему-то, хотя для меня это не проблема - могу R, SciPy, Octave, SciLab). Поскольку буферная зона заокруглена, то в местах пересечения линий получаются треугольники, которые как бы налазят на линию и пробивают её тонкой линией насквозь. В этом случае они накрывают большие значения растра меньшими, что несовсем корректно даже на растрах с большими разрешениями. Скорее всего я когда-то в ArcMap нашел функцию которая делает растр по алгоритмуgamm писал(а):1) создать бланковку (растр с нулями)
2) цикл по всем линиям
2.1) растеризовать линию
2.2) сделать дистанционное преобразование...
(растр 2 в начале поста), но, к сожалению не имею возможности сейчас проверить.gamm писал(а):
Код: Выделить всё
r.grow --overwrite input=dp_car@ruska output=dp_car_growth radius=100Сейчас этот форум просматривают: нет зарегистрированных пользователей и 0 гостей
© GIS-Lab и авторы, 2002-2017. При использовании материалов сайта, ссылка на GIS-Lab и авторов обязательна. Содержание материалов - ответственность авторов (подробнее).