2016-04-01 10 views
0

Ich suche seit einigen Tagen nach einer Lösung, wie ich große georeferenzierte Raster in ArcMap in kleinere zerlegen kann. Natürlich gibt es das Werkzeug "Teilungsraster", das Raster nach Features vom Polygon trennt, aber die Rasternamen werden nach dem Zufallsprinzip zugewiesen. Stattdessen muss ich die Namen der neuen Raster definieren - ich teile das Raster mit einem Shapefile aus dem Koordinatenraster und die Namen der Raster sollten mit den Rasterseitenzahlen identisch sein. Ich habe überhaupt nicht mit Python gearbeitet, also gibt es eine Chance, dass es dort gemacht werden kann?Geteiltes Raster, Namen neuer Raster definieren

Antwort

0

Was passiert, wenn Sie Folgendes verwenden: Datenverwaltungstools -> Raster -> Rasterverarbeitung -> Clip? Enthält das Ergebnis keine Feature-ID im Namen des erstellten Rasters oder etwas ähnliches? Ich denke, das kann nützlich sein für Hinweise: https://gis.stackexchange.com/questions/84969/clip-raster-using-shapefile-into-multiple-rasters-using-arcgis-desktop Kann mich nicht genau an die Methode erinnern, die Feature-ID im Dateinamen bereitgestellt, wenn Raster mit Fischnetz in shp-Datei abgeschnitten wurde.

+0

Dank die Methode in dem Link, den Sie einen Modellbauern geteilt beschrieben Sie- und gelöst hat die Verwendung von „clip-Raster“ Tool scheint mein Problem ! – User6819

+0

Gut zu wissen. Obwohl das Spaßige daran ist, dass die meisten Leute das Erzeugen von Fischernetz schwieriger finden, als das Clipping selbst. –

0

Hier ist ein Beispiel dafür, wie Sie mit tun können, dass 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)