Я работаю над ЦМР в районе, где мало рельефа. Я хотел бы «сжечь» речную сеть в ЦМР, чтобы я мог точно рассчитать накопление и длину потока. Я использую ArcGIS Desktop...
Я работаю над ЦМР в районе, где мало рельефа. Я хотел бы «сжечь» речную сеть в ЦМР, чтобы я мог точно рассчитать накопление и длину потока. Я использую ArcGIS Desktop...
Можно ли разделить шейп-файл на элемент в Python? (Лучше всего было бы решение, где я могу временно сохранить результирующие векторные объекты в памяти, а не на диск). Причина: я хочу использовать функцию gdal rasterizeLayer с несколькими различными подмножествами шейп-файлов. Для этой функции...
У меня есть линейный шейп-файл, представляющий дорожную сеть. Я хочу растеризовать эти данные, в результате чего полученные значения в растре показывают общую длину линий, которые попадают в растровую ячейку. Данные представлены в проекции Британской национальной сети, поэтому единицы измерения...
Я пытаюсь записать шейп-файл в растр, используя RasterizeLayer из GDAL. Я предварительно создаю растр области интереса из другого шейп-файла, учитывая определенный размер пикселя. Этот AOI затем служит базой для всех последующих растеризаций (одинаковое количество столбцов и строк, одинаковая...
Я хочу создать растр с разрешением 25 х 25 м, в котором каждая ячейка содержит расстояние до ближайшей береговой линии, рассчитанное по центру ячейки. Для этого у меня есть только шейп-файл побережья Новой Зеландии . Я пытался следовать учебному пособию Доминика Роя для того, чтобы сделать это в R,...
В моей рабочей области 100 пространственных точек данных. Я загрузил их таким образом: filenames <- list.files(path="", pattern="XYhectareTravelTimes_ez+.*shp") for(i in filenames){ filepath <- file.path("/",i) assign(i, readShapePoints(filepath)) Они называются так:...
Я конвертирую вектор в растр в R. Однако процесс был слишком долгим. Есть ли возможность поместить скрипт в многопоточную или графическую обработку, чтобы сделать это быстрее? Мой скрипт в растеризованный вектор. r.raster = raster() extent(r.raster) = extent(setor) #definindo o extent do raster...