My guess is that it's the logical column. Try
gpsori1$UTM.Ch <- as.character(gpsori1$UTM.Ch) writeOGR(gpsori1, dsn="MATA20090729", layer="MATA20090729gps1",driver="ESRI Shapefile") On Thu, Nov 5, 2009 at 9:35 AM, Agustin Lobo <[email protected]> wrote: > Thanks, > No POSIXct fields: >> summary(gpsori1) > Object of class SpatialPointsDataFrame > Coordinates: > min max > coords.x1 453020 457145 > coords.x2 4623500 4625866 > Is projected: NA > proj4string : [NA] > Number of points: 2096 > Data attributes: > No UTM.Zone UTM.Ch UTM.East > Min. : 1.0 Min. :31 Mode:logical Min. :453020 > 1st Qu.: 524.8 1st Qu.:31 TRUE:2096 1st Qu.:453601 > Median :1048.5 Median :31 NA's:0 Median :454140 > Mean :1048.5 Mean :31 Mean :454418 > 3rd Qu.:1572.2 3rd Qu.:31 3rd Qu.:455120 > Max. :2096.0 Max. :31 Max. :457145 > > UTM.North Altitude Date Time > Min. :4623500 Min. :1648 2009/07/29:2095 08:49:04: 3 > 1st Qu.:4624886 1st Qu.:1795 2009/07/30: 1 08:49:05: 2 > Median :4625170 Median :1821 08:49:06: 2 > Mean :4625068 Mean :1835 08:49:07: 2 > 3rd Qu.:4625338 3rd Qu.:1885 08:49:08: 2 > Max. :4625866 Max. :1952 08:49:09: 2 > (Other) :2083 >> str(gpsori1) > Formal class 'SpatialPointsDataFrame' [package "sp"] with 5 slots > ..@ data :'data.frame': 2096 obs. of 8 variables: > .. ..$ No : int [1:2096] 1 2 3 4 5 6 7 8 9 10 ... > .. ..$ UTM.Zone : int [1:2096] 31 31 31 31 31 31 31 31 31 31 ... > .. ..$ UTM.Ch : logi [1:2096] TRUE TRUE TRUE TRUE TRUE TRUE ... > .. ..$ UTM.East : int [1:2096] 453245 453245 453245 453245 453245 453245 > 453245 453245 453245 453245 ... > .. ..$ UTM.North: int [1:2096] 4625338 4625338 4625338 4625338 4625338 > 4625338 4625338 4625338 4625338 4625338 ... > .. ..$ Altitude : num [1:2096] 1648 1648 1648 1648 1648 ... > .. ..$ Date : Factor w/ 2 levels "2009/07/29","2009/07/30": 1 1 1 1 1 1 > 1 1 1 1 ... > .. ..$ Time : Factor w/ 1452 levels "08:23:05","08:23:07",..: 1 2 3 4 5 > 6 7 8 9 10 ... > ..@ coords.nrs : num(0) > ..@ coords : num [1:2096, 1:2] 453245 453245 453245 453245 453245 ... > .. ..- attr(*, "dimnames")=List of 2 > .. .. ..$ : NULL > .. .. ..$ : chr [1:2] "coords.x1" "coords.x2" > ..@ bbox : num [1:2, 1:2] 453020 4623500 457145 4625866 > .. ..- attr(*, "dimnames")=List of 2 > .. .. ..$ : chr [1:2] "coords.x1" "coords.x2" > .. .. ..$ : chr [1:2] "min" "max" > ..@ proj4string:Formal class 'CRS' [package "sp"] with 1 slots > .. .. ..@ projargs: chr NA >> class(gpsori1$Date) > [1] "factor" >> class(gpsori1$Time) > [1] "factor" > > After your message I thought it could be because of the factors, but the > following does not work either: > >> gpsori1 <- >> read.table("MATA20090729/MATA20090729_1.csv",sep=",",header=T,stringsAsFactors >> =F)> coordinates(gpsori1) <- cbind(gpsori1$UTM.East, gpsori1$UTM.North) >> writeOGR(gpsori1, dsn="MATA20090729", >> layer="MATA20090729gps1",driver="ESRI Shapefile") > Error in writeOGR(gpsori1, dsn = "MATA20090729", layer = "MATA20090729gps1", > : > unknown data type > >> summary(gpsori1) > Object of class SpatialPointsDataFrame > Coordinates: > min max > coords.x1 453020 457145 > coords.x2 4623500 4625866 > Is projected: NA > proj4string : [NA] > Number of points: 2096 > Data attributes: > No UTM.Zone UTM.Ch UTM.East > Min. : 1.0 Min. :31 Mode:logical Min. :453020 > 1st Qu.: 524.8 1st Qu.:31 TRUE:2096 1st Qu.:453601 > Median :1048.5 Median :31 NA's:0 Median :454140 > Mean :1048.5 Mean :31 Mean :454418 > 3rd Qu.:1572.2 3rd Qu.:31 3rd Qu.:455120 > Max. :2096.0 Max. :31 Max. :457145 > UTM.North Altitude Date Time > Min. :4623500 Min. :1648 Length:2096 Length:2096 > 1st Qu.:4624886 1st Qu.:1795 Class :character Class :character > Median :4625170 Median :1821 Mode :character Mode :character > Mean :4625068 Mean :1835 > 3rd Qu.:4625338 3rd Qu.:1885 > Max. :4625866 Max. :1952 > > Perhaps the NA for proj4string ? > > Agus > > > Michael Sumner wrote: >> >> Hi Agustin, you will need to determine which columns in the data frame >> are causing problems. >> >> Can you run this and let us know what you see? >> >> summary(gpsori1) >> >> For example, you cannot write POSIXct columns to SHP and would need to >> convert to text, or numeri: >> >> library(rgdal) >> d <- data.frame(x = 1:10, y = 1:10, z = 1:10) >> coordinates(d) <- ~x+y >> >> ## OK >> writeOGR(d, ".", "test", "ESRI Shapefile") >> unlink("test.shp") >> >> d$z <- Sys.time() + 1:10 >> ## not OK >> writeOGR(d, ".", "test", "ESRI Shapefile") >> >> >> Regards, Mike. >> >> >> >> On Thu, Nov 5, 2009 at 7:11 AM, Agustin Lobo <[email protected]> >> wrote: >>> >>> I'm puzzled with this: >>> >>>> class(gpsori1) >>> >>> [1] "SpatialPointsDataFrame" >>> attr(,"package") >>> [1] "sp" >>>> >>>> writeOGR(gpsori1, dsn="CAST20090907", >>>> layer="CAST20090907gps1",driver="ESRI Shapefile") >>> >>> Error in writeOGR(gpsori1, dsn = "CAST20090907", layer = >>> "CAST20090907gps1", >>> : >>> unknown data type >>> >>> Am I doing something wrong or is there a recent change that I'm not aware >>> of? >>> >>> Thanks >>> >>> Agus >>> >>> >>> _______________________________________________ >>> R-sig-Geo mailing list >>> [email protected] >>> https://stat.ethz.ch/mailman/listinfo/r-sig-geo >>> >>> >> > > -- > Dr. Agustin Lobo > Institut de Ciencies de la Terra "Jaume Almera" (CSIC) > LLuis Sole Sabaris s/n > 08028 Barcelona > Spain > Tel. 34 934095410 > Fax. 34 934110012 > email: [email protected] > http://www.ija.csic.es/gt/obster > _______________________________________________ R-sig-Geo mailing list [email protected] https://stat.ethz.ch/mailman/listinfo/r-sig-geo
