Zeichnen einer Heatmap in R nur auf Basis von Postleitzahlen

Das Paket spatstat ist dein Freund!

Ссылка

Wenn es sich bei Ihren Daten um eine Menge von Ereignissen handelt (wie "ein Verbrechen, das bei (x, y)" für viele (x, y) aufgetreten ist), können Sie eine Kerndichte-Schätzung verwenden, um Ihre Heatmap zu generieren. Sie können hier ein Beispiel sehen:

Ссылка

insbesondere Zeile 72.

Wenn Ihre z-Werte echte Werte sind (was wahrscheinlich der Fall ist, wenn Sie Ihre Frage lesen), können Sie die smooth -Funktion von spatstat verwenden, die einen Gaußschen Kernel verwendet, um die Interpolation durchzuführen und eine Reihe von Pixeln zurück zu bekommen interpolieren Sie Ihre Daten.

Alternativ können Sie das akima -Paket (wie von den spatstat authors empfohlen) für die Interpolation an von Ihnen angegebenen Stellen verwenden. Es verwendet entweder lineare oder spline-basierte Interpolation und scheint ziemlich einfach zu sein (obwohl ich keine Erfahrung damit habe!).

Ganz allgemein wird das, was Sie versuchen zu tun, oft "Kriging" genannt, und Sie erhalten also viele Ergebnisse, wenn Sie bei Google nachschlagen. Siehe z.B. die Pakete gstat und geoR .

Schließlich (und FTW) können Sie Gaußsche Prozesse verwenden, um die gleiche Sache zu tun. Dies gibt Ihnen eine Verteilung über mögliche Interpolationen, die Ihre Daten ergeben. Das R-Paket kernlab hat eine Implementierung, obwohl ich keine Ahnung habe, wie man es benutzt.

Jeffrey Breen hat diese Art von Handlung mit den Paketen ggplot2 und zipcode gemacht: Ссылка

Damit solltest du anfangen

Sie können sich den Code des Gewinners der ggplot-Fallstudien ansehen: hier

Der Autor hat mit geom_tile eine Heatmap erstellt, der Titel des Diagramms lautet: Violent Crime Wetterkarte von Downtown Houston