Räumliches Korrelogramm mit dem Raster-Paket
Liebe Menschenmenge
Problem
Ich habe versucht, ein räumliches Korrelogramm mit den Paketen nfc, pgirmess, SpatialPack und spdep zu berechnen. Ich war jedoch beunruhigt, den Anfang und den Endpunkt der Entfernung zu definieren. Ich interessiere mich nur für die räumliche Autokorrelation in kleineren Entfernungen, aber dort für kleinere Bins. Da das Raster ziemlich groß ist (1,8 Megapixel), stoße ich bei diesen Paketen auf Speicherprobleme, aber auf das SpatialPack.
Also habe ich versucht, meinen eigenen Code mit der Funktion Moran aus dem Paket-Raster zu erzeugen. Aber ich muss einen Fehler haben, da sich das Ergebnis für den gesamten Datensatz etwas von dem der anderen Pakete unterscheidet. Wenn es keinen Fehler in meinem Code gibt, könnte es zumindest anderen mit ähnlichen Problemen helfen.
Frage
Ich bin mir nicht sicher, ob meine Fokusmatrix fehlerhaft ist. Können Sie mir bitte sagen, ob das zentrale Pixel eingefügt werden muss? Unter Verwendung der Testdaten kann ich die Unterschiede zwischen den Methoden nicht anzeigen, aber auf meinem vollständigen Datensatz sind Unterschiede sichtbar, wie in dem Bild unten gezeigt. Die Bins sind jedoch nicht genau gleich (50m vs. 69m), das könnte Teile der Unterschiede erklären. Im ersten Fach scheint mir diese Erklärung jedoch nicht plausibel zu sein. Oder könnten die unregelmäßige Form meines Rasters und die verschiedenen Arten, mit NAs umzugehen, den Unterschied verursachen?
Vergleich der Own-Methode mit der von SpatialPack
Ausführbares Beispiel
Testdaten
Der Code zur Berechnung der Testdaten stammt aus Ссылка
%Vor%Eigener Code
%Vor%Andere Methoden zum Berechnen des räumlichen Korrelogramms
%Vor%1 Antwort
Um die Ergebnisse des Korrelogramms zu vergleichen, sollten in Ihrem Fall zwei Dinge berücksichtigt werden. (i) Ihr Code funktioniert nur für Bins proportional zur Auflösung Ihres Rasters. In diesem Fall könnte ein kleiner Unterschied in den Bins dazu führen, dass eine wichtige Menge an Paaren eingeschlossen oder ausgeschlossen wird. (ii) Die unregelmäßige Form des Rasters hat einen starken Einfluss auf die Paare, von denen angenommen wird, dass sie die Korrelation für bestimmte Abstandsintervalle berechnen. Ihr Code sollte also beides behandeln, einen beliebigen Wert für die Länge des Bin zulassen und die unregelmäßige Form des Rasters berücksichtigen. Eine kleine Modifikation Ihres Codes zur Lösung dieser Probleme finden Sie weiter unten.
%Vor%Das Bild zeigt, dass die Ergebnisse der Verwendung des SpatialPack-Pakets und des eigenen Codes identisch sind.
Tags und Links r spatial r-raster spatialpack