Ich habe zwei Raster: ein digitales Höhenmodell (DEM) und ein Landsat-Band. Die Rasterauflösung von beiden ist 30 m, aber die Raster fallen nicht zusammen (siehe Bild: Das dunkle Raster oben ist das DEM, das helle Raster darunter ist das Landsat-Band).
Ich möchte eine Transformation des Landsat-Bandes durchführen, damit sein Gitter perfekt mit dem DEM-Gitter übereinstimmt. Ich habe versucht, das Landsat-Band in ArcGIS unter Verwendung des DEM als Fangraster in den Umgebungseinstellungen auf 30 m neu abzutasten, dies hat jedoch keine Auswirkungen. Ich habe auch versucht, das Landsat-Band auf das Maß eines Vektorpolygons zu beschneiden (auf das das DEM ebenfalls beschnitten ist), während das DEM als Fangraster festgelegt wurde, aber ohne Erfolg.
Ich vermute, dass dies passiert, weil beide Raster bereits die gleiche Auflösung haben (nur auf unterschiedlichen Rastern), aber ich muss ihre Raster wirklich ausrichten, um weitere Analysen durchführen zu können.
Kann ich dies in ArcGIS Desktop tun?
quelle
Antworten:
Nach weiteren Untersuchungen stellte ich fest, dass die Fehlfunktion Fang - Raster Einstellung tatsächlich durch einen Fehler in ArcGIS 10. Die Lösung Hintergrund Geoprocessing (auszuschalten zu sein scheint verursacht wurde Geoprocessing \ Geoverarbeitungsoptionen ... dann deaktivieren Aktivieren der Hintergrundverarbeitung Abschnitt), und führen Sie den Beschneidungsvorgang mit angegebenem Fangraster in ArcCatalog aus.
In den ArcGIS-Foren wird hier diskutiert: http://forums.arcgis.com/threads/50808-Another-Snap-to-Raster-problem
Anscheinend ist der Fehler in ArcGIS SP 4 behoben. Ich habe das Service Pack installiert, aber noch nicht getestet.
Interessanterweise (oder eher ärgerlich) funktioniert der Snap-Raster, wenn ich den Clip mit dem angegebenen Snap-Raster (Snap-Raster ist DEM) gemäß der Lösung im obigen Forumbeitrag in meinem Interessenbereich-Polygon durchführe, aber mit Landsat Das Band verfügt über eine Spalte und eine Zeile mehr als das DEM-Raster, das auch in ArcGIS auf das interessierende Polygon gekürzt wurde. Dies ist ärgerlich, da für die weitere Analyse die beiden Raster genau die gleiche Anzahl von Zeilen und Spalten haben müssen. Eine einfache Problemumgehung besteht darin, eine Berechnung mit der
Con()
Funktion durchzuführen, mit der ein neues Raster mit der gleichen Anzahl von Zeilen und Spalten wie das DEM erstellt wird, das jedoch die Werte des Landsat-Bandes aufweist.quelle
Tatsächlich hatte ich lange Zeit mit diesem Problem zu kämpfen, und ich fand, was ich für eine sehr effektive und ziemlich einfache Abhilfe halte. Wählen Sie eines Ihrer Raster aus, von dem Sie glauben, dass es sich an der "richtigen" Stelle befindet, zoomen Sie in ArcGIS hinein und erstellen Sie ein neues Punkt-Shapefile direkt über dem Schnittpunkt von vier Rasterzellen. Der einzige Trick dabei ist, sicherzustellen, dass Sie so weit wie möglich hineingezoomt sind und dass Sie den neuen Punkt direkt über dem "Kreuz" aus vier Rasterzellen erstellen. Notieren Sie dann die Koordinaten des Punktes, den Sie gerade erstellt haben. Projizieren Sie abschließend Ihr zweites Raster neu (lassen Sie jedoch die Projektions- und Ausgabezellengröße unverändert) und verwenden Sie die Koordinaten Ihres neu erstellten Punkts als "Referenzpunkt". Dadurch werden Ihre beiden Raster perfekt aufeinander abgestimmt. Keine Fehler,
quelle
Sie können versuchen, mit dem Maskierungswerkzeug das DEM mit Landsat zu überlagern.
Möglicherweise müssen Sie ein Untersuchungsfeld erstellen, um Pixel zu erfassen, die vollständig übereinander liegen.
quelle
Messen Sie den Abstand zwischen einem Raster und dem anderen und verschieben Sie dann ein Raster entsprechend.
Sie können das Verschiebungstool aus der Datenverwaltungs-Toolbox verwenden.
quelle
Ich hatte das gleiche problem Ich habe die Lösung von Louis Longchamps ausprobiert , die bei einer kleinen Testprobe wunderbar funktioniert hat, aber mein Datensatz war einfach zu groß, um diese Methode praktisch anzuwenden. Nach einer Weile und vielen Möglichkeiten, dies zu tun, habe ich schließlich aus einer Laune heraus das "Extrahieren durch Rechteck" nach einem Resample verwendet, das genau so funktioniert hat, wie die Raster-Einrastfunktion im Clip-Tool funktionieren sollte. Ich habe einen Modellbauer eingerichtet und versucht, ihn zu verwenden, aber es gab einen Fehler, so dass ich diesen Schritt nach dem anderen ausführen musste, aber schließlich stimmten alle meine Raster überein.
quelle
Ich war mit der Idee, das Raster zu verschieben, nicht zufrieden, weil Sie meiner Meinung nach die Informationen ändern. Um dies zu vermeiden, habe ich eine Folge von Konvertierungen verwendet.
Hier sind die Schritte:
Properties/Source
dasCellsize
auf einem Blatt Papier mit allen verfügbaren Dezimalstellen.Conversion Tools/Raster to Point
.Spatial Analyst Tools/Extraction/Extract Multi Values to Points
, die Sie verwenden möchten, und geben Sie das soeben erstellte Punkt-Feature (in Schritt 2) als Eingabepunkt-Feature und alle Raster an, die Sie in den Eingabe-Rastern aufnehmen möchten.Conversion Tools/To Raster/Point to Raster
, indem Sie angeben, welche Ebene Sie im Feld "Wert" erstellen, und einen entsprechenden Dateinamen UND denCellsize
Wert verwenden, den Sie auf dem Blatt Papier notiert haben.Sie sollten perfekt ausgerichtet sein, und Sie können diese beispielsweise auf R übertragen und Matrix und Vektoren mit exakt den gleichen Abmessungen haben.
quelle
Wir waren so besorgt über unser DEM-Image, dass es auch nicht auf das überwachte klassifizierte Rasterbild abgestimmt war. Wir haben auch die Methode von "Matthew Long" verwendet .
Wir haben die Schnappschüsse nicht ausgerichteter und ausgerichteter Bilder hinzugefügt.
Nach der Methode von Matthew Long wurde unser DEM-Bild vollständig mit unserem überwachten klassifizierten Rasterbild abgeglichen.
Dies ist das Image des Tools, das wir zum Ausrichten unseres DEM-Images in ARCGIS 10.2.1 in Data Management-Tools verwendet haben:
Dies ist das Bild vor dem Ausrichten, Schwarzbild ist DEM und Farbbild wird klassifiziertes Rasterbild überwacht.
Dies ist das Bild nach dem Ausrichten, Schwarzbild ist DEM und Farbbild wird klassifiziertes Rasterbild überwacht.
quelle
Ich löse mein Problem mit der Methode von "Matthew Long". Um es jedoch noch einfacher und 100% genauer zu machen. Verwenden Sie ohne Zoomen einfach die linke obere Ausdehnungskoordinate des gut ausgerichteten Rasters.
Gehen Sie dazu wie folgt vor: Klicken Sie mit der rechten Maustaste auf die Eigenschaft im Bereich Quelle unter der Eigenschaft Umfang
Mithilfe der linken oberen Ausdehnungskoordinate können Sie auch ohne manuelle Eingriffe eine Python-Lösung erstellen.
quelle
Ich hatte ein ähnliches Problem mit zwei Dateien aus verschiedenen Quellen. Zunächst eine Rasterdatei für das Referenzniveau der Entwaldung auf Länderebene (Peru). Zweitens wird eine konsolidierte Datei auf der GEE-Plattform heruntergeladen und auf das erste Elementsystem projiziert. Zusammenfassend zeigte die Überlappung von beiden, dass es eine Verschiebung zwischen Zellen (30 × 30 Zellen) von nicht mehr als 5 Metern gab. Um die zweite Datei an der ersten auszurichten und anzupassen, musste ich eine * erweiterte Basis der ersten Datei erstellen, die beide Dateien * abdeckt und die Pixelausrichtung der ersten Datei mit ihrem Projektionssystem und anderen Funktionen beibehält . Zum Zeitpunkt der Neuprojektion meiner zweiten Datei oder eines Nebenprodukts mit einem ähnlichen Offset
quelle