2016-07-07 8 views
-1

Ich versuche, einen räumlichen Polygon-Datenrahmen aus einer Reihe von Breiten- und Längengraden und einem entsprechenden Wert zu erstellen. Jede Zeile muss ein einzelnes Polygon sein. Ich kann auch eine Shape-Datei erstellen und sie dann wieder lesen, sondern muß ich es tun, indem Sie einfachR, Räumliches Polygon

Hier ein räumliches Polygondatenrahmen zu schaffen ist ein Beispiel für meine Daten:

head(Mydata) 

    NA. count... zipcode Latitude Longitude 
1  1  1  0 44.200797 24.50230 
2  3  67 560001 12.976594 77.59927 
3  4  28 560002 12.963521 77.58211 
4  5  90 560003 13.001914 77.57134 
5  6  26 560004 12.943751 77.57376 
6  7  34 560005 12.997989 77.62265 

Ich habe den Grafen schildern entsprechend jedem Ort auf einer interaktiven Karte. Ich habe viele andere Beiträge durchgesehen, konnte aber keine Lösung finden. Ich bin neu in R, also bitte, ertragen Sie mit mir.

+1

Bitte lesen Sie die Informationen über [wie eine gute Frage] (http://stackoverflow.com/help/how-to-ask) und wie eine [reproduzierbar geben Beispiel] (http://stackoverflow.com/questions/5963269). Dies wird es anderen sehr erleichtern, Ihnen zu helfen. – zx8754

+0

Es tut mir leid, aber das ist mein erster Ques on Stack-Überlauf. Ich werde es im Hinterkopf behalten und meine Frage bearbeiten, um es klarer zu machen. –

+0

Was ist zählen? Sie haben für jeden Punkt einen Vount-Wert. Ich dachte, du willst einen Polygondatenrahmen. Sind Sie sicher, dass Sie keinen Punktdatenrahmen möchten? – Alex

Antwort

0

Vielleicht so etwas wie

library(dismo) 
v <- voronoi(Mydata[, c("Longitude", "Latitude")] 
Mydata$id <- 1:nrow(Mydata) 

v <- merge(v, Mydata, by='id') 
spplot(v, 'count') 
+0

Ich muss eine interaktive Choroplethemap erstellen, die eine Basis eines Ortes hat, der durch die Postleitzahl in den Daten definiert ist, wie zum Beispiel: http://openbangalore.org/maps/ , aber Bereiche, die anhand der Farben differenziert wurden. –

+0

Dies ist natürlich nicht möglich, basierend auf Ihren Daten allein (und einer anderen Frage). Aber es ist einfach zu tun: Holen Sie sich eine Datenquelle mit den Zip-Grenzen und dann 'z <- fusionieren (zip, Mydata, by = 'Postleitzahl') – RobertH