Dieses Forum scheint ja noch nicht so stark frequentiert zu sein, aber vielleicht kann mir ja trotzdem jemand helfen
Ich habe ein Problem beim speichern eines SpatialPolygonDataframes mit dem Befehl "writeOGR" des Paketes rgdal.
Das merkwürdige an der ganzen Sache ist, dass das Speichern des Objektes noch vor ein paar Wochen funktioniert hat und ich weder das Projekt, noch die Quelldatei aus der der Dataframe stammt geändert habe.
meine Arbeitsgänge bestehen darin, eine MapInfo Datei per "readOGR" einzulesen. Anschließend will ich diese (zunächst um weitere Fehlerquellen auszuschließen), sofort wieder per "writeOGR" als Map Info Datei abspeichern. Hierbei kommt es zur Fehlermeldung "Failed to create feature". Die Dateien werden dann zwar im Verzeichnis abgelegt, sind aber leer (lediglich die Namen der Spalten sind in der Datei zu finden).
Ich habe auch schon versucht, die Daten als ESRI Shape abzuspeichern, dann werden die Dateien wesentlich größer, doch wenn ich diese versuche zu öffnen geht etwas schief, also stimmt mit den Shapes auch irgendetwas nicht.
Damit ihr das ganze nachvollziehen könnt, hier noch mal ein Beispiel:
Code: Alles auswählen
library(rgdal)
setwd("G:/Löschen")
(dat <- data.frame(Number=c(1,2,3,4,5,6),
East=c(371008,371124,371346,371276,371051,371344),
North=c(5720763,5720567,5720761,5720493,5720506,5720490),
BW=c(320,350,300,380,275,390)))
coordinates(dat) <- ~East+North
proj4string(dat) <- "+proj=utm +zone=32 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs"
writeOGR(dat, dsn = "Zones.tab", layer = "Zones", driver = "MapInfo File")