Ich verwende die GDAL ReadAsArray-Methode, um mit Rasterdaten unter Verwendung von numpy (insbesondere Neuklassifizierung) zu arbeiten. Da meine Raster groß sind, verarbeite ich die Arrays in Blöcken, durchlaufe jeden Block und verarbeite sie mit einer ähnlichen Methode wie im GeoExamples-...
Ich habe zwei Ebenen in ArcGIS: Eine ist ein Vektorpunkt und die andere ist eine Rasterebene. Ich möchte den Wert des Rasterdatengitters an den Vektorpunkt anhängen. Das Problem ist, dass der Vektorpunkt von allen vorhandenen Rastergittern mit Wert entfernt ist, sodass der Befehl "Beispiel" keinen...
Ich bin an einer Möglichkeit interessiert, Raster (z. B. ESRI GRID-Format) oder Shapefiles zu identifizieren, die beschädigt sind und deren Verwendung in einem Analysetool verhindern würde. Ich habe eine große Anzahl in einem strukturierten Ordnerbaum und möchte in der Lage sein, Skripte für den...
Ich möchte alle mehreren Mosaik-Raster (jedes Mosaik besteht aus 7 Ebenen) aus dem angegebenen Ordner in mein R importieren . Dann als einzelne Multilayer-Raster darauf zugreifen. So tat ich: # read all mosaics named "mos....img" in R raster_data<-list.files(path=getwd(), pattern="mos.*.img$") #...
Ich habe ein DEM einer Bergregion. Ich möchte die Bereiche des DEM extrahieren, über denen 60% der gesamten DEM-Fläche liegen. Im Wesentlichen möchte ich also die Bereiche extrahieren, die die oberen Teile der Berge darstellen. Ich habe davon gehört, dass dies als hypsometrische Linie (H60-Linie)...
Wie kann ich unter Bezugnahme auf meine vorherige Frage ein generiertes Rasterbild mit einem generierten Rasterbild in ein Format "geocodieren" / konvertieren, das als Layer auf einer Karte angezeigt werden kann? Kurz gesagt, ich möchte die Ergebnisse der kommerziellen SpatialKey- Software (...
Ich versuche, ein Raster mit den Python-GDAL-Bibliotheken zu erstellen, und bin an dem Punkt angelangt, an dem Daten ausgegeben werden, aber die Ausgabedaten werden auf der x-Achse des Ursprungspunkts gespiegelt . Ich weiß, dass ich etwas übersehen muss, aber ich kann nicht herausfinden, wo ich...
Ich habe zwei Rasterbilder, Boolesch klassifiziert (1 oder 2). Sie überlappen sich in bestimmten Bereichen. Ich brauche ein Raster, das nur aus überlappenden Abschnitten besteht. Verwenden von Arcinfo 10. Answer kann nur Rastermethoden verwenden, keine Konvertierung in Polygone, Punkte...
Ich möchte von einem binären Raster aus Wald / Nicht-Wald mit einer Auflösung von 30 m zu einem Raster mit einer Auflösung von 240 m wechseln, wobei der Wert jeder Zelle der Anteil dieser bewaldeten Zelle ist - dh die Datenwerte für die 240 m-Zellen würden von 0 gehen (wo alle Zellen im 30-m-Raster...
Mein Hintergrund: erfahrener Programmierer, Amateurhistoriker, absoluter Anfänger in GIS (also benutze ich QGIS). Ich habe jedoch gelernt, das QGIS-Georeferencer-Tool zu verwenden. Die Eingabe: Angenommen, ich habe ein Bild wie dieses Wikipedia: Datei: VoyagesOfRabbanBarSauma.jpg . Angenommen, es...
Ich habe ein GeoTIFF-Raster in ArcGIS Desktop, das über die folgenden zwei Links heruntergeladen werden kann: http://www.rtwilson.com/downloads/MODIS_MOD04_16-06-2006_Subset.tfw http://www.rtwilson.com/downloads/MODIS_MOD04_16-06-2006_Subset.tif Wenn Sie diese in ArcGIS Desktop laden, werden große...
Ich versuche ein Raster zu projizieren. In R gibt es dazu die projectRaster()Funktion (unter einem vollständig reproduzierbaren Beispiel): # example Raster require(raster) r <- raster(xmn=-110, xmx=-90, ymn=40, ymx=60, ncols=40, nrows=40) r <- setValues(r, 1:ncell(r)) projection(r) # project...
Ich möchte ADF-Dateien in R laden. Die Daten stammen von dieser Seite: http://www.fao.org/geonetwork/srv/en/metadata.show?id=14057 Ich habe den folgenden Code ausprobiert, den ich nach einigen Recherchen im Internet gefunden habe. Das Problem ist, dass ich in der Klasse RasterLayer negative Werte...
Ich möchte einem Raster, das in einer Postgis-Datenbank gespeichert ist, Attribute hinzufügen. Ich benutze den folgenden Befehl, um das Raster in die Datenbank einzufügen: raster2pgsql -s 4236 -I -C -M /home/user/mysrc/outlabeled.tif public.test | psql -d template_postgis Das Raster...
Arbeiten mit QGIS 2.4.0 Chugiak Ich habe eine Frage zum Ausschneiden von Rasterbildern. Ich habe ein Luftbild, das ich mit den Ausmaßen eines Shapefiles ausschneiden möchte. So weit so gut, ich habe das Raster / Extraktions / Clipper-Tool verwendet, um dies zu tun. Das Problem ist, dass die...
Ich versuche, eine große ECW zu schneiden (Details unten), aber die Rasterdatei ist zu groß, um vollständig verarbeitet zu werden. Einige Details der ECW unten Treiber: ECW / ERDAS Compressed Wavelets (SDK 5.0) Dateigröße: 50 GB Größe ist 450000, 565081 Pixelgröße: 0,15 0,15 COLORSPACE = RGB...
Ich habe ein Rasterbild von Höhen. Ich habe ein Vektor-Shapefile hinzugefügt, das lokale politische Spaltungen enthält. Ich möchte die Höhenfarben beibehalten, aber ich möchte auch die schwarzen Linien beibehalten, die die politischen Spaltungen darstellen. Ich habe die Transparenz des Shapefiles...
Ich versuche, einen derzeit äußerst umständlichen Vektor / Python-Prozess für ein Naturgefahrenmodell zu verbessern. Im Moment haben wir ein langes Skript, das Entfernungs- / Peilungslinien von einem bestimmten Punkt generiert, um Folgendes zu bestimmen: die Art des Polygons, das es schneidet (z....
Ich habe den aktuellen QGIS-Master wie gewohnt kompiliert und installiert und einige Rasterdaten mit der Erweiterung postgis 2.0.2 in meine postgres 9.3-Datenbank importiert. Immer wenn ich versuche, die Rasterdaten mit dem Datenbankmanager in qgis zu importieren, wird folgende Meldung angezeigt:...
Ich habe zwei überlappende .ecw-Schichten desselben Bereichs, obwohl eine in der Mitte mit einem weißen Rand abgeschnitten ist. Wenn ich den Hintergrundwert auf 255.255.255 setze und ihn als hohl definiere, wird nur dieser bestimmte RGB-Wert ausgehöhlt. In Wirklichkeit gibt es in meiner Ebene...