Ich versuche herauszufinden, warum der Slot @ data @ values beim Einlesen eines NDVI-Rasters erst dann die tatsächlichen Werte enthält, wenn ich sie manuell einstelle. Beispielsweise: NDVI <- raster("./filename.tif", crs="+proj=longlat +datum=WGS84") NDVI@data@values ## returns: logical(0)...
Ich muss es in eine UTM umwandeln, um den Puffer funktionsfähig zu machen. wets<-readOGR(dsn=".",layer="shapefile") r.raster <- raster() extent(r.raster) <- extent(wets) res(r.raster) <- 100 wets.r <- rasterize(wet,r.raster) plot(wets.r) wetsbuf<-buffer(wets.r,width=500) Während...
Ich versuche, Rasterdaten in Polygondaten zu konvertieren. Bisher verwendete ich shapefile.py, um Punktdaten für einzelne Lat / Lon im Raster zu erstellen. Aber das wird riesige Daten erzeugen und das Abfragen wird ein großes Problem sein. Es wäre also besser, wenn ich die Rasterdaten in...
Ich versuche, ein Shapefile mit GDALs RasterizeLayer in ein Raster zu brennen. Ich erstelle ein Interessenbereich-Raster aus einem anderen Shapefile bei einer bestimmten Pixelgröße. Dieser AOI dient dann als Basis für alle folgenden Rasterisierungen (gleiche Anzahl von Spalten und Zeilen, gleiche...
Ich habe mehrere SRTM DEM-Rasterkacheln von https://dwtkns.com/srtm30m/ heruntergeladen und aus den Kacheln Hillshade-Versionen erstellt. Bei näherer Betrachtung haben die Hillshade-Kacheln jetzt Grenzlinien, die ich nicht ausblenden oder entfernen kann (wie in der Abbildung unten, in der die...
Wie wichtig ist es, Landsat für Atmosphäre zu korrigieren, Dunst zu erzeugen und DN in Reflexionsvermögen umzuwandeln? Welchen Einfluss hat es auf die Ergebnisse, etwa für viele
Ich habe eine Rasterdatei, in der alle Zellen Werte zwischen -3 und 7 haben, kein Datenwert ist -9999. Wie kann ich die Anzahl der Zellen mit einem bestimmten Wert berechnen, z. B. 6? Kann ich den Taschenrechner
Vor allem; Ich habe versucht, eine ähnliche Frage zu finden, ohne Erfolg. Vielleicht liegt es daran, dass ich ziemlich neu in GIS bin und nicht genau weiß, wonach ich suche. Wenn mich jemand auf ein ähnliches Problem hinweist, würde ich diesen Beitrag gerne entfernen. Ich muss eine...
Ich erstelle einige Gleitkomma-Raster in GDAL mit dem GeoTIFF-Treiber. Wenn ich die resultierenden Bilder in QGIS oder Arc lade, setzt die Standardsymbologie die Min auf -3.40282e38 und die Max auf 3.40282e38, sodass das Raster genau grau aussieht. Gibt es eine Möglichkeit, den tatsächlichen...
Ich habe eine Reihe von Rastern (.vrt) mit täglichen Bodenfeuchtedaten. Ich möchte die Pixel aller Raster summieren, um ein Maß pro Monat zu erhalten. Es wird jedoch nicht jeden Tag die ganze Welt abgedeckt, was zu Knotenwerten an Orten führt, an denen für die anderen Tage Maßnahmen existieren. Was...
Ich schreibe eine Python-Methode, die zwei Raster hinzufügt und eine einzelne Rasterausgabe generiert. Aus Gründen, die außerhalb meiner Kontrolle liegen, sind die Ausmaße der Eingabe-Raster unterschiedlich, sie überschneiden sich jedoch. Ist es möglich, ausschließlich die Eingabe-Rasterpixel zu...
Für eine Studie zum Ausbau des Straßennetzes im Regenwald versuche ich, Straßen aus Landsat-Bildern zu extrahieren. Wir haben bereits scharfe und wolkenfreie Verbundwerkstoffe, auf denen die Straßen mit dem Auge gut sichtbar sind, aber das Extrahieren dieser in Linienmerkmale erweist sich als...
Ich versuche, mit Python einen Histogrammabgleich durchzuführen, um den Mosaikprozess mehrerer überlappender Raster zu verbessern. Ich stütze meinen Code auf den folgenden: http://www.idlcoyote.com/ip_tips/histomatch.html Bisher habe ich es geschafft, den überlappenden Bereich zweier benachbarter...
In meiner täglichen Arbeit werde ich ständig gebeten, Bereiche globaler Raster-Datensätze in geografischer Projektion mit einer Auflösung von 30 Bogensekunden zu berechnen. Diese Datensätze sind normalerweise das Ergebnis einer Kombinationsoperation (ein typisches Beispiel sind Vegetationsklassen...
Ich habe versucht, Windspeed-Dateien ohne Erfolg in Raster zu konvertieren. Es wurde in einem Arcgis-Forum diskutiert, dem ich folgte. Dies beinhaltete das Ersetzen von Semikolons durch Leerzeichen und das Hinzufügen einer Kopfzeile. Ich hatte gehofft, eine der ursprünglichen .asc-Dateien zusammen...
Ich muss eine Analyse mit mehreren Kriterien durchführen, um die Frage zu beantworten: "Welches ist das beste Los, das entwickelt werden kann?". Einige der Kriterien sind: Entfernung zur nächsten Bushaltestelle (Punktschicht mit Bushaltestellen) Entfernung zum nächsten Geschäft (Punktschicht mit...
Ich habe mit diesem Befehl ein Raster-TIF-Bild in PostgreSQL importiert: raster2pgsql -s 32643 -I -M filepath.tif -F -t 100x100 public.databassename > filepath.sql Und importierte die Ausgabe-SQL-Datei in die PostgreSQL-Datenbank, in der Folgendes ausgeführt wird: psql -U postgres -d...
Ich habe ein Mosaik eines 1 m LiDar-abgeleiteten DEM. Ich muss eine Teilmenge der Daten als 10m DEM ausgeben. Ich verwende derzeit das Aggregat-Tool in ARCGIS 10, um einen Mittelwert für jedes neue 10-m-Pixel zu erstellen. Irgendwelche Ratschläge, ob dies die beste Technik für eine solche Aufgabe...
Dies sind 4 Bildabschnitte (NE, SE, NW, SW). Sie haben einige Überlappungen zwischen ihnen. Arbeiten mit 1 m 3-Band-Orthobildern. Verschiedene Teile werden manchmal zu unterschiedlichen Zeiten / Quellen geflogen, daher der Unterschied in den Jahreszeiten. Die Bilder sind im TIF-Format. Die...
... anstatt die Raster als lokal gespeicherte GEOTIFFs zu behalten? Ich habe ein großes Projekt, bei dem ungefähr 50.000 GEOTIFFs mit 10 MB als WMS für eine OpenLayers-basierte Webkarte verfügbar sein müssen. Es wird immer nur 1 WMS angezeigt, sodass die Leistung kein Problem darstellt....