Разделить растровые, определить имена новых растров

Question

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

Answer 1

Что произойдет, если вы используете: Инструменты управления данными -> Растр -> Растровая обработка -> Клип? Не содержит ли идентификатор функции имя созданного растра или что-то в этом роде? Я думаю, это может быть полезно для подсказок: https://gis.stackexchange.com/questions/84969/clip-raster-using-shapefile-into-multiple-rasters-using-arcgis-desktop Не помню точно метод, который предоставил идентификатор функции в имени файла, когда растр был обрезан с помощью fishnet в shp-файле.

Answer 2

Вот пример того, как вы можете сделать это с R

library(raster) 
# example data 
r <- raster() 
values(r) <- 1:ncell(r) 

# create blocks 
a <- aggregate(raster(r), 60) 
b <- as(a, 'SpatialPolygons') 

f <- function(i) { 
    fname <- paste0('raster_', i, '.tif') 
    crop(r, extent(b[i,]), filename=fname) 
} 

x <- sapply(1:length(b), f)