Создайте растр высоты сцены, используя путь наименьшей стоимости и r

У меня есть точечный шейп-файл с идентификаторами станций и высотами сцены. Я хотел бы создать растр, в котором каждая ячейка имеет значение высоты сцены (в метрах) ближайшей к этой ячейке станции на месте.

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

Наборы данных:

1) Шейп-файл точки с высотами сцены дельты реки

2) Шейп-файл протяженности дельты реки

3) Растр дельты, где NA обозначают сушу (при необходимости также могут быть нули), а 1 - вода. Два набора данных с разрешением 10 метров и разрешением 30 метров.

Одна концептуальная проблема, с которой я столкнулся, связана с количеством имеющихся у меня небольших потоков. Например (изображено на изображении ниже), станция 1 (обведена синим) технически ближе к черной области x, чем станция 2 (обведена красным), но значение высоты сцены, выделенное красным, более репрезентативно для точки x. Между двумя потоками есть NA, означает ли это, что значение не будет перепрыгивать через потоки?

Концептуальное изображение проблемы

Как я могу переназначить значения в моем растре (все единицы) на высоту сцены ближайшей станции и убедиться, что эти значения не переходят от потока к потоку? Нужно ли мне использовать путь с наименьшими затратами? Как лучше всего это сделать?

Я хотел бы использовать R, но при необходимости могу использовать ArcMap.


r gis raster arcmap spatial-interpolation
person Christiana A    schedule 03.01.2017    source источник
comment
Хорошо, я перефразирую то, что, по-моему, ваша проблема, чтобы увидеть, правильно ли я понял. У вас есть шейп-файл, представляющий данные, которые вы хотите назначить ячейкам растра. Вы хотите классифицировать ячейки некоторого растра на основе их линейной близости к значениям в вашем шейп-файле точек; вы хотите, чтобы линейная близость определялась данными в другом шейп-файле.   -  person SeldomSeenSlim    schedule 04.01.2017
comment
Я хочу, чтобы линейная близость определялась шейп-файлом точки. Извините, я не заметил этого, когда опубликовал его, но набор данных 2 и набор данных 3 почти одинаковы, за исключением того, что один является растром, а другой - вектором.   -  person Christiana A    schedule 04.01.2017
comment
Если вы можете привести воспроизводимый пример, я, вероятно, смог бы сделать шаг за шагом.   -  person SeldomSeenSlim    schedule 04.01.2017
comment
Извините, последние несколько дней я сильно заболел и не был в сети. Как мне лучше всего предоставить вам свой набор растровых данных? Могу ли я просто связать вас с учетной записью github, или это запрещено в stackoverflow?   -  person Christiana A    schedule 07.01.2017
comment
Вы можете начать с ответа на другой вопрос ниже. Похоже, что для вашей проблемы может быть готовое решение Arcmap, даже если вам придется переформатировать данные. Если вы поработали над этим, но все еще не получили искомого решения, попробуйте создать произвольный пример с некоторыми фиктивными данными. Опубликуйте эти данные, и мы сможем решить, как решить вашу проблему.   -  person SeldomSeenSlim    schedule 10.01.2017

Ответы (1)


arrow_upward
0
arrow_downward

Поэтому я не уверен, какие инструменты у вас есть, но я думаю, что этот ответ может быть полезен:

Вычисление атрибута сетевого расстояния между несколькими точками в ArcGIS Desktop?

Здесь спрашивающий пытался рассчитать расстояния на дорогах до некоторых точек, но ваша проблема кажется похожей. Я думаю, что главное, что я хотел бы здесь сказать, это то, что вам следует провести классификацию расстояния до сети, прежде чем беспокоиться о растровом слое. Возможно, вам придется преобразовать многоугольник в линии или какой-либо обходной путь, чтобы преобразовать ваши данные в формат, который работает, но это именно та работа, для которой предназначен инструмент.

После того, как вы переклассифицируете шейп-файл реки на основе их сетевого расстояния до заданной точки, преобразуйте полигоны в растр и используйте его для классификации исходного растра. Вы можете сделать это в R или Arcmap. Возможно, Arcmap будет быстрее.

person SeldomSeenSlim    schedule 04.01.2017